فهرست مطالبفصل اول1مقدمه2فصل دوم5مروری برکارهای گذشته62-1- کاهش تولید در مخازن گاز میعانی62-2- روش های برطرف کردن انسداد میعانی62-3- بررسی رفتار فازی سیالات گاز میعانی9فصل سوم14مدل سازی153-1- مقدمه153-2- معادلات حاکم بر جریان گازی تک فازی153-3- معادلات حاکم بر جریان دو فازی163-4- معادلات حالت و مشتق های آن173-4-1- مقدمه:173-4-2- معادله ی حالت مکعبی معمولی SRK173-4-3- معادله ی حالت مکعبی همراه با همبستگی (Cubic Plus Association)203-5- حل معادلات با روش عددی کولوکیشن [31]293-6-تست کاهش شار در حجم ثابت323-7- فرآیند حل معادلات مخزن با کمک روش ریاضی کولوکیشن343-7-1- تحلیل روش عددی کولوکیشن برای سیستم های تک فازی در مختصات استوانه ای در جهت محوری343-7-2- تحلیل روش عددی کولوکیشن برای سیستم های دو فازی امتزاج ناپذیر پایا در مختصات استوانه ای در جهت محوری383-7-3- تحلیل روش عددی کولوکیشن برای سیستم های دو فازی گاز میعانی در مختصات استوانه ای در جهت محوری41فصل چهارم51نتایج و بحث مدل سازی524-1- مقدمه524-2- بررسی رفتار فازی با استفاده از معادله ی حالت CPA524-2-1- تست مدل CPA برای حالت خالص534-2-2- تست مدل CPA برای حالت چند جزئی574-3-تست های رفتار فازی724-3-1-مقدمه724-3-2-تست تخلیه در حجم ثابت724-3-3- تست Flash744-4- مدل سازی جریان در محیط متخلخل794-4-1- مقدمه794-4-2- مدل سازی جریان سیال در هندسه ی مغزه804-4-2-1- جریان تک فازی در مختصات استوانه ای و جریان محوری804-4-2-2- جریان دو فازی امتزاج ناپذیر پایا در مختصات استوانه ای و جریان محوری824-4-2-3- جریان دو فازی گاز میعانی پایا در مختصات استوانه ای و جریان محوری844-4-2-4- جریان دو فازی گاز میعانی پایا به همراه حلال در مختصات استوانه ای و جریان محوری864-4-2-5- جریان دو فازی گاز میعانی ناپایا در مختصات استوانه ای و جریان محوری894-4-2-6- جریان دو فازی در مختصات استوانه ای و جریان محوری بدون تزریق حلال و با تزریق آن94فصل پنجم106تست آزمایشگاهی و نتایج و بحث آن1075-1- مقدمه1075-2- دستگاه سیلاب زنی مغزه و اجزای آن1075-2-1 مخزن نگه دارنده ی سیال1085-2-2- مغزه نگه دار1095-2-3- پمپ1095-2-4- محفظه ی گرم کننده1095-2-5- ریگلاتور1105-2-6- سیستم نمایشگر اختلاف دما و فشار1105-3- انجام آزمایش1105-3-1- آماده سازی دستگاه1105-3-2- مراحل انجام آزمایش1125-4- نتایج حاصل از انجام آزمایش و بحث روی آن113فصل ششم116نتیجه گیری و پیشنهاد ها1176-1- نتیجه گیری1176-2- پیشنهاد ها117فهرست منابع118 فهرست جدول ها عنوان صفحهجدول3-1- ثابت های معادله ی SRK18جدول3-2- معادلات جهت بدست آوردن XAها (هوانگ و رادوز 1990) [28] برای سیستم های خود- همبسته و خالص24جدول4-1- پارامترهای مربوط به معادله ی CPA53جدول4-2- ضریب برهم کنش دوجزئی مخلوط متانول و هیدروکربن [29]58جدول4-3- ضریب برهم کنش دوجزئی مخلوط مونو اتیلن گلیکول و هیدروکربن [14]58جدول4-4- ضریب برهم کنش دوجزئی مخلوط دی اتیلن گلیکول و هیدروکربن [14]58جدول4-5- ضریب برهم کنش دوجزئی مخلوط تری اتیلن گلیکول و هیدروکربن [14]59جدول4-6- ضریب برهم کنش دوجزئی مخلوط دو متوکسی اتانول و هیدروکربن [20 و 21]60جدول4-7- ضریب برهم کنش دوجزئی مخلوط دو اتوکسی اتانول و هیدروکربن [20]61جدول4-8- ضریب برهم کنش دوجزئی مخلوط دو بوتوکسی اتانول و هیدروکربن [20 و23]61جدول4-9- جزء مولی سیال گاز میعانی ساختگی دو و همکارانش(2000)[3]70جدول4-10- غلظت مخلوط گاز میعانی ساختگی شماره ی 1 مورد نظر این پروژه79جدول4-11- غلظت مخلوط گاز میعانی ساختگی شماره ی 2 مورد نظر این پروژه80جدول4-12- داده های لازم برای خواص سنگ مخزنی80جدول4-13- نتایج محاسبه ی خطای APD برای نقاط محاسبه ای مختلف81جدول4-14- داده های مربوط به سنگ و سیالهای امتزاج ناپذیر82جدول4-15- داده های مربوط به سنگ و سیالهای امتزاج پذیر84جدول4-16- داده های مربوط به سنگ و سیالهای امتزاج پذیر در جریان پایای دو فازی87جدول4-17- داده های مربوط به سنگ و سیالهای امتزاج پذیر در جریان ناپایای دو فازی89جدول4-18- داده های لازم برای خواص سنگ مخزنی94جدول4-19- داده های لازم برای خواص سنگ مخزنی101جدول4-20- نتایج حاصل از تزریق حلال های گلیکول اتری105جدول5-1- جرم لازم از هپتان نرمال برای رسیدن به غلظت معلوم دما و فشار تعیین شده111جدول5-2- نتایج حاصل از تزریق حلال های گلیکول اتری115 فهرست شکل ها عنوان صفحهشکل3-1- فلوچارت محاسبه ی ضریب فوگاسیته با کمک معادله ی حالت CPA22شکل3-2- ساختار مولکولی آب و دو متوکسی اتانول25شکل3-3- فلوچارت مربوط به محاسبه ی X ها28شکل3-4- شماتیکی از تست CVD [24]33شکل3-5- نمونه ی ساده ای از الگوریتم فرآیند CVD33شکل3-6- شکل مربوط به جریان محوری و جزء دیفرانسیلی آن34شکل4-1- منحنی ضریب دوم ویریال ناشی از مدل برای مونو اتیلن گلیکول54شکل4-2- منحنی ضریب دوم ویریال ناشی از مدل برای پروپیلن گلیکول54شکل4-3- منحنی فشار بخار برای مونو اتیلن گلیکول، داده های تجربی گرفته شده از[37]55شکل4-4- منحنی فشار بخار برای پروپیلن گلیکول، داده های تجربی گرفته شده از [37]55شکل4-5- منحنی فشار بخار برای اتیلن گلیکول مونو متیل اتر ،داده های تجربی گرفته شده از [38]56شکل4-6- منحنی فشار بخار برای اتیلن گلیکول بوتیل اتر ،داده های تجربی گرفته شده از [38]56شکل4-7- محاسبه ی چگالی اتیلن گلیکول مونو متیل اتر برای دو فاز مایع و گاز با CPAو مقایسه ی آن با داده های تجربی[38]57شکل4-8- محاسبه ی چگالی اتیلن گلیکول مونو بوتیل اتر برای دو فاز مایع و گاز با CPA و مقایسه ی آن با داده های تجربی[38]57شکل4-9- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر62شکل4-10- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین62شکل4-11- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط دی اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین63شکل4-12- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط تری اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین63شکل4-13- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط پروپیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین64شکل4-14- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از [39]64شکل4-15- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط دی اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از [39]65شکل4-16- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط تری اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از[39]65شکل4-17- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط پروپیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر[داده های تجربی از مرجع 39]66شکل4-18- مقایسه ی حلالیت های گلیکول ها در نرمال هپتان در دمای 351 درجه ی کلوین66شکل4-19- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از[21]67شکل4-20- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از[21]67شکل4-21- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال هگزا دکان در فشار 1 اتمسفر ، داده های تجربی از [21]68شکل4-22- منحنی لگاریتمی غلظت نرمال هپتان بر حسب دما برای مخلوط دو متوکسی اتانول و نرمال اکتان دکان در فشار 1 اتمسفر، داده های تجربی از[21]68شکل4-23- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو اتوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از [21]69شکل4-23- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو بوتوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از[23]69شکل4-24- منحنی فشار بر حسب غلظت دومتوکسی اتانول در دمای 135 درجه ی فارنهایت70شکل4-25- منحنی فشار بر حسب غلظت دو اتوکسی اتانول در دمای 145 درجه ی فارنهایت71شکل4-26- منحنی فشار بر حسب غلظت دو بوتوکسی اتانول در دمای 145 درجه ی فارنهایت71شکل4-27– منحنی غلظت اجزای هیدروکربی سیال برحسب فشار در تستCVD73شکل4-28– منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو متوکسی اتانول73شکل4-29- منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو اتوکسی اتانول74شکل4-30- منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو بوتوکسی اتانول74شکل4-31- منحنی اشباع مایع بر حسب غلظت کلی دو متوکسی اتانول75شکل4-32- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو متوکسی اتانول75شکل4-33- منحنی غلظت دو متوکسی اتانول در فاز گاز بر حسب غلظت کلی دو متوکسی اتانول76شکل4-34- منحنی اشباع مایع بر حسب غلظت کلی دو اتوکسی اتانول76شکل4-35- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو اتوکسی اتانول77شکل4-36- منحنی غلظت دو اتوکسی اتانول در فاز گاز بر حسب غلظت کلی دو اتوکسی اتانول77شکل4-37- منحنی اشباع مایع بر حسب غلظت کلی دو بوتوکسی اتانول78شکل4-38- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو اتوکسی اتانول78شکل4-39- منحنی غلظت دو اتوکسی اتانول در فاز گاز بر حسب غلظت کلی دو اتوکسی اتانول79شکل4-40– منحنی افت فشار بر حسب طول مغزه برای حالت ناپایای حالت تک فازی و جریان محوری81شکل4-41- مقایسه ی اثر تعداد نقاط محاسبه ای برای افت فشار کل با زمان در سیستم تک فازی و جریان محوری82شکل4-42- منحنی فشار بر حسب طول مغزه83شکل4-43- منحنی اشباع گاز بر حسب طول مغزه83شکل4-44- منحنی فشار بر حسب طول مغزه (N=6)83شکل4-45- منحنی اشباع گاز بر حسب طول مغزه (N=6)83شکل4-46- اثر تعداد نقاط محاسبه ای بر روی منحنی اشباع مایع و فشار بر حسب طول85شکل4-47- بررسی اثر تعداد نقاط محاسبه ای بر منحنی غلظت متان در مخلوط بر حسب طول در دو فاز گاز و مایع85شکل4-48- منحنی قابلیت های حرکت فاز ها بر حسب طول86شکل4-49- بررسی همراه بودن دو متوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز87شکل4-50- بررسی همراه بودن دواتوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز88شکل4-51- بررسی همراه بودن دو بوتوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز88شکل4-52- منحنی فشار بر حسب طول برای بازه ی 5 تا 1500 ثانیه. (برای بازه های 50 ثانیه ای)90شکل4-53- منحنی فشار بر حسب زمان برای چند مقطع مغزه90شکل4-54- اشباع گاز بر حسب طول برای زمان 5 تا 1500 ثانیه (در بازه های 100 ثانیه ای)91شکل4-55- اشباع گاز بر حسب زمان برای مقاطع مختلف از مغزه91شکل4-56- غلظت متان در فاز مایع بر حسب طول از زمان 5 ثانیه تا 1500 ثانیه ( در بازه های 100 ثانیه ای)92شکل4-57- غلظت متان بر حسب زمان برای مقاطع مختلف92شکل4-58- غلظت متان در فاز گاز بر حسب طول برای زمان های بین 5 ثانیه تا 1500 ثانیه (در بازه های زمانی 100 ثانیه)93شکل4-59- غلظت متان در فاز گاز بر حسب زمان برای چهار مقطع93شکل4-60- منحنی افت فشار برای جریان تک فازی بر حسب زمان95شکل4-61- منحنی های افت فشار و اشباع ورودی سنگ برای جریان دو فازی تعریف شده95شکل4-62- منحنی های تراوایی نسبی دو فاز گاز و مایع برای جریان دو فازی مورد توجه96شکل4-63- مقایسه ی افت فشار برای جریان تک فازی و دو فازی97شکل4-64- منحنی افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دومتوکسی اتانول98شکل4-65- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو متوکسی اتانول98شکل4-66- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دواتوکسی اتانول99شکل4-67- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو اتوکسی اتانول100شکل4-68- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دوبوتوکسی اتانول100شکل4-69- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو بوتوکسی اتانول101شکل4-70- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو متوکسی اتانول102شکل4-71- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو متوکسی اتانول102شکل4-72- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو اتوکسی اتانول103شکل4-73- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو اتوکسی اتانول103شکل4-74- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو بوتوکسی اتانول104شکل4-75- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو بوتوکسی اتانول104شکل5-1- شمای کلی از دستگاه سیلاب زنی108شکل5-2- منحنی رفتار فازی مخلوط متان و هپتان نرمال (85 % متان)111شکل5-3- منحنی تراوایی مطلق بر حسب زمان114شکل5-4- منحنی تراوایی نسبی بر حسب زمان قبل و بعد از تزریق دو متوکسی اتانول114شکل5-5- منحنی تراوایی نسبی بر حسب زمان قبل و بعد از تزریق دو اتوکسی اتانول115فهرست نشانه های اختصاری aترم انرژی در معادله ی حالتzکسر مولی کلیbترم حجمی درمعادله ی حالتJماتریس مشتقاتعلایم یونانیkتراوایی∆اختلافLطول مغزه∆ABشدت همبستگیMwجرم مولیδ(x)دلتای دیراکncتعداد اجزای مخلوطμگرانرویPcفشار بحرانیρچگالیRثابت جهانی گاز هاφتخلخل سنگSاشباعضریب فوگاسیته جزء در مخلوطTدماTcدمای بحرانیزیر نویس هاtزمانiشمارندهvحجم مولیjشمارندهwتابع آزمونkشمارندهxکسر مولی در فاز مایعAسایت Ayکسر مولی در فاز گازgگازZضریب تراکم پذیریlمایع مقدمه کشور ایران از نظر منابع گازی در جهان جایگاه خوبی قرار دارد. این کشور دارای ذخیره ی 981 تريليون مترمكعب به صورت در جا و تولید روزانه ی 9/111 بيليون متر مكعب است. با توجه به اینکه در ناحیه جنوبی کشور بسیاری از این مخازن از نوع گاز میعانی است، تولید با بهره وری بهتر لازم به نظر می رسد.از آنجایی که میادین گاز میعانی در اقتصاد کشور نقش بسزایی دارد، بنابراین پژوهش بر روی عملکرد این منابع ارزشمند جنبه ی حیاتی دارد. این مهم در داخل کشور چنان مورد توجه واقع شده است که در شرکت ملی نفت ایران زیر نظر مدیریت پژوهشی و فناوری آن شرکت واحدی با عنوان گروه پژوهش مخازن گاز میعانی در پژوهشکده ی ازدیاد برداشت ایجاد شده است.مخازن گازی به سه دسته تقسیم می شوند. درصد بيشتري از مخازن گازي، جزء مخازن گاز خشكهستند. این مخازن ضمن توليد گاز در سر چاههيچگونه ميعاناتي با خود توليد نمي کنند.دسته دوم مخازن گازي، مخازن گاز مرطوبهستند. در اين مخازنمقداري ميعانات ضمن توليد در سر چاه بدست مي آيد که معمولادر هر چاه بازاء توليد 1 ميليون فوت مكعبگاز در شرايط استاندارد، بین 3 تا20 بشكه ميعانات در سطح توليد مي گردد. این میعانات به علت سبکی می توانند در تولید محصولات متنوعی نظيرحلال ها، بنزين و خوراك واحد هاي پتروشيمي استفاده شوند.دسته سوم مخازنگازي، مخازن گاز ميعاني هستند. در این میادین میزان تولید میعانات در بیشترین حد قرار دارد. در مواردي حتي ممكن است بازاء توليد 1 ميليون فوت مكعب گاز در شرايطاستاندارد، 300 بشكه از اين ميعانات در سطح توليد گردد. تفاوت مهم ديگر اين مخازنبا مخازن گاز مرطوب اين است كه در مخازن گاز ميعاني بر خلاف مخازن گاز مرطوب، ممكناست مقدار زيادي از ميعانات در درون مخزن تشكيل گردد. دو نقطه ضعف اساسي تشكيل اينميعانات در درون مخزن به قرار زير است:1- تشكيل اين ميعانات ارزشمند در درونمخزن باعث مي شود تا به روش هاي متداول بازيافت نتوان آنها را توليد نمود و به اينترتيب بخش مهمي از سرمايه هيدروكربني غير قابل برداشت مي گردد. به طور كلي فضايدرون مخزن شبيه محيط اسفنجي مي باشد كه خلل و فرج آن توسط رابط هاي مويينه بهيكديگر وصل شده اند و براي شروع حركت يك سيال در اين محيط بسيار متراكم اسفنجي لازماست تا سيال به يك ميزان مشخص در آن محيط تشكيل شود و درصد معيني از حجم محيط رااشغال نمايد. اين بدان معناست كه پس از توليد اولين قطرات مايع در مخزن ميعاناتبدون حركت مانده و تا رسيدن به درصد اشباع معين در درون مخزن باقي خواهندماند.2- تشكيل اين ميعانات در درون محيط مخزن باعث انسداد مسير درون حفرهها ولولههاي مويينه محيط متخلخل مخزن شده و يك مقاومت اضافي در برابر حركت جريان گازاز درون مخزن به سمت چاه و نهايتاً توليد گاز از مخزن ايجاد ميكند و بهرهبردارياز مخزن را كاهش ميدهد.در فرا سوی مرز های ایران اسلامی، مراکز تحقیقاتی کشور های دیگر هم به اهمیت این مهم پی برده اند و سالانه مقالات و گزارشهای فراوانی را در مجلات و همایش های معتبر به چاپ می رسانند و ارائه می دهند.در بررسی مخازن هیدروکربنی زیر زمینی تا زمانی که با نفت سر و کار داریم، مطابق با انتظار مان با کاهش فشار در دمای ثابت با تشکیل گاز روبرو هستیم. اما این وضعیت برای مخازن گازی تر[1] و گاز میعانی[2] جور دیگری است. با کاهش فشار در دمای ثابت با تشکیل مایع برخورد می کنیم. این پدیده را میعان معکوس[3] نام می گذارند. با وجود چنین پدیده ای در سنگ مخزنی، با برداشت از آن به مرور زمان با افزایش مقدار مایع در حفره های مخزن روبرو خواهیم بود. چنین پدیده ای با گذشت زمان باعث ایجاد مانعی بر سر راه تولید گاز از چاه های آن منطقه می گردد. چنین پدیده ای را انسداد میعانی[4] می گویند.روش های مقابله با این پدیده تنوع زیادی دارد. در این پایان نامه هدف استفاده از حلال های مایع گلیکول و گلیکول اتری است. حلال از راه چاه تولید به مخزن تزریق می شود. حوزه ی چاه مدتی در همان شرایط می ماند تا حلال با میعانات و گاز داخل حفره های سنگ مخزن به تعادل برسد. با شروع دوباره ی تولید ابتدا میعانات به بیرون هدایت می شوند و گاز در این شرایط قابلیت حرکت بیشتری دارد و تولید آن افزایش می یابد. حلال های مورد نظر تاثیری بر روی سنگ مخزنی ندارند و تنها اثر آنها بر تعادل فازی است. مقداری از حلال که در حفره ها باقی می ماند، به میعانات اجازه ی ورود به فاز مایع را نمی دهد. بنابراین تا مدتی که حلال تبخیر نشود. سطح مایع پایین می ماند و استحصال گاز با افت فشار کمتری روبرو است.هدف از این پایان نامه مدل سازی اثر حلال های گلیکول و گلیکول اتری بر رفتار فازی سیال مخزنی در شرایط مخزن در مقیاس مغزه همراه با حل معادلات جریان دو فازی در مختصات استوانه ای در جهت شعاعی و مقایسه اثر آن در مقیاس آزمایشگاهی است. نوآوری های کار استفاده از یک روش عددی در حل معادلات مخزن و استفاده از یک معادله ی حالت موثر برای بررسی رفتار فازی حلال ها در کنار سیال هیدروکربنی است.فصل دوم مروری بر کار های گذشته در این فصل خلاصه ی مقالاتی که به استفاده ی حلال های مختلف در برطرف کردن پدیده ی انسداد میعانی پرداخته است ارائه می شود که در چند قسمت آورده شده است. 2-1- کاهش تولید در مخازن گاز میعانی افیدیک و همکارانش (1994) [1] به بررسی کاهش تولید در مخزن گاز میعانی آرون در اثر انباشتگی مععانات پرداختند و به این نتیجه رسیدند که با کاهش فشار به زیر فشار نقطه ی شبنم تولید چاه با فاکتوری حدود 2 کاهش می یابد که این در اثر تجمع میعانات در دهانه چاه است.بارنوم و همکارانش(1995) [2] به این نتیجه رسیدند که کاهش تولید برای مخازن با تولید کم با Kh کمتر از 1000 md-ft شدید است. 2-2- روش های برطرف کردن انسداد میعانی روش های متفاوتی برای مقابله با پدیده ی انسداد میعانی وجود دارد. دو روش عمده برای این برطرف کردن این پدیده به قرار زیر است:- تغییر رفتار فازی گاز میعانی- کاهش افت فشار و نگه داشتن فشار بالای نقطه ی شبنمیکی از روش ها استفاده از حلال ها است.حلال های مورد نظر از خانواده ی الکل ها، گیلکول ها و گلیکول اتر ها است که در زیر به شرح کارهایی که در این زمینه شده است می پردازیم:دو و همکارانش(2000)[3] سعی کردند که کاهش نفوذ پذیری نسبی توسط میعانات را با استفاده از یک حلال که آن حلال متانول است بهبود بخشند.دو و همکارانش احتمال می دانند که با تزریق متانول کشش بین سطحی گاز و مایع کم می شود ولی نتایج اندازه گیری این فرض را رد کرد. احتمال دیگری که مورد توجه آن ها بود اثر متانول روی رفتار فازی هیدروکربن بود که با کمک معادله حالت این فرض را مورد آزمایش قرار داد و مشاهده کردند که اگر غلظت متانول زیاد باشد در شرایط عبوری از مغزه یک فاز مایع غنی از متانول تشکیل می شود که مانع تشکیل فاز میعانی دیگری می گردد که باعث بهبود تراوایی نسبی می شود.نتایجی که از تزریق متانول بدست می آید به قرار زیر است:1- تزریق متانول به درون مغزه که نفوذ پذیری آن توسط انسداد میعانی کم شده باعث می شود که نفوذ پذیری نسبی بهتر شود.2- دلیل افزایش نفوذ پذیری نسبی این است که متانول باعث جابجایی فاز آب و میعانات می شود که این عمل توسط حل شدن متانول صورت می گیرد .3- با تزریق متانول آب کاملا حذف می شود.4- استفاده از متانول ودیگر حلال های الکلی روش های ارزان تری نسبت دیگر روش ها مانند شکست هیدرولیکی می باشند
بررسی اثر گلیکول ها و گلیکول اتر ها بر پدیده ی انسداد میعانی مدل سازی word
فهرست مطالبفصل اول1مقدمه2فصل دوم5مروری برکارهای گذشته62-1- کاهش تولید در مخازن گاز میعانی62-2- روش های برطرف کردن انسداد میعانی62-3- بررسی رفتار فازی سیالات گاز میعانی9فصل سوم14مدل سازی153-1- مقدمه153-2- معادلات حاکم بر جریان گازی تک فازی153-3- معادلات حاکم بر جریان دو فازی163-4- معادلات حالت و مشتق های آن173-4-1- مقدمه:173-4-2- معادله ی حالت مکعبی معمولی SRK173-4-3- معادله ی حالت مکعبی همراه با همبستگی (Cubic Plus Association)203-5- حل معادلات با روش عددی کولوکیشن [31]293-6-تست کاهش شار در حجم ثابت323-7- فرآیند حل معادلات مخزن با کمک روش ریاضی کولوکیشن343-7-1- تحلیل روش عددی کولوکیشن برای سیستم های تک فازی در مختصات استوانه ای در جهت محوری343-7-2- تحلیل روش عددی کولوکیشن برای سیستم های دو فازی امتزاج ناپذیر پایا در مختصات استوانه ای در جهت محوری383-7-3- تحلیل روش عددی کولوکیشن برای سیستم های دو فازی گاز میعانی در مختصات استوانه ای در جهت محوری41فصل چهارم51نتایج و بحث مدل سازی524-1- مقدمه524-2- بررسی رفتار فازی با استفاده از معادله ی حالت CPA524-2-1- تست مدل CPA برای حالت خالص534-2-2- تست مدل CPA برای حالت چند جزئی574-3-تست های رفتار فازی724-3-1-مقدمه724-3-2-تست تخلیه در حجم ثابت724-3-3- تست Flash744-4- مدل سازی جریان در محیط متخلخل794-4-1- مقدمه794-4-2- مدل سازی جریان سیال در هندسه ی مغزه804-4-2-1- جریان تک فازی در مختصات استوانه ای و جریان محوری804-4-2-2- جریان دو فازی امتزاج ناپذیر پایا در مختصات استوانه ای و جریان محوری824-4-2-3- جریان دو فازی گاز میعانی پایا در مختصات استوانه ای و جریان محوری844-4-2-4- جریان دو فازی گاز میعانی پایا به همراه حلال در مختصات استوانه ای و جریان محوری864-4-2-5- جریان دو فازی گاز میعانی ناپایا در مختصات استوانه ای و جریان محوری894-4-2-6- جریان دو فازی در مختصات استوانه ای و جریان محوری بدون تزریق حلال و با تزریق آن94فصل پنجم106تست آزمایشگاهی و نتایج و بحث آن1075-1- مقدمه1075-2- دستگاه سیلاب زنی مغزه و اجزای آن1075-2-1 مخزن نگه دارنده ی سیال1085-2-2- مغزه نگه دار1095-2-3- پمپ1095-2-4- محفظه ی گرم کننده1095-2-5- ریگلاتور1105-2-6- سیستم نمایشگر اختلاف دما و فشار1105-3- انجام آزمایش1105-3-1- آماده سازی دستگاه1105-3-2- مراحل انجام آزمایش1125-4- نتایج حاصل از انجام آزمایش و بحث روی آن113فصل ششم116نتیجه گیری و پیشنهاد ها1176-1- نتیجه گیری1176-2- پیشنهاد ها117فهرست منابع118 فهرست جدول ها عنوان صفحهجدول3-1- ثابت های معادله ی SRK18جدول3-2- معادلات جهت بدست آوردن XAها (هوانگ و رادوز 1990) [28] برای سیستم های خود- همبسته و خالص24جدول4-1- پارامترهای مربوط به معادله ی CPA53جدول4-2- ضریب برهم کنش دوجزئی مخلوط متانول و هیدروکربن [29]58جدول4-3- ضریب برهم کنش دوجزئی مخلوط مونو اتیلن گلیکول و هیدروکربن [14]58جدول4-4- ضریب برهم کنش دوجزئی مخلوط دی اتیلن گلیکول و هیدروکربن [14]58جدول4-5- ضریب برهم کنش دوجزئی مخلوط تری اتیلن گلیکول و هیدروکربن [14]59جدول4-6- ضریب برهم کنش دوجزئی مخلوط دو متوکسی اتانول و هیدروکربن [20 و 21]60جدول4-7- ضریب برهم کنش دوجزئی مخلوط دو اتوکسی اتانول و هیدروکربن [20]61جدول4-8- ضریب برهم کنش دوجزئی مخلوط دو بوتوکسی اتانول و هیدروکربن [20 و23]61جدول4-9- جزء مولی سیال گاز میعانی ساختگی دو و همکارانش(2000)[3]70جدول4-10- غلظت مخلوط گاز میعانی ساختگی شماره ی 1 مورد نظر این پروژه79جدول4-11- غلظت مخلوط گاز میعانی ساختگی شماره ی 2 مورد نظر این پروژه80جدول4-12- داده های لازم برای خواص سنگ مخزنی80جدول4-13- نتایج محاسبه ی خطای APD برای نقاط محاسبه ای مختلف81جدول4-14- داده های مربوط به سنگ و سیالهای امتزاج ناپذیر82جدول4-15- داده های مربوط به سنگ و سیالهای امتزاج پذیر84جدول4-16- داده های مربوط به سنگ و سیالهای امتزاج پذیر در جریان پایای دو فازی87جدول4-17- داده های مربوط به سنگ و سیالهای امتزاج پذیر در جریان ناپایای دو فازی89جدول4-18- داده های لازم برای خواص سنگ مخزنی94جدول4-19- داده های لازم برای خواص سنگ مخزنی101جدول4-20- نتایج حاصل از تزریق حلال های گلیکول اتری105جدول5-1- جرم لازم از هپتان نرمال برای رسیدن به غلظت معلوم دما و فشار تعیین شده111جدول5-2- نتایج حاصل از تزریق حلال های گلیکول اتری115 فهرست شکل ها عنوان صفحهشکل3-1- فلوچارت محاسبه ی ضریب فوگاسیته با کمک معادله ی حالت CPA22شکل3-2- ساختار مولکولی آب و دو متوکسی اتانول25شکل3-3- فلوچارت مربوط به محاسبه ی X ها28شکل3-4- شماتیکی از تست CVD [24]33شکل3-5- نمونه ی ساده ای از الگوریتم فرآیند CVD33شکل3-6- شکل مربوط به جریان محوری و جزء دیفرانسیلی آن34شکل4-1- منحنی ضریب دوم ویریال ناشی از مدل برای مونو اتیلن گلیکول54شکل4-2- منحنی ضریب دوم ویریال ناشی از مدل برای پروپیلن گلیکول54شکل4-3- منحنی فشار بخار برای مونو اتیلن گلیکول، داده های تجربی گرفته شده از[37]55شکل4-4- منحنی فشار بخار برای پروپیلن گلیکول، داده های تجربی گرفته شده از [37]55شکل4-5- منحنی فشار بخار برای اتیلن گلیکول مونو متیل اتر ،داده های تجربی گرفته شده از [38]56شکل4-6- منحنی فشار بخار برای اتیلن گلیکول بوتیل اتر ،داده های تجربی گرفته شده از [38]56شکل4-7- محاسبه ی چگالی اتیلن گلیکول مونو متیل اتر برای دو فاز مایع و گاز با CPAو مقایسه ی آن با داده های تجربی[38]57شکل4-8- محاسبه ی چگالی اتیلن گلیکول مونو بوتیل اتر برای دو فاز مایع و گاز با CPA و مقایسه ی آن با داده های تجربی[38]57شکل4-9- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر62شکل4-10- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین62شکل4-11- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط دی اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین63شکل4-12- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط تری اتیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین63شکل4-13- منحنی فشار بر حسب غلظت نرمال هپتان برای مخلوط پروپیلن گلیکول و نرمال هپتان در دمای 351 درجه ی کلوین64شکل4-14- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط مونو اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از [39]64شکل4-15- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط دی اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از [39]65شکل4-16- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط تری اتیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از[39]65شکل4-17- منحنی مایع- مایع دما بر حسب غلظت نرمال هپتان برای مخلوط پروپیلن گلیکول و نرمال هپتان در فشار 1 اتمسفر[داده های تجربی از مرجع 39]66شکل4-18- مقایسه ی حلالیت های گلیکول ها در نرمال هپتان در دمای 351 درجه ی کلوین66شکل4-19- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال هپتان در فشار 1 اتمسفر، داده های تجربی از[21]67شکل4-20- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از[21]67شکل4-21- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو متوکسی اتانول و نرمال هگزا دکان در فشار 1 اتمسفر ، داده های تجربی از [21]68شکل4-22- منحنی لگاریتمی غلظت نرمال هپتان بر حسب دما برای مخلوط دو متوکسی اتانول و نرمال اکتان دکان در فشار 1 اتمسفر، داده های تجربی از[21]68شکل4-23- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو اتوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از [21]69شکل4-23- منحنی دما بر حسب غلظت نرمال هپتان برای مخلوط دو بوتوکسی اتانول و نرمال اکتان در فشار 1 اتمسفر ، داده های تجربی از[23]69شکل4-24- منحنی فشار بر حسب غلظت دومتوکسی اتانول در دمای 135 درجه ی فارنهایت70شکل4-25- منحنی فشار بر حسب غلظت دو اتوکسی اتانول در دمای 145 درجه ی فارنهایت71شکل4-26- منحنی فشار بر حسب غلظت دو بوتوکسی اتانول در دمای 145 درجه ی فارنهایت71شکل4-27– منحنی غلظت اجزای هیدروکربی سیال برحسب فشار در تستCVD73شکل4-28– منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو متوکسی اتانول73شکل4-29- منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو اتوکسی اتانول74شکل4-30- منحنی مقدار مایع خروجیبر حسب فشار در دمای 360 درجه ی کلوین برای غلظت های مختلف دو بوتوکسی اتانول74شکل4-31- منحنی اشباع مایع بر حسب غلظت کلی دو متوکسی اتانول75شکل4-32- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو متوکسی اتانول75شکل4-33- منحنی غلظت دو متوکسی اتانول در فاز گاز بر حسب غلظت کلی دو متوکسی اتانول76شکل4-34- منحنی اشباع مایع بر حسب غلظت کلی دو اتوکسی اتانول76شکل4-35- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو اتوکسی اتانول77شکل4-36- منحنی غلظت دو اتوکسی اتانول در فاز گاز بر حسب غلظت کلی دو اتوکسی اتانول77شکل4-37- منحنی اشباع مایع بر حسب غلظت کلی دو بوتوکسی اتانول78شکل4-38- منحنی غلظت دو متوکسی اتانول در فاز مایع بر حسب غلظت کلی دو اتوکسی اتانول78شکل4-39- منحنی غلظت دو اتوکسی اتانول در فاز گاز بر حسب غلظت کلی دو اتوکسی اتانول79شکل4-40– منحنی افت فشار بر حسب طول مغزه برای حالت ناپایای حالت تک فازی و جریان محوری81شکل4-41- مقایسه ی اثر تعداد نقاط محاسبه ای برای افت فشار کل با زمان در سیستم تک فازی و جریان محوری82شکل4-42- منحنی فشار بر حسب طول مغزه83شکل4-43- منحنی اشباع گاز بر حسب طول مغزه83شکل4-44- منحنی فشار بر حسب طول مغزه (N=6)83شکل4-45- منحنی اشباع گاز بر حسب طول مغزه (N=6)83شکل4-46- اثر تعداد نقاط محاسبه ای بر روی منحنی اشباع مایع و فشار بر حسب طول85شکل4-47- بررسی اثر تعداد نقاط محاسبه ای بر منحنی غلظت متان در مخلوط بر حسب طول در دو فاز گاز و مایع85شکل4-48- منحنی قابلیت های حرکت فاز ها بر حسب طول86شکل4-49- بررسی همراه بودن دو متوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز87شکل4-50- بررسی همراه بودن دواتوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز88شکل4-51- بررسی همراه بودن دو بوتوکسی اتانول همراه با جریان گازمیعانی برروی فشار و اشباع گاز88شکل4-52- منحنی فشار بر حسب طول برای بازه ی 5 تا 1500 ثانیه. (برای بازه های 50 ثانیه ای)90شکل4-53- منحنی فشار بر حسب زمان برای چند مقطع مغزه90شکل4-54- اشباع گاز بر حسب طول برای زمان 5 تا 1500 ثانیه (در بازه های 100 ثانیه ای)91شکل4-55- اشباع گاز بر حسب زمان برای مقاطع مختلف از مغزه91شکل4-56- غلظت متان در فاز مایع بر حسب طول از زمان 5 ثانیه تا 1500 ثانیه ( در بازه های 100 ثانیه ای)92شکل4-57- غلظت متان بر حسب زمان برای مقاطع مختلف92شکل4-58- غلظت متان در فاز گاز بر حسب طول برای زمان های بین 5 ثانیه تا 1500 ثانیه (در بازه های زمانی 100 ثانیه)93شکل4-59- غلظت متان در فاز گاز بر حسب زمان برای چهار مقطع93شکل4-60- منحنی افت فشار برای جریان تک فازی بر حسب زمان95شکل4-61- منحنی های افت فشار و اشباع ورودی سنگ برای جریان دو فازی تعریف شده95شکل4-62- منحنی های تراوایی نسبی دو فاز گاز و مایع برای جریان دو فازی مورد توجه96شکل4-63- مقایسه ی افت فشار برای جریان تک فازی و دو فازی97شکل4-64- منحنی افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دومتوکسی اتانول98شکل4-65- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو متوکسی اتانول98شکل4-66- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دواتوکسی اتانول99شکل4-67- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو اتوکسی اتانول100شکل4-68- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دوبوتوکسی اتانول100شکل4-69- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو بوتوکسی اتانول101شکل4-70- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو متوکسی اتانول102شکل4-71- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو متوکسی اتانول102شکل4-72- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو اتوکسی اتانول103شکل4-73- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو اتوکسی اتانول103شکل4-74- افت فشار و اشباع مایع ورودی بر حسب زمان برای جریان دو فازی قبل و بعد از تزریق دو بوتوکسی اتانول104شکل4-75- منحنی تغییرات تراوایی نسبی گاز و مایع متوسط مغزه با زمان برای سه مرحله ی فرآیند در تزریق دو بوتوکسی اتانول104شکل5-1- شمای کلی از دستگاه سیلاب زنی108شکل5-2- منحنی رفتار فازی مخلوط متان و هپتان نرمال (85 % متان)111شکل5-3- منحنی تراوایی مطلق بر حسب زمان114شکل5-4- منحنی تراوایی نسبی بر حسب زمان قبل و بعد از تزریق دو متوکسی اتانول114شکل5-5- منحنی تراوایی نسبی بر حسب زمان قبل و بعد از تزریق دو اتوکسی اتانول115فهرست نشانه های اختصاری aترم انرژی در معادله ی حالتzکسر مولی کلیbترم حجمی درمعادله ی حالتJماتریس مشتقاتعلایم یونانیkتراوایی∆اختلافLطول مغزه∆ABشدت همبستگیMwجرم مولیδ(x)دلتای دیراکncتعداد اجزای مخلوطμگرانرویPcفشار بحرانیρچگالیRثابت جهانی گاز هاφتخلخل سنگSاشباعضریب فوگاسیته جزء در مخلوطTدماTcدمای بحرانیزیر نویس هاtزمانiشمارندهvحجم مولیjشمارندهwتابع آزمونkشمارندهxکسر مولی در فاز مایعAسایت Ayکسر مولی در فاز گازgگازZضریب تراکم پذیریlمایع مقدمه کشور ایران از نظر منابع گازی در جهان جایگاه خوبی قرار دارد. این کشور دارای ذخیره ی 981 تريليون مترمكعب به صورت در جا و تولید روزانه ی 9/111 بيليون متر مكعب است. با توجه به اینکه در ناحیه جنوبی کشور بسیاری از این مخازن از نوع گاز میعانی است، تولید با بهره وری بهتر لازم به نظر می رسد.از آنجایی که میادین گاز میعانی در اقتصاد کشور نقش بسزایی دارد، بنابراین پژوهش بر روی عملکرد این منابع ارزشمند جنبه ی حیاتی دارد. این مهم در داخل کشور چنان مورد توجه واقع شده است که در شرکت ملی نفت ایران زیر نظر مدیریت پژوهشی و فناوری آن شرکت واحدی با عنوان گروه پژوهش مخازن گاز میعانی در پژوهشکده ی ازدیاد برداشت ایجاد شده است.مخازن گازی به سه دسته تقسیم می شوند. درصد بيشتري از مخازن گازي، جزء مخازن گاز خشكهستند. این مخازن ضمن توليد گاز در سر چاههيچگونه ميعاناتي با خود توليد نمي کنند.دسته دوم مخازن گازي، مخازن گاز مرطوبهستند. در اين مخازنمقداري ميعانات ضمن توليد در سر چاه بدست مي آيد که معمولادر هر چاه بازاء توليد 1 ميليون فوت مكعبگاز در شرايط استاندارد، بین 3 تا20 بشكه ميعانات در سطح توليد مي گردد. این میعانات به علت سبکی می توانند در تولید محصولات متنوعی نظيرحلال ها، بنزين و خوراك واحد هاي پتروشيمي استفاده شوند.دسته سوم مخازنگازي، مخازن گاز ميعاني هستند. در این میادین میزان تولید میعانات در بیشترین حد قرار دارد. در مواردي حتي ممكن است بازاء توليد 1 ميليون فوت مكعب گاز در شرايطاستاندارد، 300 بشكه از اين ميعانات در سطح توليد گردد. تفاوت مهم ديگر اين مخازنبا مخازن گاز مرطوب اين است كه در مخازن گاز ميعاني بر خلاف مخازن گاز مرطوب، ممكناست مقدار زيادي از ميعانات در درون مخزن تشكيل گردد. دو نقطه ضعف اساسي تشكيل اينميعانات در درون مخزن به قرار زير است:1- تشكيل اين ميعانات ارزشمند در درونمخزن باعث مي شود تا به روش هاي متداول بازيافت نتوان آنها را توليد نمود و به اينترتيب بخش مهمي از سرمايه هيدروكربني غير قابل برداشت مي گردد. به طور كلي فضايدرون مخزن شبيه محيط اسفنجي مي باشد كه خلل و فرج آن توسط رابط هاي مويينه بهيكديگر وصل شده اند و براي شروع حركت يك سيال در اين محيط بسيار متراكم اسفنجي لازماست تا سيال به يك ميزان مشخص در آن محيط تشكيل شود و درصد معيني از حجم محيط رااشغال نمايد. اين بدان معناست كه پس از توليد اولين قطرات مايع در مخزن ميعاناتبدون حركت مانده و تا رسيدن به درصد اشباع معين در درون مخزن باقي خواهندماند.2- تشكيل اين ميعانات در درون محيط مخزن باعث انسداد مسير درون حفرهها ولولههاي مويينه محيط متخلخل مخزن شده و يك مقاومت اضافي در برابر حركت جريان گازاز درون مخزن به سمت چاه و نهايتاً توليد گاز از مخزن ايجاد ميكند و بهرهبردارياز مخزن را كاهش ميدهد.در فرا سوی مرز های ایران اسلامی، مراکز تحقیقاتی کشور های دیگر هم به اهمیت این مهم پی برده اند و سالانه مقالات و گزارشهای فراوانی را در مجلات و همایش های معتبر به چاپ می رسانند و ارائه می دهند.در بررسی مخازن هیدروکربنی زیر زمینی تا زمانی که با نفت سر و کار داریم، مطابق با انتظار مان با کاهش فشار در دمای ثابت با تشکیل گاز روبرو هستیم. اما این وضعیت برای مخازن گازی تر[1] و گاز میعانی[2] جور دیگری است. با کاهش فشار در دمای ثابت با تشکیل مایع برخورد می کنیم. این پدیده را میعان معکوس[3] نام می گذارند. با وجود چنین پدیده ای در سنگ مخزنی، با برداشت از آن به مرور زمان با افزایش مقدار مایع در حفره های مخزن روبرو خواهیم بود. چنین پدیده ای با گذشت زمان باعث ایجاد مانعی بر سر راه تولید گاز از چاه های آن منطقه می گردد. چنین پدیده ای را انسداد میعانی[4] می گویند.روش های مقابله با این پدیده تنوع زیادی دارد. در این پایان نامه هدف استفاده از حلال های مایع گلیکول و گلیکول اتری است. حلال از راه چاه تولید به مخزن تزریق می شود. حوزه ی چاه مدتی در همان شرایط می ماند تا حلال با میعانات و گاز داخل حفره های سنگ مخزن به تعادل برسد. با شروع دوباره ی تولید ابتدا میعانات به بیرون هدایت می شوند و گاز در این شرایط قابلیت حرکت بیشتری دارد و تولید آن افزایش می یابد. حلال های مورد نظر تاثیری بر روی سنگ مخزنی ندارند و تنها اثر آنها بر تعادل فازی است. مقداری از حلال که در حفره ها باقی می ماند، به میعانات اجازه ی ورود به فاز مایع را نمی دهد. بنابراین تا مدتی که حلال تبخیر نشود. سطح مایع پایین می ماند و استحصال گاز با افت فشار کمتری روبرو است.هدف از این پایان نامه مدل سازی اثر حلال های گلیکول و گلیکول اتری بر رفتار فازی سیال مخزنی در شرایط مخزن در مقیاس مغزه همراه با حل معادلات جریان دو فازی در مختصات استوانه ای در جهت شعاعی و مقایسه اثر آن در مقیاس آزمایشگاهی است. نوآوری های کار استفاده از یک روش عددی در حل معادلات مخزن و استفاده از یک معادله ی حالت موثر برای بررسی رفتار فازی حلال ها در کنار سیال هیدروکربنی است.فصل دوم مروری بر کار های گذشته در این فصل خلاصه ی مقالاتی که به استفاده ی حلال های مختلف در برطرف کردن پدیده ی انسداد میعانی پرداخته است ارائه می شود که در چند قسمت آورده شده است. 2-1- کاهش تولید در مخازن گاز میعانی افیدیک و همکارانش (1994) [1] به بررسی کاهش تولید در مخزن گاز میعانی آرون در اثر انباشتگی مععانات پرداختند و به این نتیجه رسیدند که با کاهش فشار به زیر فشار نقطه ی شبنم تولید چاه با فاکتوری حدود 2 کاهش می یابد که این در اثر تجمع میعانات در دهانه چاه است.بارنوم و همکارانش(1995) [2] به این نتیجه رسیدند که کاهش تولید برای مخازن با تولید کم با Kh کمتر از 1000 md-ft شدید است. 2-2- روش های برطرف کردن انسداد میعانی روش های متفاوتی برای مقابله با پدیده ی انسداد میعانی وجود دارد. دو روش عمده برای این برطرف کردن این پدیده به قرار زیر است:- تغییر رفتار فازی گاز میعانی- کاهش افت فشار و نگه داشتن فشار بالای نقطه ی شبنمیکی از روش ها استفاده از حلال ها است.حلال های مورد نظر از خانواده ی الکل ها، گیلکول ها و گلیکول اتر ها است که در زیر به شرح کارهایی که در این زمینه شده است می پردازیم:دو و همکارانش(2000)[3] سعی کردند که کاهش نفوذ پذیری نسبی توسط میعانات را با استفاده از یک حلال که آن حلال متانول است بهبود بخشند.دو و همکارانش احتمال می دانند که با تزریق متانول کشش بین سطحی گاز و مایع کم می شود ولی نتایج اندازه گیری این فرض را رد کرد. احتمال دیگری که مورد توجه آن ها بود اثر متانول روی رفتار فازی هیدروکربن بود که با کمک معادله حالت این فرض را مورد آزمایش قرار داد و مشاهده کردند که اگر غلظت متانول زیاد باشد در شرایط عبوری از مغزه یک فاز مایع غنی از متانول تشکیل می شود که مانع تشکیل فاز میعانی دیگری می گردد که باعث بهبود تراوایی نسبی می شود.نتایجی که از تزریق متانول بدست می آید به قرار زیر است:1- تزریق متانول به درون مغزه که نفوذ پذیری آن توسط انسداد میعانی کم شده باعث می شود که نفوذ پذیری نسبی بهتر شود.2- دلیل افزایش نفوذ پذیری نسبی این است که متانول باعث جابجایی فاز آب و میعانات می شود که این عمل توسط حل شدن متانول صورت می گیرد .3- با تزریق متانول آب کاملا حذف می شود.4- استفاده از متانول ودیگر حلال های الکلی روش های ارزان تری نسبت دیگر روش ها مانند شکست هیدرولیکی می باشند