SA01210708A - تحديد مواقع آبار مثلى في نموذج مكمن ثلاثي الابعاد - Google Patents
تحديد مواقع آبار مثلى في نموذج مكمن ثلاثي الابعاد Download PDFInfo
- Publication number
- SA01210708A SA01210708A SA01210708A SA01210708A SA01210708A SA 01210708 A SA01210708 A SA 01210708A SA 01210708 A SA01210708 A SA 01210708A SA 01210708 A SA01210708 A SA 01210708A SA 01210708 A SA01210708 A SA 01210708A
- Authority
- SA
- Saudi Arabia
- Prior art keywords
- value
- points
- reservoir
- wells
- quality
- Prior art date
Links
- 238000000034 method Methods 0.000 claims abstract description 82
- 239000000203 mixture Substances 0.000 claims abstract description 11
- 239000012530 fluid Substances 0.000 claims abstract description 8
- 239000011148 porous material Substances 0.000 claims abstract description 8
- 238000007405 data analysis Methods 0.000 claims abstract description 4
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims description 50
- 230000035699 permeability Effects 0.000 claims description 10
- 238000005553 drilling Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000013500 data storage Methods 0.000 claims description 6
- 238000001125 extrusion Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 206010000496 acne Diseases 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 241000234435 Lilium Species 0.000 claims description 2
- 239000003795 chemical substances by application Substances 0.000 claims description 2
- 230000037361 pathway Effects 0.000 claims description 2
- VAYOSLLFUXYJDT-RDTXWAMCSA-N Lysergic acid diethylamide Chemical compound C1=CC(C=2[C@H](N(C)C[C@@H](C=2)C(=O)N(CC)CC)C2)=C3C2=CNC3=C1 VAYOSLLFUXYJDT-RDTXWAMCSA-N 0.000 claims 1
- 238000007664 blowing Methods 0.000 claims 1
- 230000003247 decreasing effect Effects 0.000 claims 1
- 238000010899 nucleation Methods 0.000 claims 1
- 239000003208 petroleum Substances 0.000 claims 1
- 238000002834 transmittance Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 26
- 238000012360 testing method Methods 0.000 abstract description 14
- 238000005457 optimization Methods 0.000 abstract description 10
- 230000003068 static effect Effects 0.000 abstract description 9
- 238000009472 formulation Methods 0.000 abstract description 7
- 230000009897 systematic effect Effects 0.000 abstract description 5
- 238000010561 standard procedure Methods 0.000 abstract description 3
- 238000013459 approach Methods 0.000 abstract description 2
- 238000012856 packing Methods 0.000 abstract 1
- 210000004027 cell Anatomy 0.000 description 67
- 230000006870 function Effects 0.000 description 12
- 229930195733 hydrocarbon Natural products 0.000 description 8
- 150000002430 hydrocarbons Chemical class 0.000 description 8
- 238000011160 research Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 6
- 230000018109 developmental process Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 206010037888 Rash pustular Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 208000029561 pustule Diseases 0.000 description 3
- 241000282412 Homo Species 0.000 description 2
- 235000015076 Shorea robusta Nutrition 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000029052 metamorphosis Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 241001136792 Alle Species 0.000 description 1
- 241000272194 Ciconiiformes Species 0.000 description 1
- 241000079420 Fluda Species 0.000 description 1
- 244000166071 Shorea robusta Species 0.000 description 1
- 244000269722 Thea sinensis Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 235000019658 bitter taste Nutrition 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 108090001015 cancer procoagulant Proteins 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000007159 enucleation Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- HEDOODBJFVUQMS-UHFFFAOYSA-N n-[2-(5-methoxy-1h-indol-3-yl)ethyl]-n-methylpropan-2-amine Chemical group COC1=CC=C2NC=C(CCN(C)C(C)C)C2=C1 HEDOODBJFVUQMS-UHFFFAOYSA-N 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 238000010298 pulverizing process Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000002271 resection Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010025 steaming Methods 0.000 description 1
- 230000000699 topical effect Effects 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
- 238000012108 two-stage analysis Methods 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
- 238000013022 venting Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
- Image Processing (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Special Spraying Apparatus (AREA)
- Supply Devices, Intensifiers, Converters, And Telemotors (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
يتعلق الاختراع الراهن بطريقة منهجية فعالة من الناحية الإحصائية ثنائية المرحلة لتحديد مواقع آبار في نموذج مكمن ثلاثي الأبعاد مع الإيفاء بقيود مختلفة تشمل: أدنى مسافة بين الآبار، أقصى امتداد للبئر، حدود زاوية للآبار المنجزة المنحرفة وأدنى مسافة من حدود المكمن والمائع. وفي المرحلة الأولى، يحدد موقع الآبار بافتراض أن الآبار لا تكون إلا عمودية. وفي المرحلة الثانية، تختبر هذه الآبار العمودية بالنسبة للآبار المنجزة الأفقية والمنحرفة المستمثلة. ويعتبر هذا الحل ناجعا، مع أنه منهجي، ويزود مجموعة تمرير أولى جيدة لمواقع وأشكال الآبار. وتصيغ طريقة المرحلة الأولى مشكلة تحديد مواقع الآبار في صورة مشكلة برمجة باعداد صحيحة ثنائية (BIP) تستخدم طريقة تخزين البيانات التي تستغل بنية المشكلة وتعزز الصيغة المستمثلة وتقلل حجم المشكلة. وتتوفر حزم برامجية تجارية بسهولة لحل مشاكل BIP. ومن ثم تاخذ المرحلة الثانية بعين الاعتبار الآبار المنجزة العمودية المختارة لتحديد مسارات الآبار المتصلة بقيم مردود المكمن القصوى مع الإيفاء بقيود الشكل بما في ذلك: قيود المسافة بين الآبار المنجزة، قيود الانحراف الزاوي وقيود أقصى امتداد. ويتمثل المتغير الذي ينبغي استمثاله في المرحلتين في جودة المكمن المعدلة وفقا للتعرج. ويفضل أن تكون الجودة مقياسا استاتيا يعتمد على قيمة تقريبية مثل المسامية، المردود الصافي، الإنفاذية، سماكة الإنفاذية أو حجم المسام. وتنتج وحدات تخزين الخواص هذه عن طريق تقنيات قياسية من تحليل بيانات زلزالية وترجمتها وتفسير ورسم خرائط جيولوجية وطبيعة الصخور وفحص الآبار من آبار موجودة. وتوصف خوارزمية لحساب قيم الجودة المعدلة وفقا للتعرج.
Description
Y
في Ata تحديد مواقع آبار نموذج مكمن ثلاثي الأبعاد الوصف الكامل خلفية الاختراع: بطرق لتقليل تكاليف استخراج البترول من Ta gee يتعلق الاختراع الراهن التحديد يتعلق الاختراع الراهن بتحديد موقع dag مكامن جوفية إلى الحد الأدنى. وعلى بثر أمثل في نموذج مكمن جوفي ثلاثي الأبعاد. وتتمثل وظيفة حرجة لفريق إدارة المكمن في تشكيل خطة تطوير للمكمن مع اختيار ° تشكيل fag مجموعة من مواقع حفر الآبار ومواقع إنجاز تزيد الإنتاجية إلى الحد الأقصى.
Jong بمجموعة من خرائط خواص المكمن ومجموعة من قيود أساسية. Sale الخطة جيولوجيين وجيوفيزيائيين ومهندسين يختارون مواقع الآأبار باستخدام نماذج Bale الفريق مكامن. وتحدد مواقع الآبار لاستمثال بعض الخواص المنشودة للمكمن المتعلقة بإنتاجية |ّ الهيدروكربونات. وفي التطوير المبكر لحقل؛ قد تتكون هذه النماذج من خرائط المسامية - ٠ أو خصائص الصخور التي تعتمد بصفة أساسية على التفسيرات الزلزالية المقصورة على افتراض بسرعة Tle بعض الأبار التقييمية. وبمجرد تحديد النموذج؛ يطلب من الفريق مجموعة من مواقع تزيد الإنتاجية إلى الحد الأقصى. وتتعقد هذه المحاولة عن طريق الحاجة لأن تمتثل المواقع المختارة لمجموعة من القيود؛ مثل أدنى مسافة بين الآبار وأقصسى امتداد للبثر وأدنى مسافة من السطوح الملامسة للمائع أو حدود المكمن وقيود شكل vo البثر. وتكون المشكلة الناتجة توافيقية بدرجة عالية ومن ثم فقد تستغرق وق ت طويل لحلها. ويعتبر هذا صحيحاً بصفة خاصة لمكامن متغايرة الخواص ذات مناطق مغلة منفصلة. تقييم مجموعة جزئية صغيرة لتوليفات مواقع آبار Sole وتتضمن حلول عملية لهذه المشكلة ممكنة كدراسات لحالات ومن ثم اختيار تلك الحالات التي لها أعكلى قيمة لإنتاجية (asa المترية المنشودة؛ مثل المردود الصافي أو سماكة الإنفاذية (المتمثلة في Ye
وعندما يطور مكمن بآبار منتجة؛ JTS نموذج مكمن أكثر شضمولاً مع خرائط مفصلة للاستراتيجرافية والمناطق المغلة. وقد تتوفر كذلك خرائط توزيع الضغط أو خرائط الإشباع بالمائع من حالة تاريخية موائمة. ومن ثم يتطلب افتراض آبار تجاوزية أو ردمية اعتبار آخر لقيود تفرض عن طريق أداء الآبار الموجودة. وهكذا قد يصير ٠ اختيار مواقع الآبار أثناء تطوير مكمن أكشر تعقيداً. وثانية يعتبر هذا صحيحاً بصفة خاصة لمكامن متغايرة الخواص ذات مناطق مغلة منفصلة. وقد يكون إيجاد حلول لمتشكلة تحديد مواقع HUY) متزايدة التعقيد مهمة مملة ومتكررة. وثمة عدة دراسات منشورة حاولت استخدام قوانين خاصة ونماذج رياضية لتحديد مواقع و/أو أشكال JU جديدة في حقول الإنتاج. ولقد ذكرت النشرات التالية في هذا البيبان ٠ لتكون مرجعاً: -١ بحث د. سييفرت؛ جيه. جيه. إم. لويس؛ سي. واي. فهيرن و إن. سسي. تي. ستيل بعنوان "استمثال تحديد مواقع الآبار والمخاطر الناجمة عن ذلك باستخدام تقنيات تشكيل مكمن عشوائي ثلاثي الأبعاد"؛ 35520 (SPE قدم في مؤتمر تشكيل المكامن الأوروبي 7 ستافانجرء أبريل 149% -١ - بحث بي. أيه. جوتريدج و د. ئي. جاويث بعنوان "معايرة الأحجام المتصسلة اتصنيف ممرات الأآبار"» 35503 SPE قدّم في مؤتمر تشكيل مكامن ثلاثية الأبعاد الأوروبي؛ ستافانجرء ١7-١١ أبريلء A497 *- بحث جي. دبليو. روزنوالد و د. دبليو. جرين بعنوان "طريقة لتحديد المواقع المقلى لآبار في مكمن باستخدام برمجة بأعداد صحيحة .)١977( SPE J. "A ye : © 4- بحث لارس كجلسفيك و جيير جوهانسن بعنوان 'تحليل الريبة لإمكاتية إنتاج بثرء؛ Tale] على محاكاة مسار انسيابي للحصول على مكمن متعدد"؛ 8 ad في ندوة احصائيات جغرافية بترولية Ty ل «SPE/EAGE مدينة تولوزء أبريل 13499 TO بحث بي. إل. بكتر و اكس. سونج بعنوان "خطة تطوير حقل باستخدام Oli محاكى-جدولة وتحديد موقع بثر اقتصادي أمثقل"؛ 30650 (SPE قدّم في مؤتمر vo ومعرض ثتقنيات SPE السنوي؛ مدينة دالاس؛ YOY اككوبرء A390
H . >- بحث اس. فاسانثراجان و أيه. اس. سوليك بعنوان "اختيار موقع بثر باستخدام استمشال ببرمجة بأعداد صحيحة مثلى"؛ قدّم في اجتماع 14346 السنوي؛ مدينة بارسلوناء سبتمبر 1991 -١ بحث ام. جي. إرابتريتو؛ سي. أيه. فلوداس؛ اس. فاستثراجان؛ و أيه. اس. سوليك بعنوان ٠ "طريقة تعتمد على التحليل لتحديد موقع أمثل لآبار عمودية"؛ مجلة (AICHE مجلد £0( أبريل 1439 ص 44م-كعدي. 4- بحث كيه. بي. هيرد و أو. دوبرول بعنوان "تحديد مقدار ارتباطية المكمن بالنسبة لتطبيقات وصف المكمن"؛ 30571 (SPE قدّم في مؤتمر ومعرض تقنيات SPE السنويء تقييم JT SES وجيولوجية المكمن؛ V0 مدينة دالاس؛ ولاية تكساس. 0 4- بحث سي. في. دوتيش بعنوان "برامج فورتران لحساب ارتباطية نماذج عددية ثلاثية الأبعاد ولتصنيف إدراكات متعددة"؛ في مجلة علم الكمبيوتر وعلوم الأرض"؛ مجلد YE العدد ص 2-١4 ض -٠ بحث د. إل. تشوك وسي. سي. تشيين بعنوان "طريقة لتحديد مواقع الآبار وتوجيمها بتشكل أمثل للتعدين بمحلول"؛ البراءة الأمريكية رقم 5776 5,7 ٠١ فبراير AA) Vo : -١١ بحث تي. اس. لو و جيه. تشو. بعنوان "أداة تحديد ارتباطية مكمن هيدروكربوني باستخدام خلايا وكواشف للمردود"؛ البراءة الأمريكية رقم 5,707,177 YT مارسء مكح وقدم سييفرت ومن معه في المرجع )١( طريقة تستخدم نماذج مكمن إحصائية © جغرافية. ولقد أجروا تنتقيب شامل بتخميد وتدي لعدد كبير من مسارات مرشحة من مواقع منصات معينة ذات نصف قطر وزاوية ميلان وامتداد للبثر وسمت مضبوطة مسبقاً. وحلل JS مسار بئر إحصائياً بالنسبة للمردود الصافي المتداخل أو خصائص الصخور. ولم يكن موقع الآبار المرشحة متغيراً وبذلك أوجد الإجراء أقصى عدد من المواقع في مكان واحد من الناحية الإحصائية ولم يصمم للإيفاء بقيود البثر المتعددة.
واستخدم جوتريدج و جاويث في المرجع Tage (Y) الحجم المتصل لترتيب المواقف ع في بعدين ولكن لم يصفا الخوارزمية. ومن ثم بطريقة يدوية أعادا تحديد المواقع وتصميم الآبار في نموذج مكمن ثلاثي الأبعاد. وتعتبر هذه الطريقة "مكلفة" ولا تتوافق مع القيود على مواقع الآبار ويجرى اختيار مواقع الآبار في بعدين. ويعتبر هذا المرجع والمرجع السلبق
٠ حلول خاصة للمشكلة.
وقدم روزنوالد و جرين في المرجع (7) صيغة برمجة بأعداد صحيحة (IP) لتحديد الموقع JAY) لعدد قليل من الآبار. ولقد افترضا أن علاقة الإنتاج المحدد مقابل الزمسن معروفة للمكمن وأن المواقع المحتملة للآبار الجديدة محددة مسبقاً. ومن ثم تختار الخوارزمية عدد معين من الآبار من المواقع المرشحة وتحدد التسلسل الملائم لمعدلات
٠ التدفق من الأآبار. ّ( وصنف كجلسفيك و جوهانسن في المرجع )£( أحجام يمكن صرفها من الآبار باستخدام مسارات انسيابية لمواقع مختارة Tie وتزود المسارات الانسيابية دليلآ على مقدرة الصرف يعتمد على التدفق؛ ومع أن محاكاة مسار انسيابي أسرع بدرجة كبيرة من محاكاة Ap الفروق المنتهية على نطاق واسع إلا أن عدد العمليات المطلوبة في مخطط الاستمثال؛ > مثلاآً تلدين محاكى أو خوارزمية متكونة؛ لا يزال O(N?) حيث N يمثل عدد مواقع الخلايا السمتية الفعالة في النموذج. ويكون زمن الحساب dash جدآ بالمقارنة مع استخدام مقياس إستاتي. واستخدم بكنر وسونج في المرجع )0( كذلك محاكاة مسارات انسيابية متسقة مع طريقة استمثال عامة ولكنهما لم يستطيعا إجراء الاستمثال إلا على وحدات تخزين بيانات صغيرة Jan Ye وقدم فاساتثراجان و سوليك في المرجع )1( حلا لمشكلة اختيار مواقع الآبار لخرائط مكمن ثنائية الأبعاد Jie صيغة برمجة خطية بأعداد صحيحة فعالة من الناحية الإحصائية حيث استخدمت متغيرات LE لصياغة مواقع الآبار المحتملة. وهذه الصيغة غير ملائمة لوحدات تخزين بيانات ثلاثية الأبعاد. وقدّم إرابتريتو في المرجع (V) طريقة التحليل لمشاكل بيانات أكبر في خرائط ثلاثية الأبعاد.
: واستخدم هيرد و دوبرول في المرجع (A) طريقة محاكاة مسارات انسيابية في نماذج مكمن ثنائية الأبعاد لتحديد الارتباطية بين موقعين لبثرين. وكان ذلك بالنسبة لنماذج صغيرة نسبياً في بعدين وحددت الارتباطية بين نقطتين محددتين فقط. وقدم سي. في. دوتيش في المرجع (1) خوارزمية لتحديد الارتباطية تحل المتشكلة بإجراء خطوات ٠ تنقيب متداخلة عن "طبقات رقيقة" متشكّلة. وتعتبر هذه الخوارزمية بطيئفة بحيث يتعذر استخدامها. وقدم تشوك وتشيين في المرجع )٠١( طريقة خاصة لتحديد موقع dog ana مرتبة من الآبار تختار النسق الخلوي للمجموعة المرتبة من الآبار بحيث تعدل المنطقة الخلوية حسب الحاجة ويكون المحور الرئيسي للخلايا موازياً للمحور الرئيسي AY حقل ٠ > الآبار. ولا تحدد هذه الطريقة المواقع المثلى لآبار منفردة. وقدم لو. وتشو في المرجع (VY) طريقة لحساب الحجم الكلي الذي يمكن إنتاجه لبشر من موقع ثقب بئثر مختار. ولم ينشد استمثال الحجم الكلي الذي يمكن إنتاجه في هذا المرجع. ولقد فشلت النشرات المذكورة أعلاه في تزويد طريقة يمكن Lads عملياً ١ ا لاختيار مواقع إنجاز آبار متلى أو قريبة من المثلى في نموذج مكمن ثلاثفي الأبعاد لعدة أسباب»؛ ليس أقلها حجم حيز المشكلة. وتشمل نماذج زلزالية DE الأبعاد نموذجية *٠١-”٠١ نقطة صغيرة حجمية (تعرف كذلك بصور خلايا) . ولا تجد الطرق الموصوفة في النشرات أعلاه Sa بشكل فعال. ووفقآً لذلك. Lo حاجة لطريقة منهجية لتحديد مواقع آبار متلى أو قريبة من المشلى في نموذج © مكمن ثلاثي الأبعاد. ويفضل أن تكون الطريقة فعالة من الناحية الإحصائية وتفسر تقنية حفر الآبار المتطورة المتوفرة في الوقت الحاضر التي تكفل آبار منجزة أفقية و/أو منحرفة بدرجة Alle ذات امتدادات مختلفة قد تربط مواقع عالية المردود متعددة.
الوصف العام للاختراع: يتعلق الاختراع الراهن بطريقة منهجية فعالة من الناحية الإحصائية ثنائتية المرحلة لتحديد مواقع آبار في نموذج مكمن ثلاثي الأبعاد مع الإيفاء بقيود مختلفة تشمل: أدنى مسافة بين الآبارء أقصى امتداد للبشر؛ حدود زاوية للآبار المنجزة المنحرفة وأدتى م مسافة من حدود المكمن والمائع. وفي المرحلة الأولى؛ يحدد موقع الآبار بافقتراض أن الآبار لا تكون إلا عمودية. وفي المرحلة الثانية؛ تختبر هذه الآبار العمودية بالنسبة للآبار المنجزة الأفقية والمنحرفة المستمثلة. ويعتبر هذا الحل ناجعا؛ مع أنه منهجي؛ ويزود مجموعة تمرير أولى جيدة لمواقع وأشكال الآبار. وتصيغ طريقة المرحلة الأولى مشكلة تحديد مواقع الآبار في صورة مشكلة برمجبة ٠ بأعداد صحيحة ثنائية (BIP) تستخدم طريقة تخزين البيانات التي تستغل بنية المشكلة وتزز الصيغة المستمثلة وتقلل حجم المشكلة. وتتوفر حزم برامجية تجارية بسهولة لحل مشاكل . ومن ثم تأخذ المرحلة الثانية بعين الاعتبار الآبار المنجزة العمودية المختارة لتحديد مسارات الآبار المتصلة بقيم مردود المكمن القصوى مع الإيفاء بقيود الشكل بما في ذلك: قيود المسافة بين الآبار المنجسزة؛ قيود الانحراف الزاوي وقيود أقصى امتداد. ve ويتمثل المتغير الذي ينبغي استمثاله في المرحلتين في جودة المكمن المعدلة وفقا للتعرج. ويفضل أن تكون الجودة مقياسآً استاتياً يعتمد على قيمة تقريبية مل المسامية؛ المردود الصافيء الإنفاذية؛ سماكة الإنفاذية أو حجم المسام. وتنتج وحدات تخزين الخواص هذه عن طريق تقنيات قياسية من تحليل بيانات زلزالية وترجمتها وتفسير ورسم خرائط جيولوجية وطبيعة الصخور وفحص الأبار من آبار موجودة. وتوصف خوارزمية لحساب قيم الجودة Ye المعدلة وفقاً للتعرج.
A
شرح مختصر للرسوم:
Sl يمكن فهم الاختراع الراهن بشكل أفضل بالرجوع للوصف التفصيلي للتجسيد المفضل بالاقتران مع الرسوم المرفقة؛ حيث: القكلان al رسمان تخطيطيان يبينان الخطوات المتتالية لطريقة تمييز : [YI و]١[ ٠ صخرية؛ يبين وحدة تخزين بيانات مسامية ثلاثية الأبعاد تمثيلية؛ [VY] الشكل الكل [4] : يمثل وحدة تخزين بيانات تبين الكتل الصخرية المميسزة؛ الشكل ]0[ : رسم تخطيطي يبين الخطوات المتتالية لطريقة حساب جودة مكمن؛ Ye الشكل ]1[ : رسم توضيحي تخطيطي لبثر منحرفة؛ رسم تخطيطي يبين الخطوات المتتالية لطريقة اختيار ممر بثفشر : [V] JS أفقي/متنحرف. ومع أنه يمكن إجراء تعديلات مختلفة على الطريقة وفقا للاختراع وبذلك قد يكون له أشكال بديلة؛ إلا أن تجسيدات معينة تبين على سبيل المثال في الروم ye وستوصف بشكل أوفى في الوصف التفصيلي التالي. بيد أنه يتبغغي إدراك أن الرسوم والوصف التفصيلي المزودة في هذا البيان لا يقصد بها تحديد الاختراع بالتشكل الخاص والأشكال cll Sly الموصوف؛ ولكن على العكس يقصد بها أن تغطي كل التعديلات
البديلة التي تقع ضمن نطاق ومبداً الاختراع الرامن كما تحددها عناصر الحماية المرفقة. : الوصف التفصيلي: ولأغراض التوضيح؛ يركز الوصف التفصيلي التالي على مسائل اختيار مواقع : ٠ الآبار التي واجهها فريق إدارة منجم أثناء المراحل الأولى من تطوير مشروع؛ حيث تحدد مواقع الآبار لزيادة الإنتاجية إلى الحد الأقصى مع الإيفاء بالقيود. ومن المدرك أنه يمكن تطبيق الطريقة والتقنيات الموصوفة على مجموعة أوسع بكثير من المشاكل ولا يقصد بالوصف التفصيلي التالي أن يحدد نطاق الاختراع المطالب بحمايته. مقياس متري إستاتي لإنتاجية المكمن ١ يختار مقياس إنتاجية المكمن أثناء مرحلة المشروع الأولى عادة بحيث يكون مقياسا متريا إستاتيا لإنتاجية المكمن؛ مثل المردود الصافي (يعرف ب المسامية »* السماكة ” المساحة X نسبة الصافي إلى الإجمالي x درجة gS الهيدروكربوني)؛ أو سماكة الإنفاذية أو توليفة منهما. وبعبارة أخرى؛ لا تأخذ حركة المائع تحت الأرض -. غالبا بعين الاعتبار في تحديد موقع البثر في مرحلة تطوير الحقل هذه. ويركز على صياغة 46 القيود الفراغية والشكلية؛ Jie أدنى مسافة بين الآبارء أقصى امتداد cp Al الحدود الزاوية للآبار المنجزة المنحرفة؛ رأس المال ASH المتوفر أو أقصى عدد من الآبار وأدنى مسافة من حدود المكمن والمائع والمسافة من منصات بحرية أو منصات حفر التي ينبغي تحليلها إلى عوامل في اختيار هذه المواقع. ومن ثم يمكن إجراء محاكاة لمسارات انسيابية مفصلة لتحديد منهاج إنتاج ملائم من الآبار المرشحة هذه للايفاء بقيم الإنتاجية - المنشودة. وفي التجسيد المفضل؛ يتمثل المقياس الاستاتي في جودة المكمن أو يفضل في جودة مكمن مضبوطة وفقا للتعرج. ويعتمد حساب جودة المكمن على بعض الخواص المقاسة التي قد تعمل كقيمة تقريبية لمقدار الهيدروكربونات المتوفر أو المقدار الذي يمكن إنتاجه من الهيدروكربونات للاستخلاص عن طريق بثر. وتشمل أمثلة مقاييس تقريبية لإنتاجية vo بثر ملائمة: المسامية؛ المردود الصافيء الإنفاذية؛ سماكة الإنفاذية وحجم المسام. وتتوفر
تقنيات قياسية في مجالات تحليل بيانات زلزالية وترجمتها وتفسير ورسم خرائط جيولوجية وطبيعة الصخور وفحص الأآبار لتحديد هذه القيم لكل خلية حجمية (يطلق عليها في هذا البيان فيما بعد "نقطة صغيرة حجمية") في نموذج مكمن ثلاثي الأبعاد. وتحسب جودة المكمن لنتقطة صغيرة حجمية محددة بجمع القيم المقاسة ٠ التقريبية المتصلة ضمن نصف قطر تصريف مقدّر لبثر مستكشف للنقطة الصغيرة الحجمية المحددة. ويجوز ضرب القيم المقاسة التقريبية بأحجام النقاط الصغيرة المقترنة قبل الجمع. فعلى سبيل المثال؛ إذا كانت القيمة التقريبية هي عبارة عن المسامية فإن الجودة تمل حجم المسام المتصلة المجبموعة ضمن نصف قطر التصريف المفترض. وإذا كانت القيمة التقريبية هي عبارة عن المردود الصافي (يعرف ٠ - بناتج ضرب المسامية ودرجة الإشباع الهيدروكربوني والحجم ونسبة الصافي إلى الإجمالي) فإن الجودة تكافئ الحجم الهيدروكربوني الذي يمكن إنتاجه في الحجبسم المتصل بالنقطة الصغيرة الحجمية المحددة. وقد تكون الجودة قيمة تقريبيسة للإنتاجية أفضل من المسامية لوحدها لأن المسامية تعتبر مقياس موضعي ّ| بشكل كامل في حين تحدد الجودة حجم المسام المرتبطة. وقد تكيتف طريقة لو ١ وتشو في المرجع (VY) وفقآ للطلب الراهن ولكن توصف أدناه طريقة مفضلة بشكل أكبر لحساب الجودة. : وتتمثل إحدى المسائل التي تعرّضت لها الطريقة المفضلة لحساب الجودة في التعرج. وفي مكامن ذات حدود عديدة أو قنوات متعرجة أو مردود متداخل مع صخر طفلي أو متغير ما (ng قد تكون المسارات الإنسيابية الحقيقية في حجم متعرجة. ويؤدي © - التخلص من التعديل Gig للتعرج المقترن بالمقاييس التقريبية إلى تحسين وثوقية المتياس : الاستاتي. ويحسب التجسيد المفضل للطريقة الموصوفة جودة المكمن أولا بضبط القيم المقاسة التقريبية تحت قيمة قطع مختارة. ويمكن إجراء ذلك بتحديد قيم مقاسة تقريبية تساوي صفر للنقاط الصغيرة الحجمية التي لها قيم أقل من قيمة القطع أو vo بدلا من ذلك بالإشارة إلى نقاط صغيرة حجمية من هذا القبيل بنقاط غير فعالة. ومن
١١ ى ثم تنفذ خوارزمية تحديد الارتباطية لتعيين مجموعات من نقاط صغيرة حجمية فعالة متصلة (لا تساوي صفر). ويطلق على هذه المجموعات في هذا البيان فيما بعد وتنتج القيم المقاسة التقريبية من "وحدات تخزين بيانات" لخواص lie ٠ (مثل السعة؛ المعاوقة؛ المسامية وكثافة المسامية) قد تحتوي على قيم بيانات تتراوح من ٠١ مليون إلى ٠٠١ مليون. ويكون fale تقييم ارتباطية المكمن مملا. وفي الوقت السابق؛ وفّر العلماء الجيولوجيين أداة لتمييز كتلة مترابطة مفردة بمعرفة نقطة بزرية مثل موقع عند حفرة البثر. وينبغي تمييز كل كتلة وجعلها مرئية كل منها على حدة. ولوحدات تخزين كبيرة تشتمل على كتل عديدة؛ Dla حوالي ١7٠١ تستغرق هذه العملية ٠ ساعات عديدة وحتى أيام أو أسابيع. ولقد جرّبت خوارزميات آلية سابقة لتمييز الككل الصخرية. وتمتلت المشكلة في حسابها البطيء لوحدات تخزين بيانات كبيرة الحجم. فعلى سبيل المثال؛ أجرى جوتريدج وجاويث في المرجع () خطوة تمييز الكثتلة الصخرية لنماذج ثلاثية الأبعاد في طبقات رقيقة ثنائية الأبعاد لإجراء حساب عملي. وأنتجبت خوارزميات دوتيش في المرجع )3( أزمان الحساب التالية (زاد زمن الحساب بحوالي v6 ثلاثة أضعاف مقدار الزيادة في عدد الخلايا السمتية). ١ > 5". Ye “Va : ١ 1 va ¥ eo حوالي '٠١ (محسوبة بالاستقراء) وبالمقارنة؛ تشتمل خوارزمية تحديد الارتباطية الموصوفة في هذا البيان على زيادة خطية تقريبآً مع حجم وحدة تخزين البيانات. ويعتمد زمن الحساب على عدد الخاتآيا السمتية الفعالة وعدد الكل الصخرية المنفصلة. وتبين أمثلة قليلة في الجدول التالي.
ب YY Wax og Tes “To xy Moxy 76 وتحدد الخوارزمية بشكل سريع الارتباطية الداخلية ضمن وحدة تخزين بيانات ثلاثية الأبعاد كبيرة. وتصنف الكتل المترابطة؛ التي يشار إليها Gy tea J TE للحجم مما يتيح اختيارها بشكل منفرد أو في مجموعات لتصير مرئية. ويوضح الشكلان [1] و [7] خوارزمية مفضلة لتحديد الارتباطية. وبدءا من ° الخطوة )¥ + )(¢ توجه الخوارزمية تعليمات للكمبيوتر بتحميل مجموعة مرتبة ثلاثية الأبعاد من الخواص المقاسة. وفي الخطوة (؛١٠)؛ تعالج المجموعة المرتبة ثلاثية الأبعاد لتحديد الخلايا الصالحة. وتكون الخلايا صالحة )13 كانت الخواص المقترنة ضمن مدى قياس معين Sli) عندما تكون قيمة الخاصية المقاسة أكبر من قيمة قطع محددة). وإذا لم توجد خلايا صالحة؛ تنتهي الخوارزمية في الخطوة (V+) ومن ناحية ٠ ثاتية؛ في الخطوة ali (V0 A) المجموعة المرتبة من أرقام الكتل الصخرية التي لها نفس أبعاد مجموعة مرتبة ثلاثية الأبعاد إلى "٠" في الخلايا الصالحة و 'صفر" في كل الخلايا الأخرى. وفي الخطوة (١١١)؛ a gay رقم الكتل الصخرية (NGEO) إلى "٠ وفي الخطوة )١١"( يضبط دليل موقع (LOC) عند نقطة في خلية أولى. وفي الخطوة (VY £) يزيد دليل الموقع خلال كل الخلايا في المجموعة المرتبة ثلاثية الأبعاد. وفي الخطوة ١ - (١١١)؛ يجرى اختبار لمعرفة Lad )13 عولجت جميع الخلايا. وإذا تحقق ذلك؛ تعالج المجموعة المرتبة من أرقام الككل الصخرية في الخطوة )١١8( لتحديد حجم كل كتلة صخرية؛ وفي الخطوة JT aud (VY) الكتل الصخرية لتصنتف حسب الحجسم (بحيث تكون الكتلة الصخرية الأولى الأكبر). ومن ثم تنتهي الخوارزمية بعد الخطوة (VY)
\Y يجرى اختبار لتحديد فيما إذا كانت (VYY) في الخطوة Al ومن ناحية dds الخلية في المجموعة المرتبة من أرقام الكتثل الصخرية المشار إليها عن طريق الموقع صالحة وفيما إذا لم يحدد لها بعد رقم كتلة صخرية. وإذا لم يحدد لها رقم كتلة (011) صخرية؛ فإن دليل الموقع يزيد في الخطوة )£ )0( ويرجع الضبط إلى الخطوة ويحدد للخلية رقم (VY) هومن ناحية ثانية؛ يزيد رقم الكتل الصخرية في الخطوة وتمهّد قائمة خلايا صالحة معينة (VY) الصخرية الحالي في الخطوة A) ومن تم ١ عدادان لتلك القائمة إلى aT gays )١78( إلى صسفر في الخطوة (VVC) وترجع حلقات الضبط اللاحقة إلى الخطوة )١77( تنفّذ حلقة تمييز الكتل الصخرية .)١١( shall وفي .)١37( حلقة تمييز الككل الصخرية ]١[ ويبين الشكل ١ (LOC) بحيث يساوي دليل الموقع (VVC) يضبط العنصر الأول من قائمة (YY) إذا عولجت كافة عناصسر قائمة Lad يجرى اختبار لمعرفة ¢(Y + £) وفي الخطوة ومن ناحية ثانية؛ يضبط دليل (V1 £) وإذا تحقق ذلك؛ يرجع الضبط إلى الخطوة VVC عند موقع العنصر الحالي من قائمة 770 في الخطوة (CLOC) موقع الخلية الحالية |ّ بحيث يساوي دليل خلية مجاورة أولى (NCELL) ويضبط دليل الخلية المجاورة .)٠١١( ٠ في CLOC خلال كافة المواقع المتجاورة في NCELL ومن ثم يفهرس -(¥ + A) في الخطوة وقد يختلف تعريف "الخلايا المتجاورة” ولكن يفضل أن تكون الخلايا .(Y1 7) الخطوة المتجاورة هي الخلايا الستة التي تتقاسم سطح مع الخلية التي لها ©01.0. وفي الخطوة يجرى اختبار لمعرفة فيما إذا أخذت كافة الخلايا المتجاورة بعين الاعتبار. وإذا تحقق ؛)7٠١( ويرجع الضبط إلى الخطوة (؟١٠). ومن ناحية (YVY) في الخطوة (Y) ذلك؛ يزيد العداد vy. اختبار لتحديد فيما إذا كانت الخلية المجباورة (YY £) ثانية؛ يجرى في الخطوة إذا لم يحدد لها بعد رقم كتلة صخرية. وإذا لم يحدد لها بعد رقم كتلة صخرية؛ Lady صالحة li يزيد في الخطوة (116). وإذا حدث ذلك؛ يحدد للخلية المجاورة رقم الكتلة NCELL فإن | الخلية المجاورة (YYY) وتضيف الخطوتان (70؟) و (VIA) الصخرية الحالي في الخطوة وقد تعرف خلايا .)1١١( ومن ثم يزيد دليل 1108117 في الخطوة -VCC All] vo
VE
= بالخلايا الستة التي تشترك في سطح أو الخلايا ((Y +A) متجاورة بديلة (الخطوة لإضافية مشتركة الحافة أو الخلايا التسعة الإضافية مشتركة الزاوية وكل التوليفات I عشر نقطة لكافة الخلايا المتجاورة عندما يكون مردود YY منها. ويفضل التتقيب الذي يشتمل على ومنحدراً بالنسبة لتوجيه الخلايا. ويفضل التنقيب سداسي النقاط لخلايا تشترك في SU المكمن سطح عندما يكون مرود المكمن أكبر من سماكة الخلية مع تأرجح قليل بالنسبة لتوجيه ٠ للحالات المتوسطة. YA الخلية. ويفضل التنقيب بخلايا متجاورة عددها كتل صخرية باستخدام خوارزمية تحديد YS ولحساب جودة المكمن؛ تشكّل مجموعة مرتبة لخواص مقاسة ثلاثية الأبعاد [Y] الارتباطية الموصوفة. ويبين الشكل مليون خلية تقريباً. وهذه المجموعة المرتبة هي المسامية الحجمية Vo تشتمل على هي المسامية). وتشتمل المجموعة المرتبسة على Auli A al AM (أي أن ٠ متر تقريباً. WX خلية وتكون أبعاد كل خلية 9 متر 7 79 متر YE) X YOO X Yo) وتحول بيانات سعة الزلازل الأصلية إلى مقاومية حجمية وجزء من حجم الصخر الطقلي سجل البثر: وتعتبر المسامية lily باستخدام شبكات ظهرية معايرة تستخدم Vi, باستخدام قيم القطلع Vie ى| الحجمية مقدار محسوب يعتمد على توليفة من المقاومية و المحددة. وكانت قيمة القطع للمسامية 9617. وبإدراك المسامية الحجمية تنتج معلومات vo عن طريق AZ ial عن ارتباطية المسامية. ويبين الشكل [4] الككل الصخرية ALE خوارزمية تحديد الارتباطية. وتحسب قيمة جودة المكمن لكل نقطة صغيرة حجمية في النموذج بجع القيم المقاسة التقريبية ضمن حجم تصريف حول كل نقطة صغيرة حجمية المتوفرة في نفس الكتلة الصخرية بالإضافة إلى النقطة الصغيرة الحجمية وضرب الناتج بأحجام التقاط © الصغيرة الحجمية. ولضبط التعرج للمسارات الانسيابية الفهية تستخدم خوارزمية للتعرج تركيبة جوالة بشكل عشضوائي Gy للتعديل وفقآ للتعرج. وتستخدم خوارزمية التعديل لتحديد درجة عدم وجود حدود انسيابية في حجم التصريف. وتحدد التركيبات الجوالة بشسكل عشوائي بصفة أساسية امتدادات المسارات من كل موقع خلية إلى كل الحدود ضمن حجم
التصريف وتقلل مساهمة الخواص التي تقع Tae عن النتقطة الصغيرة الحجمية موضوع البحث. ويبين الشكل ]0[ تطبيق لطريقة التركيبة الجوالة بشكل عشوائي لحساب قيم
جودة المكمن المعدلة Gy للتعرج. وبدءآ من الخطوات (YY) إلى (١7)؛ توجه البرامج
تعليمات إلى الكمبيوتر بتحميل المجموعة المرتبة من الخواص المقاسة ثلاثية الأبعاد
وتحميل المجموعة المرتبة من الكتل الصخرية ثلاثية الأبعاد من الخوارزمية السابقة
وتمهيد المجموعة المرتبة لقيم الجودة ثلاثية الأبعاد إلى الصفر. وتتقاسم
هذه المجموعات المرتبة بعاد مشتركة. ويمهد دليل الموقع 1,02 في الخلية
الأولى في هذه المجموعات المرتبة في الخطوة (Y +A) وتجرى زيادته بشكل متتابع
© خلال كل الخلايا في الخطوة (YY) وفي الخطوة (١٠7)؛ يجرى اختبار لمعرفة
Lad إذا أجريت زيادة للدليل خلال كل الخلايا. وإذا تحقق ذلك»؛ تنتهي البرامج. ومن
ناحية ثانية؛ في الخطوة )1١"( يحدد مدى الخلايا التي يمكن Logi yaa من
الموقع الحالي. وفي تجسيد مفضل؛ هذا الحجم هو عبارة عن الحجم المتعامد
| للخلايا المحدد من ضرب نصف قطر Cig peal بنسبة باعية في كل اتجاه.
ve ويحسب أقصسى عدد من الحواف في الخطوة )£ (Y) ويفضل أن يساوي عدد
الحواف هذا عدد سطوح الخلايا في المساحة السطحية لحجم التصريف. بيد
أنه يختار بحيث يكون هذا العدد أقصى عدد لممرات تدفق عشوائية تشكلت من
الموقع الحالي. ويمهد عداد ممر إلى ١ في الخطوة (YT) وفي الخطوة )7١8(
يجرى اختبار لتحديد Lad إذا كان العداد أقل من أو مساو لأقصى عدد من الحواف.
© - وإذا لم يتحقق ذلك؛ تنتقل البرامج إلى موقع الخلية التالي في الخطوة (YY) ومن
تاحية ثانية؛ تبدأ تركيبة جوالة جديدة عند الموقع الحالي في الخطوة (77؟). وفي
الخطوة )£ (YY تتحرك التركيبة الجوالة بمقدار خلية واحدة باتجاه عضوائي. وفي
الخطوات (YY) إلى (١7)؛ تجرى سلسلة من الاختبارات لمعرفة فيما إذا انتققت
التركيبة الجوالة إلى خارج المجموعة المرتبة ثلاثية الأبعاد أو إلى خارج حجم
ve التصريف أو إلى خارج الكتلة الصخرية الراهنة. وإذا كان إحدى هذه الحالات صحيحاء
تزيد البرامج عداد الممرات في الخطوة (77). وقبل البدء بتشغيل تركيبة جوالة جديدة؛ تجري البرامج اختبارا لتحديد Lad إذا وصل مقياس الجودة إلسى القيمة القصسوى في الخطوة (TYE) وفي تجسيد؛ يجرى الاختبار لتحديد Lad إذا تغيرت قيمة الجودة للموقع الحالي بأكثر من مقدار التفاوت المسموح به المحدد مسبقا على عدد محدد مسبقا من م الممرات. فعلى سبيل المثالء إذا لم تتغير الجودة بأكثر من 961 في الممرات Sal الأخيرة؛ تفترض البرامج أن مقياس الجودة قد بلغ الذروة وتنتقل إلى الموقع التالي في الخطوة (770). وإذا لم يصل مقياس الجودة إلى الذروة؛ فإن البرامج تعود إلى الخطوة .)7١( وإذا بيبتت الاختبارات في الخطوات من (YY) إلى (YF) أن التركيبة الجوالة ٠ لا تزال في الحجم القابل للتصريفء فإنه يجرى اختبار في الخطوة (١776؟) لتحديد فيما إذا عين الموقع Mad) للتركيبة الجوالة. وإذا تحقق ذلك؛ تعود البرامج إلى الخطوة (4 77) للبدء بالخطوة التالية للتركيبة الجوالة. ومن ناحية ثانية؛ تضاف قيمة الخاصسية المقاسة لموقع التركيبة الجوالة الحالي إلى الجودة لموقع الخلية الحالي قبل البدء بخطوة التركيبة الجوالة التالية. وتقلل هذه الطريقة لتحديد قيمة جودة المكمن لخلية بشكل فعال 4s مساهمة قيم الخواص المقاسة للخلايا الأقل احتمالا للوصول إليها عن طريق التركيبة الجوالة بشكل عشوائي. وهذه الخلايا هي تلك الخلايا البعيدة عن موقع الخلية الحالي وتلك الخلايا المرتبطة بالخلية Adal عن طريق "شق صغير ضيق" أي مسار متعرج. وتضبط الجودة وفقا لتجسيد بديل عن طريق مقاومة التدفق في الممر كما تزودها قيم الإنفاذية في الخلايا. وينبغي ان تميز Aad الانتاجية التقريبية لجودة مضبوطة وفقا للتعرج؛ مواقع الابار القريبة من مركز يشتمل على حجم متداخل بدرجة عالية من تلك المواقع القريبة من حده الخارجي. تحديد موقع بئر ثنائي الأبعاد وبتحديد الآن مقياس استاتي يتعلق بإنتاجية المائع للمكمن؛ تكون الخطوة التالية في إدارة المنجم هي تحديد موقع وشكل الآبار. وينبغي أن تزيد الدالة المتشودة ys لاختيار مواقع الآبار ضبط إنتاجية كل الآبار إلى الحد الأقصسى مع الإيفاء بقيود
لال معينة. وفي التطبيق العملي؛ تختار مواقع SLY غالبا بمحاولة زيادة العلاقة مع المقياس الاستاتي إلى الحد الأقصى. ويكون من الصعب Tan صياغة النموذج الرياضي للحصول على مسافة بين آبار منجزة مستخدمة من هذا القبيل ويؤدي إلى تضخم كبير في حجم المشكلة لا يمكن حله عن ٠ | طريق أهلية أجهزة الكمبيوتر والخوارزميات العددية المتوفرة في الوقت الراهن. ومن نم تتمثل الطريقة المفضلة في إجراء تحليل ثنائي المرحلة يحل أولاً مشكلة تحديد مواقع إنجاز لآبار عمودية Lal ضمن وحدة تخزين بيانات مكمن ثلاثية الأبعاد. وفي المرحلة الثانية؛ تصير الآبار العمودية المختارة هي المواقع المرشحة التي ينبغي اعتبارها بسبب التدريج العالي في آبار أفقية أو متحرفة بدرجة عالية. وتحدد هذه ٠ الطريقة بشكل منهجي مسارات منحرفة بدرجة عالية قد تصل إلى مناطق Ale المردود منفصلة في حجم ثلاثي الأبعاد معين مع الإيفاء بقيود أقصى امتداد للبشر وأقصى زوايا انحراف. ويستخدم نموذج المرحلة الثانية قواعد نظرية بيانية لتزويد هيكل مدمج جديد لتحديد امتداد المسار النموذجي وسمت بئر أفقي أو متحرف لزيادة الإنتاجية إلى الحد الأقصى. Vo وبسبب الإجراء ثنائي المرحلة والطبيعة الناتجة للإجراء عالي التدرج: لا يمكن اعتبار الضبط النهائي لأشكال ومواقع الآبار المختارة مثالي بشكل تام. ومع ذلك تزود الطريقة المقترحة إجراء آلي لتحديد بسرعة مجموعة جيدة من مواقع Sad عمودية ومنحرفة بدرجة عالية متقاطعة مع مواقع خواص مكمن عالية الجودة مع Jia) لقيود المسافة بين الآبار وقيود فراغية أخرى.
ل : وفي الطريقة المفضلة؛ تصاغ مواقع الآبار في صورة برمجة بأعداد صحيحة ثنائية (BIP) حيث يكون موقع نقطة معتبرة عند موقع معين في المكمن صفر/١ لقرار القطع والوصل. ويمكن حل BIP عن طريق التعداد فقط. وهكذاء ثئمة قيود صارمة عن طريق الخوارزميات العددية المتوفرة وعن طريق القدرة الحسابية المتوفرة لحل مشاكل BIP م٠ المعقدة على نطاق كبير. وينبغي إعطاء انتباه كبيير لصياغة النموذج لتمييز بيات و/أو سمات Digna يمكن أن تستغل من قبل الخوارزميات العددية لحل مشاكل عملية. ويمكن وصف المشكلة بالكيفية التالية: لتمثل de pana 0000701 08 كل مواقع الآبار المحتملة ولتنتمي الرموز السفلية و ز إلى 1. dads متغير ثنائي :ا يختار من المجموعة (صفر؛ 4١ وجود/عدم وجود موقع ٠ بثرء وليمثل Q قيمة جودة المكمن المقترنة به. ويقترن مع كل موقع بر تكلفة معروفة للحفر والإنجازء :©. ويمكن أن يعبر عن المشكلة العامة لتحديد مواقع حفر Sy LLY نوعي في صورة: TRE اس ض i=] اد ve مع الخضوع للقيود التي تشمل: تحديد مواقع الأبارء المسافة بين الآبارء؛ شسكل الآبار ورأس المال المتوفر. وتصف الأجزاء التالية صيغ رياضية تصيغ بشسكل كمي مجموعة القتيود المذكورة أعلاه. ومع أن هذه المناقشات تركز على تطوير صيغ فعالة لوصف القيود من نوع شسكل البثرء إلا أنه يمكن إدراك أن نفس التقنيات قد تطبق لتمييز الأنواع © - الأخرى من القيود. وتعتبر كل نماذج الاستمثال المطورة قابلة للتكيف ويمكن تمثيلها بشكل بياني نسبي؛ ويمكن أن تتوافق بسهولة مع هذه القيود وقيود أخرى. وفي المرحلة الأولى؛ تستخدم وحدة تخزين بيانات جودة مكمن ثلاثية الأبعاد لتشكيل خريطة جودة ثنائية الأبعاد ٠. وتحدد خريطة الجودة ثنائية الأبعاد بضبط قيمة الجودة لخلية إلى قيمة الجودة القصوى في العمود المناظر للخلايا في وحدة تخزين البيانات ثلاثية ve الأبعاد. ويمكن اعتبار كل خلية في المجموعة المرتبة ثنائية الأبعاد كموقع محتمل حيث
to 4 OS حفر بئثر. وتكون كل خريطة ثنائية الأبعاد برتبة حوالي بضع عشرات الآلاف من الخلايا. وتتمشل المهمة في اختيار مجموعة جزئية من المواقع المحتملة هذه التي تزيد القيمة التراكمية للخاصية إلى الحد الأقصى مع ضمان أن تكون المسافة المستوية بين المواقع المختارة فوق أدنى قيمة محددة لتفادي تداخل الآبار. ° وتعرف Lad يلي المصطلحات التالية: (xi ,y:) Sad الإحداثيات المعروفة لهذه المواقع على شبكة متسامتة متعامدة؛ ليمشل زن المسافة الإقليدية بين أي موقعين لبئرين (i) حيث Gen)? + زوجو =V راص - ليمثل ...0 أدنى مسافة متشودة بين الأبار (بوحدات تسامتية)؛ ليمثل Now أقصى عدد من الآبار ينبغي اختياره. ويمكن التعبير عن صيغة BIP لاختيار موقع بثر في خرائط مكامن ثنائية الأبعاد كما يلي: TF 0) تدوج ّم \o مع الخضوع للقيود: : (7) :لا تنتمي ل (صفرء )١ Dj Din h 0 ازا م1 > الاجر )4( سند 234[ Y. وتمثل المعادلة )1( الربح والتكلفة الكليين لتحديد مواقع آبار عمودية. وتبين المعادلة (7) أن Yi هو عبارة عن متغير ثنائي. وتفرض المعادلة )7( قيد المسافة بين HUY وتحدد المعادلة (؛) عدد الآبار بحد أقصى. ولأن المعادلة (FY) تكون متكاقفقة عندما يتبادل gi فإنه ينبغي أخذ الحذر لتفادي مضاعفة معادلات القيود غير ve الضرورية.
AE . : ومن الملاحظ أن المعادلة (7) JA بشكل فعلي عدد كبير من معادلات القيود (بشكل تقريبي 07,14/2)؛ مما يؤدي إلى جعل تمييز مواقع آبار عمودية في خرائط مكامن ثنائية الأبعاد نمونجية مشكلة كبيرة يصعب معالجتها. ويمكن إعادة صياغة
المعادلة J | في صورة أخرى: = )© | منو لا > Dy ار {J | i= ,1 > +1 i (5) قز lie { 17 +5 ض J وبالإضافة إلى تقليل عدد معادلات القيود بشكل كبير؛ تضع هذه الصياغة العديد ٠ .من معادلات القيود في شكل cad بيانات" بحيث تستغل من قبل أجهزة تتاظرية لبرامج تجارية لتقليل حيز المشكلة. وعلى وجه التحديد قد تستغل أجهزة تناظرية من نوع 10 تجارية Je سبليكس (علامة تجارية مسجلة) و OSL (علامة تجارية مسجلة) صورة المعادلة )7( بإجراء تشعب على المتغيرات الثنائية المستخدمة في صسورة "مجموعة مرتبة بشكل خاص". تحديد موقع بثر شلاثي الأبعاد : عند تشكيل خرائط مكامن ثنائية الأبعاد؛ كان التركيز على ضمان أن تكون المسافة المستوية بين مواقع الآبار المختارة أكبر من أدنى قيمة محددة. وفي وحدات تخزين بيانات مكمن ثلاثية الأبعاد تظهر الخواص الطبقية للمكمن تغيرات SSK في الاتجاه العمودي أو الاتجاه )2( وإذا كان هنالك تغير كاف في خواص المكمن في © > الاتجاه ()؛ يمكن التقرير ob ينجز البشر في مناطق متعددة بأعماق مختلفة. وهكذا؛ بالنسبة لوحدات تخزين بيانات ثلاثية الأبعاد لا يكون ضمان أن نقفقي مواقع حفر SLY) بقيود المسافة في المستوى (x, y) فقط كافياً. وبالإضافة إلى ذلك؛ ينبغي ضمان أن ثفي كذلك مواقع إنجاز الآبار التي تحدد على امتداد الاتجاه )2( بهذه القيود. وكذلك؛ بالنسبة LY أفقية أو منحرفة ينبغي ضمان أن تكون هذه القيود ve موفية بالغرض على امتداد الطول الكلي لمسارات الآبار.
AR
وتوضح الأهداف المشفرة الملونة في الشكل ]£[ ككل صخرية متفصلة. وتعرّف الجودة لبثر منجز في كتلة صخرية في هذا البيان بالجبودة القصسوى المواجهة في كل النقاط الحجمية الصغيرة العمودية الموجودة في نفس ALS) الصخرية على خريطة ذلك الموقع (أي الجودة القصوى في عمود LI صخرية). ٠ وينبغي أن يكون بين الآبار أدني مسافة (Doin) إذا أنجزت ضمن نفس الكتلة الصخرية. وإذا وجدت وحدات انسيابية منفصلة في المكمن؛ أي كتل صخرية مختلفة؛ قد تبعد الأبار عن بعضها بمسافة تقل عن ئ.0. وإذا وجدت وحدات انسيابية فوقية يمكن إنجازها عن طريق حفرة بثر واحدة؛ فإنه ينبغي أن توجد تكلفة لإنجازات متعددة تشتمل عليها الدالة المتشودة. ١ وتصيغ عملية اختيار موقع البئثر وحدة تخزين البيانات ثلاثية الأبعاد في صورة طبقات ثنائية الأبعاد متكدسة. وتناظر الخلايا في الطبقة العليا الموزعة في المجال (x,y) مواقع الآبار المحتملة؛ كما هو في Ala وحدة تخزين البيانات ثنائية الأبعاد. ولتمثل 17 هذه المجموعة من مواقع الآبار المحتملة. والآن. من كل موقع من هذه المواقع عندما تجتاز الطبقات بشكل مستعرض نحو الأسفل بخط مستقيم في الاتجاه o(z) تواجه النقاط الصغيرة ١ الحجمية في الكتلة الصخرية. وتوجد مواقع إنجاز صالحة محتملة عديدة لكل موقع بشسر في المجال (x, y) بمقدار مماثل لعدد المواقع الموجودة في الاتجاه « التي تقطع ككل صخرية مختلفة (أي طبقات منفصلة طبقاية). Jagd 6 مجموعة التقاط الصغيرة الحجمية في الكتلة الصخرية. وتمثل التوليفة من هذه المجموعات؛ أي W تقاطع JSG الإنجازات الصالحة. وتثقترن قيمة جودة بكل إنجاز صالح من هذا القبيل. وتعرّف © | الصيغة مجموعة المتغيرات الثنائية؛ (Y(WAG) بمجموعة مرتبة من متغيرات ثنائية لها قيم صفر/١ للإشارة إلى وجود أو عدم وجود إنجاز. ويمشل QWNG) المجموعة المرتبة من قيم الجودة المقترنة. ومن ثم ينبغي تطبيق قيود المسافة بين الآبار على مواقع إنجاز آبار مختلفة ضمن كتلة صخرية AK Jah) صخرية). ومن الملاحظ أن الإنجازات داخل الكتلة الصخرية »| غير مقيدة. ومن الملاحظ كذلك أنه يمكن تحديد هذه القيود بالأخذ بعين الاعتبار ES
YY
واحدة في كل مرة وتسجيل مجموعة قيود المسافة بين الآبار كما هو مبين في ةيرخص٠ (1) المعادلتين (*) و يتشد Cua ويتمشل وجه جدير بالاهتمام لهذه المشكلة في صياغة الدالة المتشودة؛ مقايضة زيادة الجودة الكلية لمواقع الآبار المختارة إلى الحد الأعلى مقابل تكلفة حفر وإنجاز الآبار. ويعمل الطرف الأول في الدالة المنشودة على زيادة الجودة التراكمية إلى ٠ الحد الأعلى للمواقع المختارة: أكبر قيمة ل 1 الاجم )7( ww G : وتشمل الحدود المالية ما يلي: إذا أنجز بثر بشكل منفرد فإنه يسبب في وتشكل إنجازات إضافية جزء قليل من هذه التكلفة .» Ae تكلفة معينة؛ » 1/١ لكل منها. ولصياغة بنية التكلفة هذه يعرف حد تكلفة ابت يساوي » 1/١ متلا ينتج عندما ينجز بثر. ويمكن الملاحظة بسهولة أن هذه الصيغة تمل بنتية التكلفة المتشودة. بيد أنه لتمثيل ذلك بشكل كمي؛ يستلزم متغير إضافي لصياغة اختيار موقع إنجاز بثر في ككلة صخرية dle البشر. (ويذكر أن المتغير لا في هذه الحالة : لتشير إلى X(W) وليس اختيار موقع بثشر). ومن ثم تعرف المجموعة المرتبة الثنائية ١ أي مجال الخريطة W وجود أو عدم وجود بثر في مجموعة المواقع المستوية وبما أن كل الآبار المنجزة عمودية تمامآ؛ فإنه يضاف متغير موقع .)». y) واحد فقط لكل المتغيرات المناظرة )2 ,7 ,»)7. ويمكن دمج بنية التكلفة X(x, v) المقترحة في الدالة المتشودة في صورة: 277 60م - 52X07) - EPWRLCEL) E(B بعلاقة يمكن تمثيلها كما يلي: FD ¢ وترتبط المجموعتان من المتغيرات الثنائية ,أ XP) < 220,76 (3)
A وتكفل المجموعة المبينة أعلاه من المعادلات أنه إذا أنجز بثر في فإن ١ ؛ مساوياا ل Y(WAG) صخرينة: أي أنه إذا كان أي من المتغيرات الثنائية؛ vo
١ موقع حفر البشر المقترن + XW) يكون مساوياً ل ١ كذلك. ويكفل عكس هذه العبارة؛ أي "إذا لم تختار كل الإنجازات المقترنة بموقع بقرء أي أن Y (WAG) يساوي pia فإن المتغير الثنائي المقترن X(W) يساوي صسفر"؛ عن طريق الدالة المتشسودة في المعادلة (A) حيث أن XW) يعتبر Te a من حد التكلفة السلبي فقي ٠ دالة منشودة وجدت أكبر قيسمة لها. وفي الحقيقة. يمكن ملاحظة أنه لا يلزم التصريح بشكل واضح عن المتغيرات»؛ (XW) بحيث تكون من النوع الثتقائي ؛ ولكن قد تعتبر كمتغير مستمر محدد بين صفر و a وتكفل صيغة الدالة المنشودة وتمثيسل القيود المبينة أعلاه إمكانية استخدام XW) على قيم صحيحة ملائمة فقط. Lad يلي النموذج النهائي لتحديد المجموعة المثتلى لمواقع الآبار وإنجازات عمودية تمامآً في نموذج مكمن ثلاثي الأبعاد: (٠)أكبرقيمةل ET rwac - ضر - TOF AIF AG) re 2 wo . 2 WG مع الخضوع للقيود التالية: Dg 0 A . 1h Je FG) مو 2ك <Dy كك زع F nD = Ly ji تج ه LF \e )1( 5 Y,+TY. <1, [i i= J, Dy 1 i,je@ nG) rg J \ (OY) pI ACOET I VY ¢ Y. XW)2YW 0G) 70 (4) (Ve) }01{ > 70770 0c Xs] ("1 ve ولا تزال المشكلة المعيقة في الصياغة أعلاه هي حساب وتعيين القيود لضمان أن ض تكون الآبار التي تتجبز ضمن نفس الكتلة الصخرية مفصولة عن بعضها البعض بمسافة ,ا على الأقل. وتتعلق هذه المحاولة مباشرة بعدد النقاط الصغيرة الحجمية؛ أي الإنجازات المحتملة؛ في كتلة صخرية حيث ينبغي تحديد القيود لكل زوج من التوليفات وهكذا قد تستغرق Dippy م الأبار منجزة من هذا القبيل تبعد عن بعضها البعض بمسافة تقل عن الخرائط ثلاثية الأبعاد المترابطة بدرجة عالية؛ أي تتألف من بعض الكتل إنجاز محتمل لكل كتلة صخرية) زمنا طويلا ٠١ الصخرية كثيفة السكان (حوالي لتحديدها وتفسيرها. بيد أنه يمكن حل بشكل فعال مشكلة المكامن الضخمة غير المتجانسة التي تشتمل على مناطق مغلة منفصلة بسبب عدم تنفيذ القيود داخل الكتلة الصخرية. ولتوضيح مزايا الطريقة الموصوفة أعلاه؛ قوبل أداؤها مع طريقة "مكلفة". وتختار الطريقة المكلفة بشكل متعاقب مواقع الآبار بترتيب تنازلي لجودة المكمن مع الإيفاء بقيود المسافة بين الآبار. وفيما يلي خطوات طريقة من هذا القبيل: تحدد الجودة القصوى في عمود النقاط الصغيرة الحجمية (W عند كل موقع مستو -١ بصفتها قيمة الجودة التمثيلية له. 6 تستقصى من المواقع المأخوذة بعين الاعتبار المواقع التي لها قيم جودة تقل عن -” قيمة القطع الدنياء جودة من المواقع المتبقية؛ led يختار موقع إنجاز البثر الذي له -* ؛- تستقصى من المواقع المعتبرة الآتية كل المواقع المتبقية في نفس الكتلة للبثر المنجز المختارء Dy الصخرية التي تكون ضمن +. إذا كان عدد المواقع المختارة أقل من العدد الأقصسى المسموح به يرجبسع 0 للخطوة (؟)؛ تحسب الجودة التراكمية وتكلفة المواقع المختارة لتحديد قيمة الدالة المنتشودة = النهائية.
Yo
وقد تكون مجموعة مواقع SLY) المختارة باستخدام الخوارزمية من انوع المكلف قريبة من المجموعة المثلى؛ حيث لا يوجد طريقة منهجية لتقدير القرار والتراجع لتصحيح القرارات دون المثلى التي اتخذت مبكراً. وفي مقارنة بين الطريقتين؛ حصل على الحل الأمثل لعشرة آبار ب ase VA إنجاز في كتل صخرية متعددة م وكانت الجودة الكلية أكبر بنسبة 9647 من الطريقة المكلفة. وكان للحل الأمثل تكلفة أكبر بنسبة 96179 بافتراض أن إنجاز ثان في بثشر يستلزم نصف تكلفة البثر.
شكل a تتضمن المرحلة الثانية في إجراء تحديد موقع وشكل البثسر تحديد JIC الآبار التي حدد موقعها في المرحلة الأولى. وتتضمن هذه المرحلة صيغة رياضية جديدة تعيتن ممر بثر أفقي و/أو منتحرف بدرجة عالية باستخدام مجموعة الإنجازات العمودية المحددة مسبقآً كنقطة بدء. ويتمثل الهدف في زيادة إنتاجية الهيدروكربون الكلية؛ وفي تحقيق ذلك؛ لتحديد إذا كان من الممكن استغلال مناطق مغلة منفصلة يتطلب كل منها إنتاج آبار منجزة بشكل عمودي ومنفرد باستخدام عدد قليل من الآبار.
Vo ويبين الشكل ]1[ بئرآً منحرفاً يتصل بمواقع عالية الجودة في المكمن. ومن المدرك أن المشكلة تتمثل في تعيين مسار إنجاز منحرف بمعرفة توزيع فراغي ثلاثي الأبعاد لنقاط سمتية ذات قيم جودة مقترنة؛ أي في حيز مكعب الشكل حول موقع إنجاز عمودي مختار مسبقاً. وتشمل قيود المشكلة أقضى امتداد fll وأقصضى زاوية انحناء وأدنى مسلفة بين الإنجازات داخل الكتلة الصخرية.
0 وتزود النظرية البيانية نماذج مفيدة لهذه المشكلة. ويتكون رسم بياني (VE) =G من مجموعة منتهية غير خالية من الذروات 7“- (me oF OV) ومجموعة من الحواف 8- fo, «cep coy} التي تعتبر عناصرها مجموعة جزئية من V بالحجم 7ء أي أن =e (ز, ) Cus و ز ينتميان إلى . ويطلق عادة على عناصر V "عقد". وبذلك تزود الرسوم البيانية تقنية ملائمة لتحديد أزواج معينة من المجموعات. وتتمثل مميزة مهمة
vo لرسم بياني في "مسلك"”؛ وهو عبارة عن تسلسل متصسل من الحواف. والتعريف التقليدي
للمسلك هو تسلسل عقدء ١ > 1 7260 Vp eV حيث (Vig, Vi) تنتمسي ل 3 بالنسبة ل 1< ke) ويشار إلى المسلك بممر إذا لم يوجد تكرارات للعقد. ويشار إلى العقدة 7 بالعقدة الأصلية والعقدة ,© بالعقدة Ai" all والعقد (Vig ¢ eer Vy) بالعقد المتوسطة.
° ويمكن تصور النقاط السمتية لخريطة ثلاثية الأبعاد محددة كعقد للرسم البياني. pi مع كل عقدة قيمة معينة لخاصية المكمن المنشودة. وقد يكون مسار بثر أفقي Jf منتحرف عبارة عن ممر يصل بين مجموعة جزئية من هذه العقد. وتمثل العقدة الأصلية في هذا الممر بداية إنجاز وتمثل العقدة المعينة نهايته. وتتاظر العقد المتوسطة المناطق المغلة المتلامسة عن طريق مسار «All وتمثل الحواف المناظرة أجؤاء
الإنجاز للبثر. وفي هذه الحالة؛ تناظر مهمة تحديد ممر إنجاز متحرف أمثل حل مشكلة الاستمثال الذي يختار أفضل ممر؛ أي أفضل مجموعة جزئية من العقد التي تسهم خواص مكمتنها بأعلى قيمة محتملة للدالة المنشودة. ويمثل هذا التسلسل للعقد الامتداد والمسار والسمت النموذجية لبثر أفقي أو منحرف بدرجة عالية له أقصى منطقة تماس أو إنتاجية ضمن وحدة تخزين البيانات ثلاثية الأبعاد المحددة.
Yo وبالإضافة إلى ذلك؛ ينبغي أن يضمن أن يكون JZ S35 البثر ممككاآً من الناحية العملية. وتشمل الأنواع الثلاثة للقيود المحتملة في التطبيق العملي المأخوذة بعين الاعتبار: أن تكون المسافة بين الآبار اكبر من Din وأن يكون سمت ممر الإنجاز ض ضمن انتحراف معينن عن الأفق؛ وأن يكون الامتداد الكلي لممر الإنجاز Gada الحدود الفيزيائية لتقنيات الحفر الراهنة. ويبين الشكل ]1[ رسمآً تخطيطياً لمكونات
,+ الصيغة. وفيما يلي سيتطرق لهذه المكونات واحدآً واحداً. وللمحافظة على تعققّد المشسكلة ضمن حدود ممكنة عملياً؛ تأخذ الآبار المنحرفة بعين الاعتبار واحدآً Jas وتفرض قيود المسافة بين الآبار المنحرفة بعد استمثال المسار بحذف كل النقاط السمتية ضمن حيز مكعب الشكل في Din dg
Yv
حول مسارات آبار سابقة من اعتبارات أخرى. ويعتمد الإجبراء المتتابع هذا على
الترتيب الذي تشكّل به الآبارء وقد يؤدي إلى حلول دون المثلى. ض ولضمان أن يصمم إنجاز البثر بتطبيق فعلي؛ ينبغي التأكيد على أن يكون سمت المسار ضمن زاوية الانحراف المسموح بها عن الزاوية VAY وبعبارة أخرى. ه ينبغي أن تكون زاوية الانحناء بين حواف الرسم البياني أقل من قيمة محددة
مسبقآء Jo Sta ومن الملاحظ أن طريقة لصياغة هذه القيود تبدأ بتحديد متغيرات GE تمثل وجود أو عدم وجود النقاط السمتية (العقد) في المسار النهائي. بيد أنه يفضل ٠ تحديد متغيرات JS LE حواف الرسم البياني. ومن الملاحظ كذلك أن الرسم - البياني غير موجه أي أن الحواف (iL) (:, ز) AL Aa ونتيجة لذلك؛ af بياني cally من Made ينبغي الأخذ بعين الاعتبار الحواف الواضحة MC, ولصياغة القيودء؛ تحدد أولا الزاوية بين كل زوج من الحواف في as Jl البياني. وفي هذه المرحلة؛ يلجأ للصيغ من الهندسة التحليلية المجسمة لتحديد ve جيب التمام للزاوية. ويأخذ بعين الاعتبار أي حافتين (أو بشكل مكافئ ثلاث عقد) في رسم بياني. وتعرف الإحداثيات (x,y, z) للعقد ومن ثم يمكن حساب المسافة المستقيمة بين العقد (امتداد الحواف). ومن ثم يمكن تحديد قيم جيب ED الاتجباءه للخطوط التي تضم هذه التقاط (الحواف)؛ وأخيرآ باستخدام قيم جيب تمام الاتجاه هذه يمكن حساب جيب تمام الزاوية بين الحافتين. ويمكن كذلك استخدام طرق أخرى © لحساب الزاوية. ويمكن اختبار الزاوية المحسوبة مقابل التفاوت المسموح به المحدد. وإذا تجاوزت الزاوية التفاوت المسموح به؛ فإن زوج الحواف المقترن يكون عبارة عن توليفة غير ممكنة في التطبيق العملي.
ولهدف تمثيل زوج قيود يتعذر تنفيذه رياضياً؛ لتمثشل كل من المجموعتين (W) و (W') نقاط إنجازات ممكنة في حيز حول dy عمودي Jie وليمثشل do gana (WW) vo الأزواج المرتبة للمجموعتين (7) و (177) التي JAS كل العلاقات بين
YA
للمجموعة ١ مجموعة متغيرات ثنائية قيمتّها Y(W,W') إنجازات ممكنة. ويمثل als” المختارة من العلاقات بين نقاط إنجازات ممكنة وصفر في العلاقات الأخرى. ومن ثم عقدي": (py AT من نوع JAE هذا القيد رياضياً dela يمكن 1 77ر07 جز + 707,7 (1 غير قابلين للتنفيذ معاً. وقد يتطلب استخدام هذه المعادلة V(W,W') و Vi(W,W') حيث يكون من قيود من هذا القبيل لضمان صياغة جيدة. وكذلك تتمثل محاولة تحديد Tan عدد كبير يمثقل عدد العقد في رسم بياني. وبما أن التكلقة M تقريباً؛ حيث MP هذه القيود ب فقد يفضل تحديد Md الحسابية لتحديد كل القيود قد تستغرق زمنآً طويلآ لقيم معقولة عدد العقد المأخوذة بعين الاعتبار في وحدة تخزين بيانات ثلاثية الأبعاد لكل مشكلة ٠ تحديد مسار أفقي إلى مجموعة جزئية من العدد الكامل للعقد. ويعتمصد حجم المجموعة . الجزئية هذه على سرعة الكمبيوتر المتوفرة ولكنها تكون غالباً برتبة عدة مثشات. ولصياغة القيود التي تفرض حد أعلى على الامتداد الكلي لإنجاز منحرف؛ يلاحظ |, يمتقل الامتداد LW, W') أنه يمكن حساب امتدادات كل الحواف بشكل مسبق بفرض أن و (1)17,777. وباستخدام نفس الرموز السابقة؛ يمكن صياغة هذا القيد (WW) للعلاقات vo رياضياً في صورة:
SS ا ww > L, (A aD و ...مآ كميات معروفة. وهكذاء؛ إذا وجدت حافة في المسار L(W,W') حيث يساوي ١؛ فإن امتداد تلك الحافة Y(W,W') الأمثلء أي أن متغيرها الثنائي المقترن © سيسهم في الامتداد الكلي للبثر المنجز. ولضمان أن يمل ترتيب العقد المختار بالاستمثال ممرآً في الرسم البياني؛ يوضع قيد للتحقق من أنه لا يوجد تكرار للعقد. ويمكن تحقيق ذلك بفرض قيود بأن درجة العقدة قوس واحد على أبعد تقدير يكون ساق ا على )١( تساوي واحد في الحل النهائي؛ أي أن
Ya قوس واحد على أبعد تقدير يكون موجها بعيدا عن عقدة. ويمكن تمثيل هذه (Y) و date القيود رياضيا في صورة: 2 70777041 ود 3 YEH 5.1 (14) ولزيادة الجودة الكلية لمسار بثر مقدر إلى الحد الأقصسى؛ يفضل التعبير عن o
JA) الدالة المتشودة في صورة مجموع قيم الجودة للعقد المختارة بخوارزمية تمثل مجموعة العقد 7 في (X(W) وهكذا تزود مجموعة إضافية من المتغيرات الثنائية؛ عن طريق Ble BY و X الرسم البياني. وترتبط المجموعتان من المتغيرات الثنائية وصل" إذا وإذا فقط كان القوس المقترن Ala" في X(W) المبرهنة المنطقية: تكون العقدة عند نقاط الإنجاز ١ قيمة X(W) في "حالة وصسل". وهكذا يكون ل Y(W,W) أو 77)9/,77( | ٠ الجودة المققرنة المحددة Q(W) Jails المحتملة المختارة وصفر في أماكن أخرى. مسبقا لنقاط الإنجاز هذه. وقد تكون العبارة "إذا" في المبرهنة أعلاه مكافئة رياضيا للمجموعتين التاليتين من ا المعادلات:
XT YH HW) د xy rw wy (11) oe w 7 ولصياغة العبارة الفرعية "إذا فقط" في المبرهنة؛ يتطلب التأكد فيما إذا كانت
SYW,W) مجموعة الحواف الساقطة أو الموجهة بعيدا عن العقدة 177 غير مختارة؛ أي أن صفر كذلك. ولضمان أن X(W) جميعها صفر ومن ثم تمثل العقدة المقترنة Y(W,W) صفر بالضبط في هذه الحالة؛ تذكر المبرهنة التالية: يزيد عدد العقد في X(W) تمثل © بالضبط عن عدد الحواف. daly ممر ويعتبر هذا صحيحا لكل مسار بثر محدد بالاستمثال. وببلجبراء امتدادء يمكسن ملاحظة أنه عند تشكيل آبار متعددة في نفس الوقت يكون عدد العقد المختارة ناقص عدد صفر للحالة X(W) Jia الحواف المختارة مساويا لعدد الآبار. وتكفل المبرهنة أعلاه أن الموصوفة مسبقا. vo
: وباستخدام هذه الصيغة؛ لا ينبغي التصريح بوضوح أن تكون المتغيرات X(W) من إ: التوع الثنائي وقد يصرح بها في صورة متغير مستمر محدد بين صفر و .١ وتكفشل القيود المبينة أعلاه والمبرهنة المذكورة أعلاه أن تطبق X(W) على قيم صحيحة ° والنموذج النهائي لتحديد مسار بثشر أفققي/منحرف Jl في نموذج مكمن ثلاتي الأبعاد هو: Somyx@w) did Si) 7 مع الخضوع للقيود 2 ) 3 7777 wv SS 70777 > 1 (Y7) we. الا د SS YW FIL FF) Vo ا ٍ ا )1°( ومجيدده زرب 24 77ر7 ل LWW) raw) (MN) د رس C. : : بل 70 وروم ردص Yr =n, (PM) - 2 امك 7 ص w e ob (YY) رحج (Fv) re 0S X()<1
ويبين الشكل [V] طريقة مفضلة لتحديد مواقع إنجاز آبار أفقية/منحرففة مثلى. وفي الخطوات (Yr) إلى )£ oY تسترجع المجموعة المرتبة لجبودة مكمن ثلاثي الأبعاد والمجموعة المرتبة من الكتل الصخرية. وتسترجع مواقع الآبار العمودية من مرحلة تحديد مواقع الآبار العمودية في الخطوة (Vo) وتحمسّل القيود في الخطوة (Vo A) وتشمل
م القيود أقصى امتداد all وأقصى عدد للآبار الأفقية/المتحرفة وأقصى زاوية انحناء. وتشمل أمثلة قيود أخرى يمكن تطبيقها AIK أدنى مسافة من سطح التلامس مع الماء أو الغاز والتنفيس العمودي الكلي المسموح به وتقييد البثر لينحدر بشكل دائم إلى الأمسفل أو إلى الأعلى من موقع بدء والمسافة من المنصة والمسافة من صدع ورأس المال الكلي المتوفر.
0 وفي الخطوة (١٠7)؛ تجد الطريقة نقطة الإنجاز العمودي غير المستخدمة التي لها Je] جودة. وقد تختار أي خلية كتلة صخرية في عمود الخلايا حيث يحدد موقع بشسر عمودي كنقطة إنجاز عمودية. وتعتبر تلك الخلية غير مستخدمة إذا لم تسهم بجودة نقطة الإنجاز المختارة مسبقاً.
,| وفي الخطوة (VY) يحدد حيز حجمي حول الخلية غير المستخدمة التي لها
vo أعلى جودة. وللحيز الحجمي نصف قطر يحدد بأقصى امتداد All محدد. وفي
الخطوة )£ (TY تختار مجموعة من نقاط إنجاز محتملة من الحيز الحجمي هذا. ويزال
من النقاط المرشحة كنقاط إنجاز خلايا غير كتلية صخرية وخلايا مستخدمة. وتختار نقاط
الإنجاز المحتملة عشوائياً ويحدد عدد النقاط بعدد أقضصسى إلى حدما (مثل )٠٠١
للمحافظة على إمكانية تدبر التعقيد. ويحدد العدد الأقصسى عن طريق ذاكرة الكمبيوتر
.+ - وسرعة المعالج. ويزيد عدد الحسابات المحلولة مسبقآً بمقدار “« ويزيد عدد المتغيرات الثناثئية
بمقدار 2« ويزيد عدد معادلات القيود بمقدار an’ حيث « يمثل عدد نقاط الإنجاز المحتملة المختارة.
وفي الخطوة (FY) تحسب امتدادات جميع الأقواس بين نقاط الإنجاز المحتملة في المجموعة؛ وتلغى تلك الأقواس التي لها امتدادات أكبر من أقصى امتداد محدد coll Yo وتحسب الزوايا بين كل أزواج الأقواس وتصنف تلك الأزواج التي لها زوايا انحتاء أكبر
YY ؤ من الزاوية المحددة كأزواج غير صالحة. وفي الخطوة (VIA) يمكن إيجاد الحل ض الأمثل للمعادلات من )7١( إلى )"١( باستخدام برمجة خطية بأعداد صحيحة مركبة (MILP) . ويحفظ ممر البثر المنحرف الأمثل. وفي الخطوة )77١( يجرى اختبار لتحديد فيما )13 أمكن الوصول إلى أقصى عدد للآبار الأفقية/المنحرفة. وفي الخطوة (YY) يجرى اختبار لتحديد فيما إذا تبقتت أي نقاط إنجاز عمودية غير مستخدمة. وإذا Ge بثر آخر وبقيت نقطة إنجاز واحدة على الأقل؛ فإن الطريقة تعود للخطوة )0 (FY وفيما عدا ذلك تنتهي الطريقة. وعبّر عن الصيغ بقواعد لغوية وفقاً ل GAMS (نظام صياغة جبري معمم). وحلت الصيغ باستخدام نسخة مماثلة للجهاز التناظري من نوع CPLEX MIP (علامة تجارية مسجلة) على أي رسوم بيانية من نوع سيتيكون اس جي آي اونيكس وباستخدام نسخة مماثلة من الجهاز التناظري .087 (علامة تجارية مسجيلة) على 2 . ويفضل تزويد المستخدم بجهاز بيني بياني لمعالجة وحدات تخزين البيانات وإجراء مكونات تمييز الكتل الصخرية وحساب جودة المكمن وتحديد مواقع آبار عمودية وتحديد مواقع آبار أفقية بشكل منفصل إذا اقتضى الأمر. ويتيح الجهاز البيني للمستخدم أن يختار معايير قطع Ale ومنخفضة وعمليات تنقيب سداسية النقاط أو باستخدام VA أو 7١ نقطة؛ ومتغيرات أخرى مثل قيود نصف قطر التصريف للآبار المقترحة والمسافة بين الآبار وامتداد البثر الأفقي والزاوية السمتية. وستتضح تغييرات وتعديلات عديدة للأُشخاص المتمرسين في التقنية بمجرد إذراك الوصف أعلاه بشكل كامل. فعلى سبيل المثال؛ قد تشكثّل زاوية الانحناء القصوى ؛" بحيث تكون دالة في امتداد القوس؛ مثلاآً ١7 لكل ٠0 متر. ومن المنتشود تفسير عناصر الحماية التالية بحيث تشتمل على كل التغييرات والتعديلات من هذا القبيل.
Claims (1)
- ض . -١ طريقة لتحديد مواقع لعدة آبار » حيث تشتمل الطريقة على: استلام قيمة تقريبية لإنتاجية بثر لكل نقطة صغيرة حجمية من وحدة تخزين ° بيانات خواص مشتقة من تحليل بيانات زلزالية؛ معالجة القيم التقريبية لإنتاجية البثر لتمييز كل صخرية؛ حساب dad جودة مكمن لكل نقطة صغيرة Asean في Js الصخرية ؛و استخدام برمجة بأعداد صحيحة لتحديد موقع نقاط الإنجاز الحجمية الصغيرة التي تزيد مجموع قيم جودة المكمن المقترنة إلى الحد الأعلى مع الخضوع لقيود معينة. ١ ٠ —Y طريقة Gy لمطلب الحماية (١)؛ حيث وحدة تخزين lily الخواص المشتقة من tla | بيانات زلزالية هي عبارة عن وحدة تخزين بيانات ثلاثية الأبعاد لتكوين جيولوجي بترولي له خواص جيولوجية متغايرة وتوزيع مائع غير متجانس. 7أ- طريقة Gy لمطلب الحماية (١)؛ حيث وحدة تخزين البيانات ثلاثية الأبعاد هي عبارة عن وحدة تخزين خواص مشتقة من رسم الخرائط أو تشكيل نموذج إحصائي جيولوجي من بيانات آبار متوفرة. “*- طريقة وفقآً لمطلب الحماية (١)؛ حيث القهيمة التقريبية لإنتاجية البثر هي قيمة 0 من مجموعة قيم تقريبية؛ وتشمل المجموعة المسامية؛ المردود الصافي؛ dala) سماكة الإنفاذية وحجم المسام. ؛- طريقة 84( لمطلب الحماية ) ١ ( « حيث تشمل المعالجة المذكورة للقيم التقريبية لإنتاجية البثر:Yo1و إعادة تحديد كل القيم التقريبية لإنتاجية البئر الأقل من قيمة قطع دنيا مختارة ¢ بصفر؛ ° تحديد أحجام J الصخرية بجمع أحجام النقاط الصغيرة الحجمية المترابطة التي : لها قيم تقريبية لإنتاجية البثر لا تساوي صفر؛ و ل تحديد قيم دليلية للكتل الصخرية لتقليل حجم الكتل الصخرية. )= طريقة Gig لمطلب الحماية (١)؛ حيث تشمل المعالجة المذكورة للقيم التقريبية لإنتاجية A البثر: 1 تعيين كل النقاط الصغيرة الحجمية التي لها قيم تقريبية لإنتاجية بئر dF من م قيمة قطع دنيا مختارة كنقاط غير Aled وتعيين كل النقاط التي لها قيمة تقريبية لإنتاجية Jy ° مساوية ل أو أكبر من قيمة قطع دنيا مختارة كنقاط فعالة؛ ل" تحديد أحجام الكتل الصخرية بجمع أحجام النقاط الصغيرة الحجمية الفعالة المترابطة؛ v A تحديد قيم دليلية للكتل الصخرية لتقليل حجم الكتل الصخرية.١ +- طريقة وفقا لمطلب الحماية (ه)؛ حيث تشمل خطوة حساب قيمة جودة مكمن لنقطة \ صغيرة حجمية محددة:roo جمع القيم التقريبية لإنتاجية بر لكل النقاط الصغيرة الحجمية الفعالة المرتبططة Lally ¢ 4 الصغيرة الحجمية المحددة والتي تقع ضمن نصف قطر تصريف iy ° للنقطة الصغيرة الحجمية المحددة. -١ ١ طريقة وفقآ لمطلب الحماية (١)؛ حيث تشمل خطوة حساب قيمة جودة مكمن \ لنقطة صغيرة حجمية:a 7+ محاكاة مسارات ثلاثية الأبعاد لتركيبة جوالة بشكل عشوائي من النقطة الصغيرة 1 الحجمية المحددة إلى حد؛ حيث يحدد الحد عن طريق أي قيمة في مجموعة Jodi 5 نصف قطر تصريف وحدود كتلة صخرية وحدود غير انسيابية؛ و 1 جمع القيم التقريبية لإنتاجية بثر لكل النقاط الصغيرة الحجمية التي تمس v عن طريق مسار تركيبة جوالة بشكل عشوائي واحد على الأقل. ~A ١ طريقة وفقا لمطلب الحماية ) ١ ؛ حيث يتضمن استخدام برمجة بأعداد صحيحة Y مجموعة من القيسود تشتمل على: أقصى عدد «JL أدنى مسافة بين آبار منجزة v في كتلة صخرية AS ida أقصى مسافة من منصسة بحرية؛ أقصى تكلفة رأسمالية $ للحفر وأدنى مسافة من سطوح grok الماء <u ¢ سطوح التلامس dull للغاز eu Hl 0 الصدوع وحدود تكوين المكمن الأخرى. \ 9- طريقة وفقا لمطلب الحماية ) ٠» ١ حيث يشمل استخدام برمجة بأعداد صحيحة لتحديد Y مواقع نقاط إنجاز صغيرة حجمية:إيجاد أكبر قيمة ل r بيجاد ابر © ] S00 00707 AG) - aX XW) - 3TH NG) [WG 7 WG مع الخضوع للقيود التالية: ¢ , مدلا 1ب مسرا Lj Bd sD كوه مه 50 ض 1 : . | ل 2 م ٍ ّ 7 +7 14717 جا لح نر WG)J’ 2> XP) No ١ : AXF)2YW NG) 7 6: rv 0707,6( « fo} 00s X()<l ٠١ يمثل مجموعة من نقاط G سطحية ممكنة oul مجموعة من مواقع Sa WwW Cua ا QWNG) يمثل كل الإنجازات الصالحة WAG ب صغيرة حجمية في الكتلة الصخرية؛ يمثل قيمة جودة مقترنة بكل إنجاز صالح من هذا القبيل» (7)770 يثشل aay متغيراً Bla XW) متغيرآً ثنائيآً له قيم تشير إلى وجود أو عدم وجود إنجاز؛ ٠ للإشارة إلى وجود أو عدم وجود بر في مجموعة مواقع | لآبار السطحية الممكنة . ‘eo ا 77 » يمثل تكلفة البثر و 8 يمثل تكلفة إنجاز.: تشتمل كذلك على ١ ) لمطلب الحماية a, طريقة -١ ١\f إيجاد نقطة صغيرة حجمية غير مستغلة لها قيمة جودة قصوى؛1 اختيار بشكل عشوائي عدد محدد مسبقآً من النقاط الصغيرة الحجمية ضمن نصف t قطر محدد مسبقآً للنقطة الصغيرة الحجمية غير المستغلة؛° حساب امتدادات الأقواس بين كل أزواج النقاط الصغيرة الحجمية المختارة؛1 حساب الزوايا بين كل أزواج الأقواس المترابطة؛ و7 استخدام برمجة بأعداد صحيحة لتحديد ممر إنجاز بثر متنحرف.| طريقة وفقآ لمطلب الحماية (١٠)؛ تشتمل كذلك على: ى| -١١ voتكرار خطوات الإيجاد والاختيار والحساب والبرمجة بأعداد صحيحة المذكورة )13 7 بقيت نقاط صغيرة حجمية غير مستغلة وإذا لم يتجاوز أقصسى عدد للآبار المنحرفة.-١ ١ طريقة GE لمطلب الحماية (١٠)؛ Gua يشتمل استخدام برمجة بأعداد صحيحة لتحديد ل ممر إنجاز بشر منحرف:YA>. 0 XW) ا إيجاد أكبر قيمة ل # ص مع الخضوع للقيود التالية: 3 YYW) 1 ° z : YW) 17 . لا TT 707777 FF) 5 Los 7 وز. A YW RY+Y,W HL {4 > 180+120i} xX < رج YW.H) 1 2 : 0 Xm YW.) . Fd ١ ١١ TX) - 2717 ساد w wv wr . : 7 \YYW. Ww ) 5 oo ض 0s X)s1 Tv كل منهما مجموعة نقاط إنجاز محتملة في حيز حول بشر Sig WW و W حيث ٠ de gana يمثل XW) lad) قيمة جودة مقترنة بكل نقطة Flay QW) عمودية منجزة؛ Vo كل Jia (wW,w' ) ¢ متغيرات مرتبة تحدد للإشارة إلى وجود أو عدم وجود كل إنجاز ١ Ql jie يمثل مجموعة Y(WW) W' العلاقات بين نقاط الإنجاز المحتملة في 177 و 7YaL(W,W') ثنائية مرتبة تشير إلى العلاقات المختارة بين نقاط الإنجاز الممكنة VA Jia tol أقصى امتداد محدد مسبقا و Slag Lig امتداد مقترن بكل علاقة؛ Jig ١ التفاوت الزاوي المسموح به المحدد مسبقا. Y.-١" ١ طريقة لحساب قيمة جودة مكمن لخلية في وحدة تخزين بيانات زلزالية ثلاثية Gus ala \ تشتمل الطريقة على:1 محاكاة عدد محدد مسبقا بتركيبة جوالة بشكل عشوائي ثلاثية الأبعاد من الخلية إلى aa ¢ حيث يحدد الحد عن طريق قيود تشمل نصف قطر تصريف وحوود كتلة ° صخرية؛ و1 جمع القيم التقريبية لإنتاجية 58 لكل الخلايا الموجودة في مسار تركيبة جوالة Sy 7 عشوائي واحد على الأقل.VE ١ طريقة Gag لمطلب الحماية (OF) حيث القيمة التقريبية لإنتاجية بثر هي قية Y من مجموعة af تقريبية؛ وتشمل المجموعة المسامية؛ المردود الصافيء الإنفاتية؛ 1 سماكة الإتفاذية وحجم المسام.Yo ١ طريقة لتمييز كتل صخرية من وحدة تخزين بيانات؛ حيث تشتمل الطريقة على:v اختيار من وحدة تخزين البيانات خاصية كقيمة تقريبية لإنتاجية بثرءv تشكيل مجموعة مرتبة لأرقام الكثل الصخرية ذات عناصر مناظرة للخلايا في وحدة ¢ تخزين البيانات» حيث يحدد للعناصر المناظرة LOAN في وحدة تخزين البيانات التي لها ° قيم خواص تقل عن قيمة قطع مختارة قيمة مؤشرة أو لى ويحدد لكل الخلايا 1 المتبقية قيمة مؤشرة ثانية؛ل بحث بشكل منهجي في المجموعة المرتبة لأرقام الكتل الصخرية عن عناصسر لها A القيمة المؤشرة الثانية وعن أي عنصر راهن وجد أن له القيمة المؤشرة الثانية؛3 زيادة alae الكتل الصخرية؛١ تحديد قيمة عداد الكتل الصخرية للعنصر الراهن؛ و ١ تنفيذ حلقة لتحديد قيمة عداد الكتل الصخرية لكل العناصر المرتبطة بالعنصر الراهن. -١١ ١ طريقة Ty لمطلب الحماية o10) حيث تشتمل خطوة تنفيذ الحلقة المذكورة على: ل تمهيد مجموعة مرتبة من عناصر معينة إلى الصفر؛ و تمهيد عداد أول لعناصر معينة وعداد ثان لعناصر معينة؛ ¢ تحديد لعضو أول في المجموعة المرتبة من عناصر معينة موقع العنصر الراهن ‘ ِ ضبط موقع راهن بحيث يناظر عضو في المجموعة المرتبة من عناصر معينة مشار 1 إليه عن طريق العداد الثاني لعناصر معينة؛ v ولكل عنصر مجاور للموقع الراهن له القيمة المؤشرة الثانية: A تحدد قيمة عداد Ji الصخرية للعنصر المجاور ¢ 4 تجرى زيادة للعداد الأول لعناصر معينة ١ يعين موقع العنصر المجاور لعضو في المجموعة المرتبة من عناصر معينة مار إليه ١ عن طريق العداد ١ لأول لعناصر معينة؛ و أ تجرى زيادة للعداد الثاني لعناصر معينة. -١7 ٠ طريقة Bay لمطلب الحماية (17)؛ Cum تشمل العناصر المجاورة كل العناصر التي تتقاسم سطح مع العنصر عند الموقع الراهن. =A \ طريقة وفقاً لمطلب الحماية (ل7 ١ ل حيث تشمل العناصر المجاورة كذلك كل العخاصر \ التي تشترك في حافة مع العنصر عند الموقع الراهن. -١9 ١ طريقة وفقاً لمطلب الحماية ١ A) 1 حيث تشمل العناصر المجاورة كذلك كل العخاصر 7 التي تشترك في نقطة صغيرة حجمية مع العنصسر عند الموقع الرا هن.١ ١ طريقة وفقا لمطلب الحماية (١١)؛ تشتمل كذلك على: —Y ١ ل تحديد حجم لكل كتلة صخرية؛ و تصنيف الكتل الصخرية بترتيب حجمي متناقص. 1 الخاصية هي واحدة من مجموعة خواص Cus ؛)١١( طريقة وفقا لمطلب الحماية -7١ ١ تشمل المسامية؛ المردود الصافيء الإنفاذية؛ سماكة الإنفاذية وحجم المسام. Y طريقة لتحديد ممر لبثر منحرف؛ حيث تشمل الطريقة: —YY ١ استلام قيمة تقريبية لإنتاجية بثر لكل نقطة صغيرة حجمية من وحدة تخزين Y بيانات زلزالية؛ Y معالجة القيم التقريبية لإنتاجية البثر لتمييز الكتل الصخرية؛ حساب قيمة جودة المكمن لكل نقطة صغيرة حجمية في الكتل الصخرية؛ و ° إيجاد نقطة صغيرة حجمية غير مستغلة لها قيمة جودة قصوى تقل عن قيمة الجودة 1 : آلا لموقع البثشر المختار؛ من نقاط صغيرة حجمية ضمن نصسف Lie اختيار بشكل عشوائي عدد محدد A قطر محدد مسبقا للنقطة الصغيرة الحجمية غير المستغلة؛ 4 حساب امتدادات الأقواس بين كل أزواج النقاط الصغيرة الحجمية المختارة؛ ye حساب الزوايا بين كل أزواج الأقواس المترابطة؛ و ١ استخدام برمجة بأعداد صحيحة لتحديد ممر إنجاز بثر منحرف يزيد مجموع قيم " الجودة إلى الحد الأقصى. 9 حيث يتضمن استخدام برمجبة بأعداد صحيحة YY) طريقة وفقا لمطلب الحماية YY ١ لتحديد ممر إنجاز بثشر منحرف مجموعة من القيود تشمل: \ أدنى مسافة بين الإنجازات في كتلة صخرية مشتركة؛ أقصى انحراف عن الزاوية ٍٍ : الخطية على مسافة محددة؛ أقصى امتداد للبثر وأدنى مسافة من سطوح تلامس 6Ly is يمثل كل منهما مجموعة نقاط إنجاز محتملة في حيز حول W' و W حيث eo يمثل مجموعة X(W) يمثل قيمة جودة مقترنة بكل نقطة إنجازء QW) ب عمودي منجزء؛ كل Jay (WW) متغيرات مرتبة تحدد لتشير إلى وجود أو عدم وجود كل إنجاز» 9 يمثل مجموعة متغيرات Y(WW') W's 17 العلاقات بين نقاط الإنجاز المحتملة في 7 LW, W') ثنائية مرتبة تشير إلى العلاقات المختارة بين نقاط الإنجاز الممكنة؛ VA tol أقصى امتداد محدد مسبقاً و Jia 1... يمثل امتداد مقترن بكل علاقة من العلاقات؛ ١ve يمتل التفاوت الزاوي المسموح به المحدد مسبقاً.£Y الماء والزيت؛ سطوح التلامس البينية للغاز والزيت؛ الصدوع وحدود تكوين المكمن : 0 الأخرى. 1 يشمل استخدام برمجة بأعداد صحيحة لتحديد Cus (V1) طريقة وفقآ لمطلب الحماية —Y ١ منحرف: py ممر إنجاز Y Som x@) إيجاد أكبر قيمة ل r Ww مع الخضوع للقيود التالية: 3 Y YW HF) > 1 ١ z د 707,77( 6 1 1 ٍ SY 7070+ 1007,77( S Lo ١ ١ wow LOW WY+Y, (FW) ااا BE (Cs F5 0% (CA 0 TN Z XW 2Y rw wy ل w SX@) - 22070777 = No " 122 ww Yo < fol) " ى| : 0s X()<1 ل
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/399,857 US6549879B1 (en) | 1999-09-21 | 1999-09-21 | Determining optimal well locations from a 3D reservoir model |
Publications (1)
Publication Number | Publication Date |
---|---|
SA01210708A true SA01210708A (ar) | 2005-12-03 |
Family
ID=23581250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SA01210708A SA01210708A (ar) | 1999-09-21 | 2001-02-06 | تحديد مواقع آبار مثلى في نموذج مكمن ثلاثي الابعاد |
Country Status (13)
Country | Link |
---|---|
US (1) | US6549879B1 (ar) |
EP (1) | EP1389298B1 (ar) |
CN (1) | CN1421009A (ar) |
AT (1) | ATE500486T1 (ar) |
AU (1) | AU777657B2 (ar) |
BR (1) | BR0014186A (ar) |
CA (1) | CA2384810C (ar) |
DE (1) | DE60045693D1 (ar) |
EA (1) | EA004217B1 (ar) |
MX (1) | MXPA02003097A (ar) |
NO (1) | NO326435B1 (ar) |
SA (1) | SA01210708A (ar) |
WO (1) | WO2001023829A2 (ar) |
Families Citing this family (148)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6853921B2 (en) | 1999-07-20 | 2005-02-08 | Halliburton Energy Services, Inc. | System and method for real time reservoir management |
US6980940B1 (en) * | 2000-02-22 | 2005-12-27 | Schlumberger Technology Corp. | Intergrated reservoir optimization |
US6560501B1 (en) * | 2000-03-07 | 2003-05-06 | I2 Technologies Us, Inc. | System and method for collaborative batch aggregation and scheduling |
US20020082811A1 (en) * | 2000-03-17 | 2002-06-27 | Honjas William A. | Optimization apparatus, system, and method of use and doing business |
US6978210B1 (en) * | 2000-10-26 | 2005-12-20 | Conocophillips Company | Method for automated management of hydrocarbon gathering systems |
KR100624298B1 (ko) * | 2000-12-22 | 2006-09-13 | 주식회사 하이닉스반도체 | 플래쉬 메모리 셀의 센싱 회로 |
US7761270B2 (en) * | 2000-12-29 | 2010-07-20 | Exxonmobil Upstream Research Co. | Computer system and method having a facility management logic architecture |
US7277836B2 (en) * | 2000-12-29 | 2007-10-02 | Exxonmobil Upstream Research Company | Computer system and method having a facility network architecture |
FR2842321B1 (fr) * | 2002-07-11 | 2008-12-05 | Inst Francais Du Petrole | Methode pour contraindre un champ de permeabilite heterogene representant un reservoir souterrain par des donnees dynamiques |
US7317989B2 (en) * | 2001-05-15 | 2008-01-08 | Baker Hughes Incorporated | Method and apparatus for chemometric estimations of fluid density, viscosity, dielectric constant, and resistivity from mechanical resonator data |
FR2831917B1 (fr) * | 2001-11-08 | 2004-01-02 | Schlumberger Services Petrol | Procede de determination de la variation de la permeabilite relative a au moins un fluide d'un reservoir contenant des fluides en fonction de la saturation en l'un d'entre eux |
US7283941B2 (en) * | 2001-11-13 | 2007-10-16 | Swanson Consulting Services, Inc. | Computer system and method for modeling fluid depletion |
US7584165B2 (en) * | 2003-01-30 | 2009-09-01 | Landmark Graphics Corporation | Support apparatus, method and system for real time operations and maintenance |
US7200540B2 (en) * | 2003-01-31 | 2007-04-03 | Landmark Graphics Corporation | System and method for automated platform generation |
US7096172B2 (en) * | 2003-01-31 | 2006-08-22 | Landmark Graphics Corporation, A Division Of Halliburton Energy Services, Inc. | System and method for automated reservoir targeting |
US6810332B2 (en) * | 2003-01-31 | 2004-10-26 | Chevron U.S.A. Inc. | Method for computing complexity, confidence and technical maturity indices for reservoir evaluations |
FR2852710B1 (fr) * | 2003-03-18 | 2005-04-29 | Inst Francais Du Petrole | Methode pour former rapidement un modele stochastique representatif de la distribution d'une grandeur physique dans un milieu heterogene par une selection appropriee de realisations geostatistiques |
CA2526576A1 (en) * | 2003-05-22 | 2004-12-02 | Schlumberger Canada Limited | Method for prospect identification in asset evaluation |
US7584086B2 (en) | 2003-09-30 | 2009-09-01 | Exxonmobil Upstream Research Company | Characterizing connectivity in reservoir models using paths of least resistance |
NL1024444C2 (nl) * | 2003-10-03 | 2005-04-08 | J O A Beheer B V | Werkwijze, inrichting, computerprogramma en gegevensdrager voor het met een digitale verwerkingseenheid modelleren van een meerdimensionale heterogene structuur. |
US7725302B2 (en) * | 2003-12-02 | 2010-05-25 | Schlumberger Technology Corporation | Method and system and program storage device for generating an SWPM-MDT workflow in response to a user objective and executing the workflow to produce a reservoir response model |
GB2410550B8 (en) * | 2003-12-04 | 2008-10-01 | Schlumberger Holdings | Fluids chain-of-custody |
US7672818B2 (en) * | 2004-06-07 | 2010-03-02 | Exxonmobil Upstream Research Company | Method for solving implicit reservoir simulation matrix equation |
WO2006036389A2 (en) * | 2004-09-10 | 2006-04-06 | Exxonmobil Upstream Research Company | Numerical modeling evaluation of basin sedimentation properities |
WO2006065915A2 (en) * | 2004-12-14 | 2006-06-22 | Services Petroliers Schlumberger | Geometrical optimization of multi-well trajectories |
US7526930B2 (en) * | 2005-04-22 | 2009-05-05 | Schlumberger Technology Corporation | Method system and program storage device for synchronizing displays relative to a point in time |
US7565243B2 (en) | 2005-05-26 | 2009-07-21 | Exxonmobil Upstream Research Company | Rapid method for reservoir connectivity analysis using a fast marching method |
EA200800434A1 (ru) * | 2005-07-27 | 2008-10-30 | Эксонмобил Апстрим Рисерч Компани | Моделирование скважины, связанное с добычей углеводородов из подземных формаций |
EP1922669A2 (en) * | 2005-07-27 | 2008-05-21 | ExxonMobil Upstream Research Company | Well modeling associated with extraction of hydrocarbons from subsurface formations |
MX2007016574A (es) * | 2005-07-27 | 2008-03-04 | Exxonmobil Upstream Res Co | Modelaje de pozo asociado con extraccion de hidrocarburos a partir de yacimientos subterraneos. |
EP1999492A4 (en) * | 2006-01-20 | 2011-05-18 | Landmark Graphics Corp | DYNAMIC PRODUCTION MANAGEMENT SYSTEM |
WO2008028122A2 (en) * | 2006-09-01 | 2008-03-06 | Chevron U.S.A. Inc. | History matching and forecasting in the production of hydrocarbons |
US9043188B2 (en) * | 2006-09-01 | 2015-05-26 | Chevron U.S.A. Inc. | System and method for forecasting production from a hydrocarbon reservoir |
WO2008083004A2 (en) * | 2006-12-28 | 2008-07-10 | Chevron U.S.A. Inc. | History matching and forecasting of hydrocarbon-bearing reservoirs utilizing proxies for likelihood functions |
US8005658B2 (en) * | 2007-05-31 | 2011-08-23 | Schlumberger Technology Corporation | Automated field development planning of well and drainage locations |
US8046314B2 (en) * | 2007-07-20 | 2011-10-25 | Schlumberger Technology Corporation | Apparatus, method and system for stochastic workflow in oilfield operations |
WO2009032416A1 (en) * | 2007-09-07 | 2009-03-12 | Exxonmobill Upstream Research Company | Well performance modeling in a collaborative well planning environment |
US8139062B2 (en) * | 2007-09-12 | 2012-03-20 | Schlumberger Technology Corporation | Method and system for displaying a map using a projected coordinate system |
WO2009075945A1 (en) | 2007-12-13 | 2009-06-18 | Exxonmobil Upstream Research Company | Parallel adaptive data partitioning on a reservoir simulation using an unstructured grid |
EP2223157A4 (en) | 2007-12-13 | 2016-12-07 | Exxonmobil Upstream Res Co | ITERATIVE TANK SURVEILLANCE |
CA2705277C (en) * | 2007-12-18 | 2017-01-17 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-d and 3-d heterogeneous data |
US8751208B2 (en) * | 2007-12-20 | 2014-06-10 | Shell Oil Company | Method for producing hydrocarbons through a well or well cluster of which the trajectory is optimized by a trajectory optimisation algorithm |
AU2008340399B2 (en) | 2007-12-21 | 2013-09-26 | Exxonmobil Upstream Research Company | Method and apparatus for analyzing three-dimensional data |
WO2009085395A1 (en) * | 2007-12-31 | 2009-07-09 | Exxonmobil Upstream Research Company | Methods and systems for determining near-wellbore characteristics and reservoir properties |
US7894991B2 (en) * | 2008-02-01 | 2011-02-22 | Schlumberger Technology Corp. | Statistical determination of historical oilfield data |
CA2710809C (en) * | 2008-03-10 | 2017-06-13 | Exxonmobil Upstream Research Company | Method for determining distinct alternative paths between two object sets in 2-d and 3-d heterogeneous data |
CN101266299B (zh) * | 2008-04-14 | 2011-03-30 | 林昌荣 | 利用地震数据体结构特征预测油气的方法 |
AU2009238481B2 (en) * | 2008-04-22 | 2014-01-30 | Exxonmobil Upstream Research Company | Functional-based knowledge analysis in a 2D and 3D visual environment |
CA2717514C (en) | 2008-05-05 | 2016-07-26 | Exxonmobil Upstream Research Company | Systems and methods for connectivity analysis using functional objects |
US20100082509A1 (en) * | 2008-09-30 | 2010-04-01 | Ilya Mishev | Self-Adapting Iterative Solver |
CN102138146A (zh) * | 2008-09-30 | 2011-07-27 | 埃克森美孚上游研究公司 | 使用并行多级不完全因式分解求解储层模拟矩阵方程的方法 |
WO2010039317A1 (en) * | 2008-10-01 | 2010-04-08 | Exxonmobil Upstream Research Company | Robust well trajectory planning |
CA2733729C (en) | 2008-10-24 | 2017-01-17 | Exxonmobil Upstream Research Company | Tracking geologic object and detecting geologic anomalies in exploration seismic data volume |
AU2009311619B2 (en) * | 2008-11-06 | 2015-10-01 | Exxonmobil Upstream Research Company | System and method for planning a drilling operation |
AU2015268702B2 (en) * | 2008-11-17 | 2016-02-11 | Landmark Graphics Corporation | Systems and methods for dynamically developing wellbore plans with a reservoir simulator |
US8301426B2 (en) | 2008-11-17 | 2012-10-30 | Landmark Graphics Corporation | Systems and methods for dynamically developing wellbore plans with a reservoir simulator |
BRPI1006862B1 (pt) * | 2009-01-13 | 2020-03-17 | Exxonmobil Upstream Research Company | Método e sistema para otimização de tomada de decisão para um poço de hidrocarbonetos, e, sistema associado com a produção de hidrocarbonetos |
US8301382B2 (en) * | 2009-03-27 | 2012-10-30 | Schlumberger Technology Corporation | Continuous geomechanically stable wellbore trajectories |
US10332219B2 (en) * | 2009-03-30 | 2019-06-25 | Landmark Graphics Corporation | Systems and methods for determining optimum platform count and position |
US20100299123A1 (en) * | 2009-05-21 | 2010-11-25 | Schlumberger Technology Corporation | Well placement in a volume |
AU2010290068B2 (en) | 2009-09-01 | 2015-04-30 | Exxonmobil Upstream Research Company | Method of using human physiological responses as inputs to hydrocarbon management decisions |
US8949173B2 (en) * | 2009-10-28 | 2015-02-03 | Schlumberger Technology Corporation | Pay zone prediction |
US8931580B2 (en) | 2010-02-03 | 2015-01-13 | Exxonmobil Upstream Research Company | Method for using dynamic target region for well path/drill center optimization |
US8731872B2 (en) | 2010-03-08 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
EA201201285A1 (ru) | 2010-03-15 | 2013-05-30 | Лэндмарк Грэфикс Корпорейшн | Системы и способы для определения местоположения горизонтальных скважин в пределах заданных границ |
US8731887B2 (en) | 2010-04-12 | 2014-05-20 | Exxonmobile Upstream Research Company | System and method for obtaining a model of data describing a physical structure |
US8727017B2 (en) | 2010-04-22 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for obtaining data on an unstructured grid |
US8731873B2 (en) | 2010-04-26 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
US8532968B2 (en) * | 2010-06-16 | 2013-09-10 | Foroil | Method of improving the production of a mature gas or oil field |
CA2801382C (en) | 2010-06-29 | 2018-12-18 | Exxonmobil Upstream Research Company | Method and system for parallel simulation models |
US8731875B2 (en) | 2010-08-13 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
US9593558B2 (en) | 2010-08-24 | 2017-03-14 | Exxonmobil Upstream Research Company | System and method for planning a well path |
US9229129B2 (en) * | 2010-12-10 | 2016-01-05 | Conocophillips Company | Reservoir geobody calculation |
CA2823017A1 (en) | 2011-01-26 | 2012-08-02 | Exxonmobil Upstream Research Company | Method of reservoir compartment analysis using topological structure in 3d earth model |
AU2011360212B2 (en) | 2011-02-21 | 2017-02-02 | Exxonmobil Upstream Research Company | Reservoir connectivity analysis in a 3D earth model |
CA2823263A1 (en) * | 2011-05-09 | 2012-11-15 | Chevron U.S.A. Inc. | System and method for determining properties of a hydrocarbon reservoir based on production data |
FR2979724B1 (fr) | 2011-09-06 | 2018-11-23 | Ifp Energies Now | Procede d'exploitation d'un gisement petrolier a partir d'une technique de selection des positions des puits a forer |
US9223594B2 (en) | 2011-07-01 | 2015-12-29 | Exxonmobil Upstream Research Company | Plug-in installer framework |
US8731891B2 (en) * | 2011-07-28 | 2014-05-20 | Saudi Arabian Oil Company | Cluster 3D petrophysical uncertainty modeling |
US20130231901A1 (en) * | 2011-09-15 | 2013-09-05 | Zhengang Lu | Well pad placement |
FR2987149B1 (fr) * | 2012-02-16 | 2014-10-31 | IFP Energies Nouvelles | Procede d'exploitation d'un gisement a partir d'une technique de selection des positions de puits a forer |
FR2989200B1 (fr) * | 2012-04-10 | 2020-07-17 | IFP Energies Nouvelles | Procede de selection des positions de puits a forer pour l'exploitation d'un gisement petrolier |
WO2013169429A1 (en) | 2012-05-08 | 2013-11-14 | Exxonmobile Upstream Research Company | Canvas control for 3d data volume processing |
EP2856387B1 (en) | 2012-05-30 | 2023-03-08 | Landmark Graphics Corporation | System and method for reservoir simulation optimization |
MX2014014438A (es) * | 2012-05-31 | 2015-05-11 | Landmark Graphics Corp | Sistemas y metodos para el optimo posicionado de plataformas de perforacion. |
BR112014031922B1 (pt) | 2012-06-18 | 2022-03-15 | Technological Resources Pty. Limited | Sistemas e métodos para processar dados geofísicos |
US20140005996A1 (en) * | 2012-06-28 | 2014-01-02 | Schlumberger Technology Corporation | Interactive and three-dimensional well path design |
US9970284B2 (en) | 2012-08-14 | 2018-05-15 | Schlumberger Technology Corporation | Downlink path finding for controlling the trajectory while drilling a well |
US9183182B2 (en) * | 2012-08-31 | 2015-11-10 | Chevron U.S.A. Inc. | System and method for determining a probability of well success using stochastic inversion |
US9229910B2 (en) * | 2012-10-26 | 2016-01-05 | Schlumberger Technology Corporation | Predicting three dimensional distribution of reservoir production capacity |
US20140129296A1 (en) * | 2012-11-07 | 2014-05-08 | Schlumberger Technology Corporation | Method and system for offering and procuring well services |
US9417256B2 (en) * | 2012-12-12 | 2016-08-16 | Repsol, S. A. | System, method and program product for automatically matching new members of a population with analogous members |
EP2904530B1 (en) * | 2012-12-13 | 2018-10-10 | Landmark Graphics Corporation | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties |
US10429545B2 (en) | 2012-12-13 | 2019-10-01 | Landmark Graphics Corporation | System, method and computer program product for evaluating and ranking geobodies using a euler characteristic |
US20140214387A1 (en) * | 2013-01-25 | 2014-07-31 | Schlumberger Technology Corporation | Constrained optimization for well placement planning |
US10048396B2 (en) | 2013-03-14 | 2018-08-14 | Exxonmobil Upstream Research Company | Method for region delineation and optimal rendering transform of seismic attributes |
ES2660432T3 (es) | 2013-06-06 | 2018-03-22 | Repsol, S.A. | Método para evaluar planes de estrategia de producción |
WO2014200685A2 (en) | 2013-06-10 | 2014-12-18 | Exxonmobil Upstream Research Company | Interactively planning a well site |
SG11201510115UA (en) | 2013-07-02 | 2016-01-28 | Landmark Graphics Corp | 3d stadia algorithm for discrete network meshing |
US10689965B2 (en) * | 2013-08-26 | 2020-06-23 | Repsol, S.A. | Field development plan selection system, method and program product |
US11181662B2 (en) | 2013-08-28 | 2021-11-23 | Landmark Graphics Corporation | Static earth model grid cell scaling and property re-sampling methods and systems |
MX2016002396A (es) * | 2013-08-29 | 2016-12-14 | Landmark Graphics Corp | Metodos y sistemas de calibración de modelo terrestre estático mediante el uso de evaluaciones de turtuosidad. |
US9864098B2 (en) | 2013-09-30 | 2018-01-09 | Exxonmobil Upstream Research Company | Method and system of interactive drill center and well planning evaluation and optimization |
US9958571B2 (en) | 2013-12-30 | 2018-05-01 | Saudi Arabian Oil Company | Machines for reservoir simulation with automated well completions and reservoir grid data quality assurance |
RU2669948C2 (ru) * | 2014-01-06 | 2018-10-17 | Геоквест Системз Б.В. | Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности |
SG11201606390QA (en) * | 2014-03-10 | 2016-09-29 | Landmark Graphics Corp | Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence |
US9957781B2 (en) | 2014-03-31 | 2018-05-01 | Hitachi, Ltd. | Oil and gas rig data aggregation and modeling system |
US10062044B2 (en) * | 2014-04-12 | 2018-08-28 | Schlumberger Technology Corporation | Method and system for prioritizing and allocating well operating tasks |
CN105093308B (zh) * | 2014-05-19 | 2017-11-21 | 中国石油化工股份有限公司 | 一种井轨迹设计方法和系统 |
US10192182B2 (en) * | 2014-06-10 | 2019-01-29 | Wellaware Holdings, Inc. | Aerial drone for well-site and signal survey |
FR3023316B1 (fr) * | 2014-07-04 | 2016-08-19 | Ifp Energies Now | Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer |
US9816366B2 (en) | 2014-07-14 | 2017-11-14 | Saudi Arabian Oil Company | Methods, systems, and computer medium having computer programs stored thereon to optimize reservoir management decisions |
US11414975B2 (en) | 2014-07-14 | 2022-08-16 | Saudi Arabian Oil Company | Quantifying well productivity and near wellbore flow conditions in gas reservoirs |
CN104331537A (zh) * | 2014-09-28 | 2015-02-04 | 长江大学 | 基于储层静态因子的井位优化设计方法 |
US10221659B2 (en) * | 2014-10-08 | 2019-03-05 | Chevron U.S.A. Inc. | Automated well placement for reservoir evaluation |
CN104360412B (zh) * | 2014-11-14 | 2017-12-26 | 中国石油大学(北京) | 致密深盆气成藏预测方法和装置 |
US9858484B2 (en) * | 2014-12-30 | 2018-01-02 | Facebook, Inc. | Systems and methods for determining video feature descriptors based on convolutional neural networks |
CN104895550B (zh) * | 2015-06-04 | 2018-03-13 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
CN105257252A (zh) * | 2015-06-08 | 2016-01-20 | 中国石油集团川庆钻探工程有限公司 | 利用测井资料优选页岩气水平井分簇射孔井段的方法 |
US10502047B2 (en) * | 2015-06-30 | 2019-12-10 | Magnetic Variation Services LLC | Reservoir recovery simulation process and system |
WO2017015317A1 (en) * | 2015-07-21 | 2017-01-26 | Schlumberger Technology Corporation | Well placement plan optimization |
US9754351B2 (en) * | 2015-11-05 | 2017-09-05 | Facebook, Inc. | Systems and methods for processing content using convolutional neural networks |
CN105545275B (zh) * | 2015-12-31 | 2018-04-06 | 中国石油天然气股份有限公司 | 气田气井部署方法和装置 |
WO2017120447A1 (en) * | 2016-01-08 | 2017-07-13 | Nature Conservancy, The | Techniques for positioning energy infrastructure |
FR3046810B1 (fr) * | 2016-01-15 | 2018-01-26 | IFP Energies Nouvelles | Procede de production d'hydrocarbures comportant un index de productivite des puits sous effet thermique |
US10167703B2 (en) * | 2016-03-31 | 2019-01-01 | Saudi Arabian Oil Company | Optimal well placement under constraints |
US10941635B1 (en) | 2016-06-27 | 2021-03-09 | East Daley Capital Advisors, Inc | Optimization computer program and method |
US10482202B2 (en) | 2016-06-30 | 2019-11-19 | The Procter & Gamble Company | Method for modeling a manufacturing process for a product |
US10060227B2 (en) | 2016-08-02 | 2018-08-28 | Saudi Arabian Oil Company | Systems and methods for developing hydrocarbon reservoirs |
US10605055B2 (en) * | 2016-09-15 | 2020-03-31 | Baker Hughes, A Ge Company, Llc | Integrated hydrocarbon fluid distribution modeling |
US10678967B2 (en) * | 2016-10-21 | 2020-06-09 | International Business Machines Corporation | Adaptive resource reservoir development |
US11740384B2 (en) | 2016-12-09 | 2023-08-29 | Schlumberger Technology Corporation | Field operations neural network heuristics |
WO2018213483A1 (en) * | 2017-05-18 | 2018-11-22 | Conocophillips Company | Resource density screening tool |
WO2019064037A1 (en) | 2017-09-26 | 2019-04-04 | Total Sa | METHOD FOR DEFINING LOCATIONS OF A PLURALITY OF WELLS IN A FIELD, ASSOCIATED SYSTEM AND COMPUTER PROGRAM PRODUCT |
WO2019086938A1 (en) * | 2017-11-06 | 2019-05-09 | Abu Dhabi National Oil Company | Method and system for determining permeability of a porous medium |
US11299964B2 (en) * | 2018-07-03 | 2022-04-12 | Baker Hughes, A Ge Company, Llc | Drilling productive wells |
FR3101660B1 (fr) | 2019-10-03 | 2021-10-08 | Ifp Energies Now | Procédé pour déterminer une trajectoire d'un puits dans un réservoir pétrolier |
US11715034B2 (en) | 2020-01-16 | 2023-08-01 | Saudi Arabian Oil Company | Training of machine learning algorithms for generating a reservoir digital twin |
US11586790B2 (en) | 2020-05-06 | 2023-02-21 | Saudi Arabian Oil Company | Determining hydrocarbon production sweet spots |
US11608734B2 (en) | 2020-05-11 | 2023-03-21 | Saudi Arabian Oil Company | Systems and methods for creating hydrocarbon wells |
US11708754B2 (en) | 2020-05-11 | 2023-07-25 | Saudi Arabian Oil Company | Systems and methods for generating a drainage radius log |
WO2021236877A1 (en) * | 2020-05-20 | 2021-11-25 | Schlumberger Technology Corporation | Drilling trajectory and steering design optimization based on predicted tool performance |
RU2747019C1 (ru) * | 2020-06-18 | 2021-04-23 | Общество с ограниченной ответственностью "Тюменский нефтяной научный центр" (ООО "ТННЦ") | Способ обоснования технологического режима промысла |
CN112580851B (zh) * | 2020-11-17 | 2024-06-18 | 西安中控天地科技开发有限公司 | 丛式井场抽油机井群错峰开井间抽运行调度方法 |
US11859472B2 (en) | 2021-03-22 | 2024-01-02 | Saudi Arabian Oil Company | Apparatus and method for milling openings in an uncemented blank pipe |
CN113269879B (zh) * | 2021-05-27 | 2021-11-26 | 广东省地震局 | 一种地质体三维模型自动构建方法及装置 |
US11788377B2 (en) | 2021-11-08 | 2023-10-17 | Saudi Arabian Oil Company | Downhole inflow control |
WO2024064077A1 (en) * | 2022-09-19 | 2024-03-28 | Schlumberger Technology Corporation | Training of machine learning models for well target recommendation |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4249776A (en) | 1979-05-29 | 1981-02-10 | Wyoming Mineral Corporation | Method for optimal placement and orientation of wells for solution mining |
US4916616A (en) * | 1986-12-08 | 1990-04-10 | Bp Exploration, Inc. | Self-consistent log interpretation method |
US5012675A (en) | 1989-07-25 | 1991-05-07 | Amoco Corporation | Integrating multiple mappable variables for oil and gas exploration |
US5706194A (en) * | 1995-06-01 | 1998-01-06 | Phillips Petroleum Company | Non-unique seismic lithologic inversion for subterranean modeling |
US5757663A (en) | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
US6035255A (en) | 1997-12-01 | 2000-03-07 | Schlumberger Technology Corporation | Article of manufacturing for creating, testing, and modifying geological subsurface models |
US6070125A (en) | 1997-12-01 | 2000-05-30 | Schlumberger Technology Corporation | Apparatus for creating, testing, and modifying geological subsurface models |
US6044328A (en) * | 1997-12-01 | 2000-03-28 | Schlumberger Technology Corporation | Method for creating, testing, and modifying geological subsurface models |
US6266619B1 (en) * | 1999-07-20 | 2001-07-24 | Halliburton Energy Services, Inc. | System and method for real time reservoir management |
-
1999
- 1999-09-21 US US09/399,857 patent/US6549879B1/en not_active Expired - Lifetime
-
2000
- 2000-09-20 MX MXPA02003097A patent/MXPA02003097A/es unknown
- 2000-09-20 AU AU77061/00A patent/AU777657B2/en not_active Ceased
- 2000-09-20 CA CA002384810A patent/CA2384810C/en not_active Expired - Lifetime
- 2000-09-20 CN CN00814550A patent/CN1421009A/zh active Pending
- 2000-09-20 EP EP00966771A patent/EP1389298B1/en not_active Expired - Lifetime
- 2000-09-20 DE DE60045693T patent/DE60045693D1/de not_active Expired - Lifetime
- 2000-09-20 AT AT00966771T patent/ATE500486T1/de not_active IP Right Cessation
- 2000-09-20 WO PCT/US2000/025804 patent/WO2001023829A2/en active IP Right Grant
- 2000-09-20 EA EA200200393A patent/EA004217B1/ru not_active IP Right Cessation
- 2000-09-20 BR BR0014186-0A patent/BR0014186A/pt not_active Application Discontinuation
-
2001
- 2001-02-06 SA SA01210708A patent/SA01210708A/ar unknown
-
2002
- 2002-03-20 NO NO20021383A patent/NO326435B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
ATE500486T1 (de) | 2011-03-15 |
CA2384810A1 (en) | 2001-04-05 |
NO326435B1 (no) | 2008-12-08 |
WO2001023829A2 (en) | 2001-04-05 |
CA2384810C (en) | 2008-12-02 |
EP1389298A4 (en) | 2005-02-09 |
MXPA02003097A (es) | 2002-12-16 |
CN1421009A (zh) | 2003-05-28 |
EP1389298B1 (en) | 2011-03-02 |
BR0014186A (pt) | 2003-07-29 |
AU777657B2 (en) | 2004-10-28 |
EA004217B1 (ru) | 2004-02-26 |
NO20021383L (no) | 2002-05-15 |
US6549879B1 (en) | 2003-04-15 |
EA200200393A1 (ru) | 2002-12-26 |
EP1389298A2 (en) | 2004-02-18 |
NO20021383D0 (no) | 2002-03-20 |
WO2001023829A3 (en) | 2003-12-04 |
DE60045693D1 (de) | 2011-04-14 |
AU7706100A (en) | 2001-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
SA01210708A (ar) | تحديد مواقع آبار مثلى في نموذج مكمن ثلاثي الابعاد | |
US7783462B2 (en) | Reservoir evaluation methods | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
CN105701319B (zh) | 一种水平井规则开发井网下的沉积微相建模方法 | |
Tolstukhin et al. | Ekofisk 4D seismic-seismic history matching workflow | |
Kamali et al. | 3D geostatistical modeling and uncertainty analysis in a carbonate reservoir, SW Iran | |
AlRassas et al. | CO2 storage capacity estimation under geological uncertainty using 3-D geological modeling of unconventional reservoir rocks in Shahejie Formation, block Nv32, China | |
RU2685005C1 (ru) | Способ и компьютерная система для проектирования размещения кустовых площадок на месторождении | |
Martin et al. | An integrated approach to the modeling of permeability barrier distribution in a sedimentologically complex reservoir | |
Gomes et al. | Quality control of 3D GeoCellular models: examples from UAE carbonate reservoirs | |
Chen et al. | 3D Geomechanical Modeling for an Extra Deep Fractured Carbonate Reservoir, Northwest China | |
Guerreiro et al. | Integrated reservoir characterisation of a fractured carbonate reservoir | |
Gurpinar et al. | Numerical modeling of a triple porosity reservoir | |
Shevko et al. | Application of hybrid DPDP-dFN modeling of fractured carbonate reservoirs | |
CN105259577B (zh) | 一种确定地层界面的角度信息的方法及装置 | |
AU2012396846B2 (en) | System, method and computer program product for evaluating and ranking geobodies using a Euler Characteristic | |
Salahuddin et al. | Static and Dynamic Uncertainty Management for Probabilistic Volumetric and Production Forecast: A Case Study from Onshore Abu Dhabi | |
Budilin et al. | Integrated uncertainty quantification for development planning of a large field | |
Karpov et al. | Local Stress Shadow Effect Analysis in Multistage Hydraulic Fracturing Design Considering Small Drillhole Spacing | |
Chambers et al. | Geologic modeling, upscaling and simulation of faulted reservoirs using faulted stratigraphic grids | |
Feazel et al. | Carbonate reservoir characterization and simulation: From facies to flow units: Report from the March 2004 Hedberg Research Symposium | |
Utama et al. | Case Study: Assessment of Predictive Capability of Reservoir Simulators for Waterfloods in Carbonates: How Realistic is My Simulation Model? | |
Anyanwu et al. | Integrated Reservoir Modeling and Uncertainty Assessment of a Reservoir Complex in the Niger Delta | |
Benton et al. | Integrated Ensemble-Modeling Study Utilizing Well Surveillance Data and Interpreted Contacts for a Mature Malaysian Oilfield | |
Chongrueanglap et al. | Challenges on Building Representative 3D Static Models under Subsurface Uncertainties for a Giant Carbonate Field in Central Luconia, Offshore Sarawak |