دوره 15، شماره 1 - ( فصلنامه علمی تخصصی طب کار یزد 1402 )                   جلد 15 شماره 1 صفحات 57-37 | برگشت به فهرست نسخه ها

Ethics code: 56839


XML English Abstract Print


Download citation:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:

Torfi Alivi A, Shahraki F, Sardashti Birjandi M R, Khalilipour M M. Consequences modeling and determining of safe distance in the natural gas pressure reduction station using PHAST software (Case study: Borumi station in Ahvaz). tkj 2023; 15 (1) :37-57
URL: http://tkj.ssu.ac.ir/article-1-1221-fa.html
طرفی علیوی احمد، شهرکی فرهاد، سردشتی بیرجندی محمد رضا، خلیلی پور میر محمد. ارزیابی پیامد و تعیین حریم ایمن در ایستگاه تقلیل فشار گاز طبیعی با استفاده از نرم‌افزار PHAST (مطالعه موردی: ایستگاه برومی شهر اهواز). فصلنامه علمی تخصصی طب کار. 1402; 15 (1) :37-57

URL: http://tkj.ssu.ac.ir/article-1-1221-fa.html


دانشگاه سیستان و بلوچستان ، m.r_sardashti@eng.usb.ac.ir
متن کامل [PDF 1545 kb]   (348 دریافت)     |   چکیده (HTML)  (545 مشاهده)
متن کامل:   (308 مشاهده)
ارزیابی پیامد و تعیین حریم ایمن در ایستگاه تقلیل فشار گاز طبیعی با استفاده از نرم‌افزار PHAST
(مطالعه موردی: ایستگاه برومی شهر اهواز)


احمد طرفی علیوی[1]، فرهاد شهرکی[2]، محمدرضا سردشتی بیرجندی[3]*، میر محمد خلیلی پور لنگرودی[4]

چکیده
مقدمه: ایستگاه تقلیل فشار گاز یکی از مهم‌ترین بخش شبکه‌های توزیع گاز طبیعی هستند. ارزیابی ریسک در واحدهای عملیاتی روش مناسبی جهت ارزیابی خطرات بوده و نتایج آن‌ها را می‌توان جهت مدیریت و تصمیم‌گیری در خصوص کنترل و کاهش پیامدهای آن بدون نگرانی به کاربرد. به‌عبارتی‌دیگر ارزیابی ریسک شامل فرآیندی است که به مشخص کردن مخاطرات، ارزیابی زیان‌ها و تعیین مشخصات ریسک می‌پردازد.  هدف از مطالعه حاضر مدل‌سازی حوادث محتمل و متعاقب آن ارزیابی پیامد و تعیین حریم ایمن در ایستگاه تقلیل فشار گاز طبیعی برون‌شهری واقع در منطقه برومی شهر اهواز با استفاده از نرم‌افزار PHAST بوده است.
روش بررسی: در این مطالعه با استفاده از نرم‌افزار PHAST نسخه 22/8 و با استفاده از داده‌های فرآیندی ایستگاه تقلیل فشار گاز موردمطالعه نظیر دما، فشار گاز، آنالیز گاز، و داده‌های هواشناسی و استفاده از طبقه‌بندی شرایط جوی پاسکوئیل با تعیین سناریوهای محتمل بر اساس پیشینه موردبررسی تعیین حریم ایمن در این ایستگاه انجام گرفت. برای این مهم 6 سناریو محتمل در کلاس‌های جوی انتخاب‌شده به همراه بدترین شرایط ممکن بررسی گردید. ارزیابی پیامدها نیز بر اساس حالت‌های مختلف تخلیه گاز، انتشار گاز، آتش فورانی، آتش ناگهانی، انفجار موردبررسی قرار گرفت.
نتایج: بر اساس نتایج به‌دست‌آمده محدوده‌های با پتانسیل خطر بالا به ترتیب مربوط به نتایج انتشار و انفجار گاز بوده است. همچنین بیشترین محدوده خطر متعلق به سناریوی شماره پنج برای انفجار گاز باکلاس جوی F و شدت موج bar 02068/0 و معادل با 5/358 متر گزارش گردید. در این محدوده با توجه به جدول خسارت ناشی از انفجار با احتمال 95% بدون خسارت جدی 50% شیشه پنجره‌ها شکسته می‌شود و بدین‌جهت می‌تواند سبب صدمات انسانی گردد. بر این اساس حریم ایمن در این ایستگاه برابر 5/358 متر تعیین گردید.
نتیجه‌گیری: به سبب عدم مطالعه ارزیابی ریسک در فازهای اولیه و در حین بهره‌برداری مطالعه اخیر نشان می‌دهد تعیین نقاط ایمن و محدوده با ریسک بالا به همراه اولویت کاهش ریسک بسیار ضروری بوده که در برخی موارد در واحدهای صنعتی مقفول مانده است لذا می‌توان نتیجه گرفت که مدل‌سازی حوادث و ارزیابی پیامد ناشی از آن به همراه تعیین حریم ایمن یکی از فاکتورهای اساسی در سلامت محیط کار پرسنل و عموم مردم محسوب می‌شود.
واژه­های کلیدی: انتشار گاز سمی، مدل‌سازی انفجار، مدل‌سازی حوادث فرایندی، ارزیابی کمی ریسک، نرم‌افزار PHAST

مقدمه
 
 با توسعه صنعت نفت و گاز به دلیل کاربرد تکنولوژی پیچیده و غیرقابل ‌انعطاف و هزینه‌بر، چالش‌های ایمنی، بهداشت و محیط‌زیستی این صنعت نیز رو به افزایش است و با پیشرفت علم و فناوری در این‌گونه صنایع همیشه افزایش مخاطرات بر خواسته از کار وجود دارد و به‌منظور نظارت بر این ریسک‌ها احتیاج به سیستم مدیریتی است که باعث کم‌کربن ریسک و به‌دست‌آمده اطمینان از بهبود ایمنی، رفاه کارکنان و نیز نگهداری از محیط‌زیست می‌شود (1). ایجاد حوادث در فرآیندهای شیمیایی که باعث به ‌وجود آمدن فجایع انسانی و محیطی می‌شود متخصصان را بر آن داشته که برای برآورد اثرات این حوادث به روش‌های پیش‌بینی کننده روی‌آورند. با استفاده از این روش می­توان پیش از به وجود آمدن حادثه به‌منظور کنترل آن برنامه‌ریزی نمود. خصوصیات صنایع بزرگ فناوری  مانند صنعت نفت و گاز این است که حجم خیلی بزرگی از مواد و کالاهای خطرناک در یک بخش قرار می‌گیرند که به‌وسیله اپراتورها اداره می‌شوند لذا اتفاقات در این بخش‌ها کارکنان را تهدید می‌کنند (2). با توجه ‌به توسعه صنایع به‌ویژه صنایع نفت و گاز در دنیا و در رابطه ‌با حمل‌ونقل مواد ریسک‌ها و پیامدهای آن‌ها روزبه‌روز بیشتر می‌شود. شناخت درست ریسک‌ها و عوامل آسیب‌زا در این صنعت و مدیریت ریسک ازنظر کاهش زیان‌ها و آسیب‌های فیزیکی و اقتصادی ضروری است (3). انتشار مواد دارای پتانسیل اشتعال و انفجار و یا ویژگی سمی بودن در این صنعت یکی از تهدیدکننده‌های اصلی کارکنان، ساکنین مجاور این صنایع و محیط‌زیست می‌باشد.  اتفاقات در صنعت نفت و گاز همیشه با آغاز یک واقعه به وقوع می‌پیوندد که این اتفاقات عبارتند از ترکیدن یا شکستن خط لوله، سوراخ در مخزن و یا ایجاد واکنش‌هایی که خارج از کنترل هستند، می­باشند. بر اساس این وقایع مقدار زیادی از مواد انتشار می‌یابد و با نشت این مواد امکان به وجود آمدن حوادثی دیگر از قبیل آتش‌سوزی، انفجار و مسمومیت پیش می‌آید. از طریق بررسی این حوادث به این نتیجه می‌رسیم که امکان ایجاد حادثه و زیان‌های مرتبط به آن همواره مشابه نمی‌باشد و در چارچوب کشوری یا منطقه‌ای خاص قرار نمی‌گیرد. از سویی توسعه فناوری نیز از عدم رخداد این وقایع نمی‌تواند جلوگیری کند و همیشه امکان به‌وجودآمدن این‌گونه حوادث وجود دارد (4). PHAST (نرم‌افزار تحلیل خطر فرآیند) یک ابزار جامع تحلیل پیامدها است. این نرم‌افزار یک حادثه احتمالی را از نشت اولیه تا پراکندگی در محیط به‌صورت مدل‌سازی تبخیر و تبخیر استخر و اثرات قابل اشتعال و سمی بررسی می‌کند. PHAST قادر به شبیه‌سازی سناریوهای انتشار مختلف مانند نشت، پارگی خط، آزاد شدن طولانی خط لوله و فروریختن سقف مخزن در مخازن یا لوله‌های تحت‌فشار و بدون فشار می‌باشد.
در صنایع فرایندی و ایستگاه‌های تقلیل فشار گاز به سبب ماهیت وجود مواد قابل اشتعال و سمی ممکن است ماده‌ای از یک واحد فرآیندی رهاشده، در هوا منتشر شود و چنان خطری ایجاد کند که ایمنی کارکنان و حتی افرادی که در فاصله زیادی از تأسیسات فرآیندی هستند را در معرض خطر قرار دهد. همین عوامل منجر می‌شود صنایع با گسترش اندازه و تعداد تجهیزات، روزبه‌روز در مورد خسارات بالقوه انسانی و اقتصادی نگرانی بیشتری پیدا کنند. برخی حوادث داخلی و خارجی که مربوط به ایستگاه گاز و بیشتر در خطوط لوله اتفاق افتاده است کد در بخش‌های بعدی بیشتر به آن اشاره‌شده است. از مهم‌ترین این حوادث می‌توان به انفجار گاز شهری با پنج مصدوم و دو نفر کشته در قیام دشت اشاره نمود. همچنین انفجار خطوط لوله در شهریار تهران، غرب کشور، خیابان لاهور شهر اصفهان، آتش‌سوزی خط لوله منتهی به ایستگاه شهری در شهر گناوه در داخل کشور و انفجار خطوط لوله ایستگاه انتقال گاز در شهر رنشوفوجیا  با مرگ 10 نفر و آسیب به 40 نفر، انفجار ایستگاه گاز شهری در بیجینگ چین با 47 نفر کشته و انفجار ایستگاه گاز شهری در اتریش با مجروح شدن بیش از 60 نفر اشاره نمود (5). با توجه به توضیحات فوق لزوم بررسی و ارزیابی پیامد و تعیین حریم ایمن برای این مناطق از اهمیت بالایی برخوردار است. ایستگاه‌های تقلیل فشار گاز برون‌شهری از واحدهای عملیاتی و فرآیندی با فشار زیاد می­باشند که درنتیجه مخاطرات زیادی در پی دارند. در ادامه به توضیح مختصری از تحقیقات و مطالعات داخلی و خارجی انجام‌شده مرتبط با موضوع مقاله توضیح داده‌شده است. ستاره‌شناس صیقلانی (1389) به بررسی ارزیابی پیامد پخش گاز کلر در تصفیه‌خانه آب رشت و زمین‌های اطراف آن پرداخت. در این واحد کلرین که ماده سمی است به مقدار زیاد برای گندزدایی آب به کار می‌رود که رهایش آن اثرات مخرب زیادی دارد. این مطالعه مکان‌های پرخطر در اطراف تصفیه‌خانه را مشخص می‌کند. نتایج نشان می‌دهد که زمین‌های اطراف این واحد تا فاصله 100 متر باید خالی از سکنه باشد (6). پورقدیری (1392) به ایمنی، بهداشت و محیط‌زیست در ایستگاه‌های تقلیل فشار و خطوط انتقال گاز پرداخت صنعت گاز از مرحله اکتشاف شروع می‌شود و تا مصرف ادامه دارد. اطمینان از صحت عملکرد تجهیزات موجود در ایستگاه‌ها (فیلترها، شیرهای اطمینان و قطع‌کننده جریان و غیره)، اطمینان از عدم نشتی خطوط انتقال و تجهیزات، اطمینان از عدم تشکیل هیدرات و به حداقل رساندن آلودگی‌هایی نظیر آلودگی صوتی و غیره ازجمله نکات مهم HSE (Health, Safety and Environment) هستند که در این مقاله پیرامون آن‌ها بحث شده و راهکارهایی برای آن‌ها ارائه‌شده است (7). موسوی (1393) به بررسی ارزیابی کمی ریسک ناشی از حوادث فرآیندی در ایستگاه تقویت فشار گاز شهرستان رامسر پرداخت. در این تحقیق، ارزیابی کمی ریسک ناشی از حوادث فرآیندی، نظیر نشتی تجهیزات فرآیندی ازجمله فلنج ها و شیرها که ممکن است در ایستگاه تقویت فشار گاز شهرستان رامسر به وقوع بپیوندد موردبررسی قرارگرفته‌اند. این ارزیابی با توجه به اطلاعات جمع‌آوری‌شده نظیر فشار ماده موجود در فرآیند، دما ماده موجود در فرآیند، فاز ماده موجود در فرآیند، ترکیب درصد ماده موجود در فرآیند، شرایط آب و هوایی، نوع سناریو، محل سناریو و جانمایی به‌وسیله نرم‌افزار PHAST مدل گردید. در این پروژه سعی بر آن است تا با مدل‌سازی دو نوع آتش محتمل در ایستگاه تقویت فشار گاز رامسر یعنی آتش ناگهانی و آتش فورانی ناشی از نشتی فلنج ها و شیرهای موردبررسی در سناریوهای تعریف‌شده و با رسم نمودارهای مربوط به آن، محدوده اثر این اتفاقات تعیین گردد (8). قربانی (1393) به بررسی روش‌های مهار سولفید هیدروژن در صنعت نفت و گاز پرداخت. این مقاله شامل گردآوری و بررسی مجموعه‌ای از مطالعات است که در سال‌های اخیر برای حذف و یا به حداقل رساندن مشکلات ناشی از هیدروژن سولفید در عملیات‌های نفت و گاز با بهره‌گیری از برخی مواد شیمیایی تحت عنوان پالاینده هیدروژن سولفید صورت پذیرفته است. چالش بر سر تحقیق و توسعه پالاینده‌ای است که اهداف مدنظر ازجمله سلامتی پرسنل شرکت نفت و گاز، و نگهداری مناسب خطوط انتقال را فراهم نماید، علاوه بر این به لحاظ زیست‌محیطی نیز قابل‌قبول باشد. در این کار به‌طور مقایسه‌ای به بررسی مزایا و معایب استفاده از مواد شیمیایی مختلف برای مهار پالایش هیدروژن سولفید در سیالات حفاری و تولیدی پرداخته‌شده است (9). ابراهیم‌زاده (1394) به آنالیز حوادث فرآیندی ناشی از نشتی گاز هیدروژن سولفید از واحد لخته گیری یکی از پالایشگاه‌های پارس جنوبی با رویکرد تعیین حریم ایمن پرداخت. مطالعه صورت گرفته به‌صورت توصیفی تحلیلی و به‌طور مقطعی در پالایشگاه‌های گازی 2 و 3 پارس جنوبی صورت گرفت جهت ارزیابی چگونگی پخش و تأثیرگذاری سناریوی نشتی از نرم‌افزار PHAST نسخه 6.54 استفاده شد (10). نوریان و امیری (1395) به بررسی ارزیابی ریسک ایمنی ایستگاه تقلیل فشارگازطبیعی با نرم‌افزار PHAST پرداختند. حوادثی از قبیل نشت گاز و آتش‌سوزی به دلیل جدایی اتصالات شیرها و دیگر تجهیزات، خوردگی، سایش و همچنین آسیب دیدن و یا از تنظیم خارج شدن تجهیزات به دلیل افزایش فشار در کلیه تأسیسات انتقال و شبکه توزیع گاز ازجمله ایستگاه‌های تقلیل فشار امکان وقوع دارند. اهمیت به ارزیابی دقیق و پیش‌بینی احتمال انفجار، حریق و تعیین شعاع ایمن برای ایستگاه‌های تقلیل فشار گاز طبیعی که دارای استانداردهای یکسانی ازنظر حجم گاز ورودی و خروجی با توجه به ظرفیت طراحی‌شده آن‌ها می‌باشند، ضروری است. لذا سه سناریو برای انفجار احتمالی تعیین گردیدند. در این سناریوها اندازه نشتی گاز قبل و بعد از رگلاتور و همچنین پارگی مخزن برابر قطر تجهیزات به ترتیب 12، 6 و 20 اینچ تعیین شدند. پس از مدل‌سازی ریاضی در نرم‌افزار PHAST فاصله اطمینان 250 متر برای بدترین نوع نشتی تعیین شد و دیگر فاصله‌های ایمن ایستگاه برابر 105 و 70 متر می‌باشند. در صورت نشت مرکاپتان با غلظت 44000 پی.پی.ام، فاصله‌های 170 و 75 متر به‌عنوان نقاط ایمن تعیین گردیدند (11). مینایی (1396) به شبیه‌سازی غلظت و نحوه پراکنش گاز سولفید هیدروژن (H2S) ناشی از کوره‌های زباله‌سوز واحد بازیافت گوگرد در یک پالایشگاه گازی در عسلویه پرداخت. در این مطالعه غلظت آلاینده‌های منتشره از دو دودکش کوره‌های زباله‌سوز و همچنین غلظت آلاینده H2S در پنج ایستگاه پایش محیطی در محدوده پالایشگاه و اطراف آن طی چهارفصل از تابستان 1393 تا پایان بهار 1394 اندازه‌گیری شد (12). موحد و پروینی (1398) به مدل‌سازی پیامد نشت میعانات گازی در مخزن ۶۵۰ یک پالایشگاه گازی که خطرناک‌ترین مخزن واحد است، به‌منظور تدوین برنامه واکنش در شرایط اضطراری پرداخت. در این مطالعه برای بررسی نحوه انتشار میعانات گازی از مخزن ذخیره‌سازی پالایشگاه گازی موردمطالعه و پیامد آن از نرم‌افزار PHAST   استفاده و  مدل‌سازی برای چهار سناریوی آتش فورانی، آتش استخری، آتش ناگهانی و انفجار ابر بخار انجام شد. چهار سناریو یادشده می‌توانند پیامدهای مرگ باری را در پی داشته باشند و حادثه انفجار ابر بخار و پس‌ازآن آتش استخری دارای بیشترین پیامد هستند. بر اساس گلباد منطقه و مدل‌سازی انجام‌شده برخی تغییرات در برنامه واکنش در شرایط اضطراری، ازجمله تغییر مسیرهای دسترسی و همچنین تغییر فاصله دو (Point Muster) موردنیاز است (13). نعمتی (1398) بررسی آسیب های انتشار گاز از مخزن پروپان در یک پالایشگاه گازی به کمک نرم‌افزارPHAST مطالعه ی موردی: پالایشگاه هفتم پارس جنوبی، در این مطالعه به‌منظور تعیین حداکثر ریسک محتمل در شرایط ذخیره پروپان در مخزنی با حجم 45000 مترمکعب، دمای 46- درجه سانتی‌گراد و فشار 01/1 بار دو پدیده گسست مخزن و تخلیه کل موجودی مخزن در یک‌زمان معین بررسی شد. در اثر گسست مخزن پدیده‌های انتشار ابر ناشی از غلظت پروپان، آتش استخری تأخیر یافته، آتش ناگهانی و تبخیر حاصل از استخر ایجادشده و در اثر تخلیه کل موجودی مخزن در یک‌زمان معین پدیده‌های انتشار ابر ناشی از غلظت پروپان، آتش استخری تأخیر یافته و سریع، انفجار، آتش ناگهانی و تبخیر حاصل از استخر ایجادشده رخ خواهد داد (14). تیان و همکاران (2021)، به بررسی "تشخیص نشت در شبکه‌های توزیع گاز کم‌فشار با روش‌های احتمالی" پرداختند. رویدادهای نشتی با استفاده از داده‌های فشار و جریان به‌دست‌آمده از مدل‌سازی حالت پایا شبکه گاز برآورد می‌شوند. استحکام روش با تجزیه‌وتحلیل شبکه‌های گاز در حضور خطاهای اندازه‌گیری نشان داده‌شده است که نویز حسگر اجتناب‌ناپذیر در داده‌های جریان و فشار را به حساب می‌آورد. امکان‌سنجی روش پیشنهادی در یک شبکه گاز مصنوعی کوچک نشان داده‌شده است. علاوه بر این، این روش در بخشی از شبکه توزیع گاز سنگاپور برای سناریوهای نشت منفرد و همچنین چندگانه اعمال می‌شود. همچنین به‌طور تجربی نشان داده‌شده است که شدت نشت و مکان برای یک سناریوی نشت واحد را می‌توان به ترتیب با دقت 95% و 80% حتی در حضور نویز قوی تعیین کرد (15). ژو و لیو ( 2013)، به بررسی "چارچوب بهینه‌سازی مبتنی بر شبکه عصبی عمیق برای مسیر تخلیه ایمنی در هنگام حوادث نشت گاز سمی" پرداختند. برنامه‌ریزی تخلیه برای کاهش تلفات در حوادث نشت گاز سمی مهم است. بااین‌حال، طرح‌های تخلیه که در موقعیت‌های عملی غیرمنتظره اعمال می‌شوند، بیشتر کیفی هستند. این مطالعه سیستم پیشنهادی مسیر تخلیه را بر اساس ارزیابی کمی ریسک ارائه می‌کند که با پیش‌بینی پراکندگی گاز دینامیکی بادقت بالا و زمان محاسبه کوتاه، ایمن‌ترین مسیر را برای افراد تخلیه‌شده فراهم می‌کند. سناریوهای تخلیه دقیق، ازجمله شرایط آب‌وهوایی، شدت نشت، و اطلاعات تخلیه، در نظر گرفته شد. سیستم پیشنهادی ریسک کمی را در منطقه آسیب‌دیده با استفاده از یک مدل جایگزین شبکه عصبی عمیق برای تعیین مسیرهای تخلیه بهینه با برنامه‌ریزی عددی و ریاضی صحیح ارزیابی می‌کند. مدل جایگزین با استفاده از داده‌های شبیه‌سازی دینامیک سیالات محاسباتی آموزش داده شد. یک رمزگذار خودکار متغیر برای استخراج ویژگی‌های هندسی ناحیه آسیب‌دیده استفاده شد. ریسک پیش‌بینی‌شده با برنامه‌ریزی اعداد صحیح خطی برای تعیین مسیر بهینه در یک شبکه جاده‌ای از پیش تعریف‌شده ترکیب شد. برای مطالعه موردی از سناریوی نشت خط لوله گاز آمونیاک در مجتمع پتروشیمی استفاده شد. نتایج نشان می‌دهد که مدل توسعه‌یافته ایمن‌ترین مسیر را در عرض چند ثانیه با حداقل ریسک ارائه می‌دهد. مدل توسعه‌یافته برای تحلیل حساسیت جهت تعیین تأثیرات متغیر و مکان‌های پناهگاه امن اعمال شد (16). سانگ و همکاران (2021)، پیرامون "ریسک پیامدهای نشت گاز و انفجار در داخل ساختمان‌های مسکونی" ارزیابی کمی انجام دادند. در این مطالعه، یک روش ارزیابی کمی ریسک برای پیامدهای نشت گاز و حوادث انفجار در داخل ساختمان‌های مسکونی توسعه داده‌شده است. در همین حال، اقدامات مؤثر کاهش ریسک نیز پیشنهادشده است. رفتار نشت و پراکندگی گاز ابتدا برای ارائه یک سناریوی اولیه واقعی برای شبیه‌سازی انفجار بعدی با استفاده از کد CFD FLACS تجزیه‌وتحلیل می‌شود و سپس، عواقب انفجار با یک ابر گاز قابل اشتعال ناهمگن واقعی به‌دقت پیش‌بینی می‌شود. بر اساس مدل پروبیت، احتمال اثرات نامطلوب مختلف (به‌عنوان‌مثال، پارگی پرده گوش، آسیب ریه و ضربه سر) ناشی از حوادث انفجار گاز به ترتیب برآورد شده است. با ترکیب احتمال با نمرات شدت پیامد، یک مفهوم مبتنی بر ریسک برای محاسبه شاخص خطر انفجار در هر منطقه از ساختمان‌های مسکونی اتخاذ شده است. در همین حال، شاخص ریسک بالاتر از بین این اثرات نامطلوب در هر نقطه شبکه انتخاب می‌شود. علاوه بر این، یک مفهوم مبتنی بر شبکه و ایزولاین بیشتر برای نشان‌دادن این شاخص خطر انفجار استفاده می‌شود. علاوه بر این، اثرات نسبت سطح دریچه و فشار شکستگی شیشه‌ها برای پیشنهاد اقدامات مؤثر کاهش خطر بررسی می‌شود. این روش برای یک ساختمان مسکونی معمولی، واقع در هاربین، چین اعمال می‌شود. نتایج نشان می‌دهد که شدت پیامد انفجار در هر منطقه از ساختمان‌های مسکونی می‌تواند به‌صورت کمی ارائه شود. علاوه بر این، توزیع منطقه خطر انفجار گاز در داخل ساختمان‌های مسکونی نیز قابل شناسایی است. علاوه بر این، حداقل نسبت سطح دریچه و حداکثر فشار شکستگی شیشه‌ها را می‌توان برای کاهش شاخص خطر انفجار گاز توصیه کرد. این مطالعه یک روش تحلیل مؤثر ارائه می‌کند و به توسعه اقدامات ایمنی و بهبود طرح‌های واکنش اضطراری برای حوادث نشت گاز و انفجار در داخل ساختمان‌های مسکونی کمک می‌کند (17).
با توجه به مطالب بیان‌شده معین گردید رویکرد استفاده از نرم‌افزارهای باقابلیت مدل‌سازی انتشار گازهای سیمی و قابل اشتعال در تعیین حریم‌های ایمن و بررسی سناریوهای متفاوت ابزار مناسبی برای بررسی میزان ریسک و کاهش آن بوده است. لذا در این مطالعه، به‌منظور ارزیابی پیامد و تعیین حریم ایمن در ایستگاه تقلیل فشار گاز طبیعی برومی شهر اهواز به تحلیل تخلیه گاز، انتشار گاز، آتش فورانی، آتش ناگهانی و موج انفجار گاز و درنتیجه مشخص‌شدن پیامد ناشی از آن‌ها از طریق مدل‌سازی و آنالیز مخاطرات به روش ارزیابی کمی ریسک به‌وسیله نرم‌افزار PHAST پرداخته می‌شود. ازآنجایی‌که شرایط عملیاتی و موردمطالعه در تحقیقات پیشین بر روی واحد ایستگاه تقلیل فشار به کمک نرم‌افزار PHAST صورت نپذیرفته، تمرکز مطالب حاضر استفاده از تکنیک‌های ریاضی-سیالاتی به کمک این نرم‌افزار برای پیاده‌سازی مدل‌های انتشار بوده و توجه مقاله بر روی کاربرد عملی تکنیکهای علمی در دنیای واقعی استوار گردیده است.
روش بررسی
1-روش ارزیابی کمی ریسک: ارزیابی کمی ریسک مطابق شکل 1 دارای هفت مرحله می­باشد که در زیر به آن‌ها پرداخته می­شود:
مرحله اول : تعیین اهداف ارزیابی کمی ریسک
همان‌طور که از عنوان تحقیق مشخص است هدف ارزیابی پیامد ناشی از حوادث فرآیندی ازجمله آتش و انفجار گاز که در اثر نشت گاز از تجهیزات فرآیندی به وقوع می­پیوندند، به‌وسیله مدل‌سازی در نرم‌افزار PHAST و درنتیجه تعیین حریم ایمن در ایستگاه تقلیل فشار گاز منطقه برومی شهر اهواز می­باشد.
مرحله دوم : شرح واحد فرآیندی تحت بررسی
در این مرحله اطلاعات مربوط به ایستگاه تقلیل فشار گاز موردمطالعه مانند داده‌های فرآیندی (دما و فشار جریان گاز، فاز و ترکیب درصد گاز، ارتفاع از سطح زمین) از واحد ایمنی و داده­های هواشناسی منطقه (دمای محیط، سرعت باد، رطوبت هوا) از بخش آرشیو وب-سایت رسمی هواشناسی استان خوزستان جمع‌آوری گردیده است.
مرحله سوم : شناسایی مخاطرات
در این مرحله به کمک تحقیقات مشابه انجام‌شده به شناسایی مخاطرات بالقوه در ایستگاه تقلیل فشار گاز برون‌شهری برومی شهر اهواز نظیر خوردگی، برخورد فیزیکی، سهل‌انگاری کارکنان در حین تعمیرات و افزایش ناگهانی فشار جریان گاز پرداخته می­شود.
مرحله چهارم : تعیین و تحلیل سناریوها
سناریو، حادثه‌ای است که باعث ایجاد ریسک‌های فرآیندی از قبیل آتش، انفجار و آزادسازی مواد سمی می‌شود. سناریو مانند یک اتفاق جداگانه است و می‌تواند در یک بخش فرآیندی ایجاد شود و هر سناریو دارای چندین نتیجه از قبیل آتش، انفجار و انتشار مواد سمی می­باشد. سناریوها بیشتر مواقع به‌صورت نشت و یا پارگی در سیستم‌های فرآیندی به‌حساب می‌آیند به‌صورتی که باعث انتشار ماده سمی و یا ماده قابل اشتعال شوند (18).
مشخصات سناریو
برای ارزیابی یک سناریو داشتن اطلاعات زیر لازم است :
  1. مکان ایجاد سناریو (مختصات جغرافیایی و تجهیز مربوطه)
  2.  شرایط فرآیندی (نوع ماده، فاز، دما، فشار و غیره)
  3.  کل مقدار مواد تخلیه‌شده در اثر وقوع سناریو
  4.  نوع سناریو (نشتی، تخلیه ناگهانی و یا غیره)
  5. اندازه سوراخ به وجود آمده در مورد نشتی
  6.  ارتفاع انتشار مواد از سطح مبدأ
  7. جهت انتشارات (افقی، عمودی و مورب)
  8. زبری سطوح (18).
در این مرحله به کمک تحقیقات مشابه انجام‌شده به
تعیین و تحلیل سناریوها که اغلب به‌صورت نشت گاز از
تجهیزات فرآیندی نظیر فلنج ها و شیرها می­باشند، پرداخته می­شود.

مرحله پنجم : مدل‌سازی پیامد
در این مرحله سناریوهای تعیین‌شده بر اساس داده­های فرآیندی و اطلاعات هواشناسی منطقه به کمک نرم‌افزار PHAST و با استفاده از مدل‌های شناوری سیالات و حل معادلات مؤمنتم و انتقال جرم هم‌زمان مدل می­شوند.
مرحله ششم : تخمین تکرارپذیری سناریوها
پیامد * تکرارپذیری = ریسک
در این مرحله سناریوهایی که احتمال رخ دادن آن‌ها بیشتر و همچنین دارای شدت زیاد باشند موردبررسی قرار می­گیرند.
مرحله هفتم : محاسبه و ارزیابی ریسک
در این مرحله نرم‌افزار PHAST به‌صورت پیش‌فرض پیامدهای ناشی از سناریوهای تعیین‌شده را برحسب میزان تشعشع برای آتش و شدت موج برای انفجار گاز محاسبه می­کند و درنتیجه با مقایسه نتایج خروجی نرم‌افزار PHAST با جدول سطوح آسیب ناشی از تشعشع برای آتش و جدول خسارت ناشی از انفجار می­توان به محدوده خطر و درنتیجه به حریم ایمن ایستگاه تقلیل فشار گاز برون‌شهری برومی شهر اهواز دست یافت.
 


شکل 1: دیاگرام روش تحقیق (18). 

شکل 2: ایستگاه تقلیل فشار گاز برون‌شهری (19)
 
  1. ایستگاه تقلیل فشار گاز برون‌شهری مورد مطالعه
جریان گاز طبیعی با دبی حجمی m3/hr 120.000  و فشاری معادل psi 400 وارد ایستگاه تقلیل فشار گاز شده و نخست جهت فیلتر شدن (زدودن مواد معلق و غیر معلق ازجمله روغن کمپرسور، گردوغبار، آب باقیمانده) از صافی گاز رد شده و سپس در دو مرحله کاهش فشار به‌ واسطه دو عدد رگلاتور (فشارشکن) که در رگلاتور اولی فشار از psi 400 به psi 320 کاهش پیدا کرده و نهایتاً در مرحله دوم فشار گاز از psi 320 به psi 250 به‌واسطه رگلاتور (فشارشکن) دومی کاهش پیدا می‌کند. این ایستگاه چهار عدد خط دارد که یکی از آن‌ها به‌عنوان (bypass) خط کنارگذر می­باشد که در واقع این خط کنارگذر جهت اطمینان از تأمین فشار psi250 خروجی از ایستگاه در مواقع کاهش فشار ورودی ایستگاه به‌ویژه هنگام تعمیرات یکی از خطوط به کار گرفته می‌شود. در این ایستگاه از شیرهای قطع فشار (shutoff valve) جهت قطع جریان گاز هنگام بالارفتن فشار گاز از فشار تنظیمی (فشار خروجی رگلاتور) به میزان 15% استفاده شده است. همچنین در این ایستگاه جهت مقابله با افزایش ناگهانی فشار از شیرهای اطمینان (safety valve) استفاده شده است که در واقع در هنگام بروز اختلال در عملکرد رگلاتورها و افزایش فشار گاز به میزان 10% بالاتر از فشار تنظیمی جریان (فشار خروجی از رگلاتور) عمل می‌کنند (به‌صورت تخلیه بخشی از گاز به محیط اطراف). در شکل 2 نمایی کلی از ایستگاه تقلیل فشار گاز برون‌شهری را نمایش می­دهد. گاز تصفیه شده خروجی فیلتر به سمت هیتر هدایت می‌گردد، هیتر ها جهت گرم نمودن گاز ورودی عمل می‌نمایند، چرا که طبق قوانین حاکم بر گازها، در اثر کاهش فشار گاز افت دما حاصل‌شده و موجب یخ بستن و تشکیل هیدراتهای یخی می‌گردد. لذا با پیش گرم کردن گاز و افزایش دما قبل از کاهش فشار، از یخ بستن در قسمت پایین‌دستی جریان جلوگیری می­گردد. در قسمت خروجی، گاز را به‌منظور تشخیص و کشف نشت گاز قبل از ایجاد شرایط انفجار و همچنین تأمین ایمنی مصرف‌کنندگان به‌واسطه دستگاه بودار کننده (Gas odorizer) به روش تزریقی، بودار می‌کنند. مقدار ماده بودار کننده متناسب با حجم گاز و به میزانی که یک شخص با بویایی متوسط بتواند موقعی که مقدار گاز در هوا به یک‌پنجم حداقل قابلیت اشتعال (LFL) (Lower Flammable Limit) برسد به‌راحتی آن را تشخیص بدهد، تزریق می‌شود. همین‌طور در این ایستگاه جهت سنجش مقدار جریان گاز از کنتور استفاده می‌شود.
ارزیابی پیامدها: سناریوهای تعیین‌شده مذکور در جدول 1 به‌منظور تعیین میزان خطرات احتمالی، محدوده خطر و درنتیجه تعیین حریم در ایستگاه تقلیل فشار گاز برون‌شهری موردمطالعه در نرم‌افزار PHAST ارزیابی می‌گردند. پیامد ناشی از انفجار به‌صورت میزان توزیع شدت موج انفجار و پیامد ناشی از آتش به‌صورت میزان توزیع شدت تشعشع می­باشند که در نرم‌افزار PHAST بر مبنای سه شدت موج تعریف‌شده برای انفجار و سه سطح تشعشع تعریف‌شده برای آتش مورد ارزیابی قرار می‌گیرند. اصولاً انتخاب سناریوها با توجه به تعریف ریسک که تابعی از احتمال رخداد و شدت پیامد حاصله است بر مبنای این دو تابع خواهد بود. لذا در این مرحله سناریوهایی بررسی می­شوند که احتمال رخ دادن آن‌ها بالا و همچنین پیامدهای آن‌ها دارای شدت و تأثیر زیاد باشند. محاسبه و شناسایی نقاط با احتمال و شدت بالا بر اساس اطلاعات و اسناد موجود در شناسایی مخاطرات هر واحد مختص آن واحد استخراج می‌گردد. با توجه به توضیحات داده‌شده، سناریوهای محتمل برای ایستگاه تقلیل فشار گاز برون‌شهری واقع در منطقه برومی شهر اهواز تجربیات مهندسین فرآیند واحد و با توجه تحقیقات گذشته در این خصوص، به‌صورت نشت گاز از رگلاتور مرحله اول، رگلاتور مرحله دوم و فلنج شیر قطع کننده اضطراری و بر اساس داده‌های فرآیندی مذکور در جدول 1 تعیین شدند.
سناریوهای تعیین‌شده در جدول 1 بر اساس داده‌های فرآیندی نظیر فشار، دما، ترکیب درصد گاز طبیعی (جدول 2)، دبی جریان گاز، ارتفاع از سطح زمین و همچنین مطابق معیار پایداری جوی پاسکوییل (جدول 3) و داده­های هواشناسی منطقه  (جدول 4) در نرم‌افزار PHAST نسخه 8.22 مدل می‌شوند. دبی حجمی جریان گاز در ایستگاه موردمطالعه برابر 120 هزار مترمکعب در ساعت می‌باشد.  خروجی مدل‌ها برای سناریوهای تعیین‌شده در جدول 1 به‌صورت گزارش خلاصه پیامد و بر اساس سطوح تشعشع آتش و شدت‌های موج انفجار تعریف‌شده در نرم‌افزار PHAST می‌باشند. در جدول 3، در تقسیم‌بندی پاسکوییل، کلاس جوی A و B نشان‌دهنده جو ناپایدار، C و D نشان‌دهنده جو خنثی و E و F  نشان‌دهنده جو پایدار می‌باشند. کلاس جوی A/B  بر اساس تقسیم‌بندی پاسکوییل بیانگر جو ناپایدار که اغلب در روز با درنظرگرفتن میزان تابش آفتاب در منطقه و کلاس جوی F بیانگر جو پایدار که اغلب در شب با درنظرگرفتن خفیف بودن باد اختیار شدند. همچنین بر اساس دستورالعمل مدل‌سازی حوادث بدترین شرایط همواره جزء سناریوهای موردبررسی خواهد بود و پیش فرض نرم‌افزار نیز بر روی اندازه نشتی در این حالت تنظیم‌شده است.
 
جدول 1: سناریوهای تعیین‌شده
شماره سناریو سناریو محل سناریو اندازه نشتی
( in )
فشار (psi) دمای میانگین (°C) ارتفاع از سطح زمین (m) فصل
1 leak رگلاتور مرحله اول 4 400 13 1 زمستان
2 leak رگلاتور مرحله اول 4 400 30 1 تابستان
3 leak
رگلاتور مرحله دوم 4 320 12 1 زمستان
4 leak رگلاتور مرحله دوم 4 320 29 1 تابستان
5 leak فلنج شیر قطع‌کننده اضطراری 6 400 13 1 زمستان
6 leak فلنج شیر قطع‌کننده اضطراری 6 400 30 1 تابستان

جدول 2: آنالیز گاز طبیعی
ترکیب درصد (mol %)
دی‌اکسیدکربن متان اتان پروپان ایزو بوتان نرمال بوتان ایزو پنتان نرمال پنتان
23/0 56/94 36/4 63/0 06/0 13/0 02/0 01/0

جدول 3: داده‌های ترمودینامیکی گوگرد گازی
معیار پایداری پاسکوییل
(m/s)سرعت باد روز : میزان تابش خورشید شب : درجه ابری بودن آسمان هوای کاملاً ابری در هر زمان از شبانه‌روز
کم متوسط زیاد کمتر از 4/0 بین 4/0 و 8/0
کمتر از 2 B A-B A - - D
2-3 C B A-B E F D
3-5 C B-C B D E D
5-6 D C-D C D D D
بیشتر از 6 D D C D D D

جدول 4: داده‌های هواشناسی به‌کاررفته در نرم‌افزار PHAST  بر اساس داده‌های سازمان هواشناسی برای اهواز و مطابق معیار پایداری جوی پاسکوییل
فصل دمای میانگین (°C) میانگین نسبت رطوبت % سرعت باد میانگین (m/s ) معیار پایداری جوی پاسکوییل
شب روز شب روز
زمستان 9 19 41 2 F A/B
تابستان 29 46 10 2 F A/B
 
- ارزیابی آسیب‌ها
نتایج خروجی نرم‌افزار PHAST برای سناریوهای تعیین‌شده در جدول 1 بر مبنای معیار سنجش آسیب‌های ایجاده شده در اثر تشعشع آتش و شدت موج انفجار مذکور بر اساس موارد موجود در مراجع مورد ارزیابی قرار می­گیرند (20).
به‌عنوان نمونه در خصوص آسیب ایجادشده در اثر تشعشع آتش موارد شدت تشعشع با میزان  kW/m2 5/4 شروع درد و سوختگی تا مقدار  kW/m2 5/37 سبب خسارت به واحدها و تجهیزات فرآیندی، ایجاد مرگ آنی برای افراد در معرض آن می‌گردد. همچنین برای معیار آسیب ایجادشده در اثر موج انفجار حد آستانه عمومی برای شکستن پنجره‌ها در مقدار bar 01/0 تا احتمال تخریب کلی ساختمان‌ها در مقدار bar 7/0 وجود دارد. موارد بیشتر در مقادیر آستانه آسیب و تخریب در مرجع (20) در دسترس است.
نتایج و بحث
نتایج ارزیابی پیامد نیز بر اساس حالت‌های مختلف تخلیه گاز، انتشار گاز، آتش فورانی، آتش ناگهانی، انفجار موردبررسی قرار گرفت و در جداول 5 الی 12 گزارش‌شده است. نتایج خروجی مربوط به هر یک از سناریوها به‌صورت شش جدول مرتبط با تخلیه گاز، انتشار گاز، آتش فورانی، آتش ناگهانی و انفجار گاز می­باشند. تفاوت داده­های
خروجی در هر جدول یک سناریو مربوط به داده‌های هواشناسی و معیار پایداری پاسکوییل ویژه اوقات روز و شب
می­باشند.

نتایج خروجی انتشار گاز مربوط به هر یک از سناریوها بر اساس حد پایین اشتعال گاز، نتایج خروجی آتش فورانی در هر یک از سناریوها بر مبنای سه سطح تشعشع پیش فرض نرم‌افزار PHAST، نتایج خروجی آتش ناگهانی برای هر یک از سناریوها بر مبنای حد پایین اشتعال گاز و نتایج
خروجی انفجار گاز مربوط به هر یک از سناریوها بر مبنای سه شدت موج انفجار پیش فرض نرم‌افزار
PHAST
می­باشند.   
 

جدول 5: گزارش خلاصه پیامد سناریوها
شماره سناریو نوع سناریو اطلاعات هواشناسی* دبی جریان (kg/s) دما(C) جزء جرمی مایع در جریان قطر حباب μm قطر انبساط سرعت m/s زمان نهایی رهایش s
سناریو 1 نشتی Day 2/A/B 7/36 -6/109 0 0 23/0 2/670 1/19
Night 2/F 7/36 -6/109 0 0 23/0 2/670 1/19
سناریو 2 نشتی Day 2/A/B 4/35 -1/98 0 0 23/0 9/694 4/18
Night 2/F 4/35 -1/98 0 0 23/0 9/694 4/18
سناریو 3 نشتی Day 2/A/B 5/29 -1/107 0 0 21/0 6/668 19
Night 2/F 5/29 -1/107 0 0 21/0 6/668 19
سناریو 4 نشتی Day 2/A/B 4/28 -6/95 0 0 21/0 4/692 4/18
Night 2/F 4/28 -6/95 0 0 21/0 4/692 4/18
سناریو 5 نشتی Day 2/A/B 7/82 -6/109 0 0 35/0 2/670 4/8
Night 2/F 7/82 -6/109 0 0 35/0 2/670 4/8
سناریو 6 نشتی Day 2/A/B 8/79 -1/98 0 0 35/0 9/694 2/8
Night 2/F 8/79 -1/98 0 0 35/0 9/694 2/8
*اطلاعات هوا شناسایی به ترتیب شامل روز/شب، سرعت باد(m/s) و کلاس پایداری پاسکویل است.    
 
گزارش خلاصه پیامد سناریوهای تعیین‌شده: سناریوهای بررسی‌شده بر مبنای جدول 1 با موضوع زمان پایان تخلیه گاز، سرعت تخلیه گاز، قطر انبساط گاز هنگام تخلیه، دمای گاز هنگام تخلیه و اوج دبی جرمی جریان گاز هنگام تخلیه در جدول 5 آورده شده است که اوج دبی جرمی جریان گاز هنگام تخلیه برای هر دو کلاس جوی A/B  و F در تمامی سناریوها با هم برابر می‌باشند.
سناریوهای بررسی‌شده بر مبنای جدول 1 برای نتایج انتشار بر اساس فواصل محدوده در جدول 6 آورده شده است که فاصله در جهت باد تا حد پایین اشتعال گاز طبیعی، و کسر حد پایین  برای هر دو کلاس جوی A/B  و F در تمامی سناریوها قابل‌مشاهده است. بنا بر نتایج به‌دست‌آمده فاصله تا محدوده UFL در مدل‌سازی انجام‌شده حاصل نگردیده است.
 
جدول 6. فاصله در جهت باد تا غلظت‌های تعریف‌شده
شماره سناریو نوع سناریو اطلاعات هواشناسی فاصله تا محدوده UFL (m) فاصله تا محدوده LFL (m) فاصله تا کسر LFL (m)
سناریو 1 نشتی Day 2/A/B NA 2/59 6/154
Night 2/F NA 5/61 3/126
سناریو 2 نشتی Day 2/A/B NA 60 8/155
Night 2/F NA 3/61 7/125
سناریو 3 نشتی Day 2/A/B NA 6/50 6/134
Night 2/F NA 6/53 4/113
سناریو 4 نشتی Day 2/A/B NA 5/51 136
Night 2/F NA 6/53 1/113
سناریو 5 نشتی Day 2/A/B NA 1/102 197
Night 2/F NA 4/87 8/143
سناریو 6 نشتی Day 2/A/B NA 1/103 8/197
Night 2/F NA 2/87 3/143
 

نتایج بر اساس فاصله در جهت باد تا سطوح 1، 2 و 3 تشعشع آتش فورانی و همچنین طول شعله برای هر دو کلاس جوی A/B و F برای سناریوهای بررسی‌شده بر مبنای جدول 1 در جدول 7 قابل‌دسترس است. در خصوص آتش ناگهانی نتایج بر اساس فاصله در جهت باد تا حد پایین اشتعال گاز طبیعی برای هر دو کلاس جوی A/B و F در جدول 8 ارائه‌شده است بنا بر نتایج سناریوهای تعیین‌شده فواصل به‌دست‌آمده برای آتش ناگهانی با نتایج انتشار گاز یکسان می‌باشند.
 
جدول 7: فاصله در جهت باد تا سطوح تشعشع تعریف‌شده
شماره سناریو نوع سناریو اطلاعات هواشناسی طول شعله (m) فاصله تا شدت تشعشع سطح1 * در مسیر باد (m) فاصله تا شدت تشعشع سطح2* در مسیر باد (m) فاصله تا شدت تشعشع سطح3* در مسیر باد (m)
سناریو 1 نشتی Day 2/A/B 55 3/102 1/78 6/61
Night 2/F 1/54 1/103 3/78 9/61
سناریو 2 نشتی Day 2/A/B 9/55 9/100 6/77 1/61
Night 2/F 5/54 8/101 7/77 5/61
سناریو 3 نشتی Day 2/A/B 50 9/91 4/70 8/55
Night 2/F 2/49 5/92 6/70 56
سناریو 4 نشتی Day 2/A/B 8/50 7/90 1/70 4/55
Night 2/F 6/49 4/91 1/70 8/55
سناریو 5 نشتی Day 2/A/B 4/77 2/151 3/113 2/87
Night 2/F 2/76 152 3/113 5/87
سناریو 6 نشتی Day 2/A/B 8/78 9/148 4/112 5/86
Night 2/F 77 1/151 1/113 7/87
*شدت تشعشع سطح 1: kW/m2 4، سطح2: kW/m2 5/12 و سطح 3: kW/m2 5/37
جدول 8: فاصله در جهت باد تا غلظت‌های تعریف‌شده
شماره سناریو نوع سناریو اطلاعات هواشناسی* فاصله تا محدوده LFL در مسیر باد (m) فاصله تا کسر در مسیر باد LFL در مسیر باد (m)
سناریو 1 نشتی Day 2/A/B 2/59 6/154
Night 2/F 5/61 3/126
سناریو 2 نشتی Day 2/A/B 60 8/155
Night 2/F 3/61 7/125
سناریو 3 نشتی Day 2/A/B 6/50 6/134
Night 2/F 6/53 4/113
سناریو 4 نشتی Day 2/A/B 5/51 136
Night 2/F 6/53 1/113
سناریو 5 نشتی Day 2/A/B 1/102 197
Night 2/F 4/87 8/143
سناریو 6 نشتی Day 2/A/B 1/103 8/197
Night 2/F 2/87 3/143
 
نتایج به‌صورت قطر انفجار و بیشینه فاصله از محل سناریو در جهت باد تا شدت‌های موج انفجار تعریف‌شده برای هر دو کلاس جوی A/B و F برای 6 سناریو تعیین‌شده در جدول 9 ارائه‌شده است.
نتایج برای داده‌های تکمیلی برای بدترین سناریوهای انفجار در جدول 10 آورده شده است. نتایج به‌صورت مرکز انفجار، مرکز ابر، منبع احتراق (از محل سناریو و در جهت باد) و همچنین زمان احتراق و جرم قابل اشتعال انفجار برای شدت‌های موج تعریف‌شده و هر دو کلاس جوی A/B و F به دست آمد.
 
جدول 9: سناریوهای انفجار در بدترین حالت تا شدت‌های موج انفجار تعریف‌شده
شرایط جوی سطح  Overpressure [bar] بیشینه فاصله [m] قطر [m] نوع سناریو شماره
Day 2/A/B 02068/0
1379/0
2068/0
9/271
2/170
1/165
8/263
4/40
3/30
نشتی سناریو 1
Night 2/F 02068/0
1379/0
2068/0
6/265
3/148
2/141
2/291
6/56
4/42
Day 2/A/B 02068/0
1379/0
2068/0
270
3/173
4/167
1/240
6/46
9/34
نشتی سناریو 2
Night 2/F 02068/0
1379/0
2068/0
3/262
6/147
7/140
7/284
3/55
4/41
Day 2/A/B 02068/0
1379/0
2068/0
8/246
7/152
147
6/233
4/45
34
نشتی سناریو 3
Night 2/F 02068/0
1379/0
2068/0
6/238
135
7/128
2/257
50
4/37
Day 2/A/B 02068/0
1379/0
2068/0
6/244
2/152
7/146
2/229
5/44
4/33
نشتی
سناریو 4
Night 2/F 02068/0
1379/0
2068/0
9/235
4/134
3/128
8/251
9/48
7/36
Day 2/A/B 02068/0
1379/0
2068/0
7/346
5/190
4/190
5/413
1/1
8/0
نشتی
سناریو 5
Night 2/F 02068/0
1379/0
2068/0
5/358
4/174
1/164
1/457
8/88
3/48
Day 2/A/B 02068/0
1379/0
2068/0
342
191
7/190
1/404
1/2
5/1
نشتی
سناریو 6
Night 2/F 02068/0
1379/0
2068/0
5/352
2/173
163
1/445
5/86
46


جدول 10: داده‌های تکمیلی برای بدترین سناریوهای انفجار
نوع سناریو شرایط جوی سطح  Overpressure [bar] جرم قابل اشتعال/انفجار  [kg] زمان رسیدن به جرقه [s] منبع جرقه [m] مرکز ابر گازی[m] مرکز انفجار [m] شماره
نشتی Day 2/A/B 02068/0
1379/0
2068/0
4/25
4/12
4/12
4/15
8/19
8/19
140
150
150
31
7/43
7/43
140
150
150
سناریو 1
Night 2/F 02068/0
1379/0
2068/0
2/34
2/34
2/34
5/17
5/17
5/17
120
120
120
3/33
3/33
3/33
120
120
120
نشتی Day 2/A/B 02068/0
1379/0
2068/0
2/19
2/19
2/19
7/18
7/18
7/18
150
150
150
6/36
6/36
6/36
150
150
150
سناریو 2
Night 2/F 02068/0
1379/0
2068/0
32
32
32
2/17
2/17
2/17
120
120
120
4/33
4/33
4/33
120
120
120
نشتی Day 2/A/B 02068/0
1379/0
2068/0
7/17
7/17
7/17
4/16
4/16
4/16
130
130
130
2/27
2/27
2/27
130
130
130
سناریو 3
Night 2/F 02068/0
1379/0
2068/0
6/23
6/23
6/23
4/18
4/18
4/18
110
110
110
2/29
2/29
2/29
110
110
110
نشتی Day 2/A/B 02068/0
1379/0
2068/0
7/16
7/16
7/16
5/15
5/15
5/15
130
130
130
6/27
6/27
6/27
130
130
130
سناریو 4
Night 2/F 02068/0
1379/0
2068/0
1/22
1/22
1/22
18
18
18
110
110
110
4/29
4/29
4/29
110
110
110
نشتی Day 2/A/B 02068/0
1379/0
2068/0
98
0002/0
0002/0
1/8
9/16
9/16
140
190
190
8/50
86/0
86/0
140
190
190
سناریو 5
Night 2/F 02068/0
1379/0
2068/0
4/132
4/132
5/50
4/8
4/8
3/10
130
130
130
54
54
3/79
130
130
130
نشتی Day 2/A/B 02068/0
1379/0
2068/0
5/91
001/0
001/0
7/7
1/16
1/16
140
190
190
1/51
4/5
4/5
140
190
190
سناریو 6
Night 2/F 02068/0
1379/0
2068/0
3/122
3/122
6/43
1/8
1/8
10
130
130
140
8/53
8/53
1/80
130
130
140


Description: F:سناریوهاسناریو شماره 6نمودارهای سناریو 6Dispersion_Max Conc vs Distance.png
شکل 3: نمودار حداکثر غلظت در مقابل فاصله
Description: F:سناریوهاسناریو شماره 6نمودارهای سناریو 6Explosions_Worst Case Radii.png
شکل 4: نمودار شعاع انفجار در بدترین حالت


جدول 11: نتایج (محدوده) انتشار گاز
شماره سناریو کلاس جوی فاصله در جهت باد از محل سناریو تا حد پایین اشتعال گاز (m) بیشینه پوشش ابرگاز با غلظت ppm3/21663 در جهت باد (2m)
1 F 3/126 3/1072
A/B 6/154 5/1246
2 F 7/125 6/1066
A/B 8/155 1/1265
3 F 4/113 7/864
A/B 6/134 8/951
4 F 1/113 9/861
A/B 136 968
5 F 8/143 1/1608
A/B 197 9/2416
6 F 3/143 8/1599
A/B 8/197 1/2442
 
نتیجه‌گیری
با توجه ‌به یافته‌های تحقیق (گزارش پیامد سناریوهای موردمطالعه) بیشترین محدوده­های خطر مربوط به نتایج انتشار گاز و نتایج انفجار گاز می­باشند و لذا تعیین حریم ایمن در ایستگاه تقلیل فشار گاز طبیعی برون‌شهری واقع در منطقه برومی اهواز بر اساس نتایج انتشار و انفجار گاز صورت می­گیرد. برای روشن شدن بهتر نتایج علاوه بر اعداد گزارش‌شده در شکلهای 3 و4 به ترتیب حداکثر غلظت در مقابل فاصله و شعاع انفجار در بدترین حالت نمایش داده‌شده است. بر مبنای شکل 3 حداکثر غلظت گاز نشر یافته در بازه 10 تا 15 متری در مسیر جریان باد اتفاق افتاده است. به‌وضوح در شکل کاهش غلظت به سبب انتشار گاز در مسیر باد قابل‌مشاهده می‌باشد. این کاهش به کمینه مقدار خود در فاصله 140 متری و 200 متری به ترتیب در شب و روز خواهد رسید. همچنین بر مبنای سناریو بدترین حالت ممکن در مسیر باد شعاع انفجار می‌تواند تا 360 متر و در خلاف مسیر باد تا حدود 100 متر و در جهت عمود بر باد بیش از 220 متری توسعه یابد. همچنین با توجه به نتایج (محدوده) انتشار گاز (جدول 11) و نتایج (محدوده) انفجار گاز (جدول 12)، بیشترین میزان انتشار گاز برحسب فاصله در جهت باد تا حد پایین اشتعال گاز تعریف‌شده در نرم‌افزار PHAST و همچنین حداکثر فاصله در جهت باد در بدترین حالت تا شدت‌های موج انفجار تعریف‌شده در نرم‌افزار PHAST مربوط به سناریوهای شماره پنج و شش به ترتیب زیر می‌باشند :
  1. نتایج سناریوی شماره پنج برای انتشار گاز عبارتند از :
197 متر برای کلاس جوی A/B .
8/143 متر برای کلاس جوی F .
  1. نتایج سناریوی شماره شش برای انتشار گاز عبارتند از :
8/197 متر برای کلاس جوی  A/B.
3/143 متر برای کلاس جوی F .
  1. نتایج سناریوی شماره پنج برای انفجار گاز با شدت موجbar   02068/0 عبارتند از :
7/346 متر برای کلاس جوی A/B .
5/358 متر برای کلاس جوی F
  1. نتایج سناریوی شماره شش برای انفجار گاز با شدت موجbar   02068/0 عبارتند از :
342 متر برای کلاس جوی A/B
5/352 متر برای کلاس جوی F
با توجه به نتایج، بیشترین محدوده خطر متعلق به سناریوی شماره پنج برای انفجار گاز باکلاس جوی F و شدت موج bar 02068/0 و برابر 5/358 متر می‌باشد. در این محدوده با توجه به جدول خسارت ناشی از انفجار (جداول 10 و 12) به احتمال 95% خسارت جدی به وجود نمی‌آید اما 50% شیشه پنجره‌ها شکسته می‌شود و بدین‌جهت می‌تواند سبب صدمات انسانی گردد. بر این اساس حریم ایمن در ایستگاه تقلیل فشار گاز برون‌شهری منطقه برومی شهر اهواز برابر 5/358 متر تعیین میگردد.
 

جدول 12: نتایج (محدوده) انفجار گاز
شماره سناریو کلاس جوی شدت موج انفجار
(bar)
حداکثر فاصله از محل سناریو در جهت باد در بدترین حالت (m)
۱ F 02068/0
1379/0
2068/0
6/265
3/148
2/141
A/B 02068/0
1379/0
2068/0
9/271
2/170
1/165
۲ F 02068/0
1379/0
2068/0
3/262
6/147
7/140
A/B 02068/0
1379/0
2068/0
270
3/173
4/167
۳ F 02068/0
1379/0
2068/0
6/238
135
7/128
A/B 02068/0
1379/0
2068/0
8/246
7/152
147
۴ F 02068/0
1379/0
2068/0
9/235
4/134
3/128
A/B 02068/0
1379/0
2068/0
6/244
2/152
7/146
۵ F 02068/0
1379/0
2068/0
5/358
4/174
1/164
A/B 02068/0
1379/0
2068/0
7/346
5/190
4/190
۶ F 02068/0
1379/0
2068/0
5/352
2/173
163
A/B 02068/0
1379/0
2068/0
342
191
7/190
 
سپاسگزاری
این مقاله بر اساس طرح حمایت از پایان‌نامه‌های کارشناسی ارشد با پشتیبانی شرکت ملی گاز ایران اهواز انجام گردیده شده است. بدین‌وسیله از همکاری این شرکت و دانشگاه سیستان و بلوچستان کمال تشکر و قدردانی را داریم.



Refrences
1. Ebrahimi S, Agha Seyed Jafar Kashfi M. Identification and evaluation of environmental, social, safety, security and health consequences (ESHIA) of upstream projects of the oil and gas industry. Comparative Law Studies. 2015; 7(1): 1-22.
2. Anna D, Yuli W, Yi G. Development and application of screw expander in natural gas pressure energy recovery at city gas station. Applied Thermal Engineering. 2018; 142: 665-673.
3. Ebrahemzadih M, Foroghi nasab F, Mortazavi M, Soltani Gerdfaramarzi R. Analysis of Processing Accidents due to H2S in Clot Stuck Unit of one of the South Pars Refineries Leakage Using Determination of Safe Privacy Approach. Occupational Medicine Quarterly Journal. 2015; 7(3): 11-18.
4. Ghasemi M, Givehchi S, Nasrabadi M. Modeling The Outcome Of Tank Explosion Using Phast Software And Presentation Of Emergency Operations Plan (Case Study: Ethylene Tank Of Kavian Petrochemical Company). Journal of Fundamental and Applied Sciences. 2016; 8: 2007-2021.
5. Nowrozi B. Evaluating the consequences and determining the safety of gas pipelines with PHAST software in one of the provincial gas companies. Master's thesis in Chemical Engineering, Faculty of Chemical, Oil and Gas Engineering, 2016.
6. Setareshenas seghalani N. Analyzing the consequences of spreading chlorine gas in urban and industrial areas. Master's Thesis in Chemical Engineering, University of Sistan and Baluchistan, 2019.
7. Pourqadiri A, Ashrafi F, Farhadian Esfahani M. Safety, health and environment in pressure reduction stations and gas transmission lines. The first national HSE conference with the approach of upstream oil and gas industries, Abadan, 2013.
8. Mousavi SA. Quantitative assessment of risk caused by process accidents in Ramsar gas pressure boosting station, Master's thesis in chemical engineering. Islamic Azad University, Shahrood Branch, Faculty of Engineering, 2013.
9. Tabrizian SA. review of different methods of removing hydrogen sulfide. The second national conference of modern sciences and technologies of Iran, Tehran, 2015.
10. Ebrahimzadeh M, Foroghi Nesab F, Mortazavi Mehrizi M, Soltani Gardframarzi R. Analysis of process accidents caused by hydrogen sulfide gas leakage from the coagulation unit of one of the South Pars refineries with the approach of determining the safe zone. tkj. 2014; 7(3): 11-18.
11. Norian G, Amiri F. safety risk assessment of suburban natural gas pressure reduction station with FAST software. Fourth International Research Conference on Science and Technology, 2015.
12. Minabi A. Simulation of the concentration and distribution of hydrogen sulfide gas caused by the waste incinerators of the sulfur recycling unit in a gas refinery in Asalouye. Faculty of Environment and Energy, Islamic Azad University, Science and Research Unit, Tehran, 2016.
13. Movahed AR, Jahani F, Parvini M, Shakib M. Consequence Analysis of Gas Condensate Leakage in a Gas Refinery to Develop an Emergency Response Plan. J Occup Hyg Eng. 201 9 ; 6 ( 2): 1 - 8. DOI: 10.52547/johe. 6.2.1.
14. Nemati A. investigation of gas release damage from propane tank using PHAST software. Oil and gas exploration and production magazine. 2018; 168: 17-27.
15. Tian, Xinghao, Jiao, Wenling, Liu, Tianjie, Ren, Lemei, & Song, Bin (2021). Leakage detection of low-pressure gas distribution pipeline system based on linear fitting and extreme learning machine. International Journal of Pressure Vessels and Piping, v. doi:101016/jijpvp2021104553
16. Hao Xu, Hongchao Liu, Optimize Evacuation Route Considering the Operational Cost as a Constraint, International Journal of Transportation Science and Technology, Volume 2, Issue 1, 2013, Pages 65-78, ISSN 2046-0430, https://doi.org/10.1260/2046-0430.2.1.65.
17. Seung-Kwon S, Young-Gak Y, Ju-sung L, Jonggeol N, Chul-Jin L. Deep Neural Network-based Optimization Framework for Safety Evacuation Route during Toxic Gas Leak Incidents. Reliability Engineering & System Safety. 2022; 218: 108102-108118.
18. Abdolhamidzadeh B, Abbasi T, Rashtchian D, Abbasi SA. Domino effect in process-industry accidents–an inventory of past events and identification of some patterns. Journal of Loss Prevention in the Process Industries. 2011; 24(5): 575-593.
19. Haider Asadi S. Safety of natural gas in gas supply. National Gas Company, 2019.
20. Khanche Sepehruddin S. Analysis of the consequences of the explosion and fire incident in the suburban natural gas pressure reduction station number two, Urmia city. Master's thesis in chemical engineering, Urmia University of Technology. 2016.

پیوست 1. اطلاعات تکمیلی در خصوص جداول ناشی از خسارات

جدول پ-1. معیار ارزیابی آسیب ایجاده شده در اثر تشعشع آتش [20].
شدت تشعشع (kw/m2) آسیب
5/37 خسارت به واحدها و تجهیزات فرآیندی، ایجاد مرگ آنی برای افراد در معرض آن
20 آسیب جدی به افراد در معرض آن، در صورت نرسیدن تیم نجات موجب مرگ می‌شود
5/12 حداقل انرژی لازم برای ایجاد جرقه در پایلوت های چوبی و ذوب شدن مواد پلاستیکی
5/4 ایجاد درد در افرادی که حداقل 20 ثانیه در معرض آن می‌باشند، سوختگی درجه اول
6/1 در اثر تماس طولانی ایجاد عوارض نسبتاً خفیفی می کند
7/0 تابش خورشیدی

جدول پ-2. درصد تلفات ناشی از تشعشع آتش بر حسب مدت‌ زمان اعمال [20].

میزان تشعشع (kw/m2)
درصد تلفات ایجادشده بر حسب مدت‌زمان اعمال تشعشع ناشی از آتش (sec)
1% 50% %99
4 150 370 930
5/12 30 80 200
5/37 8 20 50

جدول پ-3. معیار ارزیابی آسیب ایجادشده در اثر موج انفجار [20].
آسیب شدت موج
(bar)
صدای آزار دهنده 0/0015
عدم خسارت به سازه ها، در برخی مواقع شکستن شیشه پنجره‌های بزرگ که قبلا تحت تنش کششی بوده اند 0/002
صدای بلند شبیه شکستن دیوار صوتی، در برخی مواقع شکستن شیشه ها 0/003
شکستن پنجره‌های کوچک تحت تنش کششی 007/0
حد آستانه عمومی برای شکستن پنجره‌ها 01/0
به احتمال 95% خسارت جدی در این محدوده وجود ندارد، 50% شیشه پنجره‌ها شکسته می‌شود 02/0
خسارت جزیی و محدود به سازه ها 03/0
پنجره‌ها عموما شکسته می‌شوند، در برخی مواقع صدمه به قاب پنجره‌ها 035/0-07/0
خسارت جزیی به سازه خانه ها 05/0
انهدام سقف تانک ها 07/0
تخریب قسمتی از خانه ها، غیر قابل سکونت شدن آن‌ها 08/0
آزبست های موج دار شکسته می‌شوند، بست پانل های فولادی یا آلومینیمی موج دار تخریب شده و خم خواهد شد، بست پانل های چوبی تخریب شده و پانل ها فرو می ریزند 07/0-15/0
تخریب جزیی قاب فلزی ساختمان های نماسازی شده 1/0
تخریب قسمتی از دیوارها و سقف خانه ها 15/0
تخریب بتن غیر مسلح یا دیوارهای ساخته شده با بلوک های سیمانی 15/0-2/0
حد پایین خسارت جدی به ساختمان ها، 50% تخریب قسمت های آجری خانه ها 18/0
ماشین آلات سنگین در ساختمان های صنعتی خسارت جزیی می بینند، قابهای فولادی ساختمان ها تخریب شده و از محل فونداسیون بیرون کشیده می‌شوند 2/0
ساختمان های دارای پانل فولادی قاب دار یا بدون قاب ویران می‌شوند، شکست مخازن ذخیره نفت 2/0-28/0
درخت های بزرگ سقوط می کنند 2/0-4/0
روکش فلزی ساختمان های صنعتی سبک پاره می‌شود، پانل های چوبی تخریب می‌شوند 3/0
شکست تیرهای چوبی برق و تلفن، تخریب بیشتر ساختمان ها بجز ساختمان های با دیواره های بتن مسلح، بدنه اتومبیل ها و کامیون ها به سمت داخل فرو می روند 35/0
جابجایی پل های لوله ها، شکستگی لوله کشی ها 35/0-4/0
تخریب تقریبا کامل خانه ها 35/0-5/0
فروریختن پل های لوله ها 4/0-55/0
ماشین های تانکردار پرشده، واگن های قطار و دیوارهای آجری به ضخامت 20-30 سانتیمتر واژگون و ویران می‌شوند 5/0
پانل های آجری غیر مسلح، دیوارهای آجری به ضخامت 25-35 سانتیمتر تخریب شده یا خم می‌شوند 5/0-55/0
قطارهای بارگیری شده کاملا تخریب می‌شوند 6/0
احتمال تخریب کلی ساختمان ها وجود دارد، ابزارآلات ماشینی سنگین جابجا شده و به شدت صدمه می بینند 7/0



 
[1] دانشجوی کارشناسی ارشد، گروه مهندسی شیمی، دانشکده مهندسی شهید نیکبخت، دانشگاه سیستان و بلوچستان، زاهدان، ایران
Consequences modeling and determining of safe distance in the natural gas pressure reduction station using PHAST software
(Case study: Borumi station in Ahvaz)

Torfi Alivi A 1, Shahraki F 2, Sardashti Birjandi MR 3[**], Khalilipour Langerudi MM 4
1 Masters student, Department of Chemical Engineering, Shahid Nikbakht Faculty of Engineering, University of Sistan and Baluchestan, Zahedan, Iran
2 Professor Department of Chemical Engineering, Shahid Nikbakht Faculty of Engineering, University of Sistan and Baluchestan, Zahedan, Iran
3 Assistant Professor, Department of Chemical Engineering, Shahid Nikbakht Faculty of Engineering, University of Sistan and Baluchestan, Zahedan, Iran
4 Assistant Professor, Department of Chemical Engineering, Shahid Nikbakht Faculty of Engineering, University of Sistan and Baluchestan, Zahedan, Iran
Abstract
Introduction: Gas pressure reduction stations are one of the most important parts of natural gas distribution networks. Risk assessment in operational units is a suitable method to assess risks, and their results can be used for management and decision-making regarding control and reduction of its consequences without critical concern. In other words, risk assessment is a process that identifies risks, assesses losses and determines risk characteristics. The present study aimed to model possible accidents, the subsequent evaluation of the consequences, and determination of safe boundaries in the suburban natural gas pressure reduction station located in the Bromi area of Ahvaz city using PHAST software.
Materials and Methods: In this study, possible scenarios based on the investigated background and utilizing the process data of the studied gas pressure reduction station, such as temperature, gas pressure, gas analysis, and meteorological data, and also using Pasquil's classification of atmospheric conditions the safe boundaries in the suburban natural gas pressure reduction station was investigated using PHAST software version 22.8 and. For this purpose, 6 possible scenarios were investigated in the selected atmospheric classes as well as the worst possible conditions. The consequences were evaluated based on different states of gas discharge, gas release, flash fire, sudden fire, and explosion.
Results: Based on the obtained results, the areas with high-risk potential were related to gas release and explosion, respectively. Also, the highest risk range belonging to scenario number five was reported for a gas explosion with atmospheric class F and wave intensity of 0.02068 bars and equivalent to 358.5 meters. In this range, according to the damage caused by the explosion, with a probability of 95%, without serious damage, 50% of the windows will be broken, and therefore it can cause human injuries. Based on this, the safe distance in this station was determined to be 358.5 meters.
Conclusion: Based on concerns about the lack of risk assessment studies over the initial phases and during the operation of the station, the study shows that determining the safe boundaries and high-risk areas along with the priority of risk reduction is very necessary. It can be concluded that the modeling of accidents and the evaluation of the resulting consequences along with the determination of safe boundaries is one of the basic factors in the health of the working environment of personnel and the general public.
Keywords: Toxic Gas Release, Explosion Modeling, Process Accident Modeling, Quantitative Risk Assessment, PHAST Software

This paper should be cited as:
Torfi Alivi A, Shahraki F, Sardashti Birjandi MR, Khalilipour Langerudi MM.. Consequences modeling and determining of safe distance in the natural gas pressure reduction station using PHAST software (Case study: Borumi station in Ahvaz). Occupational Medicine Quarterly Journal. 2023; 15(1): 37-57.
[2] استاد، گروه مهندسی شیمی، دانشکده مهندسی شهید نیکبخت، دانشگاه سیستان و بلوچستان، زاهدان، ایران
[3] استادیار، گروه مهندسی شیمی، دانشکده مهندسی شهید نیکبخت، دانشگاه سیستان و بلوچستان، زاهدان، ایران
[4] استادیار، گروه مهندسی شیمی، دانشکده مهندسی شهید نیکبخت، دانشگاه سیستان و بلوچستان، زاهدان، ایران
*(نویسنده مسئول)؛ تلفن 05431132440، پست الکترونیک: m.r_sardashti@eng.usb.ac.ir
تاریخ دریافت: 25/09/1401                                                                      تاریخ پذیرش:16/03/1402
* Corresponding Author:
Email: m.r_sardashti@eng.usb.ac.ir
Tel: +985431132440
Received: 16.12.2022                             Accepted: 06.06.2023
نوع مطالعه: پژوهشي | موضوع مقاله: ایمنی و حوادث ناشی از کار
دریافت: 1401/9/25 | پذیرش: 1402/3/16 | انتشار: 1402/3/29

ارسال نظر درباره این مقاله : نام کاربری یا پست الکترونیک شما:
CAPTCHA

ارسال پیام به نویسنده مسئول


بازنشر اطلاعات
Creative Commons License این مقاله تحت شرایط Creative Commons Attribution-NonCommercial 4.0 International License قابل بازنشر است.

کلیه حقوق این وب سایت متعلق به فصلنامه طب کار می باشد.

طراحی و برنامه نویسی : یکتاوب افزار شرق

© 2024 CC BY-NC 4.0 | Occupational Medicine Quarterly Journal

Designed & Developed by : Yektaweb