در این نوشتار با استفاده از مدل نفت سیاه، یک روش عددی برای شبیه سازی مخازن هیدروکربنی سه بعدی ارائه شده است. در این راستا، سیال مخزن به صورت سه فازی- سه جزئی مدل سازی می شود، از اثر تغییرات دما صرف نظر می شود، تراکم پذیری سنگ مخزن و نیز اثر فشار موئینگی لحاظ می شود. معادلات حاکم بر جریان برای هر دو حالت اشباع و زیر اشباع ارائه شده و اثر جملات مختلف آن تبیین می شود، سپس روش حل عددی اختلاف محدود که بر روی یک شبکه یا سازمان پیاده سازی شده به طور کامل تشریح می شود. روش عددی مورد استفاده از اختلاف مرکزی در مکان و روش اویلر پسرو در زمان استفاده می کند. همچنین روش خطی سازی معادلات و الگوریتم حل عددی ارائه می شود. برای ارزیابی روش حل پیشنهادی، دو مساله نمونه مربوط به مخازن واقعی مشهور به SPE-1 و SPE-8 حل، و نتایج شبیه سازی با داده های موجود در منابع علمی مقایسه می شود. نتایج عددی حاصل حاکی از تطابق خوب این نتایج با نتایج حاصل از سایر شبیه سازها است.