مدل‌سازی سازند آسماری به‌منظور تعیین میزان تولید کربنات و تأمین رسوب و پارامترهای مؤثر بر آن در میادین نفتی انتخابی دشت آبادان

اسدی مهماندوستی, الهام; لشگری, سمیه; معلمی, سید علی; دانشیان, جهانبخش

doi:10.22108/jssr.2021.25604

تعداد نشریات	43
تعداد شماره‌ها	1,682
تعداد مقالات	13,762
تعداد مشاهده مقاله	32,223,499
تعداد دریافت فایل اصل مقاله	12,752,554

مدل‌سازی سازند آسماری به‌منظور تعیین میزان تولید کربنات و تأمین رسوب و پارامترهای مؤثر بر آن در میادین نفتی انتخابی دشت آبادان

پژوهش های چینه نگاری و رسوب شناسی

مقاله 5، دوره 37، شماره 2 - شماره پیاپی 83، تیر 1400، صفحه 77-96 اصل مقاله (2.35 M)

نوع مقاله: مقاله پژوهشی

شناسه دیجیتال (DOI): 10.22108/jssr.2021.25604

نویسندگان

الهام اسدی مهماندوستی^* ¹؛ سمیه لشگری²؛ سید علی معلمی³؛ جهانبخش دانشیان⁴

¹استادیار، گروه زمینشناسی، دانشکدة علوم زمین، دانشگاه خوارزمی، تهران، ایران

²دانشجوی کارشناسی ارشد، گروه زمینشناسی، دانشکدة علوم زمین، دانشگاه خوارزمی، تهران، ایران

³استادیار، گروه اکتشاف و تولید، پژوهشگاه صنعت نفت تهران، ایران

⁴دانشیار، گروه زمین‌شناسی، دانشکدة علوم زمین، دانشگاه خوارزمی، تهران، ایران

چکیده

سازند آسماری توالی ضخیم کربناته- تخریبی در حوضة فورلند زاگرس است که در زمان الیگوسن- میوسن پیشین نهشته شده است. در این مطالعه به‌منظور تعیین میزان تأمین رسوب و پارامترهای مؤثر بر انتقال آن و میزان تولید کربنات حوضة رسوبی سازند آسماری در دشت آبادان از مدل‌سازی چینه‌ای رو به جلو استفاده شده است. نرم‌افزار مطالعه‌شده، نرم‌افزار Dionisos Flow و میادین انتخابی شامل میادین جفیر، آزادگان، یادآوران، امید و دارخوین بوده است. سازند آسماری با ترکیب سنگ‌شناسی سنگ آهک رسی، دولومیت ماسه‌ای، سنگ آهک ماسه‌ای، ماسه‌سنگ شیلی، سنگ آهک، ماسه‌سنگ و شیل در چاههای مطالعه‌شده از 4 تا 5 توالی رسوبی تشکیل شده است. محیط رسوبی سازند آسماری در چاههای مطالعه‌‌شده شامل رمپ داخلی و رمپ خارجی (رسوبات همی پلاژیک و پلاژیک) است. پارامترهای ورودی به نرم‌افزار به‌منظور مدل‌سازی شامل داده‌های لاگ گاما، تغییرات جهانی سطح آب دریا، نقشه‌های ضخامت، نقشه‌های عمق‌سنجی دیرینه و رخساره‌های رسوبی است. بررسی داده‌های محاسبه‌شده نشان می‌دهد انرژی حاصل از امواج آب و جریان آب از عوامل مؤثر در انتقال رسوبات و تغییر در روند رسوب‌گذاری منطقه است. بیشترین میزان تولید کربنات در سیستم تراکت تراز پیش‌روندة توالی دوم چاهها در محیط رمپ داخلی تخمین زده شده است.

کلیدواژه‌ها

تأمین رسوب؛ تولید کربنات؛ سازند آسماری؛ دشت آبادان؛ نرم‌افزار Dionisos Flow

اصل مقاله

مقدمه

سازند آسماری، توالی کربناته‌ای است که به دلیل توان مخزنی زیاد به‌مثابة یکی از مهم‌ترین سازندهای کربناته در ایران شناخته شده است (Agha Nabati 2006). این سازند برای نخستین‌بار در تنگ گل‌ترش واقع در کوه آسماری شناسایی شد (Richardson 1924). برش نمونة سازند آسماری با 314 متر ضخامت شامل سنگ آهک، سنگ آهک دولومیتی و سنگ آهک رسی است (Adams and Bourgeois 1967; Motie 1993; Agha Nabati 2006). به علت مقاومت و سختی خاص سنگ‌های آهک سازند آسماری، این سازند ارتفاعات و بخش‌های مقاوم منطقة زاگرس چین‌خورده را تشکیل می‌دهد (Darvishzadeh 2009). این سازند در برش‌های کامل به دو بخش عضو ماسه‌سنگ اهواز در نواحی جنوبی فروافتادگی دزفول و عضو تبخیری کلهر در جنوب غرب لرستان تقسیم می‌شود (Motie 1993) و در زمان الیگوسن تا میوسن در حوضة فورلند زاگرس تشکیل شده است (Alavi 2004). با توجه به اهمیت مخزنی سازند آسماری تاکنون مطالعات گوناگونی ازنظر ویژگی‌های مختلف زمین‌شناسی روی این سازند صورت پذیرفته است (برای نمونه Taheri et al. 2017; Farshi et al. 2019; Gharechelou et al. 2020)؛ اما تعداد مطالعات انجام‌شده با استفاده از نرم‌افزار Dionisos Flow به‌منظور مدل‌سازی چینه‌ای رو به جلو در حوضه‌های رسوبی ایران بسیار اندک است.

شبیه‌سازی حوضة رسوبی، علمی مؤثر در پاسخ به چالش‌هایی است که زمین‌شناسان نفت با آنها مواجه هستند. این علم به‌منظور اکتشاف هرچه بیشتر هیدروکربن و پارامترهای ناشناختة زمین‌شناسی استفاده می‌شود. مدل‌سازی چینه‌ای رو به جلو یکی از شاخه‌های مدل‌سازی حوضه‌های رسوبی است که در دو محیط کربناته و تخریبی یا ترکیبی از این دو به زمین‌شناسان نفت کمک می‌کند (Mirzaei et al. 2014). Dionisos یک نرم‌افزار مدل‌سازی چینه‌شناسی رو به جلو است (Stratigraphy Forward Modeling) که مؤسسة Francais du Petrole Institut در سال 1992 ساخته است (Granjeon and Joseph 1999). این نرم‌افزار میزان رسوب‌گذاری، تجمع و انتقال رسوبات، تولید کربنات و... را در طول زمان زمین‌شناسی محاسبه می‌کند و به تصویر می‌کشد (Mirzaei et al. 2014). جدیدترین مطالعة انجام‌شده با نرم‌افزار Dionisos Flow روی سازند آسماری در 8 میدان و 9 چاه نفتی در دشت آبادان و فروافتادگی دزفول صورت پذیرفته است که هندسة ساختاری حوضة رسوبی سازند آسماری را تعیین و ساختار هندسی پیش‌رونده با ماهیت پلکانی و تجمعی به سمت حوضه را برای این سازند پیش‌بینی کرده است (Asadi Mehmandosti et al. 2020).

بنابراین هدف از این مطالعه، استفاده از نتایج چینه‌شناسی توالی‌ها، مطالعات محیط رسوبی و داده‌های لاگ گاما در بررسی تحولات حوضة رسوبی سازند آسماری در میادین انتخابی دشت آبادان با استفاده از روش مدل‌سازی چینه‌ای رو به جلو است؛ به این منظور میزان تأمین رسوب و پارامترهای مؤثر بر انتقال آن و همچنین میزان تولید کربنات سازند مطالعه‌شده باید تعیین شود. نتایج بررسی تحولات حوضة رسوبی سازند آسماری با روش مدل‌سازی چینه‌ای رو به جلو می‌تواند در مطالعات اکتشافی صورت‌گرفته در دشت آبادان به کار رود.

زمین‌شناسی منطقه

در این پژوهش، سازند آسماری در 5 میدان و 6 چاه نفتی شامل میادین نفتی جفیر، آزادگان، یادآوران (چاه کوشک و حسینیه)، امید و دارخوین واقع در منطقة زاگرس برای مطالعه انتخاب شده است. موقعیت میادین نسبت به یکدیگر و میادین نفتی مجاور در شکل 1 نمایش داده شده است. میادین مطالعه‌شده ازنظر تکتونیکی در دشت آبادان قرار گرفته‌اند (Alavi 2004).

دشت آبادان، زون ساختاری واقع در انتهای جنوب غربی زاگرس است. مرز شمالی و شمال خاوری آن به جبهه چین‌های زاگرس (لبة جنوبی طاقدیس‌های سوسنگرد، آب تیمور و منصوری) محدود است و پس از عبور از جنوب میدان رگ سفید به خلیج فارس وارد می‌شود. مرز جنوبی دشت آبادان، خلیج فارس و عربستان است. دشت آبادان، بخشی از جلگة‌ میان‌رودان (بین‌النهرین) است که ازنظر زمین‌شناختی پایانة سکوی عربی به شمار می‌آید. دشت آبادان چین‌خوردگی‌های ملایم با روندی شمالی- جنوبی دارد که این روند مشابه با روند ساختارهای جنوب عراق، کویت و شمال شرقی عربستان است (Motie 1995).

شکل 1- موقعیت میادین مطالعه‌شده نسبت به یکدیگر (اقتباس همراه با تغییراتی از2015 Zeinalzadeh et al.)

Fig 1- Location of the studied oil fields relative to each other in Abadan plain (adapted with changes from Zeinalzadeh et al. 2015)

سازند آسماری در چاههای مطالعه‌شده از 4 تا 5 توالی رسوبی تشکیل شده است (گزارش‌های داخلی زمین‌شناسی شرکت ملی نفت ایران، شکل 2). توالی رسوبی اول هریک از چاههای مطالعه‌شده، رخسارة مناطق عمیق تا نیمه‌عمیق (عمدتاً آهک رسی) را شامل بوده و به بخش پایینی سازند آسماری به سن روپلین مربوط است. توالی دوم تا پنجم سازند آسماری در چاههای مطالعه‌شده عمدتاً ترکیب سنگ‌شناسی ماسه‌سنگی و ماسه‌سنگ شیلی دارد و با توجه به رخساره‌های تشکیل‌دهنده، رمپ داخلی (بخش تحتانی و فوقانی رمپ داخلی) را شامل می‌شود (شکل 2)؛ بنابراین در این مطالعه بخش‌هایی از محیط رمپ داخلی که به سمت خشکی است و عمق کمتر و انرژی بیشتری دارد، بخش فوقانی رمپ داخلی نامیده شده و درمقابل بخش‌هایی از محیط رمپ داخلی که به سمت حوضه است و عمق بیشتر و انرژی کمتری دارد، بخش تحتانی رمپ داخلی نام‌گذاری شده است. توالی دوم تا پنجم چاههای مطالعه‌شده، به بخش بالایی سازند آسماری به سن اکیتانین مربوط است. ضخامت سازند آسماری در چاههای مطالعه‌شده یکسان نبوده است. ضخامت این سازند در چاههای مطالعه‌شده و ضخامت توالی‌های پیش‌رونده و پس‌روندة هر توالی در جدول 1 نشان داده شده است (گزارش‌های داخلی زمین‌شناسی شرکت ملی نفت ایران). گفتنی است عمق دیرینة توالی اول در بعضی چاهها در سیستم تراکت HST بیشتر از TST در نظر گرفته شده است. علت این امر، فرسایش توالی اول در بعضی چاههای مطالعه‌شده یا در نظر گرفتن این توالی در بخش بالایی سازند پابده است.

جدول 1- ضخامت کل سازند آسماری و ضخامت توالی‌های پیش‌رونده و پس‌روندة هر توالی برحسب متر در چاههای مطالعه‌شده

Table 1- Total thickness of the Asmari Formation and thickness of progressive and regressive sequences of each sequence (in meters) in the studied wells

نام چاه	ضخامت کل (m)	توالی اول		توالی دوم		توالی سوم		توالی چهارم		توالی پنجم
نام چاه	ضخامت کل (m)	TST	HST	TST	HST	TST	HST	TST	HST	TST	HST
جفیر a	461	112	66	78	10	13	7	84	55	24	12
آزادگان b	429	64	85	80	16	11	11	85	57	13	7
کوشک a	367	52	41	100	20	10	7	83	36	12	6
حسینیه a	241	-	-	63	8	13	6	92	39	14	6
امید a	435	54	41	111	26	86	33	40	8	20	16
دارخوین c	216	-	-	20	9	65	32	50	5	22	13

مواد و روش پژوهش

در این مطالعه به‌منظور تعیین میزان تأمین رسوب و نحوة انتقال آن و میزان تولید کربنات در میادین مطالعه‌شده از نرم‌افزار Dionisos Flow استفاده شده است؛ بنابراین شناسایی توالی‌های رسوبی چاههای مطالعه‌شده در دشت آبادان مبنای کار قرار گرفته است (شکل 2). ژئومتری حوضة رسوبی سازند آسماری در محدودة مطالعاتی شامل Km160 X= (طول) و Km 100 Y= (عرض) است. مدل‌سازی چینه‌ای رو به جلو از گذشته تا حال حاضر در یک توالی از گام زمانی انجام می‌شود (Pellerin et al. 2018)؛ بنابراین در این مطالعه سن شروع تا پایانی شبیه‌سازی با استفاده از اطلاعات پالئولاگ هریک از چاههای بررسی‌شده و همچنین مطالعات ایزوتوپ استرانسیوم، My 4_/19-9_/33 تعیین شده است.

پارامترهای ورودی به نرم‌افزار شامل داده‌های لاگ گاما، تغییرات جهانی سطح آب دریا براساس منحنی Haq et al. 1988، نقشه‌های عمقی، نقشه‌های ضخامت، نقشه‌های عمق‌سنجی دیرینه و رخساره‌های رسوبی به‌منظور تعیین میزان تأمین رسوب، میزان تولید کربنات و پارامترهای انتقال رسوبات در حوضة رسوبی مطالعه‌شده است. داده‌های موجود به‌صورت عددی درمی‌آیند و سپس مطابق با روابط تعیین‌شده (روابط 1 تا 3) برای هریک از پارامترهای بیان‌شده، محاسبه و با نرم‌افزار Dionisos Flow پردازش می‌شوند (Granjeon and Joseph 1999).

رابطة (1) زمان/ میانگین ضخامت رسوبات × مساحت ناحیه= میزان تأمین رسوب

رابطة (2) شیب حوضه × جریان آب × ضریب انتشار آب= جریان رسوب

رابطة (3) شیب حوضه × انرژی حاصل از امواج × ضریب انتشار امواج= جریان رسوب

شکل 2- تطابق توالی‌های رسوبی سازند آسماری در دشت آبادان

Fig 2- Sedimentary sequences correlation of the Asmari Formation in the Abadan Plain

داده‌های محاسبه‌شده

میزان تأمین رسوب

منظور از میزان تأمین رسوب، مقدار و نوع رسوب (ازنظر مشخصات بافتی و ترکیب) واردشده به حوضه از خشکی‌ها (آواری) یا تشکیل‌شده در درون حوضه‌های رسوبی (شیمیایی و بیوشیمیایی) است. پارامتر اصلی در کنترل این فاکتور، منبع تأمین رسوبات در حوضه است. فرایندهای تکتونیکی، تغییر شکل‌های ساختاری، تغییرات جهانی سطح آب دریا (شکل 3)، اقلیم، پوشش گیاهی، لیتولوژی منطقة برخاستگاه و فیزیوگرافی کف حوضة آبریز از عوامل مؤثر در مقدار بار رسوبی و نوع رسوبات حمل‌شده به حوضه‌های رسوبی است (Amini 2011).

در مطالعة مدنظر میزان تأمین رسوب از تجزیه و تحلیل لاگ‌های چاهها به‌ویژه لاگ گاما برای هریک از توالی‌های تشکیل‌دهندة چاهها در واحد زمان مطابق رابطة 1 تخمین زده شده است (Granjeon and Joseph 1999)؛ (جدول 2).

شکل 3- منحنی تغییرات جهانی سطح آب دریا در بازة مطالعه‌شده (اقتباس با تغییراتی از Asadi Mehmandosti et al. 2020 Haq et al. 1988;).

Fig 3- Global sea level change curve in studied period (adapted with changes from Haq et al. 1988; Asadi Mehmandosti et al. . 2020)

بررسی میزان تأمین رسوب در توالی‌های مطالعه‌شدة سازند آسماری نشان می‌دهد میزان تأمین رسوب در واحد زمان یک روند افزایشی- کاهشی دارد؛ افزایش میزان تأمین رسوب با افزایش سطح جهانی آب دریا و فرایندهای تکتونیکی تأثیرگذار و کاهش میزان تأمین رسوب با کاهش جهانی سطح آب دریا صورت گرفته است. ارتباط میزان تأمین رسوب با تغییرات جهانی سطح آب دریا در واحد زمان براساس نظر Haq et al. 1988 در شکل 4 نشان داده شده است.

شکل 4- ارتباط میزان تأمین رسوب با تغییرات جهانی سطح آب دریا در هریک از توالی‌های چاههای مطالعه‌شده در واحد زمانی حوضة رسوبی سازند آسماری

Fig 4- Relationship between sediment supply and global sea level changes in each of the studied well sequences in the time unit of sedimentary basin of the Asmari Formation

جدول 2- مقادیر دقیق محاسبه‌شدة میزان تأمین رسوب در هریک از توالی‌های تشکیل‌دهندة چاهها در واحد زمان برحسب Km³/My

Table 2- Accurate calculated values of sediment supply in each of the sequences of wells in units of time according to Km³/My

نام چاهها	توالی اول		توالی دوم		توالی سوم		توالی چهارم		توالی پنجم
نام چاهها	TST	HST	TST	HST	TST	HST	TST	HST	TST	HST
جفیر a	3744	882	1825	360	869	273	2808	1029	1604	1123
آزادگان b	2139	1136	1872	576	735	429	2841	1067	869	655
کوشک a	1738	548	2340	720	668	273	2774	673	802	561
حسینیه a	0	0	1474	288	869	234	3075	730	936	561
امید a	1805	548	2597	936	5749	1287	1337	149	1337	1497
دارخوین c	0	0	468	324	4345	1248	1671	93	1470	1216

میزان تولید کربنات

بررسی‌های انجام‌شده در چاههای مطالعه‌شده بیان‌کنندة یک سیستم آواری- کربناته در حوضة رسوبی سازند آسماری است. در این سیستم تولید کربنات تابعی از زمان تشکیل هریک از توالی‌های تشکیل‌دهندة چاهها و عمق تشکیل رخساره در آب است؛ بنابراین در این مطالعه با بررسی رخساره‌ها مطابق شکل 5 و عمق تشکیل هریک از آنها در محدوده‌های توالی چاههای مطالعه‌شده، میزان عمق دیرینه برای هر توالی رسوبی محاسبه و در جدول 3 نمایش داده شده است.

با توجه به نوع رخساره، عمق دیرینه، مجموع ضخامت رخساره‌ها و زمان تشکیل هریک از آنها در چاههای مطالعه‌شده، اثر زمان بر میزان تولید کربنات براساس M/ My با نرم‌افزار Dionisos Flow تخمین زده شده است (شکل 7؛ جدول 4)؛ (Pellerin et al. 2018).

جدول 3- مقادیر عمق دیرینة هریک از توالی‌های تشکیل‌دهندة سازند آسماری در چاههای مطالعه‌شده برحسب متر

Table 3- Bathymetry values of each of the sequences of the Asmari Formation in the studied wells in meter

نام چاهها	توالی اول		توالی دوم		توالی سوم		توالی چهارم		توالی پنجم
نام چاهها	TST	HST	TST	HST	TST	HST	TST	HST	TST	HST
جفیر a	65	125	10	10	10	10	25_/5	25_/5	25_/5	25_/5
آزادگان b	65	65	10	10	10	10	10	10	10	10
کوشک a	65	125	10	10	10	10	10	10	10	10
حسینیه a	55 -	55 -	25_/5	10	10	10	10	10	10	10
امید a	125	125	10	10	10	10	10	10	10	10
دارخوین c	55-	55-	25_/5	25_/5	25_/5	25_/5	25_/5	25_/5	25_/5	25_/5

جدول 4- میزان تولید کربنات براساس تغییرات عمق در محیط‌های تشکیل‌دهندة سازند آسماری در دشت آبادان

Table 4- Carbonate production rate based on depth changes in the environments of the Asmari Formation in the Abadan Plain

محیط رسوبی	Pelagic		Hemi Pelagic		Bottom of Inner Ramp		Top of Inner ramp
محیط رسوبی	Min	Max	Min	Max	Min	Max	Min	Max
تولید کربنات (m/My)	0	3	8	32	51	198	25	46

شکل 5- مدل تغییرات رخساره‌ای هریک از توالی‌های سازند آسماری در چاههای مطالعه‌شده

Fig 5- Facies changes model of each sequence of the Asmari Formation in the studied wells

انتقال رسوبات

مدل‌سازی چینه‌ای رو به جلو مبتنی بر فرایندهای انتشار است؛ بنابراین بررسی انتقال رسوبات در حوضة رسوبی سازند آسماری در دشت آبادان با نرم‌افزار Dionisos Flow از طریق معادله‌های انتشار صورت گرفته است (روابط 2 و 3)؛ (Granjeon and Joseph 1999). عوامل مؤثر در انتقال رسوبات در حوضة رسوبی تحت کنترل بعضی پارامترها ازجمله نیروی کشش زمین، شیب حوضة رسوبی، انرژی موج و جریان آب است. براساس مطالعات انجام‌شده در سازند آسماری نظیر Adams and Bourgeois 1967 Wells 1967;، منشأ ماسه‌سنگ‌های بخش اهواز، محیط دلتایی با سن الیگوسن- میوسن است؛ بنابراین در این مطالعه عوامل مؤثر در انتقال رسوبات شامل جریان آب رودخانه (جدول 5) و انرژی حاصل از امواج در نظر گرفته شده است.

جدول 5- مقادیر دقیق محاسبه‌شدة میزان جریان آب در هریک از توالی‌های تشکیل‌دهندة چاهها در واحد زمان برحسب m³/s.

Table 5- Accurate calculated values of water flow in each of the sequences of wells in units of time according to m³/s

نام چاهها	توالی اول		توالی دوم		توالی سوم		توالی چهارم		توالی پنجم
نام چاهها	TST	HST	TST	HST	TST	HST	TST	HST	TST	HST
جفیر a	4003	2580	1209	632	357	372	1372	1924	1621	1464
آزادگان b	1352	1800	1240	1012	302	585	1392	1994	877	8540
کوشک a	1034	868	1550	1265	275	372	1360	1259	810	732
حسینیه a	0	0	976	506	357	319	1507	1364	945	732
امید a	1987	868	1720	1645	2367	1755	655	279	1351	1952
دارخوین c	0	0	310	569	1789	1702	819	174	1490	1584

بحث و بررسی

بررسی میزان تأمین رسوب در توالی‌های رسوبی سازند آسماری (شکل 6) نشان می‌دهد رسوبات نهشته‌شده در توالی اول چاههای مدنظر بیشتر از آهک رسی تشکیل شده است؛ بنابراین بیشترین (56%) میزان تأمین رسوب در این توالی با رسوبات پلاژیک از چاههای کوشک a (kushk- a)، آزادگان b (AZ-b)، امیدa (OD-a) و جفیر a (Jufyr-a) تأمین شده است (شکل 6- الف و 6- ب). در سیستم تراکت تراز پیش‌روندة توالی اول چاههای مطالعه‌شده با افزایش m20 سطح جهانی آب دریا (Haq et al. 1988)، میزان کل تأمین رسوب Km³/My9426 تخمین زده شده که بیشتر با چاه جفیر a (Km³/My 3744) تأمین شده است (جدول 2؛ شکل 6- الف)؛ به همین ترتیب با گذشت زمان و در سیستم تراکت تراز بالای توالی اول چاهها، با کاهش جهانی سطح آب دریا (m70) براساس شکل 6- ب، میزان تأمین رسوب به‌طور چشمگیری در این بخش از توالی چاهها کاهش یافته است (جدول 2). در بخش دستة رخسارة تراز پیش‌رونده و در بازة زمانی My29 تا 27 میزان کل تأمین رسوب، Km³/My 10576 محاسبه شده است. افزایش میزان تأمین رسوب در این قسمت از توالی به دلیل افزایش m25 سطح جهانی آب دریاست که بیشتر با رسوبات ماسه‌سنگی (شکل 2) بخش فوقانی رمپ داخلی چاههای امید a و کوشک a تأمین شده است (جدول 2؛ شکل 6- پ). به همین صورت با گذشت زمان کاهش جهانی سطح آب دریا به کاهش میزان تأمین رسوب به‌ویژه در چاههای حسینیه a و دارخوین c (بخش‌های نارنجی‌رنگ در مدل) در انتهای توالی مدنظر منجر شده است (شکل 6- ت). به طور کلی رسوبات تأمین‌کنندة توالی دوم شامل %49 رسوبات ماسه‌سنگی، %34 رسوبات کربناته و %17 متعلق به رسوبات شیلی است. رسوبات کربناتة توالی دوم به ترتیب با چاههای جفیر a، امید a و حسینیه a (Hos-a) (شکل 2) تأمین شده است. در توالی رسوبی سوم حداکثر تأمین رسوب در دوره‌های بالابودن سطح آب دریا (TST) به ترتیب با رسوبات ماسه‌سنگی، ماسه‌سنگ شیلی و آهک ماسه‌ای، Km³/My 13235 (شکل 2) تخمین زده شده است. در بازة زمانی My 7_/25- 9_/24 و در ابتدای توالی رسوبی سوم چا‌هها، سطح جهانی آب دریا m 25 بالا آمده است (Haq et al. 1988)؛ بنابراین فضا برای تجمع رسوبات در چاههای امید a و دارخوین c افزایش یافته، ولی در قسمت‌های شمالی مدل در چاههای حسینیه a، کوشک a، آزادگان b و جفیر a با افزایش سطح آب دریا، کاهش در میزان فضای رسوب‌گذاری صورت گرفته که ناشی از بالاآمدگی کف حوضه است (Asadi Mehmandosti et al. 2020)؛ بر این اساس با توجه به مطالعات انجام‌شده در توالی سوم هریک از چاهها با نرم‌افزار Dionisos Flow، بیشترین میزان تأمین رسوب در این توالی با چاههای واقع در شمال شرق مدل (چاه امید a و دارخوین c DN-c))) مطابق با افزایش تغییرات جهانی آب دریا (Haq et al. 1988) در این بخش از حوضه (شکل 6- ث) تخمین زده شده است؛ ولی به دلیل بالاآمدگی کف حوضه در چاههای حسینیه a، کوشک a، جفیر a و آزادگان b به کاهش فرونشست و میزان فضای رسوب‌گذاری (Asadi Mehmandosti et al. 2020) و درنهایت کاهش میزان تأمین رسوب در این بخش از حوضه منجر شده است (شکل 6- ث). با گذشت زمان کاهش جهانی سطح آب دریا به کاهش میزان تأمین رسوب به میزان Km³/My3744 در انتهای توالی سوم (HST) انجامیده است (شکل 6- ج)؛ به همین ترتیب در بازة زمانی My 2_/24 تا My8_/22 و در سیستم تراکت تراز پیش‌روندة توالی رسوبی چهارم چاهها، سطح جهانی آب دریا ثابت مانده و تغییری نکرده و ساختار هندسی حوضه ناشی از پایین‌افتادگی کف حوضه در چاههای حسینیه a، کوشک a، جفیر a و آزادگان b است (Asadi et al. 2020) که با افزایش فضای رسوب‌گذاری (Amini 2011) و درنهایت افزایش میزان تأمین رسوب در این بخش از حوضه همراه است؛ بنابراین میزان کل تأمین رسوب در این بخش از توالی، Km³/My 14506 تخمین زده شده است که حدود %73 از این میزان با رسوبات شیلی و ماسه‌سنگ شیلی تأمین شده است. در سیستم تراکت تراز بالای توالی بررسی‌شده با کاهش m 10 سطح جهانی آب دریا (Haq et al. 1988)، میزان تأمین رسوب تا Km3/My 3741 کاهش یافته و بیشتر با رسوبات ماسه‌سنگی تأمین شده است (شکل 6- ح)؛ به همین ترتیب با گذشت زمان در توالی پنجم با توجه به نوع و میزان رسوبات تشکیل‌دهندة توالی و تغییرات جهانی سطح آب دریا، میزان تأمین رسوب در قسمت‌های ابتدایی توالی مدنظر Km³/My 7018 و در بخش انتهایی توالی Km³/My 5613 محاسبه شده است (شکل 6- خ و 6- د).

شکل 6- مدل میزان تأمین رسوب در هریک از توالی‌های تشکیل‌دهندة چاهها در دشت آبادان در واحد زمان

Fig 6- Sediment supply rate model in each of the sequences of wells in the Abadan Plain per unit time

ارزیابی میزان تولید کربنات در حوضة رسوبی مطالعه‌شده با توجه به عمق آب (m 125 - 5) تشکیل رخساره‌های محیط‌های کربناته طی دوره‌های زمانی الیگوسن- میوسن نشان می‌دهد در بازة زمانی My 9_/33 تا My 4_/32، بیشتر رخساره‌ها و رسوبات تشکیل‌شده در محیط‌های عمیق براساس شکل 2، رسوبات پلاژیک است. با توجه به عمق دیرینه و مجموع ضخامت رخساره‌های تشکیل‌دهندة این بخش از توالی در چاهها، میزان تولید رسوبات پلاژیک در بازة زمانی یادشده، حدود m/My57 محاسبه شده (شکل 7- الف) که بیشتر با چاه امید a (شکل 2) تأمین شده است. در زمان و عمق بیان‌شده، رسوبات تخریبی (ماسه‌سنگ) مشاهده نشده و رسوبات کربناته به میزان خیلی کم (m/My 3) مشاهده شده است (جدول 4)؛ بنابراین ازنظر پارامترهای پتروفیزیکی همچون لیتولوژی و تخلخل، لایة تولیدی نفت نبوده و کیفیت مخزنی بسیار کمی دارد (Mosavi Harami 2012) (شکل 7- الف و 7- ب). با گذشت زمان، در بازة سنی My 29-4_/32 و با کاهش عمق آب (m 65)، در محیط همی‌پلاژیک، میزان رسوبات پلاژیک با رسیدن به حداکثر مقدار خود تقریباً m/My65 تخمین زده شده است.

با کاهش بیشتر عمق آب در محیط بخش تحتانی رمپ داخلی، میزان تولید رسوبات پلاژیک در حوضة رسوبی کاهش یافته است. کاهش میزان رسوبات پلاژیک به افزایش تولید کربنات در بازة زمانی توالی‌های بالاتر منجر شده است. در عمق 10متری (شکل 7- ب) از آب شامل محدودة زمانی My 29-8_/22، افزایش رسوب‌گذاری کربنات‌ها تا m/My198 در ابتدای توالی دوم چاهها در بازة زمانی My 29- 27 (شکل 7- الف) با چاه جفیر a، امید a و حسینیه a تأمین شده است.

رسوب‌گذاری کربنات‌ها در اعماق کم آب (m 5) و در محیط پرانرژی بخش فوقانی رمپ داخلی (شکل 7- ب) به حداکثر مقدار خود در ابتدای توالی پنجم رسیده است. میزان تولید کربنات در این بازه‌های زمانی براساس ضخامت رخساره‌های کربناته، عمق دیرینه و با عملکرد انرژی موج بیشتر توسط چاههای مطالعه‌شده، m/My 46 محاسبه شده است؛ به همین ترتیب با کاهش ضخامت رخساره‌های کربناته در بخش سیستم تراکت تراز بالای توالی پنجم چاهها در بازة زمانی 5_/19 - 19، میزان تولید کربنات m/My 25 محاسبه شده است (pellerin et al. 2018)؛ (جدول 4).

شکل 7- میزان تولید کربنات براساس دو پارامتر (الف) زمان و (ب) عمق رخساره‌ها در توالی‌های تشکیل‌دهندة چاههای مطالعه‌شده با نرم‌افزار Dionisos Flow.

Fig 7- Carbonate production rate based on two parameters of (a) time and (b) depth of facies in the sequences of wells studied by Dionisos Flow Software

مدل تغییرات چینه‌ای ترسیم‌شده در شکل 5 با توجه به نوع رخساره‌های تشکیل‌دهندة هریک از توالی‌ها و عمق تشکیل هریک از رخساره‌ها در واحد زمان (جدول 3) به‌خوبی بیان‌کنندة آن است که حوضة رسوبی سازند آسماری در منطقة مطالعه‌شده بیشتر تحت شرایط رسوب‌گذاری در بخش رمپ داخلی (شامل بخش تحتانی و فوقانی رمپ داخلی) قرار گرفته است؛ بنابراین ارزیابی انتقال رسوبات در حوضة رسوبی سازند آسماری در چاههای انتخابی دشت آبادان براساس رخساره‌ها و عمق تشکیل هریک از رخساره‌های تشکیل‌دهندة چاهها در واحد زمان نشان می‌دهد انرژی حاصل از امواج، یکی از عوامل مؤثر در انتقال رسوبات و تغییر در روند رسوب‌گذاری منطقه است. این عامل در نزدیکی ساحل اثر بسیار زیادی در تغییر روند رسوب‌گذاری دارد و با افزایش فاصله از ساحل و افزایش عمق آب دریا اثر آن کاهش می‌یابد (شکل 8). با توجه به اینکه توالی پنجم چاههای مطالعه‌شده براساس رخساره‌ها و محیط تشکیل‌دهندة آنها در اعماق کم قرار گرفته است، می‌توان نتیجه گرفت که انرژی حاصل از امواج در انتقال رسوبات توالی بررسی‌شده نقش مهمی دارد. انرژی محاسبه‌شده برای انتقال رسوبات در این توالی براساس محدودة عمق فعال موج، KW/m2_/48 است (رابطة 3).

شکل 8- ارتباط انرژی حاصل از امواج با عمق فعال موج در توالی پنجم چاههای انتخابی حوضة رسوبی سازند آسماری در دشت آبادان

Fig 8- Correlation of wave energy with active wave depth in the fifth sequence of selected wells of the Asmari Formation sedimentary basin in the Abadan Plain

بررسی میزان تأمین رسوب در هریک از توالی‌ها به تعیین میزان سرعت جریان آب در انتقال رسوبات در حوضة رسوبی سازند آسماری منجر شده است. در توالی رسوبی اول بیشترین میزان تأمین رسوب در چاههای واقع در بخش‌های شمالی مدل ازجمله کوشک a، آزادگان b و جفیر a صورت گرفته است؛ بنابراین سرعت جریان آب برای انتقال رسوبات در چاههای بیان‌شده افزایش یافته است (جدول 5؛ شکل 9- الف). ضریب انتشار آب نیز با توجه به رخساره‌های عمیق تشکیل‌شده در این بخش از توالی با اندازة رسوبات و محیط رسوبی کنترل می‌شود (Hawie et al. 2015)؛ بنابراین ضریب انتشار محیط عمیق پلاژیک، Km²/kyr4_/2 و میانگین شیب تشکیل‌شده در حوضة بررسی‌شده m/km35 با نرم‌افزار Dionisos Flow محاسبه شده است. با تعیین چنین پارامترهایی میزان جریان آب (Q_w) در ابتدای توالی اول m3/s 9488 تخمین زده شده است (شکل 9- الف). با گذشت زمان و با کاهش میزان تأمین رسوب، سرعت جریان آب برای انتقال رسوبات m³/s6116 محاسبه شده است (جدول 5؛ شکل 9- ب).

در دستة رخسارة تراز پیش‌روندة توالی دوم چاهها و در بازة زمانی My29 - 27 به دلیل بالابودن سطح نسبی آب دریا، میزان تأمین رسوب افزایش یافته است؛ بنابراین سرعت جریان آب در انتقال رسوبات تشکیل‌دهندة این بخش از توالی نسبت به قسمت‌های انتهایی توالی که با کاهش 10 متر از سطح جهانی آب دریا طبق منحنی Haq et al. 1988 همراه است، افزایش بارزی را نشان می‌دهد. ضریب انتشار محیط فوقانی رمپ داخلی در این توالی، Km²/Kyr3/1509 تخمین زده شده است؛ به بیان دیگر سرعت جریان آب برای انتقال رسوبات ماسه‌سنگ، آهکی، ماسه‌سنگ شیلی، شیل و آهک دولومیتی حداکثر 7005 متر مکعب بر ثانیه (m³/s) تخمین زده شده است (شکل 9- پ). با گذشت زمان در دستة رخسارة تراز بالا، کاهش میزان تأمین رسوب در این بخش از توالی به کاهش ضریب انتشار آب در محیط رمپ داخلی تا میزان Km²/Kyr 8/568 منجر شده است؛ بنابراین میزان جریان آب برای انتقال رسوبات در بخش انتهایی توالی دوم به‌ویژه در چاه جفیر a و حسینیه a کاهش یافته است (شکل 9- ت). در ابتدای توالی سوم چاهها به دلیل بالاآمدن 25 متر سطح جهانی آب دریا در چاههای دارخوین c و امید a، فضا برای تجمع رسوبات افزایش یافته و به همین ترتیب میزان تأمین رسوب در این بخش از توالی افزایش می‌یابد؛ بنابراین سرعت جریان آب برای انتقال رسوبات در بخش‌های شمال شرقی مدل افزایش یافته است. سرعت جریان آب در انتقال رسوبات در چا‌ههای واقع در شمال مدل به دلیل تأثیر فرایندهای تکتونیکی (بالاآمدگی کف حوضه) کاهش یافته است (جدول 5؛ شکل 9- ث). در سیستم تراکت تراز بالای توالی سوم و با کاهش سطح جهانی آب دریا، سرعت جریان آب برای انتقال رسوبات ماسه‌سنگی، ماسه‌سنگ شیلی، آهک ماسه‌ای و شیلی m3/s 5105 تخمین زده شده است (شکل 9- ج). با گذشت زمان در دستة رخسارة تراز پیش‌روندة توالی چهارم چاهها به دلیل پایین‌افتادگی کف حوضه در چاههای حسینیه a، کوشک a، آزادگان b و جفیر a، سرعت جریان آب برای انتقال رسوبات زیاد شده است (m³/s 7105)؛ (شکل 9- چ)؛ به همین ترتیب در قسمت‌های انتهایی توالی با کاهش سطح آب دریا، میزان تأمین رسوب کاهش یافته است و به کاهش ضرایب انتشار آب در رسوبات (Km²/Kyr 535) تشکیل‌شده منجر می‌شود (شکل 9- ح). با گذشت زمان در توالی پنجم چاهها با کم‌عمق‌شدن محیط، افزایش سرعت جریان آب برای انتقال رسوبات در چاههای امید a و دارخوین c صورت گرفته است (شکل 9- خ). در انتهای توالی پنجم چاهها، کاهش سطح جهانی آب دریا به کاهش سرعت جریان آب در انتقال رسوبات منجر شده است (شکل 9- د).

شکل 9- مدل میزان جریان آب در انتقال رسوبات در هریک از توالی‌های تشکیل‌دهندة چاهها در دشت آبادان در واحد زمان

Fig 9- Model of water flow rate in sediment transportation in each of the sequences of wells in the Abadan Plain per unit time

نتیجه‌

بررسی مدل‌سازی چینه‌ای رو به جلو در توالی‌های رسوبی سازند آسماری در 5 میدان نفتی مطالعه‌شده در حوضة زاگرس نشان می‌دهد میزان رسوب‌گذاری در هریک از توالی‌های تشکیل‌دهندة چاهها در کنترل دو پارامتر اصلی ازجمله تغییرات سطح جهانی آب دریا و فرایندهای تکتونیکی حاکم بر منطقه (بالاآمدگی و پایین‌افتادگی کف حوضه) است.

- در این مطالعه میزان تأمین رسوب در واحد زمان روندی افزایشی- کاهشی دارد. افزایش میزان تأمین رسوب با افزایش سطح جهانی آب دریا و فرایندهای تکتونیکی تأثیرگذار و کاهش میزان تأمین رسوب با کاهش جهانی سطح آب دریا همراه است.

- بیشترین میزان تأمین رسوب در بین توالی‌های رسوبی مطالعه‌شده در سیستم تراکت تراز پیش‌روندة توالی چهارم چاهها، Km³/My 14506 تخمین زده شده که با رسوبات ماسه‌سنگی، ماسه‌سنگ شیلی و شیل تأمین شده است.

- ارزیابی‌های صورت‌گرفته در این مطالعه نشان داده عوامل مؤثر بر انتقال رسوبات، جریان آب و انرژی حاصل از امواج است؛ بنابراین با تعیین میزان تأمین رسوب در هریک از توالی‌ها، انرژی حاصل از امواج و سرعت جریان آب در انتقال رسوبات در ناحیة دشت آبادان محاسبه شده است.

- انرژی حاصل از امواج آب در توالی پنجم چاهها، KW/m 2_/48 تخمین زده شده که با افزایش عمق در توالی‌های رسوبی پایین‌تر تأثیر انرژی موج کاهش یافته است.

- بررسی سرعت جریان آب در حوضة رسوبی سازند آسماری در دشت آبادان نشان داد میزان جریان آب در انتقال رسوبات تحت کنترل پارامترهایی همچون اندازة رسوبات، محیط رسوبی، تغییرات جهانی سطح آب دریا و فرایندهای تکتونیکی حاکم بر حوضة رسوبی است.

- مطالعة حاضر در حوضة رسوبی سازند آسماری در دشت آبادان نشان می‌دهد میزان تولید کربنات در توالی‌های رسوبی مطالعه‌شده تابعی از زمان تشکیل هریک از توالی‌های تشکیل‌دهندة چاهها و عمق تشکیل رخساره در آب است. در عمق 10متری از آب و در محیط بخش تحتانی رمپ داخلی که جزئی از محیط رمپ داخلی است، افزایش رسوب‌گذاری کربنات‌ها تا m/My198 در ابتدای توالی دوم چاهها در بازة زمانی My 29-27 توسط چاه جفیر a، امید a و حسینیه a تأمین شده است

مراجع

Adams T.D. and Bourgeois, F. 1967. Asmari biostratigraphy: Geological and Exploration Iranian Offshore Oil Company Report, 1074.

Agha Nabati A. 2006. Geology of Iran. Geological Survey of Iran, Second Edition, 603 p.

Alavi M. 2004. Regional stratigraphy of the Zagros fold – thrust belt of Iran and proforland evolution. American Journal of Sciences, 304: 1 - 20.

Amini A.H. 2011. Fundamentals of sequence stratigraphy. University of Tehran Publishing lnstitute, Second Edition, 332 p.

Asadi Mehmandosti E. Moallemi A. Daneshian J. and Lashgari S. 2020. Determining the structural geometry of Asmari Formation and investigating its changes in Zagros sedimentary basin using forward stratigraphic modeling method. Journal of Applied Sedimentology, 14:1-14.

Bordenave M. L. 1995. The cretaceous source rooks in the Zagros foothills of Iran. oil and gas science and Technology Revue de Institut Francais du petrole, 50: 727 – 752.

Coe, A. L. and Church, K. D. 2003. Sequence stratigraphy, in Coe, A. L. (ed.), The Sedimentary Record of Sea – level Change: Cambridge University Press, 57 – 98.

Darvishzadeh A. 2009. Geology of Iran Stratigraphy, Tectonics, Metamorphism and Magmatism. Amirkabir Publications, 434 p.

Farshi M. Moussavi-Harami R. Mahboubi A. Khanehbad M. and Golafshani T. 2019. 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.

Gawthorpe R. Hunt D. Taylor A. and Underhill J. 1994. NERC sequence stratigraphy course. Manchester University, Course notes, 172 p.

Gharechelou S. Amini A. Bohloli B. Swennen R. Nikandish A. and Farajpour V. 2020. Distribution of geomechanical units constrained by sequence stratigraphic framework: Useful data improving reservoir characterization. Marine and Petroleum Geology. 117:1048398.

Granjeon D. 1994. 3D stratigraphic modeling in complex tectonics area. In Abstract of the Assoc Petroleum Geology, Annual Convention and Exhibition, Abstract of the Assoc Petroleum Geology Search and Discovery Long Beach, California, 121–148.

Granjeon D. and Joseph P. 1999. Concepts and applications of a 3D multiple lithology, diffusive model in stratigraphic modeling. Recent Advances in Stratigraphic and Sedimentologic Computer Simulations: SEPM, Special Publication, 62:197–210.

Haq B.U. Hardenbol J. and Vail P.R. 1988. Mesozoic and Cenozoic chronostratigraphy and cycles of sea – level change.In: Wilgus C.K. Hastings B.S. Kendall C.G. Posamentier H.W. Ross C.A. Van Wagoner J.C. (Eds) Sea level changes: An integrated approach SEPM, Spes. Publication No. 42: Tulsa, 71 – 108.

Hawie N. Barrois A. Marfisi E. Murat B. Hall J. and Aillud G. 2015. Forward stratigraphic modelling, deterministic approach to improve carbonate heterogeneity prediction; Lower Cretaceous, Abu Dhabi. In Abu Dhabi International Petroleum Exhibition and Conference. Society of Petroleum Engineers, 2409 – 2446.

Madon M. 2008. A review of stratigraphic simulation techniques and their applications in sequence stratigraphy and basin analysis. Bulletin of the Geological Society of Malaysia, 10: 81-89.

Mirzaei Sh. Ghobashavi A. Ghazvini A. and Shakib S. 2014. Stratigraphic modelling using Dionisos Flow Software, Scientific – Extension Monthly of Oil and Gas Exploration and Production, 115: 38-45.

Mosavi Harami R. 2012. Sedimentology. Astan Quds Razavi Publications, 474 p.

Motie H. 1993. Geology of Iran, Zagros stratigraphy. Publications of the Geological Survey of Iran, 536 p.

Motie H. 1995. Geology of Zagros Oil. Publications of the Geological Survey of Iran, 589 p.

Pellerin M. Lanteaume C. Fournier F. and Borgomano J. 2018. Testing geologic assumptions and scenarios in carbonate exploration: Insights from integrated stratigraphic, diagenetic and seismic forward modeling. The Leading Edge, 37: 672 – 680.

Rahimpour-Bonab H. 2005. Carbonate Petrology correlates diagenesis and porosity evolution. University of Tehran Press, 487 p.

Richardson R. K. 1924. The geology and oil measures of southwest Persia. Journal of Petroleum Technology, 10 (43): 256- 283.

Taheri M.R. Vaziri-Mogaddam H. Taheri A. and Ghabeishavi A. 2017. Biostratigraphy and paleoecology of the Oligo - Miocene Asmari Formation in the Izeh zone, Zagros basin. Boletin de la Sociedad Geologica Mexicana, 69:59-85.

Wells A. J. 1967. Lithofacies and geological history of lower Tertiary sediments in southwestern Iran. (Asmari Lithofacies study). Iranian Oil Operating Companies, Geological and Exploration Division, Report No, 1108: 48 p.

Zeinalzadeh A. Moussavi-Harami R. Mahboubi A. and Sajjadian VA. 2015. Basin and petroleum system modeling of the Cretaceous and Jurassic source rock of the gas and oil reservoirs in Darquain ﬁeld, south west Iran. Journal Scientific – Extension Monthly of Oil and Gas Exploration and Production 26:419–426.

آمار

تعداد مشاهده مقاله: 853

تعداد دریافت فایل اصل مقاله: 412

سامانه مدیریت نشریات علمی. طراحی و پیاده سازی از سیناوب

پیوندهای مفید

اخبار و اعلانات

آمار

مدل‌سازی سازند آسماری به‌منظور تعیین میزان تولید کربنات و تأمین رسوب و پارامترهای مؤثر بر آن در میادین نفتی انتخابی دشت آبادان