مدلسازی سه بعدی گونه های سنگی با استفاده از ادغام داده های مغزه، نگارهای چاه و لرزه ای در یکی ازمخازن کربناته جنوب غربی ایران
الموضوعات :
مهدی خیرالهی
1
,
گلناز جوزانی کهن
2
,
رضا محبیان
3
,
علی مرادزاده
4
1 - دانشگاه تهران
2 - دانشگاه تهران
3 - دانشگاه تهران
4 - دانشگاه تهران
الکلمات المفتاحية: تعیین گونه های سنگی, خواص مخزنی, شاخص منطقه جریان, داده های مغزه,
ملخص المقالة :
تعیین گونه های سنگی، فرآیند اختصاص خواص مخزن به رخساره های زمینشناسی است و یک گونه سنگی تعیین شده دارای خواص زمین شناسی و مخزنی مشابه و حدوداً یکسان است. امروزه روش های بسیار مختلفی برای تعیین گونه های سنگی پیشنهاد شده و توسعه یافته است. از آن جمله می توان به روش رخساره حفرات و روش لوسیا اشاره کرد که از ویژگی های سنگ شناسی در کنار تخلخل و تراوایی استفاده می کند. از مهمترین این روش ها، تعیین گونه های سنگی به روش شاخص منطقه جریان می باشد. همچنین، معادلات مختلفی مانند وینلند-پیتمن و تابع جی-لورت در این زمینه وجود دارد که از روش های بسیار معتبر در تعیین گونه های سنگی می باشد. در این مطالعه با طراحی یک مسیر کار جدید و با استفاده از مهم ترین و کاراترین روش های تعیین گونه های سنگی اعم از روش های شاخص منطقه جریان، تفکیک تخلخل، شاخص وینلند-پیتمن، شاخص تغییر یافته منطقه جریان و لگاریتم های شاخص منطقه جریان و شاخص تغییر یافته منطقه جریان و اعمال روش های مذکور بر روی داده های مغزه، نگار چاه و لرزه ای، به تعیین گونه های سنگی در سه بعد و در کل میدان مورد بررسی پرداخته ایم. پس از اعتبار سنجی نهایی، به ترتیب، مقادیر ضرایب همبستگی بیش از 83%، 57%، 70%، 70%، 73% و 78% برای روش های مورد استفاده به دست آمد که بیانگر دقت و کارایی بیشتر روش شاخص منطقه جریان جهت تفکیک گونه های سنگی در این میدان می باشد. همچنین مقایسه جزئی نتایج اعتبارسنجی انجام شده پس از به کار گیری هر روش نیز تاییدکننده بیشترین دقت برای روش مبتنی بر شاخص منطقه جریان است. در نتیجه، با ادغام این روش با نشانگرهای لرزه ای، گونه های سنگی را در کل میدان و در سه بعد تفکیک کرده که چهار گونه سنگی مجزا تعیین گردید.
A. Kadkhodaie and R. Kadkhodaie, “A Review of Reservoir Rock Typing Methods in Carbonate Reservoirs: Relation between Geological, Seismic, and Reservoir Rock Types,” Pet. Eng. Iran. J. Oil Gas Sci. Technol., vol. 7, no. 4, pp. 13–35, 2018, [Online]. Available: http://ijogst.put.ac.ir.
N. Bize-forest, V. Baines, A. Boyd, A. Moss, and R. Oliveira, “Carbonate Reservoir Rock Typing and the Link Between Routine Core Analysis and Special Core Analysis,” Sca, pp. 1–6, 2014.
J. S. Gomes, M. T. Ribeiro, C. J. Strohmenger, S. Negahban, and M. Z. Kalam, “Carbonate reservoir rock typing - The link between geology and SCAL,” Soc. Pet. Eng. - 13th Abu Dhabi Int. Pet. Exhib. Conf. ADIPEC 2008, vol. 3, pp. 1643–1656, 2008, doi: 10.2118/118284-ms.
E. A. Clerke, D. Ph, and S. Aramco, “SPE 126086 Electrofacies and Geological Facies for Petrophysical Rock Typing : Khuff C,” no. May, pp. 9–11, 2009.
T. O. Services, T. O. Services, T. O. Services, and T. O. Services, “SPE 163294 Rock Typing and Characterization of Carbonate Reservoirs : A Case Study from South East Kuwait,” no. December, pp. 10–12, 2012.
P. Mukherjee, D. Singharay, S. Matar, and D. M. A. Meshari, “Rock-Typing : An Integrated Reservoir Characterization Tool for Tight Jurassic Carbonates , West Kuwait *,” vol. 70372, 2018, doi: 10.1306/70372Mukherjee2018.
Y. Liu et al., “Petrophysical static rock typing for carbonate reservoirs based on mercury injection capillary pressure curves using principal component analysis,” J. Pet. Sci. Eng., vol. 181, no. January, p. 106175, 2019, doi: 10.1016/j.petrol.2019.06.039.
S. N. A. Al-Jawad, M. A. Ahmed, and A. H. Saleh, “Integrated reservoir characterization and quality analysis of the carbonate rock types, case study, southern Iraq,” J. Pet. Explor. Prod. Technol., vol. 10, no. 8, pp. 3157–3177, 2020, doi: 10.1007/s13202-020-00982-6.
M. Skalinski, J. A. M. Kenter, M. Skalinski, and J. A. M. Kenter, “Geological Society , London , Special Publications Online First Carbonate petrophysical rock typing : integrating geological attributes and petrophysical properties while linking with dynamic behaviour Carbonate petrophysical rock typing : integrating geo,” 2014, doi: 10.1144/SP406.6.
L. Saputelli, R. Celma, D. Boyd, H. Shebl, J. Gomes, and F. Bahrini, “SPE-196704-MS Deriving Permeability and Reservoir Rock Typing Supported with Self-Organized Maps SOM and Artificial Neural Networks ANN - Optimal Workflow for Enabling Core-Log Integration Predictive Data Analytics ( PDA ) and Machine Learning ( ML ),” 2019.
M. Shabaninejad, M. Bagheripour, S. Z. Oil, and P. Company, “SPE 150819 Rock Typing and Generalization of Permeability- Porosity Relationship for an Iranian Carbonate Gas Reservoir,” 2011.
E. Aliakbardoust, “Integration of rock typing methods for carbonate reservoir characterization,” vol. 055004, doi: 10.1088/1742-2132/10/5/055004.
Z. Riazi, “Journal of Petroleum Science and Engineering Application of integrated rock typing and fl ow units identi fi cation methods for an Iranian carbonate reservoir,” J. Pet. Sci. Eng., vol. 160, no. October 2017, pp. 483–497, 2018, doi: 10.1016/j.petrol.2017.10.025.
S. Hosseinzadeh, A. Kadkhodaie, and S. Yarmohammadi, “NMR derived capillary pressure and relative permeability curves as an aid in rock typing of carbonate reservoirs,” J. Pet. Sci. Eng., vol. 184, p. 106593, 2020, doi: 10.1016/j.petrol.2019.106593.
R. Mohebian, M. A. Riahi, and A. Kadkhodaie, “Characterization of hydraulic flow units from seismic attributes and well data based on a new fuzzy procedure using ANFIS and FCM algorithms, example from an Iranian carbonate reservoir,” Carbonates and Evaporites, vol. 34, no. 2, pp. 349–358, 2019, doi: 10.1007/s13146-017-0393-y.
M. Farshi, R. Moussavi-Harami, A. Mahboubi, M. Khanehbad, and T. Golafshani, “Reservoir rock typing using integrating geological and petrophysical properties for the Asmari Formation in the Gachsaran oil field, Zagros basin,” J. Pet. Sci. Eng., vol. 176, no. May 2018, pp. 161–171, 2019, doi: 10.1016/j.petrol.2018.12.068.
A. Mirzaei-Paiaman, M. Ostadhassan, R. Rezaee, H. Saboorian-Jooybari, and Z. Chen, “A new approach in petrophysical rock typing,” J. Pet. Sci. Eng., vol. 166, pp. 445–464, 2018, doi: 10.1016/j.petrol.2018.03.075.
A. Mirzaei-paiaman et al., “Journal of Petroleum Science and Engineering A further veri fi cation of FZI * and PSRTI : Newly developed petrophysical rock typing indices,” vol. 175, no. December 2018, pp. 693–705, 2019, doi: 10.1016/j.petrol.2019.01.014.
H. Mehrabi, R. R. Karami, and M. R. Nejad, “Reservoir rock typing and zonation in sequence stratigraphic framework of the Cretaceous Dariyan Formation , Persian Gulf,” Carbonates and Evaporites, no. 0123456789, 2019, doi: 10.1007/s13146-019-00530-2.
B. Khadem, M. R. Saberi, M. Eslahati, and B. Arbab, “Integration of rock physics and seismic inversion for rock typing and flow unit analysis: A case study,” Geophys. Prospect., vol. 68, no. 5, pp. 1613–1632, 2020, doi: 10.1111/1365-2478.12952.
نشریه علمی– پژوهشی زمین شناسی نفت ایران سال یازدهم، شماره 21، بهار و تابستان 1400 ص1-15
Iranian Journal of Petroleum Geology No. 21, Spring & Summer, 2021, pp. 1-15
Dor: 20.1001.1.22518738.1400.11.22.1.1
مدلسازی سه بعدی گونههای سنگی با استفاده از ادغام دادههای مغزه، نگار و لرزهای، مطالعه موردی یکی ازمخازن کربناته جنوب ایران
مهدی خیرالهی1، گلناز جوزانی کهن2، رضا محبیان3*، علی مرادزاده4
1- کارشناسی ارشد مهندسی اکتشاف نفت، دانشکده مهندسی معدن، دانشگاه تهران، تهران، ایران
2، 3، 4- استادیار، دانشکده مهندسی معدن، دانشکدگان فنی دانشگاه تهران، دانشگاه تهران، تهران، ایران
دريافت تیر 1401، پذيرش آبان 1401
تعیین گونههای سنگی1، فرآیند اختصاص خواص مخزن به رخسارههای زمینشناسی است و یک گونه سنگی تعیین شده دارای خواص زمینشناسی و مخزنی مشابه و حدوداً یکسان است. امروزه روشهای بسیار مختلفی برای تعیین گونههای سنگی پیشنهاد شده و توسعه یافته است. برخی از این روشها عمدتاً بر جنبههای زمینشناسی (پتروگرافی) تأکید دارد که از آن جمله میتوان به روش رخساره حفرات2 و روش لوسیا اشاره کرد که از ویژگیهای سنگشناسی در کنار تخلخل و تراوایی استفاده میکند. روشهای پتروفیزیکی تعیین گونههای سنگی غالباً بر اساس دادههای تخلخل و تراوایی (K/Ø) و فشار موئینگی (PC) است. از سادهترین این روشها، نمودار تخلخل در مقابل تراوایی و تعیین حد برش3 است. از مهمترین و کاربردیترین این روشها، تعیین گونههای سنگی به روش شاخص منطقه جریان4 میباشد. همچنین معادلات مختلفی مانند وینلند-پیتمن (R35) و تابع جی-لورت در این زمینه وجود دارد که از روشهای بسیار معتبر در تعیین گونههای سنگی میباشد. در این مطالعه، با استفاده از مهمترین و کاراترین روشهای تعیین گونههای سنگی اعم از روشهای مبتنی بر تخلخل، شاخص منطقه جریان، شاخص تغییر یافته منطقه جریان، شاخص وینلند-پیتمن و لگاریتمهای آنها و اعمال روشهای مذکور بر روی دادههای مغزه، نگار چاه و لرزهای، به تعیین گونههای سنگی در سه بعد در کل میدان مورد بررسی پرداختهایم و با مقایسه ضرایب همبستگی پی بردیم که روش شاخص منطقه جریان دارای بیشترین دقت و کارایی میباشد.
واژه‌هاي كليدي: تعیین گونههای سنگی، خواص مخزنی، شاخص منطقه جریان، دادههای مغزه
1- مقدمه
مخازن هیدروکربوری در محیطهای رسوبی متفاوتی تشکیل شده و تحت فرآیندهای دیاژنزی متفاوتی قرار گرفتهاند. همچنین به دلیل عواملی مانند ناهمگونی مخازن و فرآیندهای دیاژنتیکی متفاوت، تعیین پارامترهای مخزنی و یافتن رابطه بین تخلخل و تراوایی، امری دشوار میباشد [1]. لذا، با استفاده از روشهای گوناگون که در ادامه به بررسی و بیان هر کدام پرداخته خواهد شد، سعی بر دستهبندی این خواص پتروفیزیکی در گروههایی مشابه، تحت عنوان گونههای سنگی، میگردد تا بتوان رفتار مخزن را مورد تحلیل و بررسی قرار داد. معیارهایی برای تعیین گونههای سنگی وجود دارد که از مهمترین آنها میتوان، محیط رسوبی و دیاژنزی مشابه، توزیع یکسان اندازه دهانه حفرات5، تراوایی نسبی و فشار موئینگی یکسان در یک نوع ترشوندگی و رابطه تخلخل-تراوایی مشابه را نام برد.
تعیین گونههای سنگی مخزنی فرآیندی است که به واسطه آن، رخسارههای زمینشناسی بر اساس رفتارهای دینامیکی، بررسی و دستهبندی میشوند. گونههای سنگی مخزنی از عوامل مهم در تعیین خصوصیات مخزن بوده و تعیین آنها یکی از چالش برانگیزترین امور در مخازن کربناته میباشد. عواملی همچون فرآیندهای دیاژنتیکی و رفتارهای بین سنگ و سیال باعث پیچیدگی ارتباط میان رخسارههای زمینشناسی و گونههای سنگی مخزنی میگردد. لذا شناخت درست محیط رسوبی، فرآیندهای دیاژنزی و رفتار سنگ-سیال با استفاده از بررسیهای آنالیز ویژه مغزه6، میتواند بسیار کمککننده باشد [2]. در نتیجه، پژوهشگران زیادی به انجام این فرآیند و تعیین گونههای سنگی در مخازن مختلف پرداختهاند. به عنوان مثال، برای تعیین یکپارچه گونههای سنگی مخزنی، باید روندهایی جهت استخراج خصوصیات دینامیکی مخزن از دادههای زیرسطحی مثل زمینشناسی، لرزهای و در نتیجه دادههای پتروفیزیکی به درستی طی شوند. لذا، کلرک و همکاران در سال 2009، برای بالابردن دقت مسیر کار، از یک مجموعه داده چندرشتهای7 در مخزن گازی خوفc در عربستان سعودی استفاده نمودند [3].
مخرجی و همکاران نیز در سال 2018، با ترکیب و ادغام دادههای مختلف از جمله واحدهای سنگچینهشناسی8 حاصل از نگار، آنالیزهای رسوبشناسی9 بر روی مغزه، بررسیهای کمی دادههای نگار و بهدست آوردن کانیشناسی، محتوای سیال و اشباع نمونهها، تعیین رخسارههای الکتریکی حاصل شده از دادههای نگار و مغزه و در نهایت، آنالیزهای ویژه مغزه، به تعیین گونههای سنگی موجود در یک مخزن کربناته ژوراسیک در غرب کویت پرداختند [4]. همچنین در سال 2020، الجواد و همکاران نیز با استفاده و ادغام روشهای موجود، به تعیین گونههای سنگی و بررسی کیفیت مخزنی آنها در یکی از مخازن کربناته جنوب عراق پرداختند [5].
در سال 2014، اسکالینسکی و کنتر با تعیین گونههای سنگی ترکیبی و بر اساس روشهای متداول موجود، نتایج مطلوبی را در دو میدان کربناته بزرگ مورد مطالعه بهدست آوردند. تفاوت کار آنها، چگونگی ادغام دادهها و روشها با یکدیگر بود و در نهایت با اعتبارسنجی نتایج به وسیله مدل سهبعدی زمینشناسی، به صحت کار انجام شده پی برده شد [6]. در ادامه، ساپوتلی و همکاران در سال 2019 با استفاده از نقشههای خودسازمانده10 و شبکههای عصبی مصنوعی11 یک روند انجام کار بهینه جهت ادغام دادههای مغزه و نگار برای بهدست آوردن تراوایی و تعیین گونههای سنگی مخزنی در نظر گرفتند که توانایی پیشبینی خواص پتروفیزیکی مخزن به وسیله مدل مخزنی بهدست آمده را به خوبی افزایش داد [7].
شعباننژاد و همکاران در سال 2011 با استفاده از روشهای معمول تعیین گونههای سنگی مثل روش وینلند و تعیین واحدهای جریان با استفاده از نشانگرهای مناطق جریان، به بررسی یکی از مخازن کربناته گازی در ایران پرداختند [8]. همچنین، علیاکبردوست و همکاران در سال 2013 با بهکارگیری روشهای مختلفی همچون استفاده از دادههای SCAL (منحنیهای فشار موئینگی و تراوایی نسبی)، خواص زمینشناسی (ساخت رسوبی، نوع تخلخل و محتوای سیمان) و برخی پارامترهای پتروفیزیکی مخزن مثل تخلخل و تراوایی، به تعیین گونههای سنگی مخزن کربناته مورد مطالعه پرداختند [9]. حسینزاده و همکاران نیز در سال 2019 از نگار NMR و منحنیهای تراوایی نسبی و فشار موئینگی بهدست آمده از آن برای تعیین گونههای سنگی مخازن کربناته استفاده نمودند [10].
محبیان و همکاران در سال 2017 توانستند با بهکارگیری الگوریتمهایی مثل ANFIS و FCM در یک روند فازی12 جدید در کنار دادههای چاه و نشانگرهای لرزهای، تعیین مشخصات و دستهبندی مناسبی برای واحدهای جریان هیدرولیکی و گونههای سنگی در یکی از مخازن هیدروکربوری کربناته ایران انجام دهند [11]. سپس، فرشی و همکاران در سال 2018 به تعیین گونههای سنگی مخزنی با استفاده از ترکیب مطالعات میکروسکوپی (ریزرخساره ها و فرآیندهای دیاژنزی) و پتروفیزیکی (رخسارههای الکتریکی و وحدهای جریان هیدرولیکی) پرداختند [12].
در سال 2018، میرزایی و همکاران دست به ابداع روشی جدید برای تعیین گونههای سنگی پتروفیزیکی در مخازن کربناته زدند. آنها تعیین گونههای سنگی را به دو دسته کلی تعیین گونههای سنگی پتروفیزیکی استاتیک و دینامیک تقسیم نمودند [13]. همچنین در سال 2019، میرزایی و همکاران به ادامه مطالعه در این زمینه پرداخته و پارامترهای معرفی شده را بهبود و توسعه دادند [14]. در سال 2020 نیز، خادم و همکاران با ترکیب فیزیک سنگ و وارونسازی دادههای لرزهای در قالب یک روند انجام کار مدون، سعی بر تعیین گونههای سنگی و واحدهای جریان در یک میدان نفتی ماسهسنگی در حوضه خلیج فارس نمودند [15].
با توجه به موضوع ذکر شده و مطالعات پیشین صورت گرفته در این زمینه، به عدم کارایی مشخصههای حاصله از هر بخش در سایر بخشهای مطالعات مخزن پی میبریم. به عنوان مثال، نتایج حاصل از تعیین گونههای سنگی بر اساس دادههای زمینشناسی در بسیاری از مطالعات تطابق و همخوانی مناسبی با نتایج مطالعات مخزنی ندارد چرا که روشهای مورد استفاده هیچ همبستگی خاصی با یکدیگر ندارند. لذا در این پژوهش با طراحی گام به گام مراحل مطالعه به صورت زیر، سعی بر رفع ابهامات موجود و تعیین گونههای سنگی مخزن به صورتی که مورد اطمینان و قابل قبول در همهی حوزههای مطالعاتی باشد، مینماییم. در ابتدا و با به کار گیری دادههای مغزه و اعمال روشهای مورد نظر، گونههای سنگی را در مخزن مورد مطالعه مشخص مینماییم. سپس روشها را برای دادههای نگار چاه نیز اعمال نموده و گونههای سنگی را تفکیک میکنیم. در ادامه، با وارد کردن نتایج و دادههای لرزهای در فرآیند ادغام و اعمال حد تفکیک گونههای سنگی بهدست آمده از مراحل قبل، مدل سه بعدی گونههای سنگی را برای کل میدان مورد بررسی به دست می آوریم.
میدان نفتی مورد مطالعه در بخش جنوبی ایران و شمال غربی خلیج فارس واقع شده است. این میدان در سال ۱۳۳۹ به وسیله شرکت سابق سیریپ13 با حفر چاه شمارهی یک کشف گردید. تعداد پنج مخزن هیدروکربنی شامل بخش ماسه سنگ اهواز (غار)، آسماری کربناته، سروک، کژدمی و فهلیان در این میدان دیده میشود. چاه شمارهی ده عمیقترین چاه این میدان با عمق ۴۲۵۳ متر (۱۳۹۵۴ پا) در یال غربی تا سازند گوتینا حفر گردید. تعداد ۱۰ حلقه چاه در این میدان حفاری گردیده است.
ردیفهای رسوبی موجود در میدان نفتی مورد مطالعه از بالا به پائین به قرار زیر میباشد:
رسوبات عهد حاضر، سازندهای بختیاری و آغاجاری، سازند میشان، سازند گچساران، بخش ماسه سنگی اهواز، سازند آسماری، سازند جهرم، سازند پابده، سازند گورپی، سازند ایلام، بخش لافان، سازند سروک، سازند کژدمی، سازند داریان، سازند گدون، سازند فهلیان (ياماما)، سازند گرو و سازند گوتنیا
میدان نفتی مورد مطالعه مشتمل بر پنج مخزن هیدروکربنی بوده که به ترتیب، ازدیاد عمق شامل مخزن آسماری (دو بخش ماسه سنگی (غار) و کربناته)، سروک، کژدمی (نهر امر) و فهلیان (ياماما) میباشد. تولید نفت از مخازن آسماری و سروک صورت میگیرد. تولید نفت از مخزن فهلیان برای مدت محدودی انجام گرفته، ولی بدلیل افت شدید فشار تولید از این مخزن ادامه نیافته است.
سنگهای مخازن بخش ماسه سنگ اهواز (غار) و کژدمی (نهر امر)، عمدتا از جنس ماسه سنگ و بقیه مخازن از جنس سنگهای کربناته (دولومیت و آهک) میباشد. سنگ مخزن آسماری کربناته و سروک بیشتر از جنس دولومیت بوده و تخلخل زمينه14 و حفرهای میکروسکوپی15 نقش اساسی را در تولید ایفا میکند.
بیشترین تخلخل براساس تفاسیر پتروفیزیکی در ماسه سنگ اهواز (غار) دیده میشود و کربناتهای مخزن فهلیان (ياماما) دارای کمترین تخلخل میباشد. مطالعه نتایج آزمایشات بر روی مغزهها و نتایج تفاسیر پتروفیزیکی نشان میدهد که مقادیر تخلخل و درصد اشباع آب در سنگهای مخازن بسیار متغیر است.
مساحت میدان تحت بررسی 35 کیلومتر مربع میباشد (به عرض 5 و به طول ۷ کیلومتر). این میدان مانند بیشتر میادین حوزه خلیج فارس ازیک تاقدیس نامتقارن که احتمالا از روند ساختمانهای زاگرس تبعیت میکند تشکیل شده است . در این میدان جمعا ۱۰ حلقه چاه حفر شده است که اکثرا از مخزن آسماری بهرهبرداری میشود. چاه شمارهی ۱۰ این میدان به دلیل تولید گاز زیاد، متروک گردیده است. این میدان در دو افق دارای نفت اقتصادی میباشد. یکی مخزن غار-آسماری و دیگری مخزن سروک. بر خلاف نواحی خشکی در حوضه زاگرس طبق مطالعات انجام شده، سازند آسماری در این حوضه از تراوائی و تخلخل اولیه نسبتا خوبی برخوردار است.
3- مواد و روشها
در این مطالعه، علاوه بر دادههای حاصل از آنالیزهای معمول و ویژه مغزه، از دادههای نگار چاهها نیز در عمق مخزنی میدان مورد مطالعه جهت تعیین و تفکیک گونههای سنگی استفاده شده است. علاوه بر دادههای مذکور، دادههای حاصل از برداشت لرزهای در این میدان در دسترس بوده و جهت پیشبینی گونههای سنگی مختلف در سه بعد مورد استفاده قرار گرفته است. همچنین از روشهای رایج و کارآمد تعیین گونههای سنگی از جمله تفکیک بر اساس تخلخل، شاخص منطقه جریان و شاخص وینلند-پیتمن، شاخص تغییر یافته منطقه جریان، لگاریتم شاخص منطقه جریان و لگاریتم شاخص تغییر یافته منطقه جریان در این پژوهش استفاده شده است. پس از به کار گیری روشهای ذکر شده و تعیین گونههای سنگی در محل چاهها، با اعمال وارونسازی و ادغام، نتایج را در سه بعد و در کل میدان تخمین میزنیم. این مسیر و روند کار پیشنهادی در این مطالعه همان گونه که در شکل 1 نشان داده شده است، نتایجی حاصل مینماید که در بررسیهای زمینشناسی، پتروفیزیکی و مخزنی میدان، مورد قبول و دارای همخوانی مناسب میباشد.
1-3- روش تفکیک بر اساس تخلخل
جهت اعمال روش مبتنی بر دادههای تخلخل، چه بر اساس دادههای حاصل از آنالیز مغزه و چه دادههای نگار چاه، ابتدا به ترسیم نمودار تجمعی دادهها میپردازیم. سپس با مشخص نمودن نقاط شکستگی و تغییر شیب در این نمودار، به حدود برش میان گونههای سنگی مختلف پی خواهیم برد. لازم به ذکر است که روش آماری تعیین حدود برش بر اساس تغییر شیب نمودار تجمعی برای تمامی پارامترها در روشهای مورد استفاده به کار گرفته میشود.
2-3- شاخص منطقه جریان
شاخص منطقه جریان معیاری جهت توصیف واحدهای جریان هیدرولیکی در مخزن میباشد که با توجه به شاخص کیفیت مخزنی16 و تخلخل موثر17 و به صورت زیر محاسبه میشود:
رابطه1: |
|
رابطه2: |
|
رابطه3: |
|
رابطه4: |
|
رابطه5: |
|
روش مورد استفاده | ضریب همبستگی در چاه شماره یک (%) | ضریب همبستگی در چاه شماره پنج (%) |
شاخص منطقه جریان | 54/75 | 23/63 |
شاخص وینلند-پیتمن | 05/57 | 27/61 |
شاخص تغییر یافته منطقه جریان | 46/33 | 28/43 |
لگاریتم شاخص منطقه جریان | 60/68 | 01/59 |
لگاریتم شاخص تغییر یافته منطقه جریان | 41/43 | 46/55 |
با مقایسه مقادیر ضرایب همبستگی مربوط به هر روش متوجه میشویم که کاراترین روش جهت تعیین گونههای سنگی بر اساس دادههای حاصل از آنالیز مغزه، روش شاخص منطقه جریان میباشد. از این رو، نتایج چگونگی تفکیک گونههای سنگی در چاههای شماره یک و پنج در شکل 2 نشان داده شده است.
شکل2-الف) شاخص منطقه جریان در چاه شماره یک، ب) شاخص منطقه جریان در چاه شماره پنج
2-4- تعیین گونههای سنگی بر اساس دادههای نگار در چاههای شماره یک، پنج، شش و ده
حال که تعیین گونههای سنگی بر اساس دادههای حاصل از آنالیز مغزههای موجود از چاههای شماره یک و پنج صورت پذیرفت، میتوان روشهای موجود و مورد استفاده جهت تعیین گونههای سنگی در این مطالعه را بر روی دادههای حاصل از نگار چاههای موجود در میدان مورد پژوهش، اعمال نمود. با توجه به مقادیر ضرایب همبستگی بهدست آمده، در مییابیم که بهترین و بادقتترین روش جهت تعیین گونههای سنگی در این میدان و بر اساس دادههای مغزه، روش شاخص منطقه جریان است. در ادامه و در جدول زیر، نتایج به کار گیری روشهای ذکر شده بر روی دادههای نگار چاههای مختلف به صورت ضرایب همبستگی هر کدام از روشها ذکر شده است. با توجه به این که با استفاده از دادههای نگار چاهها نیز روش شاخص منطقه جریان بهترین عملکرد را از خود نشان میدهد، شکل 3 نتایج استفاده از این روش را برای چاههای مختلف و چگونگی تعیین و تفکیک گونههای سنگی را نشان میدهد.
جدول2- ضرایب همبستگی مربوط به هر روش تعیین گونههای سنگی بر اساس دادههای نگار برای چاههای شماره یک، پنج، شش و ده
روش مورد استفاده | ضریب همبستگی در چاه شماره یک (%) | ضریب همبستگی در چاه شماره پنج (%) | ضریب همبستگی در چاه شماره شش (%) | ضریب همبستگی در چاه شماره ده (%) |
شاخص منطقه جریان | 75/90 | 29/89 | 16/87 | 94/91 |
شاخص وینلند-پیتمن | 12/57 | 15/76 | 72/40 | 56/79 |
شاخص تغییر یافته منطقه جریان | 60/39 | 10/56 | 50/31 | 09/66 |
لگاریتم شاخص منطقه جریان | 84/88 | 28/85 | 06/89 | 21/90 |
لگاریتم شاخص تغییر یافته منطقه جریان | 51/33 | 31/50 | 57/32 | 77/65 |
شکل3- الف) چهار گونه سنگی موجود در چاه شماره یک، ب) پنج گونه سنگی موجود در چاه شماره پنج، ج) چهار گونه سنگی موجود در چاه شماره شش، د) چهار گونه سنگی موجود در چاه شماره ده بر اساس روش شاخص منطقه جریان
پس از تعیین گونههای سنگی بر اساس دادههای مغزه و نگار چاهها، و در اولین گام به منظور تعیین گونههای سنگی موجود در کل میدان و در سه بعد، دادههایی از قبیل نگارهای چاهها، دادهی لرزهای و افقهای تفسیر شده موجود در میدان را به کار گرفته تا فرآیند وارون سازی پس از بر انبارش را انجام داده و مکعب سه بعدی امپدانس صوتی را برای میدان مورد مطالعه بهدست آوریم. هنگام ساخت مدل اولیه برای وارون سازی، از افقها برای درونیابی بین چاهها استفاده میکنیم و در ادامه نیز به یک موجک استخراج شدهی میانگین از چاهها نیاز داریم لذا موجک را استخراج کرده که فاز متوسط 97- درجه را نشان میدهد و در شکل زیر نشان داده شده است. پس از استفاده از این موجک میانگین استخراج شده، همبستگی متقاطع بین لرزه نگاشت مصنوعی و نگار مرکب به مقدار 71% میرسد.
شکل 4- موجک میانگین استخراج شده از چاهها
با استفاده از موجک میانگین استخراج شده و پس از درونیابی منحنیهای نگار چاه، یک فیلتر فرکانس پایینگذر اعمال میشود که به طور پیش فرض تمام فرکانسها را تا 10 هرتز عبور داده، تمام فرکانسهای بالای 15 هرتز را فیلتر کرده و فیلتر را بین این محدودهها درونیابی میکند. در نتیجه، مدل اولیه محاسبه شده به صورت شکل زیر نمایش داده میشود.
شکل5- مدل اولیه محاسبه شده با استفاده از موجک میانگین استخراج شده از تمامی چاههای موجود
شکل زیر نیز دادهی امپدانس صوتی را در کنار مدل اولیهی محاسبه شده و دادهی لرزهای برداشت شده در میدان مورد مطالعه نشان میدهد.
شکل6- الف) دادهی لرزهای میدان تحت بررسی، ب) امپدانس صوتی بهدست آمده، ج) مدل اولیه محاسبه شده
با توجه به همخوانی خوب و قابل قبول نتایج وارونسازی با دادههای لرزهای موجود در محل چاهها که در شکلهای فوق نشان داده شده است، به صحت و درستی فرآیند وارون سازی انجام شده پی میبریم. در نتیجه و پس از اتمام فرآیند وارون سازی، اکنون میتوانیم به انجام فرآیند ادغام جهت ایجاد و پیشبینی نگار ویژگیهای مورد نظر و در نهایت ایجاد آنها در سه بعد بپردازیم.
ادغام فرآیندی است که هدف آن ادغام نگارهای چاه و دادههای لرزهای است. هدف کلی این است که یک ویژگی نمودار چاه را با استفاده از نشانگرهای دادههای لرزهای پیشبینی کنیم. این ویژگی ممکن است هر نوع نگار اندازهگیری شده مانند سرعت یا تخلخل بوده، یا حتی ممکن است یک ویژگی سنگشناسی مشتق شده مانند حجم شیل و یا واحدهای جریان باشد. پس از بارگذاری دادههای مورد نیاز از جمله داده لرزهای، نتیجه وارونسازی (مکعب امپدانس صوتی)، چاهها و ...، برای انجام فرآیند ادغام آماده هستیم. این تحلیل در دو مرحله انجام میشود. در مرحله اول، مرحله آموزش19، نگار هدف و دادههای لرزهای را در مکانهای چاه تجزیه و تحلیل میشود تا یک رابطه آماری بین آنها بدست بیاید. در مرحله دوم، که همان مرحله کاربرد است، رابطه مشتق شده را به کل حجم اعمال میکنیم تا مقادیر نگار را در سراسر آن حجم ایجاد نماییم. در این پژوهش و با توجه به پارامترهای موردنظر جهت تعیین گونههای سنگی، به استخراج مکعب تخلخل، شاخص منطقه جریان، شاخص وینلند-پیتمن، شاخص تغییر یافته منطقه جریان، لگاریتم شاخص منطقه جریان و لگاریتم شاخص تغییر یافته منطقه جریان میپردازیم. با توجه به این که روش مبتنی بر شاخص منطقه جریان در بخشهای قبل نتایج بهتری از خود نشان داد، در ادامه، نتیجهی پیشبینی این شاخص در سه بعد به صورت شکل زیر نشان داده میشود.
شکل 7- شاخص منطقه جریان بهدست آمده در سه بعد در کنار دادهی لرزهای و نتیجهی وارون سازی
نمایش نهایی که با این دادهها ایجاد خواهیم کرد، یک برش داده از زمان20 مورد علاقه است. برش داده را در یک زمان 1500 میلی‌ثانیه که نزدیک به منطقه هدف است اعمال میکنیم. در حوالی آن زمان، نمونهها را در یک پنجره 10 میلی‌ثانیه به طور میانگین می‌گیریم. پس از تکمیل مراحل، برش دادهی ایجاد شده به صورت زیر نمایان میشود:
شکل 8- برش زمانی در زمان 1500 میلی‌ثانیه از نتیجهی شاخص منطقه جریان محاسبه شده در سه بعد و گونههای سنگی تفکیک شده
با توجه به شکل فوق که نشان دهنده پراکندگی شاخص منطقه جریان میباشد، میتوان دریافت که گونههای سنگی در کل مخزن و با استفاده از روش مبتنی بر شاخص منطقه جریان، به خوبی و به وسیلهی رنگهای مختلف، از هم تفکیک شدهاند به این صورت که رنگهای مشابه نشان دهنده شاخصهای منطقه جریان مربوط به یک گونه سنگی میباشند.
تمامی مراحل فوق برای سایر روشهای تعیین گونه سنگی مورد استفاده در این پژوهش انجام گردیده است و اعتبار سنجی نتایج هر کدام از روشها به صورت نمودار متقاطع در شکلهای زیر نمایش داده شده است.
شکل 9-نمودارهای متقاطع و مقادیر ضرایب همبستگی برای الف) روشهای شاخص منطقه جریان، ب) تفکیک تخلخل، ج) شاخص وینلند-پیتمن، د) شاخص تغییر یافته منطقه جریان، ه) لگاریتم شاخص منطقه جریان و و) لگاریتم شاخص تغییر یافته منطقه جریان
با توجه به نمودارهای متقاطع بهدست آمده برای هر روش که مقدار پیشبینی شده را در مقابل مقدار دادهی واقعی نشان میدهد، مقادیر ضرایب همبستگی برای روشهای شاخص منطقه جریان، تفکیک تخلخل، شاخص وینلند-پیتمن، شاخص تغییر یافته منطقه جریان، لگاریتم شاخص منطقه جریان و لگاریتم شاخص تغییر یافته منطقه جریان به ترتیب بیش از 83%، 57%، 70%، 70%، 73% و 78% میباشد. در نتیجه درمییابیم که دقیقترین روش برای تعیین گونههای سنگی در سه بعد و در کل میدان مورد مطالعه، روش شاخص منطقه جریان میباشد.
با توجه به نتایج بهدست آمده از اعمال روشهای گوناگون مطرح شده در این پژوهش جهت تعیین گونههای سنگی در میدان مورد مطالعه، نتایج زیر حاصل میشود:
-بهترین و کاراترین روش تعیین گونههای سنگی در میدان مورد مطالعه، روش شاخص منطقه جریان میباشد. این مسئله با توجه به میزان ضرایب همبستگی این روش در چاههای مختلف و با در نظر گرفتن دادههای مختلف مانند دادههای حاصل از آنالیز مغزهها و دادههای نگار چاه مشخص میگردد.
-به کار گیری روشهای پیچیدهتر مانند شاخص تغییر یافته منطقه جریان و همچنین اعمال عملگرهایی همچون لگاریتم جهت بالا بردن دقت و کیفیت تفکیک گونههای سنگی در این میدان موثر نبوده و میزان ضرایب همبستگی را نسبت به روش شاخص منطقه جریان افزایش نداد.
-پس از انجام فرآیند ادغام و بهدست آوردن مکعب پارامترهای مد نظر در سه بعد و در کل میدان مورد مطالعه همچون تخلخل، شاخص منطقه جریان، شاخص وینلند-پیتمن، شاخص تغییر یافته منطقه جریان و لگاریتمهای آنها، جهت مقایسه و بررسی دقت پیشبینی، نمودارهای متقاطع مقادیر تخمین زده شده هر پارامتر را در مقابل دادهی واقعی رسم نموده و در این سطح نیز به کارا بودن و داشتن بیشترین دقت روش شاخص منطقه جریان پی بردیم.
-با توجه به نتایج حاصل از اعمال روشهای تعیین گونههای سنگی در میدان مورد مطالعه و با استفاده از دادههای مغزه و نگار چاه در چاههای مختلف به این نتیجه رسیدیم که به طور معمول چهار الی پنج گونهی سنگی قابل تفکیک در این میدان وجود دارد.
-با توجه به شکل مکعب نهایی حاصل شده در کل میدان مورد مطالعه و بر اساس روش شاخص منطقه جریان میتوان به وجود حدودی چهار الی پنج گونهی سنگی مختلف که هر کدام با یک طیف رنگی منحصر به فرد در شکل نهایی نشان داده شده است، پی برد.
-نتایج حاصل از فرآیند ادغام که نشان دهندهی مکعب تخمین زده شده برای هر پارامتر مورد نظر و بر اساس روشهای گوناگون مورد استفاده در این مطالعه است، تطابق بسیار خوبی با نتایج حاصل از اعمال روشهای تعیین گونههای سنگی بر دادههای مغزه و نگار چاه دارد که نشان دهندهی درستی و دقت بالای نتایج حاصل از این پژوهش در سه بعد و در مقیاس بزرگ دارد.
سپاس و قدردانی
منابع
[1] KADKHODAIE, A., and KADKHODAIE, R., 2018, A Review of Reservoir Rock Typing Methods in Carbonate Reservoirs: Relation between Geological, Seismic, and Reservoir Rock Types, Pet. Eng. Iran. J. Oil Gas Sci. Technol., 7(4), 13–35.
[2] BIZE-FOREST, N., BAINES, V., BOYD, A., MOSS, A., and OLIVEIRA, R., 2014, Carbonate Reservoir Rock Typing and the Link between Routine Core Analysis and Special Core Analysis, International Symposium of the Society of Core Analysts, Avignon, France, 8-11.
[3] CLERKE, E. A., and ARAMCO, S., 2009, Electrofacies and Geological Facies for Petrophysical Rock Typing: Khuff C, SPE Saudi Arabia Section Technical Symposium, SPE 126086, 9–11.
[4] MUKHERJEE, P., SINGHARAY, D., MATAR, S., and MESHARI, D. M. A., 2018, Rock-Typing: An Integrated Reservoir Characterization Tool for Tight Jurassic Carbonates, West Kuwait, 70372.
[5] AL-JAWAD, S. N. A., AHMED, M. A., and SALEH, A. H., 2020, Integrated reservoir characterization and quality analysis of the carbonate rock types, case study, southern Iraq, J. Pet. Explor. Prod. Technol., 10(8), 3157–3177.
[6] SKALINSKI, M., and KENTER, J. A. M., 2014, Carbonate petrophysical rock typing: Integrating geological attributes and petrophysical properties while linking with dynamic behaviour, Geological Society London Special Publications, 406 (1), 229-259.
[7] SAPUTELLI, L., CELMA, R., BOYD, D., SHEBL, H., GOMES, J., and BAHRINI, F., 2019, Deriving Permeability and Reservoir Rock Typing Supported with Self-Organized Maps SOM and Artificial Neural Networks ANN - Optimal Workflow for Enabling Core-Log Integration, predictive data analytics (PDA) and machine learning (ML), In SPE Reservoir Characterisation and Simulation Conference and Exhibition. OnePetro.
[8] SHABANINEJAD, M., BAGHERIPOUR, M., 2011, Rock Typing and Generalization of Permeability- Porosity Relationship for an Iranian Carbonate Gas Reservoir, The Nigeria Annual International Conference and Exhibition, Abuja, Nigeria, SPE 150819.
[9] ALIAKBARDOUST, E., 2013, Integration of rock typing methods for carbonate reservoir characterization, Journal of Geophysics and Engineering, 10(5), 1-11.
[10] HOSSEINZADEH, S., KADKHODAIE, A., and YARMOHAMMADI, S., 2019, NMR derived capillary pressure and relative permeability curves as an aid in rock typing of carbonate reservoirs, Journal of Petroleum Science and Engineering, 184, 106593.
[11] MOHEBIAN, R., RIAHI, M., A., and KADKHODAIE, A., 2019, Characterization of hydraulic flow units from seismic attributes and well data based on a new fuzzy procedure using ANFIS and FCM algorithms, example from an Iranian carbonate reservoir, Carbonates and Evaporites, 34(2), 349–358.
[12] FARSHI, M., MOUSSAVI-HARAMI, R., MAHBOUBI, A., KHANEHBAD, M., and GOLAFSHANI, T., 2018, Reservoir rock typing using integrating geological and petrophysical properties for the Asmari Formation in the Gachsaran oil field, Zagros basin, Journal of Petroleum Science and Engineering, 176, 161–171.
[13] MIRZAEI-PAIAMAN, A., OSTADHASSAN, M., REZAEE, R., SABOORIAN-JOOYBARI, H., and CHEN, Z., 2018, A New Approach in Petrophysical Rock Typing. Journal of Petroleum Science and Engineering, 166, 445–464.
[14] MIRZAEI-PAIAMAN, A., SABBAGH, M., OSTADHASSAN, M., SHAFIEI, A., REZAEE, R., SABORIAN-JOOYBARI and CHEN, Z., H., 2019, A further verification of FZI* and PSRTI: Newly developed petrophysical rock typing indices, Journal of Petroleum Science and Engineering, 175, 693–705.
[15] KHADEM, B., SABERI, M. R., ESLAHATI, M., and ARBAB, B., 2020, Integration of rock physics and seismic inversion for rock typing and flow unit analysis: A case study, Geophys. Prospect. 68(5), 1613–1632.
[1] Rock Typing
[2] Pore Facies
[3] Cut Off
[4] Flow Zone Indicator(FZI)
[5] Pore Throat Distribution
[6] SCAL
[7] multi-disciplinary dataset
[8] Lithostratigraphical Units
[9] Sedimentological Analysis
[10] Self-Organized Maps (SOM )
[11] Artificial Neural Networks (ANN)
[12] Fuzzy
[13] SIRIP
[14] matrix porosity
[15] microvugular porosity
[16] Rock Quality Index(RQI)
[18] Emerge
[19] training
[20] time slice
Modeling of rock types using combination of core, well logs and seismic data: in a carbonate hydrocarbon reservoir in south west of Iran
Mahdi Kheirollahi1, Golnaz Jozanikohan2, Reza Mohebian3*, Ali Moradzadeh4
1-Masters in petroleum exploration engineering, School of Mining Engineering, College of Engineering, University of Tehran, Tehran, Iran
2, 3, 4-Assistant Professor, School of Mining Engineering, College of Engineering, University of Tehran, Tehran, Iran
* mohebian@ut.ac.ir
Received: July 2022, Accepted: October 2022
Abstract
Rock typing is the process of assigning reservoir properties to geological facies, and an identified rock type has similar geological and reservoir properties. Due to the importance of identifying and separating rock types in hydrocarbon reservoirs, various methods have been proposed and developed today for the determination of rock types. One of the simplest methods is the porosity chart against permeability and cut-off determination, and one of the most important and practical of these methods is to determine rock types by flow zone indicator. In this study, we examine one of the most important hydrocarbon fields in the south of Iran where core, well and seismic data are available for the field studied so that by designing a new workflow and with use of the most important and efficient methods of rock typing, including FZI, porosity, Winland-Pitman index, FZI*, and logarithms of FZI and FZI* we have identified rock types in three dimensions and through the whole field. The partial comparison of the validation results after each method employment also confirms the highest accuracy belongs to the FZI method. As a result, by integrating this method with seismic attributes, the rock types have been separated in four different groups throughout the field in three dimensions.
Keywords: Rock typing, Hydrocarbon reservoirs, Reservoir properties, Seismic attributes, Flow zone indicator
رایمگ
الروابط
المراكز ذات الصلة
دعامة
الصفحات الرسمية
حقوق هذا الموقع الإلكتروني مملوكة لنظام Rimag Press Management. © 1438-1446