تعیین فشارهای سازندی با تلفیق روش‌های فرکتالی و زمین آماری در یکی از سازندهای هیدروکربوری جنوب غرب ایران

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

نویسندگان

1 1- گروه مهندسی نفت و معدن، دانشگاه آزاد اسلامی، واحد تهران جنوب، تهران، ایران. 2- شرکت ملی نفت ایران، مدیریت اکتشاف، اداره حفاری،

2 گروه پژوهش ژئوفیزیک، پژوهشکده علوم زمین، پژوهشگاه صنعت نفت، تهران، ایران

3 گروه مهندسی نفت و معدن، دانشگاه آزاد اسلامی، واحد تهران جنوب، تهران، ایران

4 شرکت ملی نفت ایران، مدیریت اکتشاف، اداره حفاری، تهران، ایران

5 گروه علوم زمین، دانشکده علوم طبیعی، دانشگاه تبریز، تبریز، ایران

10.22034/anm.2024.19796.1590

چکیده

آگاهی دقیق از فشارهای منفذی و شکست سازند برای حفر چاه ها به صورت ایمن با وزن گل مطلوب ضروری است. این مطالعه چالشی در زمینه مطالعات فشار سازند در میدان آزادگان جنوبی است که عموماً کربناته بوده و بجز سازند کژدمی فاقد لایه شیلی هستند. جهت مطالعات فشار سازند در عمیق ترین سازند مخزنی فهلیان، نیاز به مدل‌سازی کلیه سازندهای بالایی جهت کسب نتایج دقیق تر می‌باشد. با توجه تعداد کم جفت داده‌های فشارمؤثر-سرعت در کل داده‌های سازند فهلیان و ضریب همبستگی بسیار پایین رابطه باورز، نیاز به تفکیک این سازند به دو بخش بالایی و پایینی و انجام مدل‌سازی به تفکیک سازند‌ها پس از تکمیل مکعب‌های سرعت فشاری و فشار مؤثر بوده است. مطالعه صورت گرفته بر اساس داده‌های 23 حلقه چاه و تعبیر و تفسیر داده‌های لرزه‌ای صورت پذیرفته است. همچنین مدل‌های فشار مؤثر، منفذی و شکست سازند از مدل‌های زمین آماری ترکیبی تعیین شده و با مدل فرکتالی فشار-حجم مورد صحت سنجی قرار گرفته اند. بیشترین میزان همبستگی بین مکعب نهایی فشار مؤثر و مکعب سرعت مربوط به سازندهای فهلیان پایینی با 86/0 و ایلام با 71/0 می‌باشد. بر اساس مکعب های نهایی فشار سازند، حداکثر فشار منفذی به میزان 10000 پام در سازندهای گدوان تا فهلیان بالایی و حداکثر فشار شکست سازند نیز به میزان 13000 پام در سازندهای فهلیان پایینی تا گوتنیا بدست آمده است. در این تحقیق نوآوری جدیدی برای مطالعه فشارهای سازند به روش فرکتالی فشار-حجم انجام شده است. همچنین جهت ساخت مدل نهایی مکعب فشارهای سازندی در کل وسعت میدان آزادگان جنوبی، برای اولین بار از ترکیب روش‌های زمین آماری شبیه‌سازی گوسی متوالی و کوکریجینگ با مکعب مقاومت صوتی حاصل از وارون سازی لرزه‌ای به صورت توأم استفاده شده است. بر اساس محاسبه ماتریس لوگرشیو حاصل از مدل فرکتالی مقدار-حجم، بیشترین میزان تطبیق نهایی در بازه‌های سنگ آهک غالب به میزان 74/0 مربوط به سازندهای آسماری تا سروک محاسبه شده است، که نشان از تطابق بالای مدل مکعب فشار منفذی با استفاده از ترکیب شبیه‌سازی گوسی متوالی توام با کوکریجینگ و امپدانس صوتی حاصل از وارون سازی می‌باشد.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Formation Pressures Determination Utilizing the Integration of Fractal and Geostatistical Modelling in a Hydrocarbon Formation of SW Iran

نویسندگان [English]

  • Pooria Kianoush 1
  • Nasser Keshavarz Faraj Khah 2
  • Peyman Afzal 3
  • Emad Jamshidi 4
  • Amir Hossein Bangian Tabrizi 3
  • Ali Kadkhodaie 5
1 1- Department of Petroleum and Mining Engineering, South Tehran Branch, Islamic Azad University, Tehran, Iran. 2- National Iranian Oil Company, Exploration Directorate (NIOC-EXP), Tehran, Iran.
2 Geoscience Faculty, Research Institute of Petroleum Industry (RIPI), Tehran, Iran
3 Department of Petroleum and Mining Engineering, South Tehran Branch, Islamic Azad University, Tehran, Iran
4 National Iranian Oil Company, Exploration Directorate (NIOC-EXP), Tehran, Iran
5 Earth Sciences Department, Faculty of Natural Sciences, University of Tabriz, Tabriz, Iran
چکیده [English]

An authentic understanding of formation pore and fracture pressures is essential to define a safe and optimum mud window in drilling oil and gas wells. This investigation is a challenge in formation pressure studies in the South Azadegan field, which is typically carbonated with very low traces of shale beds, except in the Kazhdumi Formation. The wells drilled in this field hinged on the geological information and pore pressure alterations, can be categorized into three classes; Gachsaran, Pabdeh, and surface formations containing marl. These parameters directly affect the selection of the casing shoe depth and the well schematic. Correspondingly, target reservoir formations, i.e., Sarvak, Kazhdumi, Gadvan, and Fahliyan, and well profiles are other parameters that can classify wells in terms of drilling. It is necessary to analyze and model all upper formation pressures to obtain more precise results to investigate the pore pressure profile in these formations. Effective pressure log data reveals an increasing trend in formation pressure with depth in all wells. Besides, there are a few effective pressure-velocity data pairs in the total data of the Fahliyan Formation of Azadegan wells. With the small number of effective pressure-velocity data pairs in the total data of the Fahliyan Formations of Azadegan wells and the very low correlation coefficient of the Bowers relation for the wells of the Fahliyan Formations, it was necessary to separate this formation into two upper and lower parts. So the modeling has been performed by constructing compressional velocity-effective pressure cubes. This research was based on the data gathered from various drilled wells in this region and the interpretation of seismic data. Also, the effective, pore and formation fracture models have been determined from the integrated geostatistical models validated with the pressure-volume fractal model. The most heightened correlation between the final velocity and effective pressure cubes corresponds to the Ilam with 0.71 and the lower Fahliyan Formations with 0.86, which signifies the model's accuracy with the original data. Based on the final pressure cubes of the formation, the maximum pore pressure of 10,000 psi in the Gadvan to Upper Fahliyan Formations and the maximum fracture pressure of 13,000 psi in the lower Fahliyan to Gotnia Formations have been obtained. In this research, an innovation has been made to study the formation pressures utilizing fractal pressure-volume (P-V) methods. Also, for the construction of the final formation pressure cube model in the entire area of the South Azadegan field, for the first time, the combination of geostatistical methods of sequential Gaussian simulation (SGS) and co-kriging with the acoustic impedance (AI) cube obtained from seismic inversion has been used together. Computation of the Logratio matrix resulting from the fractal pressure-volume model revealed the maximum overall accuracy (OA) in the dominant limestone intervals as 0.74 at the depths of 2000-3000 meters, corresponding to the Sarvak to Asmari Formations. The results exhibit the high correspondence of the pore pressure cube model, obtained by sequential Gaussian simulation (SGS) combined with co-kriging and acoustic impedance inversion.

کلیدواژه‌ها [English]

  • Seismic velocity model
  • sequential Gaussian simulation
  • formation pressure cube
  • P-V Fractal model
  • Logratio matrix

مقالات آماده انتشار، پذیرفته شده
انتشار آنلاین از تاریخ 16 اردیبهشت 1403
  • تاریخ دریافت: 09 اسفند 1401
  • تاریخ بازنگری: 23 آذر 1402
  • تاریخ پذیرش: 16 اردیبهشت 1403