این پژوهش به منظور بررسی تاثیر فرآیندهای موثر بر حرکت آلاینده‎های نفتی (pahs) در پروفیل خاکی با جریان ناماندگار در ناحی? غیر اشباع خاک و الگوسازی ریاضی آن انجام شده است. با توجه به ماندگاری و مخاطره انگیز بودن آلاینده‎های نفتی، نشت این مواد به محیط زیست تبدیل به یکی از چالش‎های بزرگ گردیده است. یکی از مناطق کشور که با حجم گسترد? نشت مواد نفتی روبرو می‎باشد، محدود? پالایشگاه تهران است که نشت مواد و فاضلاب‎های نفتی پتانسیل بالقوه‎ای جهت تخریب منابع آب زیرزمینی و خاک را ایجاد نموده است. در این پژوهش در ابتدا با استفاده از مبانی حاکم بر حرکت آلاینده در محیط و با فرض انتقال یک بعدی قائم ماده در محیط متخلخل با منظور نمودن انتقال، انتشار، جذب سطحی، تبخیر و تجزی? بیولوژیکی معادل? کلی حرکت آب در خاک بررسی گردیده و سپس با استفاده از روش تقریبی اختلاف محدود معادلات گسسته سازی و با بکارگیری زبان برنامه نویسی fortran حل گردیده است. در نهایت با توجه به نتایج بدست آمده از غلظت آلاینده‎ها در 6 چاه مشاهده‎ای منطقه مورد مطالعه و الگوسازی آن میزان تاثیر هر یک از فرآیندهای یاد شده مورد بررسی قرار گرفته؛ و مشخص شده است که با توجه به ویژگی‎های آلایند? مورد مطالعه و شرایط خاک منطقه فرآیند جذب سطحی موثرترین فرآیند کاهنده بوده و دیگر فرآیندهای کاهنده از جمله تبخیر و تجزی? بیولوژیکی عمدتاً تا عمق کمی از سطح پروفیل خاکی موثر واقع شده اند. ضریب انتشار نیز منجر به جابجایی هرچه سریع‎تر ماده گردیده است و انتقال ماده نیز با پایان یافتن نفوذ جریان خاتمه نیافته و تا برقراری حالت تعادل در میدان جریان ادامه خواهد داشت. کلید واژه الگوسازی ریاضی ، هیدروکربن‎های آروماتیک چند حلقه‎ای ، انتقال ، انتشار ، اختلاف محدود