موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-ردیابی موج با استفاده از اصل زمان همرسی و معکوس زمان سیرهای باقی مانده10862FAمحمدعلىریاحىسعیداکبرىJournal Article19700101In recent years, several methods for calculating the seismic wave’s travel time have been presented (Moser 1989, Vidale 1988). Among them, the finite difference ray tracing has given reliable results. In this study in order to calculate the travel time in a cell model, Vidal’s finite difference eikonal solver was used. To evaluate, different velocity models were tested and in all cases the results were compatible with the theoretical existing basis.
Then by the damped least square method and conjugate gradient algorithm, which was introduced by Scales (1987), the inversion of the travel time residuals for producing a velocity section of the subsurface media was performed.
In order to qualify the final section, statistics parameters such as resolution matrix, covariance matrix, uncertainty matrix and rms, were used.
This method has been applied on the real data from a weighted concrete dam. At the bottom of the final tomogram a low velocity zone is detected which was attributed to poor quality of the existing rock in that zone. The back projection method was used on the real data to evaluate the method that has been applied in this article. The result obtained from the back projection method compared to the conjugate gradient method was not satisfactoryدر سال های اخیر روش های مختلفی برای محاسبهء زمان سیر امواج لرزه ای در محیط های پیچیدهء سرعتی ارائه شده است (موزر، 1989،ویدال، 1988). در این میان ردیابی موج با روش تفاضل متناهی (finite difference) نتا یج قابل قبولی داشته است. در این مطالعه برای محاسبهء زمان سیر امواج لرزه ای دریک مدل سلولی، از تقریب تفاضل که ویدال (1988) ان را ارائه کرده برای حل معادله آیکونال استفاده شده است. به منظور ارزیابی کارایی این روش، انواع صتفا، تی از مدل های سرعتی مورد بررسی قرار ثرفت که در همهء موارد نتا- به دست آمده با مبانی نظری موجود تطابق لازم را داشت. سپس با استفاده از روش کم ترین مربعات میرا والگوریتم گرادیان بندادی معرفی شده از سوی اسکیلز (1987) مرحلهء معکوس باقی مانده های زمان سیر برای ایجاد یک مقطع سرعتی از محیط زیرسطحی، صورت گرفت. براک ارزیابی کیفیت مقطع نهایی به دست آمده، از شاخص های آماری ماتریس دقت، ماتریس کوواریانس، ماتریس عدم قطعیت، جذر میانگین مربعی خطا استفاده شده است.
این روش در مورد داده های حقیقی به دست آمده از محل زیرساخت یک سد وزنی بتونی به کار برده شده است. در توموگرام نهایی به دست آمده، در قسمت پایین، یک منطقهء کم سرعت مشاهده شد که ممکن است بتوان آن را به سست بودن سنگ های تشکیل دهنده این زون نسبت داد. به منظور ارزیابی روش مورد استفاده در این مقاله، از روش بک پروجکشن ((back projection در محدودهء مورد مطالعه استفاده شد و توموگرام نهایی به دست آمده از این روش با توجه به شاخص جذر میانگین مربعی، خطا به دست آمده در مقایسه با روش زمان همرسی (reciprocal time) از دقت قابل قبول برخوردار نبود.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-مقایسهء نتایج حاصل از روش های مختلف پیش بینی امواج ناشی از باد در خلیج فارس (منطقه بوشهر)10863FAکامرانلارىعباسعلى على اکبرىبیدختىمهدى شفیعىفرJournal Article19700101مطالعهء پارامترهای میدانی همراه با دید و درک فیزیکی مسائل، در مسیر روشن کردن برخی از پدیده های فیزیک دریا بسیار مفید است. در این تحقیق پارامترهای مختلف میدانی ثبت شده توسط بویه موج نگار سازمان هواشناسی کشور در خلیج فارس، سواحل بوشهر (28,57’ شمالی و "30، 50 شرقی) از جمله مشخصه های امواج، باد، دمای هوا، سطح آب درطی سال های 1999.1996,1995 که در حدود 0 0 39 رویداد است مورد مطالعه قرار گرفته است. فراوانی وقوع سه حالت کلی محدو د شدن طول موجگاه، محدود شدن مدت زمان وزش باد و موج کاملأ رشد یافته، همچنین فراوانی وقوع شرایط پایدار، ناپایدار، خنثی لایهء سطحی جو درطی سال 1995 محاسبه شد که نشان داد بیشتر شرایط محدود شدن مدت زمان وزش باد، هم چنین ناپایداری لایهء سطحی جو بر این منطقه حاکم است. روابط به دست آمده بین انرژی €، بسامد v بدون بعد برحسب طول موجگاه بدون بعدX در منطقه بوشهر همخوانی خوبی با نتیجهء سایر محققان دارد. همین طور با درنظر گرفتن پایداری لایهء سطحی جو (عدد ریچاردلون توده ای Rb) نشان داده شد که برای 0.091Rb0.88,v=0.056 Rb-0.03,(0.005<Rb<1) =?. است. این گونه روا بط را می توان در مدل های پیش بینی موج اعمال کرد ودقت مدل ها را افزایش داد.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-مدل سازی وارون دوبعدی داده های گرانی10864FAعلى نجاتىکلاتهوحید ابراهیم زادهاردستانىJournal Article19700101This study deals with modeling of gravity anomalies by Compact Inversion Method of Last and Kubic (1983) and Lewi et al (1994). The principle used is to minimize the volume, to maximize of the causative body. The anomalous density distribution is obtained using an iterative technique which is numerically stable and rapidly convergent. This principle can also be adapted to include modeling of gravity anomalies by single density sources.
The advantage of this method is that desirable geological characteristics are automatically incorporated into the model with minimum subjective judgments on the part of interpreter. The treatment of noise in data fits naturally into formulations of the inversion procedure
The practical effectiveness of this method is assessed by inversion of synthetic and real acquired over chromite ore in Hormozgan province of Iran.در این مقاله از روش مدل سازی خطی نوینی به نام وارو ن سازی فشرده (compact inversion) برای تفسیر داده های گرانی سنجی استفاده شده است. مبنای کار در رو ش وارون فشرده کمینه کردن حجم چشمهء بی هنجاری است. این امر معادل بیشینه کردن فشردگی (compactness) چشمهء بی هنجاری است. توزیع چگالی با استفاده از روشی بر مبنای تکرارهای متوالی با همگرایی قابل ملاحظه محاسبه می شود. این روش را می توان به آسانی برای مدل هایی که در آنها چگالی بی هنجاری را می توان به صورت مقداری معین در کل توده در نظر گرفت، به کار برد. در این روش رفتار نوفه در داده ها برای وارون سازی داده ها در نظر گرفته شده است. ابتدا به توضیح روش، سپس به کاربرد عملی آن در مورد داده های مصنوعی و وا قعی خواهیم پرداخت.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-تضعیف نوفه های تصادفی از داده های برانبارش شدهء لرزه نگاری در حوزهء F-XY10865FAحسین عالىدانشورحمیدرضاسیاهکوهىJournal Article19700101One of the main problems in geophysics data analysis is the presence of noise. This problem in seismic survey is more obvious than in the other branches.
In this study the random noise suppression is presented by a filter which is called eigenimage filter and operates on stacked 3D seismic data in frequency domain. Our tools in suppression of random noises were SVD and Lanczos. The Lanczos method works much faster than SVD, specially when we have a sparse matrix. The special feature of the Lanczos method is its high performance. The F-xy filter has its own abilities such as fantastic signal preservation and well noise suppression, and can be used even before stacking. In this study presented performance of the filter on staked synthetic seismic data.حضور نوفه ها در لرزه نگاری امری اجتناب ناپذیر است. پاره ای از این نوفه ها ماهیت تصادفی دارند. و پاره ای همدوس هستند. به منظور ارائه تصویر صحیح از ساختارهای زمین شناسی منطقه لازم است. این نوفه ها در مراحل پردازش تضعیف شوند. در این مطالعه برای تضعیف نوفه های تصادفی روش ویژه تصویر در حوزه F-XY معرفی می شود. هدف این مطالعه ‘ تضعیف نوفه های تصادفی از مکعب داده های لرزه ای بر انبارش شده در لرزه شناسی سه بعدی به روش ویژه تصویر در حوزه F-XYاست در این روش عمل فیلتر کردن با جایگزینی مقاطع بسامد ثابت به وسیله مجموع تعداد معینی از ویژه تصاویر دارای بیشترین سهم انجام می شود. این روش در مقاطع با لایه بندی افقی و مقاطع با لایه بندی شیب دار کارایی خوبی دارد. این روش برخلاف فیلتر پیشگو‘ مستقل از بسیاری از خواص وابسته به XوY است مانند ترتیب وتصحیحات استاتیک ‘ این روش در نقاط مرزی محدوده برداشت سه بعدی نیز به خوبی عمل می کند . زمان اجرای این فیلتر با استفاده از روش لنکزوس به جای SVD به مراتب کاهش یافته و گزینه مناسبی در مقایسه با فیلترهای مشابه خواهد بود. کارایی فیلتر روی داده های مصنوعی اجرا و نتایج ارائه شده است.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-بررسی نمودارهای پاشش امواج ریلی10866FAسید حمیدسیدینعبدالرحیمجواهریاننصراللهکمالیانJournal Article19700101Rayleigh waves are surface waves which are dispersive in layered media. In order to study Rayleigh wave dispersion, there are several methods. Among these methods the application of matrix methods for elastic layered media is preferred to other methods. In the matrix method for a layered medium, boundary conditions for different interfaces are satisfied to obtain a relation for calculating Rayleigh wave dispersion which should be resolved by numerical methods. A program for the computing of Rayleigh waves dispersion was developed. The resulted dispersion curves from this program were compared to similar methods which showed good agreement. According to the resulted graphs, wherever the earth structure is complex the dispersion curves also will be complex and different models will be moved oscillatory. Phase and group velocity curves tend to lower values when there are a surface sedimentary layer in the earth model. Thus, synthesis of these curves leads to clarification of different seismic phases, history of waves propagation from source to receiver and consequently imaging of earth structureامواج ریلی امواجی سطحی به شمار می روند که در محیط های لایه ای، پاششی اند. برای تحلیل پاشش امواج یادشده،روش های متعددی وجود دارند که از این بین استفاده از روش های ماتریسی برای محیط های لایه ای کشسان بر رو ش های دیگر برتری دارد. درروش ماتریسی برای محیط های لایه ای، شرایط مرزی مربوط به مرزهای مختلف به یکدیگر مرتبط می شوندو رابطه ای برای محاسبهء پاشش امواج ریلی به دست می آید که باید با استفاده از محاسبات عددی حل شود. در این مقاله با استفاده ازروش ماتریسی برنامه ای رایانه ای برای محاسبهء پاشش امواج ریلی نوشته شد. نمودارهای پاشش حاصل از برنامهء مذکور با نمودارهای منتشر شدهء نظیر مقایسه شد که همخوانی کاملی بین آنها مشاهده شد. با توجه به نمودارهای یادشده مشاهده می شود که هم زمان با پیچیده تر شدن ساختار زمین، نمودارهای پاشش نیز پیچیده تر می شوندو بر نوسانات مدهای مختلف افزوده می شود. همچنین وجود لایهء رسوبی در سطح زمین موجب گرایش نمودارهای سرعت فازو سرعت گروه به سمت مقادیر کم تر می شود. تفسیر نوسانات مذکور کمک مؤثری برای تشخیص فازهای مختلف لرزه ای، تعیین ساختار زمین و سرگذشت عبور امواج از چشمه تا گیرنده است.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-مطالعهء موردی جزیرهء گرمایی تهران و شبیه سازی عددی آن10867FAعباس رنجبرسعادتآبادىمجیدآزادىعباسعلى على اکبرىبیدختىسیدعلیرضا صادقىحسینىJournal Article19700101Temperature and wind fields over Tehran were simulated using the MM5 modeling system. To do this a nested domain with three nests with grid spacing of 27, 9 and 3 km were chosen. The model was run for two different days, one with a continental polar air mass (CPA) and the other one with a continental tropical air mass (CTA) over the region. The model outputs for the innermost nest were studied for the first 24 hours of the forecast. The model results were compared with observed 2m temperature and l0m wind speed at Mehrabad, Aghdasieh and Chitgar stations. The results show a fairly good agreement with the corresponding observations. The results also show that during the time of CPA, the Tehran heat island is well formed and has a relatively high intensity in the northern part of the city while during the time of CTA the heat island does not grow during the night time and has a lower intensity as well. Examining the model results also shows that the model is capable of capturing the micro and mesoscale characteristics features of the heat island and catabatic and anabatic winds over the city. Also, the results show that the mountain cooling intensities the urban heat island circulations at night and these circulations reduce heat island intensity.در این مطالعه جزیرهء گرمایی تهران با استفاده از یک مدل منطقه ای محدود برای پیش بینی عددی وضع هوا، موسوم به مدلMM5شبیه سازی شده است. برای این منظور سه دامنه با تفکیک افقی 27، 9و3 کیلومتر در نظر گرفته شد. مدل MM5برای 2روز مختلف، یکی در زمان استقرار تودهء هوای قاره ای قطبی (تودهء هوای سرد) و دیگری در زمان استقرار تودهء هوای قاره ای حاره ای (تودهء هوای گرم) اجرا وخروجی های دامنة با درجه تفکیک زیاد یعنی 3 کیلومتر، طی یک دورهء پیش بینی 24 ساعته بررسی شدند. پیش بینی های مدل با مقادیر دیدبانی شدهء سمت و سرعت باد در ارتفاع 10 متری و دما در ارتفاع 2 متری سطح زمین در ایستگاه های هواشناسی مهرآباد، اقدسیه و چیتگر مقایسه شدند نتایج نشان می دهند که الگوی میدان بادو دمای پیش بینی شده با مقادیر متناظر دیدبانی شده در این ایستگاه ها از سازگاری قابل قبولی برخوردار است. بنابراین در صورت پیش بینی میدان های دماو باد، امکان پیش بینی آلودگی هوای تهران نیز می تواند فراهم آورد. همچنین نتایج حاصله از شبیه سازی نشان می دهد که در زمان استقرار تودهء هوای سرد جزیرهء گرمایی به خوبی در نیمهء شمالی تهران تشکیل می شودو از شدت بیشتری نیز برخوردار است، در حالی که در زمان استقرار تودهء هوای گرم، امکان تشکیل و توسعهء جزیرهء گرمایی در ساعات شب به خوبی فراهم نمی شودو در نتیجه جزیرهء گرمایی از شدت ضعیفی برخوردار است. همچنین تغییرات آهنگ های سرمایش در شب می تواند سبب ایجاد گردش های حرارتی شود، که این گردش ها از افزایش شدت جزیرهء گرمایی جلوگیری می کنند.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-پیش بینی نفوذپذیری از روی داده های چاه نگاری با استفاده از شبکه های عصبی مصنوعی در یکی از مخازن کربناته جنوب ایران10868FAبنفشهحبیبیانمجید نبىبیدهندىعزت الله کاظمزادهJournal Article19700101The study presented here is an artificial neural network modeling for prediction of permeability in one of the limestone reservoirs in the south of Iran. The permeability ANN is a three layer network. We have used two separate ANNs to model permeability. Initially, density, gamma ray, neutron and sonic logs were applied for input. Then depth data related to these data were added to the input; the results of the two networks have then been compared. This comparison has shown that in the second state the accuracy of the model has been improved significantlyتخمین نفوذپذیری از روی داده های چاه نگاری، یک مرحلهء اساسی در مهندسی مخزن و ارزیابی سازنده است. تحقیقی که در این جا عرضه می شود، نوعی مدل سازی شبکه های عصبی مصنوعی، به منظور استفاده از نگارهای چاه برای پیش بینی نفوذپذیری، در یکی از مخازن کربناته جنوب ایران است. در این تحقیق دو شبکه پس انتشار سه لایه، برای پیش بینی نفوذپذیری، مورد استفاده قرار گرفت. شبکه او ل، نگارهای نوترون، چگالی، صوتی و پرتو گاما را به صورت ورودی به کار می برد; حال آن که در شبکهء دوم، عمق های مربوط به این داده ها نیز به ورو دی ها اضافه و نتایج در دو حالت مقایسه شده اند. بر اساس این مقایسه، مشخص شد که دقت مدل در حالت دوم بهبود قابل ملاحظه ای یافته است.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-محاسبات سرعت فاز و سرعت گروه با استفاده از نگاشت های دورهء بلند10869FAمحمدکمال قاسمالعسکرىJournal Article19700101In this investigation long-period records of earthquake data were used for separation of surface waves and calculation of phase and group velocity dispersion curves.
Kazeroun main shock (3 May, 1995) recorded by Iranian Long Period Array (ILPA) at station 1, 2, and 7; SantHellana main shock (27 April, 1986) recorded by two very long-period European stations (KON in Spain and TOL in Norway) have been processed. First, we have corrected all records by special techniques for the arm and tilting of instrumentation. Then, we processed the earthquake data for separation of surface waves (Love and Rayleigh waves) by two dimensional frequency-time algorithms. After separation of surface waves from main shocks full automotive intelligent programs have been used for detection and computation of phase and group velocity dispersion curves.
Results show that phase velocities calculated from three-components of ILPA data at stations 1, 2, and 7 are less than the group velocities for the same stations (inverse dispersion). Vertical components of the three stations show that the phase velocity with respect to frequency is increased, but for the horizontal components (N-S and E-W directions) the phase velocities with respect to frequency are decreased. This behavior shows that the phase velocity of shear waves is reduced with depth of penetration. Besides, the higher modes of surface waves are mimicked as they pass through the Zagros belt from the earthquake epicenter to ILPA recording stations. Similar computations of phase and group velocities from the TOL and the KON data stations show that the phase velocities are greater than the group velocities (forward dispersion).در این بررسی داده های زمین لرزه های دوره بلند تحلیل شد و برای محاسبه سرعت گروه و سرعت فاز مورد استفاده قرار گرفت. نگاشت های سه مولفه ای زمین لرزه سوم مه 1995 کازرون در ایران که پایگاه ایلیا (ILPA) شبکه بین المللی زلزله نگاری مؤسسه ژئوفیزیک دانشگاه تهران ثبت کرده بود. وهمچنین نگاشت های دوره بلند زمین لرزه سنت هلنا (بیست و هفت هفتم آوریل 1986 ) در دو ایستگاه زلزله نگاری KON اسپانیا و TOL نروژ مورد بررسی قرار گرفت. و برای محاسبات سرعت گروه و سرعت فاز از آنها استفاده شد. نگاشت های یاد شده که با دستگاه های زلزله نگار دوره بلند ثبت شده اند پس از تصحیحات لازم (تصیح برای دستگاه ها و تصیح انحراف محور) با روش طیفی و تحلیلی در صفحه بسامد- زمان پردازش شدند و فازهای امواج سطحی از آنها جداسازی شد. با استفاده از تحلیل طیفی محاسبات سرعت گروه و سرعت فاز بین پایگاه های یاد شده برای هر کدام از نگاشت ها با برنامه ای خودکار صورت گرفت. ومنحنی های پاشیدگی فازهای امواج سطحی از نگاشت های مربوط محاسبه شدند. نتایج به دست آمده از محاسبات سرعت گروه و سرعت فاز نگاشت های سه مؤلفه ای پایگاه ایلیا نشان می دهد . که سرعت فاز بین سه پایگاه شماره 1 و 2 و 7 از سرعت گروه محاسبه شده درهر سه پایگاه کمتر است. (پاشیدگی وارون ‘ Inverse dispersion)) ونگاشت های زلزله بیست و هفت آوریل 1986 یک نوع پاشیدگی نرمال (Normal dispersion) رانشان می دهد. به علاوه سرعت فاز محاسبه شده بین مؤلفه های قائم سه پایگاه ایلیا با بسامد امواج سطحی افزایش قابل ملاحظه ای را نشان می دهد. در صورتی که سرعت فاز محاسبه شده بین مولفه های افقی سه پایگاه (شرقی –غربی و شمالی - جنوبی) با بسامد امواج سطحی کاهش می یابد. این گونه رفتار در منحنی های پاشیدگی نشان می دهد که سرعت فاز امواج برشی بر حسب عمق کاهش قابل ملاحظه ای دارد. به علت عبور امواج سطحی از رشته کوه های زاگرس درامتداد رومرکز زمین لرزه سوم مه 1995 کازرون نسبت به پایگاه ایلیا بیشترین مقدار انرژی بسته ای مدهای بالایی امواج سطحی حذف شده است. (به علت پدیده های تکتونیکی شدید) به علاوه عبور اموج سطحی از دو ناحیه با اختلاف امپدانس کاملا متفاوت (آب به خشکی و بالعکس) مدهای بالایی امواج سطحی از نگاشت های زلزله 27 آوریل 1986 در دو ایستگاه و انرژی بسته ای خود را کاملا ازدست داده اند وضعیف شده اند و به همین دلیل روی صفحه بسامد –زمان کمتر ظاهر شده اند . لذا هر قدر فاصله پایگاهها از کانون زمین لرزه بیشتر باشد. تراکم امواج سطحی با طول موج زیاد یا تناوب بلند به نسبت زیادتر خواهد بود و اطلاعات به دست آمده از عمق های بیشتری امکان پذیر است. از طرف دیگر خطا درمحاسبات سرعت گروه وسرعت فاز از نگاشت های اصلی با پربود بلند افزایش می یابد. با کاربرد الگوریتم عرضه شده می توان نسبت این گونه خطاها راکاهش داد.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-تعیین فاکتور کیفیت امواج برشی مستقیمQ?، در منطقه البرز به کمک داده های میدان نزدیک حرکت نیرومند زمین لرزهء کجور در محدودهء بسامدی ا تا 32 هرتز10870FAهادىقاسمىنصراللهکمالیانحسین حمزهلوعلى بیتالهىJournal Article19700101Near-field data obtained from the main shock of May 28, 2004 Baladeh-Kojour earthquake, with Mw = 6.2 is used
to determine the frequency dependence of direct shear wave quality factor Q?, in the range of 1 to 32 Hz.
Based on the regression analysis, the following equation is proposed: Q? = 90f0.74
High attenuation factor and frequency dependence of quality factor for shear waves imply the high seismicity of the region and agrees well with similar studies in other seismically active regions such as Greece, Japan, India, etc.داده های میدان نزدیک حرکت نیرومند زمین مربوط به تکان اصلی زمین لرزه هشتم خرداد 1383 بلده- کجور با بزرگای گشتاوری MW=6.2(NEIC) به منظور تعیین وابستگی بسامدی ضریب کیفیت امواج برشی مستقیم Q?، در محدودهء بسامدی ا تا 34 هرتز مورد استفاده قرار گرفته است. بر اساس آنالیز رگرسیون انجام گرفته رابطهء زیر برای،وابستگی بسامدی Q? پیشنهاد شده است:
Q=90f
ضریب میرایی بالا و نیزوابستگی شدید بسامدی مشاهده شده، بیانگر لرزه خیزی بالای منطقه مورد مطالعه بوده، با نتایج مطالعات مشابهی که در سایر نقاط لرزه خیز دنیا مانند یونان، ژاپن , هندو ستان،... صورت گرفته است همخوانی دارد.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-حل عددی مدار ماهواره های نزدیک سطح زمین با گام متغیر10871FAمهدىاسحاقJournal Article19700101Investigation of the step variable numerical integration methods is the main purpose of this paper. Also it is shown how the one-dimensional numerical integration method of solving an ordinary differential equation can be employed to solve a three-dimensional differential equation numerically. In order to integrate the orbit, the step variable methods of Runge-Kutta-Fehlberg and Adams are employed and their accuracies are investigated.بررسی روش حل عددی مدار ماهواره موضوع اصلی این تحقیق است. همچنین در این مقاله نشان داده می شود که روش حل معادلات دیفرانسیل مرتبه اول یک بعدی برای حل عددی یک معادلهء دیفرانسیل سه بعدی چگونه به کار گرفته می شود. برای حل عددی مدار روش های با گام متغیر رونگ-کوتا- فلبرگ و روش آدامز به کارگرفته می شود و در نهایت دقت آن ها مورد بررسی قرار می گیرد.موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X31120050321-مقیاس بزرگی در شبکهء لرزه نگاری تبریز10872FAمهدىرضاپور.Journal Article19700101To record local earthquakes, the telemetric-digital seismic network of the Tabriz comprising of eight three component seismic stations was installed in the north-west of Iran. This array started operation at the end of 1995. Investigation of archived data in the Tabriz seismic network between 1999 to September 2004 shows that magnitude values in the database have not been determined using a single formula or a specific method. Comparison of magnitude values in the database with those values that have been published in the ISC/NEIC bulletin shows that magnitude values calculated in the Tabriz network are underestimated for events occurred in greater distances. By using magnitude values in the database and corresponding mblSC/NETC (mbISC or mbNETC) values, a scale as M log (v/4?)+2.6log (?)-2.2 is derived for determination of the magnitude in distance range of 170 to 1000 km, where v is peak-to-peak amplitude in micrometer/second and A is epicentral distance in kilometers. This formula gives a better estimate for the magnitude of events in comparison with the formula that is presently used in the Tabriz seismic network.شبکهء لرزه نگاری دیجیتالی تبریز مرکب از هشت ایستگاه لرزه نگاری سه مؤلفه ای، برای ثبت زمین لرزه های محلی در شمال غرب ایران نصب شده است. این شبکه در اوا خر سال 1995 شروع به فعالیت کرد. بررسی داده های آرشیو شده در فاصلهء زمانی 1995 تا دسامبر 2004 نشان می دهد که برای تعیین بزرگی در این شبکه از یک روش و فرمول واحد استفاده نشده است. مقایسهء مقادیر بزرگی تعیین شده برای زمین لرزه های ثبت شده در شبکهء لرزه نگاری تبریز با آنچه که دربولتن ISC/NEIC منتشر شده است نشان می دهد که در شبکهء لرزه نگاری تبریز بزرگی زمین لرزه های رخ داده در فواصل دورتر کم تر از أنچه که باید باشد برآورد می شوند. با استفاده از مقادیر بزرگی تعیین شده در شبکهء تبریز، مقادیرmb یک فرمول به صورت M=log(?/4?)+2.6log(?)-2.2کهی با دامنه (peak-to-peak) بر حسب میکرومتر در ثانیه و ? فاصله رومرکزی برحسب ثانیه است برای تعیین بزرگی رویدادهای واقع در فواصل بیشتر از 170 کیلومتر از ایستگاه لرزه نگاری ارائه شد. بررسی صورت گرفته نشان می دهد که مقادیر بزرگی برآ ورد شده با استفاده از این رابطه در مقایسه با فرمولی که اکنون در شبکه تبریز استفاده می شود، دقیق تر است.