باتوجه به اینکه راندمان بازیافت گوگرد در پالایشگاه های نفت و گاز ایران حدود 90 درصد می باشد، لذا روزانه حدود پانصد تن دی اکسید گوگرد از طریق واحدهای مذکور در محیط زیست پراکنده می شود که اثرات سوء زیست محیطی فراوانی به همراه دارد. به همین جهت، پیش بینی نحوه انتشار دی اکسید گوگرد در محیط زیست از اهمیت بسیاری برخوردار است. در این پژوهش، ابتدا نسبت به مدلسازی ریاضی گذرای دوبعدی انتشار آلاینده مورد نظر در محیط زیست پرداخته و سپس با استفاده از روش تفاضل محدود صریحforward time central space (ftcs)و تفاضل محدود ضمنی crank nicolson (c-n) به حل معادلات دیفرانسیل پاره ای مدل ریاضی فرآیند در محیط matlab اقدام می گردد. پس از تایید روش حل (verification) در خصوص شبیه سازی انجام شده، تاثیر پارامترهایی نظیر تغییرات سرعت، جهت باد و نیز تعداد دود کش های موجود در منطقه بر الگوی پراکندگی آلاینده در محدوده مورد نظر تحت بررسی قرار گرفته است. همچنین با توجه به عدم دسترسی به داده های تجربی لازم، به منظور اطمینان از صحت شبیه سازی های انجام شده از تحلیل منطقی نتایج بهره گرفته شده است. نتایج حاصل از شبیه سازی انجام شده موید آن است که در صورت عدم استفاده از تکنولوژی های نوین جهت کاهش مقدار دی اکسید گوگرد به محیط زیست، انتشار بی رویه آلاینده مذکور می تواند باعث زیان های جبران ناپذیر به محیط زیست گردد.