SA01220301B1 - طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data - Google Patents
طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data Download PDFInfo
- Publication number
- SA01220301B1 SA01220301B1 SA01220301A SA01220301A SA01220301B1 SA 01220301 B1 SA01220301 B1 SA 01220301B1 SA 01220301 A SA01220301 A SA 01220301A SA 01220301 A SA01220301 A SA 01220301A SA 01220301 B1 SA01220301 B1 SA 01220301B1
- Authority
- SA
- Saudi Arabia
- Prior art keywords
- rock
- liquid
- value
- predictive
- acoustic impedance
- Prior art date
Links
- 239000011435 rock Substances 0.000 title claims abstract description 159
- 239000007788 liquid Substances 0.000 title claims abstract description 157
- 238000000034 method Methods 0.000 title claims abstract description 82
- 230000004044 response Effects 0.000 claims abstract description 43
- 229930195733 hydrocarbon Natural products 0.000 claims description 27
- 150000002430 hydrocarbons Chemical class 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 18
- 239000004215 Carbon black (E152) Substances 0.000 claims description 14
- 238000004519 manufacturing process Methods 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 125000001183 hydrocarbyl group Chemical group 0.000 claims 4
- 238000011156 evaluation Methods 0.000 claims 1
- 230000005477 standard model Effects 0.000 claims 1
- 239000010410 layer Substances 0.000 description 9
- 230000004888 barrier function Effects 0.000 description 8
- 239000000203 mixture Substances 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000013517 stratification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241001247287 Pentalinon luteum Species 0.000 description 1
- 230000001464 adherent effect Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 239000011344 liquid material Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- JNAQOQQXXBMQKT-UHFFFAOYSA-M sodium;3,4,6-trichloro-2-[(2,3,5-trichloro-6-hydroxyphenyl)methyl]phenolate Chemical compound [Na+].OC1=C(Cl)C=C(Cl)C(Cl)=C1CC1=C([O-])C(Cl)=CC(Cl)=C1Cl JNAQOQQXXBMQKT-UHFFFAOYSA-M 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
الملخص: تقدم هنا طريقة للتنبؤ بالقيم الخاصة بخاصية معينة لصخر أو لسائل في الحجم الجيولوجي تحت الأرض subterranean geologic volume . يتم تحديد قيمة زلزالية تجريبية seismic value للمعاوقة الصوتية acoustic impedance لنموذج حجمي يرتبط بالحجم الجيولوجي تحت الأرض . وأيضا يتم تحديد قيمة أولى تنبؤية للخاصية المعينة للصخر أو السائل لنموذج حجمي. القيمة الأولى التنبؤية للمعاوقة الصوتية لنموذج الحجم يتم حسابها من خلال نموذج استجابة response model باستخدام القيمة الأولى التنبؤية للخاصية المعينة للصخر أو السائل حيث أن نموذج الاستجابة يستجيب للتغيرات الحادثة في القيم التنبؤية للخاصية المعينة للصخر أو السائل . يتم مقارنة القيمة الأولى التنبؤية للمعاوقة الصوتية بالقيمة الزلازلية للمعاوقة الصوتية لتحديد الفرق بين القيم التنبؤية والقيم الزلزالية للمعاوقة الصوتية. ويتم ضبط القيمة الأولى التنبؤية للخاصية المعينة للصخر أو السائل حسب هذا الفرق للحصول علىقيمة ثانية تنبؤية للخاصية المعينة للصخر أو السائل حيث تقلل القيمة التنبؤية الثانية من الفرق .،
Description
Y
طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان seismic data باستخدام بيانات خاصة بالزلازل الوصف الكامل : خلفبة الإختراع يرتبط الاختراع الحالي بصفة عامة بطريقة للحصول على نموذج جيولوجي متعدد الأبعاد dimensional geologic model لخزان سائل تحت الأرض subterranean fluid reservoir . وبالإخص يقدم الاختراع طريقة لملء population النموذج الجيولوجي متعدد الأبعاد بقيم كمية quantitative values لخاصية معينة لصخر أو سائل Cun تكون هذه القيم بأقل نسبة minimized tas 70 باستخدام نموذج استجابة بتروفيزيائي petrophysical response model وبيانات زلزالية seismic data . ويسهل النموذج الجيولوجي متعدد . الأبعاد_ Ad من اكتشاف exploration أو إدارة إنتاج ills «U3 production management للهيدروكربونات hydrocarbon- . bearing reservoirs ١ خلال عمليات الكشف عن الهيدروكربونات واستغلالها من البيئة تحت الأرض هناك حاجة مستمرة لتحديد الخزان تحت الأرض بدقة وبمعرفة العمق الحقيقي والمحتوي الهيدروكربوني hydrocarbon content ونفاذية السائل fluid permeability للخزان المحتوي على الهيدروكربون يكون مهما جداً تقليل خطورة الخسارة الاقتصادية وعلى ve العكس_من ذلك زيادة معدل إنتاج الهيدروكربونات في الخزان . وهذه المعلومات بخصوص الخزان تحت الأرض يمكن الحصول علىها عن طريق حفر بئر أو أكثر في الخزانات . ويعطي معدل الحفر Driling rates والعينات المأخوذة من الحفر drill 895 والتغيرات الحادثة في تركيب dub الحفر drilling mud composition
- سجلات تقييم الأداء 5 والتي تتكون عن طريق إمرار معدات تقييم الأداء خلال البئر مصدراً جيداً للمعلومات ؛ Cua تقدم هذا السجلات معلومات قيمة بخصوص خواص الصخر والسائل الموجودين في الخزان تحت الأرض Jie المسامية porosity وتحديد نوع السائل fluid identification وحجم الصخر الطيني shale volume . وتحتوي 0 هذه السجلات أيضاً على المقاومة 7881507117 وأشعة جاماً gamma ray والكثافة density وسرعة الضغط compressional velocity وسرعة القص ٠ shear velocity وسجلات النيوترون neutron logs . وحيث أن سجلات تقييم الأداء تقوم فقط بقياس خواص الصخر والسائل حتى ارتفاع قدم واحد من عمق A ومعظم الخزانات لم يتم حفر آبار فيها فلسوء الحظ لا يمكن للسجلات ٠ إلا تحديد خواص جزء صغير جداً من الخزان . وزيادة على ذلك تؤثر عملية الحفر على الصخر حول عمق البثر بذلك تغير خصائص الصخر وتؤدي إلي حدوث خطأ في القياسات المأخوذة عن طريق أجهزة قياس الأداء وتحليل قلب البثر . ويجب الانتظار لفترة طويلة حتى يمكن قياس خواص الصخر والسائل بدقة عبر كل خزان تحت الأرض وبالأخص لتحديد خواص الصخر والسائل بدقة في مناطق الخزانات التي لم يحفر فيها ve أآبار . Laie تكون بيانات all محدودة ويتم استقراء extrapolated بيانات البئر بعيداً عن عمق J لتحديد خواص الخزان كلها . وتصف تقنيات الاستقراء التقليدي للخزان تحت الأرض كمجموعة من الكتل المرتبة ثلاثية الأبعاد ومتحدة مع بعضها لتكوين نموذج ثلاثي الأبعاد للخزان . ويتم تحديد المحاور و 7 و72 لكل كتلة في كل من الارتفاع ٠ المطلق والطبقات السطحية ويتم استخدام لوغاريتم بحث لتحديد البيانات النسبية بقرب كل كتلة . وبالإضافة J) ذلك يتم تحديد خواص الصخر لكل كتلة بواسطة طرق تقدير مثل
الطرق التي تعتمد على المسافة باستخدام طرق تقدير تعتمد على القيم التقريبية للبيانات وطرق إحصاءات جيولوجية والتي تعتمد على كل من المسافة والاستمرار المكاني لخواص الصخر. وتستخدم أيضاً طرق المسح الزلزالي Seismic surveys لإعطاء معلومات زلزالية عن ٠ _المناطق في الخزانات التي لم يتم حفر آبار فيها . ويتم استخدام eal ضغط Tmpact Jie devices مصدر اهتزازات vibratory sources أو مداقع الغاز J gas guns مداقع الهواء air guns أو إلقاء أثقال weight drops عند سطح الأرض أو عند قاع SA كمصدر زلزالي لتوليد موجات ضغط وقص في الطبقات تحت الأرضية . وهذه الموجات waves تنتقل عبر الطبقات تحت الأرضية وتتعكس على هيئة تغيرات في ٠ عوائق الصوت وتسجل عادة عند سطح الأرض بواسطة أجهزة تسجيل موضوعة بترتيب معين . وهذه البيانات المسجلة يتم معالجتها بواسطة برامج مصممة لتقليل الضوضاء والحفاظ على قوة الانعكاس . وتظهر نتائج المسح الزازالي في صورة بيانات ثلاثية الأبعاد تمثل قياساً مباشراً لأسطح الصخر الذي يحدد الخزان تحت الأرض ؛ وتستخدم هذه البيانات لتقييم ولعمل خرائط ١٠ للتركيبات تحت السطح لأغراض اكتشاف واستغلال النفط والغاز والمعادن . وعلى أي حال فإن البيانات الزلزالية لا تستخدم بصفة عامة في النماذج الجيولوجية ثلاثية الأبعاد لأي غرض عبر تحديد قمة وقاع التموذج . ويعرف الاختراع الحالي الحاجة لاتحاد_البيانات الزلزالية مع النماذج الجيولوجية بصورة أكثر فاعلىة للحصول على تصور أدق للمخزن تحت الأرض . gy ٠ لذلك فإن هذا هدف الاختراع الحالي هو تقديم طريقة للتنبؤ الدقيق بالقيم الكمية al sa صخر أو سائل في الخزان تحت الأرض عن طريق دمج استخدام البيانات
الزلزالية والنماذج الجيولوجية ذات الأبعاد . وهناك هدف آخر للاختراع وهو تقديم طريقة sal بالقيم الكمية لخواص صخر أو سائل في الخزان تحت الأرض والتي يكون لها استخدام خاص لاكتشاف الهيدروكريونات وتعطي القدرة للمزاول على تحديد الكم والحدود للخزان المحتوي على الهيدروكربونات بدقة . وهناك هدف آخر للاختراع وهو تقديم > طريقة للتنبؤ بالقيم الكمية لخواص صخر أو سائل في الخزان تحت الأرض والتي يكون لها نفع خاص في إدارة الخزان المحتوي على الهيدروكربونات وتعطي القدرة للمزاول على تحقيق أعلى معدل أو بمعني آخر أمثل معدل لإنتاج الهيدروكربونات من الخزان . وهذه الأهداف وغيرها يتم تحقيقها في الاختراع الموصوف فيما بعد . الوصف العام للاختراع :
٠ الاختراع الحالي هو طريق للحصول على نموذج جيولوجي ذو أبعاد لخزان سائل تحت الأرض يحتوي على بيانات كمية دقيقة نسبياً من خواص صخر أو سائل . ويعطي النموذج الجيولوجي ذو الأبعاد الناتج تحديداً دقيقاً لخزان السائل ولذلك يسهل اكتشاف أو إدارة إنتاج الخزان المحتوي على الهيدروكربونات . وهذه الطريقة تطبق بداية عن طريق تحديد الحجم الجيولوجي ٠ ويتضمن ذلك خزان السائل موضع الاهتمام في صورة نموذج
06 حجميمقسم جزئياً إلي مجموعة من نماذج أحجام جزئية . ترتبط نماذج الأحجام الجزئية لموقع معين خلال كل الحجم الجيولوجي . ويتم تحديد قيم ذات نسبة خطاً قليلة جداً لخواص صخر أو سائل عن طريق التكرار لكل نموذج حجميجزئي . وهذه القيم ذات النسبة القليلة من Und) الناتجة لخاصية معينة لصخر أو Jibs لها استخدام خاص لاكتشاف الهيدروكربونات وتعطي القدرة للمزاول على التحديد الدقيق لكمية
A وحدود الخزان المحتوي على الهيدروكربونات في الحجم الجيولوجي . وبالإضافة ٠ فإن هذه القيم ذات النسبة القليلة من الخطاً لخاصية معينة لصخر أو سائل لها استخدام في
إدارة الخزان المحتوي على الهيدروكربونات في الحجم الجيولوجي ١ وتعطي المزاول القدرة على تحقيق أعلى وأمثل معدل لإنتاج الهيدروكربونات من الخزان . وبالتوافق مع صيغة معينة للطريقة الحالية يتم تقديم حجم جيولوجي تحت الأرض ؛ حيث يكون توزيع القيم الزلزالية التجريبية للمعاوقة الصوتية تم تحديده عملياً بالتجارب للحجم © الجيولوجي . ويتم تحديد الحجم الجيولوجي بنموذج حجمي له مجموعة من نماذج الحجم الجزئي . ويتم اختيار نموذج حجمي جزئي من مجموعة نماذج الأحجام الجزئية ؛ ويعطي قيمة زلزالية تجريبيةللمعاوقة الصوتية من التوزيع السابق . يتم تعيين خاصية لصخر أو سائل مرتبطة بالحجم الجيولوجي ؛ ويتم حساب أول قيمة تنبؤية للمعاوقة الصوتية لنموذج ٠ الحجم الجزئي عن طريق نموذج استجابة باستخدام أول قيمة تنبؤية للخاصية المعينة لصخر أو Jil حيث يستجيب نموذج الاستجابة للتغيرات الحادثة في القيم التنبؤية للخاصية المعينة لصخر أو سائل . يتم مقارنة القيمة التنبؤية الأولى للمعاوقة الصوتية بالقيمة الزلزالية للمعاوقة الصوتية لتحديد فرق اول بين القيم التنبؤية والزلزالية ؛ ويتم ضبط القيمة التنبؤية الأولى للخاصية ٠ المعينة لصخر أو سائل حسب الفرق الأول للحصول على قيمة تنبؤية ثائية للخاصية المعينة لصخر أو Bla حيث تقلل القيمة التنبؤية الثائية من الفرق الأول . وبعد ذلك يتم حساب قيمة تنبؤية ثانية للمعاوقة الصوتية لنموذج حجمي جزئي بواسطة نموذج استجابة باستخدام قيمة تنبؤية ثانية للخاصية المعينة لصخر أو سائل . يتم مقارنة القيمة التنبؤية الثائية للمعاوقة الصوتية بالقيمة الزلزالية للمعاوقة الصوتية ٠ لتحديد فرق ثاني بين القيم التنبؤية والزلزالية حيث يكون الفرق الثاني أقل من الفرق الأول .
ويتم تكرار هذه الخطوات حتى يكون الفرق لا يزيد عن حد معين تم تحديده مسبقاً ؛ وتكون القيمة التنبؤية للخاصية المعينة لصخر أو سائل عند نهاية التكرار والإعادة قيمة لها أقل نسبة (a جداً للخاصية المعينة لصخر أو سائل للحجم الجزئي المختار . ويتم تكرار الطريقة لنموذج حجمي جزئي آخر حتى يتم الحصول على قيم لها أقل نسبة خطأً © للخاصية المعينة لصخر أو سائل لكل نموذج حجمي جزئي في نموذج الحجم . وسيتم فهم الاختراع بصورة أكثر وضوحاً من خلال الرسوم والوصف المصاحب . شرح مختصر للرسومات شكل (١):عبارة عن مخطط انسياب منهجي يقدم فكرة عامة عن الاختراع الحالي . شكل (Y) : عبارة عن أطروحة منهجية schematic representation للحجم الجيولوجي ٠ تحت الأرض متضمناً خزان السائل موضع الاهتمام . ض شكل (؟) :3 ؛ 38 و 30 مخطط انسياب منهجي يعطي صيغة تفصيلية للطريقة في شكل .)١( شكل (؛) : أطروحة منهجية لنموذج جيولوجي ثلاثي الأبعاد ؛ تم بناءه طبقاً للصيغة في أشكال «3A 38 و 30 . ve شكل )0( : أطروحة منهجية لخلية من النموذج الجيولوجي ثلاثي الأبعاد في شكل (4) ؛ حيث تمتليء الخلية ببيانات أولية . شكل )1( : أطروحة رسم بياني للبيانات التنبؤية لخواص الصخر أو السائل تم تحديدها طبقاً للطريقة في شكل )١( . الوصف التفصيلي ٠ يوضح الشكل )١( فكرة عامة عن الطريقة الحالية في مخطط الانسياب وتوصف هذه الطريقة كالتالي ؛ وبالرجوع إلي شكل )١( .
A
تتكون هذه الطريقة من خمسة مراحل . المرحلة الأولى هي بتاء نموذج جيولوجي ذو أبعاد يعتمد على حجم جيولوجي فيزيائي حقيقي ؛ ويحتوي على خزان السائل موضع ويتكون نموذج الحجم الجيولوجي ذو الأبعاد من نموذج حجمي ينقسم إلي ٠ الاهتمام ترتيب معين من نماذج أحجام جزئية ؛ ويقدم هذا الترتيب من نماذج الأحجام الجزئية 0 محيط عمل ذو بعد Gee لربط القيم التجريبية والتنبؤية للخواص الجيوفيزيائية والبتروفيزيائية لمواقع معينة داخل الحجم الجيولوجي . وتم بناء الحجم الجيولوجي ذو الأبعاد باستخدام توزيع_للقيم الزلزالية ١ تم تحديده عملياً للمعاوقة الصوتية للحجم الجيولوجي ؛ وبيانات أخرى تجريبية أو مشتقة لها علاقة بالحجم الجيولوجي والذي يتم الحصول ale ابتدائياً في هذه الطريقة الحالية . ٠ المرحلة dnl) من الطريقة هو الملء الأولى لبيانات القيم في النموذج الجيولوجي ذو الأبعاد ٠ مطابقة لمرحلة ملء البيانات الأولئة ؛ يتم إعطاء قيم كمية لخواص جيوفيزيائية وبتروفيزيائية معينة لكل نموذج جيولوجي ذو أبعاد . Jes الأخص القيم الزلزالية المحددة تجريبياً للمعاوقة الصوتية من التوزيع السابق يتم إعطاؤها حجم جزئي من النموذج الجيولوجي ذو الأبعاد والقيم المقدرة لخواص الصخر أو السائل في الحجم ٠ الجيولوجي . وهذه القيم المقدرة لخواص الصخر أو السائل تتضمن القيم المقدرة الابتدائية لخاصية معينة لصخر أو سائل (وتسمي القيم التنبوية الأولى) التي لها أهمية خاصة في الاستخدام لهذه الطريقة الحالية . المرحلة الثالثة من هذه الطريقة هي تأسيس نموذج استجابة بتروفيزيائي . في مرحلة تأسيس نموذج بتروفيزيائي يتم اختيار واحدة أو أكثر من معادلات استجابة ؛ أو يتم ٠ اشتقاقها من علاقات فيزيائية بين الصخر و السائل والخواص الزلزالية في الخزان ؛ ويعرف ذلك عند ذوي المهارة في هذا المجال بفيزياء الصخر والساتل . وتعطي معادلات
الاستجابة القدرة على حساب القيمة التنبؤية للمعاوقة الصوتية بواسطة قيم مقدرة لخواص الصخر والسائل تتضمن القيمة التنبؤية لخاصية معينة لصخر أو سائل إلي معادلات استجابة ٠ وبذلك يولد نموذج الاستجابة البتروفيزيائي قيمة تنبؤية للمعاوقة الصوتية تستجيب للتغيرات في القيمة التنبؤية للخاصية المعينة لصخر أو سائل . © المرحلة الرابعة في هذه الطريقة هي تقليل الخطأً في القيم التنبؤية للخاصية المعينة لصخر أو سائل . تبدأ مرحلة تقليل Daal باستخدام نموذج استجابة بتروفيزيائي وقيم مقدرة لخواص الصخر أو السائل تتضمن القيمة التنبؤية الأولى للخاصية المعيئة لصخر أو سائل لحساب قيمة تنبؤية أولى للمعاوقة الصوتية لنموذج حجمي جزئي مختار . ويتم تحديد القيمة ذات أقل نسبة خطأً للخاصية المعينة لصخر أو سائل عن طريق مقارنة ٠ القيمة التنبؤية الأولى للمعاوقة الصوتية بالقيمة الزلزالية للمعاوقة الصوتية ؛ وتكرار ضبط القيمة التنبؤية للخاصية المعينة لصخر أو سائل حسب هذه المقارنة وإعادة حساب القيمة التنبؤية للمعاوقة الصوتية باستخدام نموذج استجابة بتروفيزيائي + والقيمة المضبوطة للخاصية المعينة لصخر أو سائل ؛ حتى يتم تحقيق توافق مقبول بين القيمة التنبؤية والقيمة الزلزالية للمعاوقة الصوتية . 0° وتكون القيمة التنبؤية للخاصية المعينة لصخر أو ile والتي تحقق هذا التوافق ؛ هي القيمة التي لها أقل نسبة خطأً للخاصية المعينة لصخر أو سائل لنموذج الحجم الجزئي المختار ٠. ويتم بعد ذلك اختيار نموذج حجمي جزئي جديد ويتم إعادة مرحلة تقليل all حتى يتم تحديد قيمة ذات أقل نسبة Tas للخاصية المعينة لصخر أو سائل لكل نموذج حجمي جزئي في نموذج الحجم .
Yo
ومن الملاحظ أن شكل نموذج الاستجابة البتروفيزيائي والقيم المقدرة لخواص الصخر
والسائل غير الخاصية المعينة لصخر أو سائل . والقيم الزلزالية للمعاوقة الصوتية يتم
الحفاظ على ثباتها خلال مرحلة تقليل الخطأ .
المرحلة الخامسة والأخيرة من هذه الطريقة هي إعادة ملء القيم ذات أقل نسبة خطأ في
© النموذج الجيولوجي ذي الأبعاد .
وطبقاً لمرحلة إعادة ملء القيم ذات أقل نسبة خطأً لخاصية معينة لصخر أو سائل يتم
استبدال القيم ذات أقل نسبة خطأً بكل قيمة تنبؤية أولى مقابلة لها للخاصية المعينة لصخر
أو سائل في المرحلة الثانية . ويكون الناتج النهائي للطريقة الحالية هو نموذج جيولوجي
ذو أبعاد ؛ مملوء بالكامل بقيم ذات نسبة خطأً للخاصية المعينة لصخر أو سائل + والتي ٠ ترتبط بالحجم الجيولوجي موضع الاهتمام .
Leds يلي وصف للصيغة المفضلة لهذه الطريقة في الاختراع الحالي بالعودة إلي أشكال
من (؟)-(0). في البداية عودة J شكل (؟) يوضح حجم جيولوجي فيزيائي ؛ ويحدد
بصفة عامة برقم )٠١( يتم تطبيق الطريقة في الصيغة الحالية على الحجم الجيولوجي
)٠١( الذي يتكون من سطح الأرض (VY) ومجموعة طبقات تحت الأرض (V8) و )١١( ٠ و (YA) تحتوي على مواد صخر وسائل . وبذلك يحتوي الحجم الجيولوجي )٠١(
على خزان سائل والذي لم يوضح بصفة خاصة في الأشكال ؛ تمتد الطبقات )0( و
. في طبقات متراصة على التوالي حسب العمق (VY) تحت سطح الأرض )١8(و (V1)
ويتم تمييز الطبقات (VE) و (VA) (V1) عن بعضها عن طريق اختلاف خصائصس
الصخر والسائل . وطبقاً لذلك يتم فصل الطبقات )١4( و (V1) و (VA) بأحداث طبقاتية ٠ تحدد Cal ga وحدود الخزان (Ve) و (VY)
IR
ومن المفهوم أن شكل (Y) ما هو إلا أطروحة لمفاهيم عن الحجم الجيولوجي تحت الأرض ؛ والتطبيق العملي لهذه الطريقة ليس مقتصراً على حجم جيولوجي معين ولكنه يمكن تطبيقه بصفة عامة لأي حجم جيولوجي يمكن توليد بيانات زلزالية منه . الأشكال 3C « 3B ١ 3A عبارة عن مخطط انسياب تفصيلي للصيغة الحالية والتي 0 تصور الطريقة على هيئة خطوات ؛ وقبل البدء في الطريقة الحالية يتم الحصول على بيانات زلزالية عملياً من الحجم الجيولوجي )٠١( موضع الاهتمام بواسطة مسح زلزالي باستخدام الطرق التقليدية المعروفة لذوي المهارة في هذا المجال . يتم تعريف البيانات الزلزالية بصفة عامة لأغراض الطريقة الحالية على أنها معلومات يتم الحصول علىها باستخدام موجات زلزالية dah حجم جيولوجي من مصادر طاقة زلزالية صناعية ٠ وملاحظة أوقات وصول وقوة هذه الموجات والتي تتكسر بسرعات عالية داخل الجسم الجيولوجي أو تنعكس من الأسطح داخل الحجم الجيولوجي استجابة للموجات الزلزالية . وهذه الأسطح تكون ناتجة عن تغيرات في سرعة الصوت أو في كثافة المحتوي . يتم تطبيق عدد من أساليب معالجة البيانات للبيانات الزلزالية لتقليل الضوضاء وجعل البيانات أكثر وضوحاً . وبالأخص تستعمل هذه الطريقة البيانات الزلزالية بعد تحويل هذه ٠ البيانات إلي عوائق صوت باستخدام أسلوب تقليدي يعرف باسم (العكس الزلزالي) . وهناك عدد من برامج العكس الزلزالي متاحة تجارياً والتي تعالج البيانات الزلزالية وتحولها إلي توزيع للقيم الزلزالية للمعاوقة الصوتية بالنسبة إلي الوقت أو العمق داخل الحجم الجيولوجي . ومن أمثلة برامج العكس الزلزالي تلك المتاحة تحت علامة تجارية (TDROV) من Americas Inc., 16430 Park Ten Place, Houston, 0060 . Tex. 77084, USA Y- y £14 vy أخرى خاصة بالبثر ترتبط بخواص الصخر Clg الزلزالية يمكن lid بالإضافة إلي متاحة للمزاول قبل بدء العملية ؛ Lad أن تكون )٠١( أو السائل في الحجم الجيولوجي وهذه البيانات تكون مطلوبة من آبار الاستكشاف أو الإنتاج (ليست موجودة بالشكل) وكما هو ملاحظ سابقاً فإن بيانات البثر تكون محدودة )٠١( وتخترق الحجم الجيولوجي وعلى . ill .في نطاق هوائي وتكون موثوقاً بها عند ارتفاع قدم واحد عمودي من عمق © يمكن الاستفادة منها في تطبيق الطريقة الحالية ةحاتم Sul) أي حال فعندما تكون بيانات . كما سيذكر فيما بعد ومن المفهوم أن وجود بيانات البثر ليس متطلباً أساسياً لتطبيق هذه الطريقة الحالية ؛ ويكون )٠١( فالمطلوب فقط هو توزيع للقيم الزلزالية للمعاوقة الصوتية للحجم الجيولوجي أو تكون البيانات الزلزالية متاحة أولاً ثم يمكن الحصول على Bled قبل بدء abe 0٠ . التوزيع منها وإذا لم تكون البيانات الزلزالية متاحة قبل بدء العملية يمكن الحصول علىها من الحجم . بواسطة عمل مسح زلزالي كخطوة عرضية في الطريقة الحالية )٠١( الجيولوجي تبداً الصيغة الحالية بالخطوة 90 ؛ حيث يعطي توزيع 30: 3B ؛ 3A بالعودة إلي أشكال من الحجم dle الزلزالية للمعاوقة الصوتية بالنسبة للزمن والذي تم الحصول A ٠5 . )٠١( الجيولوجي (AT في شكل (؟) يعرف التوزيع حجم زلزالي ثلاثي الأبعاد للمعاوقة الصوتية (حجم . و 7 للحجم لخ في وحدات طول والمحور 7 في وحدات زمن X حيث تكون المحاور و (TY) في خطوات (AT) ويتم بناء نموذج جيولوجي ثلاثي الأبعاد (3-12) من حجم ويكون للنموذج الجيولوجي ثلاثي الأبعاد نموذج حجمي ١ (£4) و (TA) د )7( و (TE) ٠ . حيث تكون المحاور الثلاثة + و ا و 7 في وحدات طول
VY
بناء النموذج الجيولوجي ثلاثي الأبعاد عن طريق فهم الحدود العامة Jas (YY) الخطوة . وعمل ربط بين حدود عوائق الصوت وحدود الخزان (AD) للمعاوقة الصوتية في حجم يجب أن (AD) إن التغيرات الحادثة في عوائق الصوت الزلزالية بالنسبة للزمن في حجم العمق في الحجم الجيولوجي (J ترتبط بالتغيرات في خواص الصخر أو السائل بالنسبة م وبعد التثبت من حدوث واجهة زمنية صحيحة للمعاوقة الصوتية في أحداث طبقاتية (Al) صحيحة ؛ يتم رسم شبكة خطوط لتحديد العلاقات الزمنية السطحية خلال حجم . باستخدام شبكة من الخطوط والخطوط المتقاطعة يتم اشتقاق الأفق الزمني من شبكة العلاقات الزمنية السطحية والتي يتم (VE) في خطوة فحصها للعثور على قيم زمنية ثابتة سالبة . ويتم تجميع الآفاق الزمنية لإكمال إطار عمل ٠ . زمني طبقي الخطوة (7؟) تحول كل شبكة علاقات زمنية سطحية إلي علاقات عمق سطحية باستخدام ويمكن إجراء عمليات ضبط . AL أفضل تقدير للسرعة المتوسطة لكل أفق داخل حجم قليلة لعلاقات العمق سطحية بعيداً عن مواقع البثر إذا كان هناك احتياج لذلك ؛ وذلك بالاعتماد على حقيقة أن التقاء السوائل يتم تحديده غالباً بتغيرات مكانية كبيرة في قيمة . ١٠ ضبط السرعات المتوسطة المقدرة حسب عمليات الضبط sale) عوائق الصوت ؛ ويمكن . في العلاقات العمق سطحية تشتق آفاق عمق من العلاقات العمق سطحية ويتم فحصها للعثور على قيم (YA) الخطوة العمق لإكمال إطار عمل طبقي . وهناك تقابل أفقي واحد GUI أيزوباك سالبة ؛ يتم تجميع . إلي واحد بين إطار العمل الطبقي للعمق وإطار العمل الطبقي للزمن ٠ ie الخطوة (40) يتم الحصول على نموذج جيولوجي ثلاثي الأبعاد عن طريق تقسيم إطار العمل الطبقي للعمق إلي مجموعة من الخلايا أو الكتل ثلاثية الأبعاد . وبالعودة إلي شكل (£6) (؟) يظهر النموذج الجيولوجي ثلاثي الأبعاد برقم )£7( والخلايات المكونة له برقم وتحتوي النماذج الجيولوجية ثلاثية الأبعاد من هذا النوع على عشرات الملايين من ٠» الخلايا ؛ وتحدد الحدود الخارجية لإطار العمل الطبقي للعمق حجم النموذج (47)ء oe ؛ بحيث لا يحدث (£Y) وتكون الاخلايا منظمة داخل النموذج الجيولوجي ثلاثي الأبعاد . تداخل بينها وبالعودة إلي شكل )0( تظهر الخلية )£2( لها بعدان في وحدات الطول على طول تمثل شكل مربع أو مستطيل . وللخلية )£8( أيضاً بعد ثالث في SY 5X المحاور . وحدات طول بطول المحور 7 والتي تمثل العمق والسمك ٠ الذي هو معطي في النموذج AT pan يتم اختيار الأبعاد الخاصة بكل خلية حسب إذا كانت متاحة . على سبيل المثال يمكن اختيار Ad الجيولوجي ثلاثي الأبعاد وبيانات أبعاد الخلية كدالة خطية عن طريق ترك مسافات بين الخطوط المتقاطعة بواسطة معدل أخذ العينات من البثر ؛ ويمكن تشكيل طبقات رقيقة من الخلايا بحيث يتم تقسيم الخلايا بالتساوي على إطاري عمل أفقيين للعمق (نسبي) موازية لقمة الأفق أو موازية لقاع الأفق ٠ يكون النموذج الجيولوجي ثلاثي الأبعاد ٠ ؛ ويتم تقديم أشكال )%( و (*) هنا للتوضيح (7؛) والخلايا المرتبطة به شكلاً محتملاً لنموذج جيولوجي ذو أبعاد داخل الاختراع الحالي . والمفهوم أن النموذج الجيولوجي ذو الأبعاد المستخدم هنا ليس محدداً بأي عدد معين من الأبعاد أو الوحدات ؛٠ ولا الخلايا المرتبطة به لها شكل هندسي معين ؛ ولذلك فإنه من صميم الاختراع الحالي بناء نموذج جيولوجي ذو أبعاد من عدد (0) أبعاد وله ٠ تتشكل حسب هذه الأبعاد. LIA yo ؛ 30 يتم ملء النموذج الجيولوجي ثلاثي 3B ٠» 3/8 وبالعودة مرة أخرى إلي أشكال و )00( ؛ ويتحدد إطار العمل الطبقي للزمن (EA) في خطوات ASN الأبعاد بالقيم وإطار العمل الطبقي للعمق لتقديم مجال سرعة ثلاثي الأبعاد والذي يعطي القدرة على عمل خرائط عمق إلي زمن ؛ ويعطي صلة بين النموذج الجيولوجي ثلاثي الأبعاد والحجم لل ء؛ وتكون القيمة الزلزالية للمعاوقة الصوتية معطاة في خط وخط متقاطع وزمن في © ووضعها في الخلية المقابلة داخل النموذج Al ويتم استخراجها من الحجم AL حجم ويتم تكرار هذه العملية حتى يتم ملء كل خلية في النموذج ٠ الجيولوجي ثلاثي الأبعاد
Al الجيولوجي ثلاثي الأبعاد بالقيمة الزلزالية المناسبة للمعاوقة الصوتية من الحجم .وبسبب الطبيعة المحدودة لنطاق الذبذبات غالباً ما يتم وضع نفس القيمة للمعاوقة الصوتية . في خلايا عديدة ؛ بمعني آخر تحليل قليل نسبياً للصوت في الحجم آل ٠ الخطوة )00( يتم تقدير قيم لخواص الصخر والسائل تتضمن قيمة أولى تنبؤية للخاصية المعينة لصخر أو سائل وتوزيع القيم المقدرة لكل خلية في النموذج الجيولوجي ثلاثي الأبعاد © ويتم اختيار خواص الصخر أو السائل بالاعتماد على استخدامها في نموذج ويتم تقدير القيم لخواص (OY) الاستجابة البتروفيزيائي والذي يأتي وصفه في خطوة الصخر أو السائل باستخدام قواعد علمية وهندسية معروفة + وكل البيانات البتروفيزيائية ٠ . المتاحة للحجم الجيولوجي على سبيل المثال يمكن أن تكون بيانات البئثر متاحة عن طريق سجلات أداء البئر والتي لخواص الصخر أو السائل بطول الطبقات عالية A تعطي قيم محددة الموقع في الإنحلال للنموذج الجيولوجي ثلاثي الأبعاد المقابل لمواقع البثر؛ وتعطي بيانات البئر قيم تحديد نوع السائل مثل : الغاز والنفط J) محددة الموقع لمسامية الصخور بالإضافة ٠ . والماء وقيم تشبع السائل
ف وهناك بيانات أخرى متاحة Jie حسابات (PVT) التي تعطي تقديرات لإنضغاطية السائل وتعطي عينات قاع JA بيانات عن كثافة الحبيبات والتركيب المعدني للصخر . ويمكن عمل توزيع القيم المقدرة لخواص الصخر أو السائل في التموذج الجيولوجي ثلاثي الأبعاد 0 بالتطابق مع أي عدد من الأساليب التقليدية Jie ؛ وزن المسافة والموقع المشترك ؛ إلي آخر هذه الأساليب .
ومن الضروري الحفاظ على علاقات معينة بين القيم الأولىة المقدرة لخاصية الصخر أو السائل والمعطاة بواسطة التوزيع في خطوات الضبط التالية في الطريقة الحالية + على سبيل المثال فإن الطبيعة النسبية لمسامية الطبقات كدالة في الموقع الرأسي
. المعطي بواسطة التوزيع يجب الحفاظ علئه في خطوات الضبط التالية ٠ تؤسس نموذج الاستجابة البتروفيزيائي والذي يعتمد على علاقات معروفة (OF) الخطوة بين خواص الصخر و السائل وعوائق الصوت ؛ يتم تمييز نموذج الاستجابة البتروفيزيائي عن النموذج الجيولوجي ثلاثي الأبعاد + حيث أنه نظام معادلات استجابة تكون فيه الخاصية المعينة التتبؤية لالصخر أو السائل متغير أولي حر . ويستخدم نموذج الاستجابة البتروفيزيائي في دفع حساب القيم التنبؤية للمعاوقة الصوتية إلي الأمام باستخدام قيم ve مقدرة لخواص الصخر أو السائل ؛ ولذلك يصف نموذج الاستجابة البتروفزيائي كيف تتغير عوائق الصوت مع واحد أو أكثر من خواص الصخر والسائل والتي تتضمن الخاصية المعينة لصخر أو سائل داخل الحجم الجيولوجي . يتم اشتقاق نموذج الاستجابة البتروفيزيائي من معادلات معاملات تمدد الصخر والسائل Ye التقليدية التنبؤية ؛ وهنا مثال لهذه المعادلات في :
VY
Mavko, G., et al, editors; Rock Physics Formulas, Rock Physics
Laboratory, Stanford University, 1993 . . والمذكور هنا كمرجع : بخواص الصخر والسائل عن طريق المعادلات التالية (AD) ترتبط عوائق الصوت : حيث AISV prhoy © سرعة الاتنضغاطية ؛ = Vv, .عدار = كثافة الكتلة ؛ عن طريق معادلة كريستوفل dilly سرعة الانضغاطية بخواص الصخر hag : Christoffle equation : حيث vp =[(K* +3 G*)/ thos] 2, معامل الكتلة (مقلوب نظام الانضغاطية) = K¥ . معامل القص = G¥
Gassmann expression : ويحسب معامل الكتلة من
K*=K.A+[(1 (Ka K on) V[(PHL/K p)+((1- PHL VE aKa Koa) : حيث Ve الإطار (مقلوب انضغاطية الصخر الصلد) alee = Koa معامل المعادن (مقلوب اتضغاطية الحبيبات) = Kop مسامية الصخر = . معامل السائل (مقلوب انضغاطية السائل) = Ko : _معامل القص يقدر كدالة خطية في مسامية الصخر تحسب كثافة الكتلة كالتالي 1 يم.مطاتدحن. مطل *(1-.phi)+.phi.(S.y * ,مط +(1-S.)* (و.10.
YA
: حيث . كثافة المادة الأولئة = thou, . كثافة المياه الجوفية = rho. . صحيحة PVT كثافة الهيدروكربونات في ظروف = 1110.12 . تشبع خزان الماء < 35 © سبق ؛ وبالعودة إلي القيم المقدرة لخواص الصخر والسائل فإن تطبيق Lad كما هو ملاحظ النموذج يحتاج إلي استعمال مضني لكل البيانات البتروفيزيائية المتاحة للحجم الجيولوجي . موضع الاهتمام لتقدير التغيرات والتبعيات والحدود لكل مقياس مطلوب في Sl وتستخدم سجلات أداء . لتقدير انضغاطية السائل PUT نموذج الاستجابة البتروفيزيائي ؛ يتم عمل حسابات ٠ لتقدير تشبع السائل ؛ وتستخدم عينات " 1 " Levert functions وتستخدم دوال ليفرت لتحديد كثافة الحبيبات والتركيب المعدني للصخر ؛ ولذلك فإن قيم الخاصية ull قاع المعينة لصخر أو سائل التي يتم التنبؤ بها في الطريقة الحالية تكون مقيدة بالقيم الفيزيائية . للخواص المتبقية للصخر والسائل المستخدمة في نموذج الاستجابة البتروفيزيائي ٠ الخطوة Tag (of حلقات اختيار تجمع خلايا ؛ حيث يتم اختيار مجموعة من الخلايا المتجاورة من النموذج الجيولوجي ثلاثي الأبعاد © Gls يكون_ لها قيمة زلزالية تجريبيةثابتة للمعاوقة الصوتية ؛ وتكون مجموعة الخلايا الموصوفة هنا عبارة عن تجمعات من الخلايا متراصة رأسياً ومتلاصقة ولكن من المفهوم أيضاً أن هذه الطريقة تطبق Lad باختيار تجمع خلايا متراصة La ومتلاصقة ؛ des اختيار مجموعة الخلايا Xo تدأ Lad Jas) ass داخل حلقة اختيار تجمع الخلايا وهذا ما يكون خطوات (VE) (17) (10) s(eA) 527) y £14
٠ تحسب الخطوة )01( القيمة التنبؤية الأولى للمعاوقة الصوتية لكل خلية في المجموعة باستخدام نموذج استجابة بتروفيزيائي والقيم المقدرة لخواص الصخر والسائل والتي . تتضمن القيمة التنبؤية الأولى للخاصية المعينة لصخر أو سائل متوسطات القيم التنبؤية الأولى التي يتم حسابها للمعاوقة الصوتية لكل (0A) تزن الخطوة . خلية في المجموعة وعلى مستوي المجموعة الكلية © تحسب الخطوة (10) الفرق بين القيمة التنبؤية الأولى المتوسطة للمعاوقة الصوتية والقيمة الزلزالية للمعاوقة الصوتية لمجموعة الخلايا بواسطة طرح القيمة التنبؤية الأولى المتوسطة_من القيمة الزلزالية . والقيمة المطلقة للفرق بين القيمة التنبؤية المتوسطة والقيمة الزلزالية للمعاوقة الصوتية هي الدالة المحددة التي يتم تقليلها بواسطة لوغاريتم وعندما تزيد القيمة المطلقة للفرق عن حد معين 961 مثلاً تقوم ٠ حلقة تقليل الخطأً ٠ بضبط القيمة التنبؤية الأولى للخاصية المعينة لصخر او سائل داخل حدود (TY) الخطوة قيمة تتبؤية ثانية . وكما هو موضح في خطوة )12( يتم ضبط القيمة التتبؤية J مقبولة التنبؤية الثانية الدالة المحددة Aad) الأولى للخاصية المعينة لصخر أو سائل بحيث تقلل للفرق بين القيمة التنبؤية الأولى المتوسطة والقيمة الزلزالية للمعاوقة الصوتية لمجموعة القيمة التنبؤية الثانية محل القيمة التتبؤية الأولى للخاصية المعينة Jad الخلايا عندما ve لصخر أو سائل في نموذج الاستجابة البتروفيزيائي. ويكون اتجاه ضبط القيمة التنبؤية الأولى للخاصية المعينة لصخر أو سائل محدداً بعلامة الفرق بين القيمة التتبؤية المتوسطة . والقيمة الزلزالية للمعاوقة الصوتئية تنبؤية ثانية ded تحسب Cus (10) للخطوة Thal الخطوة (14) تعيد حلقية تقليل اللمعاوقة الصوتية باستخدام نموذج استجابة بتروفيزيائي والقيم المقدرة لخواص الصخر أو ٠ . والتي تحل محل القيمة التنبؤية الأولى Jill)
0 ويتم تكرار خطوات (54)و(70)و(77)و(14) كلما كان ذلك ضرورياً حتى تصبح القيمة المطلقة للفرق بين القيمة التنبؤية للمعاوقة الصوتية والقيمة الزلزالية للمعاوقة الصوتية لمجموعة الخلايا المختارة أقل من أو تساوي الحد المعين الذي سبق تحديده . وتكون القيمة التتبؤية للخاصية المعينة لصخر أو سائل لكل خلية في المجموعة والتي تحقق هذه © النتيجة هي قيمة لها أقل نسبة خطأ للخاصية المعينة لصخر أو سائل لتلك الخلية .
ومن الموضح أن القيم التنبؤية للمعاوقة الصوتية ؛ والقيم التنبؤية للخاصية المعينة لصخر أو سائل هي فقط التي تتغير خلال حلقة تقليل al بينما تكون القيم المقدرة المتبقية لخواص الصخر أو السائل والقيم الزلزالية للمعاوقة الصوتية ومعادلات الاستجابة في
نموذج الاستجابة البتروفيزيائي يتم الحفاظ علىيها ثابتة .
٠ في الخطوة )11( يتم تحديد النموذج الجيولوجي ثلاثي الأبعاد بإعادة ملء الخلايا المختارة في النموذج الجيولوجي ثلاثي الأبعاد بالقيم ذات أقل نسبة خطأً للخاصية المعينة لصخر أو سائل ؛ وبالأخص القيم ذات أقل نسبة خطأ التي تم تحديدها حديثاً للخاصية المعينة لصخر أو سائل تم وضعها بدلاً من كل قيمة تنبؤية أولى مقابلة لها للخاصية المعينة لصخر أو سائل والتي كانت موجودة في النموذج الجيولوجي ثلاثي الأبعاد في خطوة )00( وتعيد
٠ الخطوة (WA) حلقة اختيار تجمع الخلايا إلي الخطوة (of) ؛» حيث يتم اختيار مجموعة أخر ى من الخلايا المتلتصقة من النموذج الجيولوجي ثلاثي الأبعاد والتي يكون لها قيمة زلزالية تجريبيةثابتة للمعاوقة الصوتية . ويتم عمل حلقة تقليل الخطاً للمجموعة الجديدة التي تم اختيارها من الخلايا المتلاصقة لتحديد القيم ذات القيم ذات أقل نسبة خطاً للخاصية المعينة لصخر أو سائل تلك الخلايا ؛ ويتم تحديث القيمة ذات أقل نسبة خطأً في النموذج
٠ الجيولوجي ثلاثي الأبعاد ؛» ويتم إعادة تكرار حلقة اختيار تجمعات الخلايا كلما كان ذلك
ف ضرورياً بالطريقة الموضحة سابقاً ؛ إلي أن يتم تحديد قيمة ذات أقل نسبة خط للخاصية المعينة لصخر أو سائل لكل خلية في النموذج الجيولوجي ثلاثي الأبعاد . وهناك ميزة إضافية لهذه الطريقة وهي اختيار الخاصية المعنية لصخر أو سائل . ويعطي اختيار خاصية معينة مناسبة لصخر أو سائل المزاول القدرة على التحديد الدقيق لخزان السائل داخل الحجم الجيولوجي عندما يكون النموذج الجيولوجي ثلاثي الأبعاد مملوء بالقيم ذات أقل نسبة Una لهذه الخاصية . ويختار المزاول الخاصة المعنية لصخر أو سائل مع أو قبل الخطوات في الشكل (4). وتكون معايير الاختيار داخل تصور ذوي المهارة في هذا المجال وتكون أيضاً متوافقة مع الاحتياجات الخاصة للمزاول ؛ وتكون أيضاً متوافقة مع التطبيق المطلوب للنموذج ٠ الجيولوجي ثلاثي الأبعاد الناتج والقيم ذات أقل نسبة Uns ؛ فعلى سبيل المثال يقوم أحد الماهرين في هذا المجال بالتعرف على المسامية أو نوع السائل أو التشبع كخاصية ذات معني للصخر أو السائل في تطبيقات إدارة الخزان ؛ ولذلك فبالرغم من أن الطريقة الحالية ليست مقتصرة على أي خاصية معينة لصخر أو سائل فإن الخواص السابقة (المسامية ؛» ونوع السائل ؛ والتشبع والنفاذية) هي أمثلة للخاصية المعينة لصخر أو سائل ٠ التي تستخدم في الطريقة الحالية . ويوضح الشكل (1) قدرة الطريقة الحالية على التنبؤ الدقيق بالقيم للخاصية المعينة لصخر أو سائل ؛ ويتم اختيار المسامية كخاصية معيئة لصخر أو سائل لهذا التطبيق ؛ ويتم تحديد قيم ذات أقل نسبة Tha للمسامية على طول ارتفاع البئر بالتطابق مع الطريقة الحالية ؛ ويتم Land تحديد قيم مقدرة لخواص الصخر أو السائل ؛ مثل التشبع والنفاذية وقيم تتبؤية ٠ للمعاوقة الصوتية على طول ارتفاع Sl بالتطابق مع الطريقة الحالية يتم رسم هذه القيم بيانيا على محاور رأسية منفصلة بخطوط منقطة ؛ ويتم رسم القيم الحقيقة التي تم قياسها
YY
عملياً للمسافة والتشبع والنفاذية وعوائق الصوت على نفس المحور السابقة بخطوط . متصلة للمقارنة بدقة بقيم المسامية في الحجم الجيولوجي ؛ وتحقق تناسباً Lim يتتضح أن الطريقة الحالية قريباً بين القيم ذات أقل نسبة خطأً والقيم المقاسة عملياً للمسامية حتى عندما تكون النسبة لنموذج الاستجابة البتروفيزيائي أقل دقة ؛ بينما Jil) المقدرة لبقية الخواص للصخر أو oo تم وصف الصيغة السابقة المفضلة لهذا الاختراع يكون من المفهوم أن البدائل والتعديلات التي تم اقتراحها وغيرها يمكن أن يتم عملها وتكون في نفس مجال وهدف الاختراع Ce . الحالي ل
Claims (1)
- YY عناصر الحماية <١ ١ طريقة لتحديد قيمة لخاصية معينة في صخر أو سائل في الحجم الجيولوجي : ؛ تشتمل الطريقة على subterranean geologic volume تحت الأرض Y acoustic impedance 43 gall للمعاوقة seismic valuesdiy ad تعيين قيمة زلزالية 1 subterranean geologic لنموذج حجمي يرتبط بالحجم الجيولوجي تحت الأرض t. volume 8 1 تعيين خاصية لصخر أو لسائل مرتبطة بالحجم الجيولوجي المذكور . pat قيمة تنبؤية أولى للخاصية المعينة لصخر أو لسائل للحجم المذكور. A حساب قيمة تنبؤية أولى للمعاوقة الصوتية لنموذج الحجم المذكور من خلال نموذج 4 استجابة باستخدام القيمة التنبؤية الأولى المذكورة للخاصية المعينة لصخر أو لسائل ٠ حيث يستجيب نموذج الاستجابة المذكور للتغيرات الحادثة في القيم التنبؤية للخاصية ١ المعينة المذكورة لصخر أو لسائل. OY مقارنة القيمة التنبؤية الأولى المذكورة للمعاوقة الصوتية بالقيمة الزلزالية المذكورة ن للمعاوقة الصوتية لتحديد الفرق بين القيم التنبؤية والزلزالية المذكورة للمعاوقة \ الصوتية . Vo ضبط ded) التنبؤية الأولى المذكورة للخاصية المعينة المذكورة لصخر أو لسائل 7 حسب الفرق المذكور للحصول على قيمة تنبؤية ثانية للخاصية المعينة المذكورة OY لصخر أو لسائل حيث JI القيمة التنبؤية الثانية المذكورة من القرق المذكور .AEA . ةيناثلا استخدام القيمة التنبؤية Lad الطريقة وفقا ” لعنصر الحماية(١) تشمل -" ١ قيمة تنبؤية تالية مضبوطة للخاصية المعينة المذكورة لصخر أو لسائل SoS hydrocarbon —bearing reservoir للكشف عن الخزان المحتوي على الهيدروكربونات 0 . في الحجم الجيولوجي المذكور 1 استخدام القيمة التنبؤية الثانية Lad الطريقة وفقا ” لعنصر الحماية(١) تشمل -" ١ ل المذكورة أو قيمة تنبؤية تالية للخاصية المعينة المذكورة لصخر أو لسائل لاستغلال من الخزان المحتوي على الهيدروكربونات hydrocarbon إنتاج الهيدروكربونات >" . في الحجم الجيولوجي المذكور 0 —bearing reservoir ¢ ؟- طريقة للتحديد المتكرر لقيمة مثل أقل نسبة خطأً للخاصية المعينة لصخر أو ١ subterranean geologic لسائل عند موقع ما في الحجم الجيولوجي تحت الأرض Y : وتشمل هذه الطريقة المذكورة على volume ¥ له توزيع محدد subterranean geologic volume إيجاد حجم جيولوجي تحت الأرض . اقيم زلزالية تجريبية للمعاوقة الصوتية ١ تحديد الحجم الجيولوجي المذكور بواسطة نموذج حجمي يحتوي على مجموعة 1 . لأحجام جزئية zie لا. من مجموعة نماذج الأحجام الجزئية المذكورة (FO اختيار نموذج حجمي A تين قيمة زلزالية تجريبية للمعاوقة الصوتية من التوزيع المذكور لنموذج الحجم . الجزئي المذكور ٠١ . تحديد خاصية لصخر أو سائل مرتبطة بالحجم الجيولوجي المذكور ١Yoyy تعيين قيمة تنبؤية أولى للخاصية المذكورة لصخر أو لسائل لنموذج الحجم الجزئي yy المذكور . > حساب القيمة التنبؤية الأولى للمعاوقة الصوتية الخاصة بنموذج الحجم الجزئي - المذكور من نموذج استجابة باستخدام القيمة التنبؤية الأولى للخاصية المعينة لصخر أو لسائل cua يستجيب نموذج الاستجابة للتغيرات الحادثة في القيم التنبؤية yy للخاصية المعينة المذكورة لصخر أو لسائل. -_مقارنة القيمة التنبؤية الأولى المذكورة للمعاوقة الصوتية بالقيمة الزلزالية المذكورة للمعاوقة الصوتية لتحديد الفرق بين القيم التنبؤية والزلزالية المذكورة للمعاوقة .+ - الصوتية . yy ضبط القيمة التنبؤية الأولى المذكورة للخاصية المعينة المذكورة لصخر أو سائل xy حسب الفرق المذكور للحصول على قيمة تنبؤية ثانية للخاصية المعينة المذكورة + الصخر أو سائل Cus تقلل القيمة التنبؤية الثانية المذكورة من الفرق الأول المذكور mo الطريقة وفقا” لعنصر الحماية(؛) تشتمل أيضاً على تكرار وإعادة الخطوات ٠ المذكورة لضبط القيم التنبوية المذكورة للخاصية المعيئة المذكورة لصخر أو سائل Clans القيم التنبؤية المذكورة للمعاوقة الصوتية المذكورة باستخدام القيم التنبؤية ؛ - المضبوطة المذكورة الخاصة بالخاصية المعينة المذكورة لصخر أو سائل ومقارنة م القيم التنبوية المحسوبة المذكورة للمعاوقة الصوتية بالقيم الزلزالية المذكورة للمعاوقة oy الصوتية GA pail بين القيم التنبؤية المحسوبة المذكورة والقيم الزلزالية + المذكورة للمعاوقة الصوتية .١ 1 الطريقة 8“ لعنصر الحماية( °( : حيث يتم إنهاء التكرار والإعادة المذكورين ب عندما لا يزيد الفرق المذكور عن حد معين تم تحديده مسبقاً . ٠ 7- الطريقة TUES لعنصر الحماية(7) : حيث تكون القيمة التنبؤية المذكورة ب للخاصية المعينة المذكورة لصخر أو ساتئل لها أقل قيمة خطأ عند إنهاء التكرار v وإعادة الخطوات المذكورة . =A ١ الطريقة وفقا لعنصر الحماية( ¢ ( : حيث يشمل نموذج الاستجابة المذكور على 0 مجموعة من الخصائص المعينة لصخر أو سائل كنماذج مقاييس تحتوي على ب الخاصية المعنية المذكورة لصخر أو سائل . ١ 4- الطريقة “ad لعنصر الحماية(8) : حيث تحتوي مجموعة الخصائص المذكورة 0 لصخر أو سائل على نماذج من صخر أو سائل مرنة (قابلة للتمدد) . -٠ الطريقة وفقا” لعنصر الحماية(9) تشتمل أيضاً على تعيين قيم مقدرة لمجموعة v خواص صخر أو سائل إلي نموذج الحجم الجزئي المذكور بالإضافة إلي القيمة ب التنبؤية الأولى المذكورة للخاصية المعينة لصخر أو سائل . -١ ١ الطريقة وفقا لعنصر الحماية( Lad Jam )٠ على الحفاظ على القيم 0 المقدرة لمجموعة خواص صخر أو سائل المذكورة غير الخاصية المعينة لصخر أوكله \ سائل التي تم تثبيتها عند ضبط القيمة التنبؤية الأولى المذكورة للخاصية المعينة 7 المذكورة لصخر أو سائل . VY ١ الطريقة وفقا ” لعنصر الحماية(؛) : تشتمل loaf على ربط نموذج الحجم َ الجزئي المذكور بالموقع في الحجم الجيولوجي المذكور . -\Y \ الطريقة Lad ’ لعنصر الحماية( ¢ ( : حيث يكون نموذج الحجم الجزئي المذكور Y الذي تم اختياره هو النموذج الأول للحجم الجزئي المتضمن في مجموعات متجاورة".من نماذج الأحجام الجزئية المتجاورة لها قيمة زلزالية مشتركة للمعاوقة الصوتية . -١40 ١ الطريقة وفقا” لعنصر الحماية(3١) حيث يتم الحصول على معدل للقيم التتبؤية OY للمعاوقة الصوتية لكل نماذج الأحجام الجزئية المذكورة في المجموعات المذكورة ويتم تحديد الفرق بين معدل القيمة ay til) المذكورة للمعاوقة الصوتية والقيمة g الزلزالية المذكورة للمعاوقة الصوتية للمجموعات المذكورة . -١#© ١ الطريقة وفقا” لعنصر الحماية( 2( حيث يتم اختيار الخاصية المعينة المذكورة ١" الصخر أو سائل من مجموعة تتكون من خواص : المسامية ونوع السائل والتشبع "| والنفاذية .YA: الطريقة وفقا” لعنصر الحماية(؛) تشمل أيضاً على - ١0 تعيين قيمة زلزالية تجريبية للمعاوقة الصوتية من التوزيع المذكور إلي النموذج ay الثاني المذكور .؛ -- تعيين قيمة تنبؤية أولى للخاصية المعينة المذكورة لصخر أو سائل إلي نموذج الحجم م الجزتي الثاني المذكور .1 حساب قيمة تنبؤية أولى للمعاوقة الصوتية لنموذج الحجم الجزئي الثاني المذكور من خلال نموذج استجابة باستخدام القيمة التنبوية الأولى للخاصية المعينة المذكورة 4 لصخر أو سائل لنموذج الحجم الجزئي الثاني المذكور .3 مقارنة القيمة التنبؤية الأولى المذكورة للمعاوقة الصوتية بالقيمة الزلزالية المذكورة ٠١ للمعاوقة الصوتية لتحديد الفرق بين القيم التنبؤية والزلزالية للمعاوقة الصوتية ا لنموذج الحجم الجزئي الثاني المذكور .VY ضبط القيمة التنبؤية الأولى المذكورة للخاصية المعينة لصخر أو سائل حسب الفرق Vy المذكور للحصول على قيمة تنبؤية ثانية للخاصية المذكورة لصخر أو سائل لنموذج Vt الحجم الجزئي الثاني المذكور حيث تقلل القيمة التنبؤية الثانية المذكورة عن الفرق الأول المذكور .SY الطريقة وفقا” لعنصر الحماية(7١) : حيث يكون نموذج الحجم الجزئي الثاني X ضمن مجموعات ثانية من نماذج الأحجام الجزئية المذكورة ولها قيمة زلزالية v تجريبية للمعاوقة الصوتية .SIA الطريقة في العنصر الحماية (17) : حيث يتم الحصول على معدل للقيم 0 التتبؤية للمعاوقة الصوتية لكل نماذج J لأحجام الجزئية المذكورة هى المجموعات + الثانية المذكورة ويتم تحديد الفرق بين معدل القيمة التنبؤية المذكورة للمعاوقة Ad gall والقيمة الزلزالية المذكورة للمعاوقة الصوتية للمجموعات الثانية المذكورة -١١ ١ الطريقة وفقا” لعنصر الحماية(؛) تشتمل Gag على حساب ded تنبؤية ثانية ٠ للمعاوقة الصوتية لنموذج الحجم الجزئي المذكور من خلال نموذج استجابة yp باستخد ام قيمة تنبؤية ثانية للخاصية المعينة المذكورة لصخر أو سائل . -”١ ٠ الطريقة وفقا” لعنصر الحماية(5١) تشتمل load على مقارنة القيمة التنبؤية ٠ الثانية للمعاوقة الصوتية بالقيمة الزلزالية للمعاوقة الصوتية لتحديد فرق ثاني بين + القيم التنبؤية والزلزالية للمعاوقة الصوتية حيث يكون الفرق الثاني المذكور أقل من ؛ الفرق الأول المذكور . -7١ ١ الطريقة “Udy لعنصر الحماية(؛): تشتمل Und على استخدام القيمة التنبؤية + المذكورة أو قيمة تنبؤية تالية مضبوطة للخاصية المعينة المذكورة لصخر أو Jit م للكشف عن خزان محتوي على هيدروكربونات hydrocarbon —bearing reservoir في ,+ الحجم الجيولوجي المذكور. ٠ 77- الطريقة “Uy لعنصر الحماية(؛) تشتمل Lad على استخدام القيمة التنبؤية ٠ الثانية المذكورة أو قيمة تنبؤية تالية مضبوطة للخاصية المعينة المذكورة لصخر أوYePliny Jl إنتاج الهيدروكربونات oe hydrocarbon خزان محتوي على +" هيدروكربونات hydrocarbon bearing reservoir في الحجم الجيولوجي المذكور . YY طريقة لتحديد قيمة لخاصية صخر أو سائل في نموذج الحجم الجزئي 0 الجيولوجي تحت الأرض subterranean geologic volume تتضمن :+ إيجاد حجم جيولوجي تحت الأرض subterranean geologic volume له توزيع محدد ,- بالتجربة العملية من القيم الزلزالية للمعاوقة الصوتية .م تخصيص الحجم الجيولوجي المذكور بنموذج حجمي يحتوي على مجموعة نماذج الأحجام جزئية .ct قيمة زلزالية تجريبيةللمعاوقة الصوتية من التوزيع المذكور لكل نموذج من م نماذج الأحجام الجزئية المذكورة .+ تعبين خاصية لصخر أو سائل ترتبط بالحجم الجيولوجي المذكور .٠ تعيين ded تنبؤية أولى للخاصية المعينة المذكورة لصخر أو سائل لكل نموذج من yy نماذج الأحجام الجزئية . اختيار مجموعة أولى من نماذج الأحجام الجزئية المتجاورة من مجموعة نماذج الأحجام الجزئية المذكورة حيث يكون كل نموذج حجمي جزئي من المجموعة yg الأولى المذكورة له قيمة زلزالية تجريبية للمعاوقة الصوتية . 6 حساب قيمة تنبؤية أولى للمعاوقة الصوئية لكل نموذج حجمي جزئي في المجموعة _الأولى المذكورة من خلال نموذج استجابة باستخدام القيمة التنبؤية الأولى المذكورة للخاصية المعينة المذكورة لصخر أو سائل لنموذج الحجم الجزئي المذكورء حيثض زه yA يستجيب نموذج الاستجابة إلي التغيرات الحادثة في القيم التنبؤية للخاصية المعينة 4 المذكورة لصخر أو سائل. .‘ يتم الحصول على معدل للقيم التنبؤية الأولى للمعاوقة الصوتية لكل نماذج الأحجام 2 الجزئية المذكورة في المجموعة الأولى المذكورة للوصول إلي قيمة متوسطة تنبؤية YY أولى لعوائق للمجموعة الأولى المذكورة . ey | --_مقارنة القيمة المتوسطة التنبؤية الأولى المذكورة للمعاوقة الصوتية بالقيمة الزلزالية Ye للمعاوقة الصوتية لتحديد الفرق بين القيمة المتوسطة الأولى المذكورة والقيمة Yo الزلزالية المذكورة للمعاوقة الصوتية . © ضبط القيمة لتتبؤية الأولى للخاصية المعينة المذكورة لصخر أو سائل لكل تموذج © ب" حجمي جزئي مذكور في المجموعة الأولى المذكورة حسب الفرق الأول المذكور YA للحصول على قيمة تنبؤية ثانية للخاصية المعينة لصخر أو سائل لكل نموذج حجمي Yq جزئي في المجموعة الأولى المذكورة حيث تقلل القيمة التنبؤية الثانية المذكورة منv. الفرق الأول . © > حساب قيمة تنبؤية متوسطة ثائية للمعاوقة الصوتية للمجموعة الأولى المذكورة من vy نموذج الاستجابة البتروفيزيائي حيث Jad القيمة التنبؤية الثانية المذكورة للخاصية vr المعينة لصخر أو سائل محل القيمة التنبؤية المذكورة للخاصية المعينة المذكورة لصخر أو سائل في نموذج الاستجابة البتروفيزيائي المذكور . vo تحديد الفرق المذكور بين القيمة التنبؤية المتوسطة الثانية المذكورة والقيمة الزلزالية v1 للمعاوقة الصوتية وتكرار وإعادة الخطوات المذكورة لضبط القيمة التنبؤية Ty للخاصية المذكورة لصخر أو سائل Bale) حساب القيمة التنبؤية المتوسطة المذكورة - V £14 |م © للمعاوقة الصوتية للمجموعة الأولى المذكورة حتى يكون الفرق المذكور لا يزيد عن vy حد معين تم تحديده مسبقاً . ِ ٠ 14- الطريقة وفقا” لعنصر الحماية(؟1) تشتمل أيضاً على : Y اختيار مجموعة ثانية من نماذج أحجام جزئية متجاورة من مجموعة نماذج ال الأحجام الجزئية المذكورة حيث يكون لكل نموذج حجمي جزئي مذكور من 1 المجموعة الثانية المذكورة له قيمة زلزالية تجريبية للمعاوقة الصوتية . ded Clues تنبؤية أولى للمعاوقة الصوتية لكل انموذج حجمي جزئي في 1 المجموعة الثانية المذكورة من خلال نموذج استجابة باستخدام قيمة تنبؤية أولى v للخاصية المذكورة لصخر أو سائل لنموذج الحجم الجزئي المذكور . A الحصول على معدل للقيم التنبؤية الأولى للمعاوقة الصوتية لكل نماذج الأحجام 9 الجزئية المذكورة للوصول إلي قيمة تنبؤية أولى متوسطة للمعاوقة الصوتبة "1 للمجموعة الثانية المذكورة "1 مقارنة القيمة التنبؤية الأولى المتوسطة المذكورة للمعاوقة الصوتية بالقيمة VY الزلزالية للمعاوقة الصوتية للمجموعة الثانية المذكورة لتحديد فرق أول بين Vv القيم المتوسطة الأولى والزلزالية . Ve ضبط القيمة التنبؤية الأولى المذكورة للخاصية المعينة لصخر أو سائل لكل نموذج حجمي جزئي من المجموعة الثانية المذكورة حسب الفرق الأول 8 المذكور للحصول على قيمة تنبؤية ثانية للخاصية المعينة المذكورة لصخر أو VY سائل لكل نموذج حجمي جزئي من المجموعة الثانية المذكورة حيث تقلل القيمة م التنبؤية الثانية من الفرق الأول . ض14 حساب قيمة تنبؤية ثانية متوسطة للمعاوقة الصوتية للمجموعة الثانية المذكورة 0 من خلال نموذج استجابة بتروفيزيائي تحل القيمة التنبؤية الثانية المذكورة أ للخاصية المذكورة لصخر أو سائل محل Anal ded الأولى للخاصية YY المذكورة لصخر أو سائل في نموذج الاستجابة البتروفيزيائي المذكور . vy تتحديد فرق بين القيمة التنبؤية الثانية المتوسطة المذكورة والقيمة 3 للمعاوقة ve الصوتية للمجوعة الثانية المذكورة وتكرار وإعادة الخطوات المذكورة من ضيط yo القيمة التنبؤية للخاصية المعينة المذكورة لصخر أو سائل وإعادة حساب القيم a المتوسطة المذكورة للمعاوقة الصوتية للمجموعة الثانية المذكورة حتى يكون الفرق yy المذكور لا يزيد عن حد معين تم تحديده مسبقاً . ٠ ©*- الطريقة وفقا" لعنصر الحماية(77): تحتوي Lid على تكرار خطوات اختيار ٠ مجموعة من نماذج الأحجام الجزئية المتجاورة والحصول على معدل للقيم التنبؤية + للمعاوقة الصوتية المذكور وتحديد فرق بين القيمة التنبؤية المتوسطة المذكورة ؛ - والقيمة الزلزالية للمعاوقة الصوتية وضبط القيمة التنبؤية الأولى المذكورة للخاصية ٠ - المعينة المذكورة لصخر أو سائل لكل نموذج حجمي جزئي من المجموعة المذكورة + لتقليل الفرق المذكور وإعادة حساب القيمة المتوسطة المذكورة للمعاوقة الصوتية من نموذج_الاستجابة_البتروفيزيائي وإعادة تحديد الفرق المذكور وإعادة وتكرار + - الخطوات المذكورة حتى يكون الفرق لا يزيد عن حد معين تم تحديده مسبقاً لكل مجموعة من نماذج الحجم الجزئي في نموذج حجمي له قيمة زلزالية تجريبية - للمعاوقة الصوتية . 2 ٍْ vt الطريقة وفقا ” لعنصر الحماية(؟ ؟) تشتمل أيضاً على استخدام قيمة تنبؤية .- ١ مضبوطة للخاصية المعينة المذكورة لصخر أو سائل للكشف عن خزان يحتوي على 7 . في الحجم الجيولوجي المذكور hydrocarbon —bearing reservoir الهيدروكربونات 0 على استخدام قيمة تنيؤية Load لعنصر الحماية(؟7) تشتمل ” Gi الطريقة -77 لاستغلال إنتاج Jl مضبوطة للخاصية المعينة المذكورة لصخر أو X المحتوي على الهيدروكربونات goal من hydrocarbon تانوبركورديهلا ‘ . في الحجم الجيولوجي المذكور hydrocarbon —bearing reservoir ¢
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/584,063 US6302221B1 (en) | 2000-05-31 | 2000-05-31 | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data |
Publications (1)
Publication Number | Publication Date |
---|---|
SA01220301B1 true SA01220301B1 (ar) | 2006-10-29 |
Family
ID=24335778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SA01220301A SA01220301B1 (ar) | 2000-05-31 | 2001-08-22 | طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data |
Country Status (16)
Country | Link |
---|---|
US (1) | US6302221B1 (ar) |
CN (1) | CN1213308C (ar) |
AU (1) | AU782432B2 (ar) |
CA (1) | CA2380463C (ar) |
CO (1) | CO5221130A1 (ar) |
DZ (1) | DZ3203A1 (ar) |
EC (1) | ECSP014091A (ar) |
GB (1) | GB2368431B (ar) |
MX (1) | MXPA02001002A (ar) |
NO (1) | NO333339B1 (ar) |
OA (1) | OA12151A (ar) |
PE (1) | PE20011346A1 (ar) |
RU (1) | RU2271023C2 (ar) |
SA (1) | SA01220301B1 (ar) |
UA (1) | UA76697C2 (ar) |
WO (1) | WO2001092915A1 (ar) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2792419B1 (fr) * | 1999-04-16 | 2001-09-07 | Inst Francais Du Petrole | Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol |
WO2002033443A2 (en) * | 2000-06-14 | 2002-04-25 | Vermeer Manufacturing Company | Utility mapping and data distribution system and method |
US7248259B2 (en) * | 2001-12-12 | 2007-07-24 | Technoguide As | Three dimensional geological model construction |
US20070121423A1 (en) * | 2001-12-20 | 2007-05-31 | Daniel Rioux | Head-mounted display apparatus for profiling system |
US6718265B2 (en) * | 2002-08-15 | 2004-04-06 | Schlumberger Technology Corporation | Petrophysical property estimation using an acoustic calibration relationship |
US6807486B2 (en) * | 2002-09-27 | 2004-10-19 | Weatherford/Lamb | Method of using underbalanced well data for seismic attribute analysis |
US7072768B1 (en) | 2003-05-02 | 2006-07-04 | Young Alan G | Method for laterally extrapolating soil property data using soil samples and seismic amplitude data within a seismic coverage area |
US6970397B2 (en) | 2003-07-09 | 2005-11-29 | Gas Technology Institute | Determination of fluid properties of earth formations using stochastic inversion |
WO2005076124A1 (en) * | 2004-01-30 | 2005-08-18 | Exxonmobil Upstream Research Company | Reservoir evaluation methods |
US7844430B2 (en) * | 2004-01-30 | 2010-11-30 | Exxonmobil Upstream Research Co. | Reservoir model building methods |
AU2005250760B2 (en) * | 2004-05-27 | 2010-09-23 | Exxonmobil Upstream Research Company | Method for predicting lithology and porosity from seismic reflection data |
US20060047429A1 (en) | 2004-08-24 | 2006-03-02 | Adams Steven L | Method of estimating geological formation depths by converting interpreted seismic horizons from the time domain to the depth domain |
US7373251B2 (en) * | 2004-12-22 | 2008-05-13 | Marathon Oil Company | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data |
WO2007001759A2 (en) * | 2005-06-24 | 2007-01-04 | Exxonmobil Upstream Research Company | Method for obtaining porosity and shale volume from seismic data |
FR2895091B1 (fr) * | 2005-12-21 | 2008-02-22 | Inst Francais Du Petrole | Methode pour mettre a jour un modele geologique par des donnees sismiques |
US9646415B2 (en) * | 2006-05-16 | 2017-05-09 | Underground Imaging Technologies, Inc. | System and method for visualizing multiple-sensor subsurface imaging data |
US8089390B2 (en) * | 2006-05-16 | 2012-01-03 | Underground Imaging Technologies, Inc. | Sensor cart positioning system and method |
ITMI20070746A1 (it) * | 2007-04-13 | 2008-10-14 | Eni Spa | Metodo di stima dei volumi di fluidi movimentati in aree compartimentate del sottosuolo |
CN101414013B (zh) * | 2007-10-17 | 2011-10-05 | 中国石油天然气股份有限公司 | 一种利用地震资料确定地下流体的方法 |
CN101446645B (zh) * | 2007-11-27 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用地震流体阻抗进行流体确定的方法 |
MY164735A (en) * | 2009-02-16 | 2018-01-30 | Maersk Olie & Gas | Finite element modelling of borehole seismic responses in layered anisotropic formations and its use for elastic inversion |
CN101604030B (zh) * | 2009-07-17 | 2012-05-09 | 中国石化集团胜利石油管理局 | 一种利用转换横波地震资料进行流体识别的方法及装置 |
US8498853B2 (en) * | 2009-07-20 | 2013-07-30 | Exxonmobil Upstream Research Company | Petrophysical method for predicting plastic mechanical properties in rock formations |
US8355872B2 (en) * | 2009-11-19 | 2013-01-15 | Chevron U.S.A. Inc. | System and method for reservoir analysis background |
CN102445708B (zh) * | 2010-10-14 | 2013-10-23 | 中国石油大学(北京) | 三维等效富泥质砂岩速度预测模型 |
CN102426390B (zh) * | 2011-10-21 | 2013-07-03 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN103116186B (zh) * | 2011-11-16 | 2015-08-19 | 中国石油天然气集团公司 | 一种小尺度非均质储集体的容积确定方法 |
EP2825901A1 (en) | 2012-03-12 | 2015-01-21 | Vermeer Manufacturing Co., Inc | Offset frequency homodyne ground penetrating radar |
US9411071B2 (en) | 2012-08-31 | 2016-08-09 | Exxonmobil Upstream Research Company | Method of estimating rock mechanical properties |
US9835017B2 (en) * | 2012-09-24 | 2017-12-05 | Schlumberger Technology Corporation | Seismic monitoring system and method |
US9739133B2 (en) | 2013-03-15 | 2017-08-22 | Vermeer Corporation | Imaging underground objects using spatial sampling customization |
US10048403B2 (en) | 2013-06-20 | 2018-08-14 | Exxonmobil Upstream Research Company | Method and system for generation of upscaled mechanical stratigraphy from petrophysical measurements |
FR3032532B1 (fr) * | 2015-02-05 | 2020-02-28 | Services Petroliers Schlumberger | Derivation d'attributs sismiques a partir d'une propriete d'age geologique relatif d'un modele base sur le volume |
US20180217286A1 (en) * | 2015-07-20 | 2018-08-02 | Halliburton Energy Services, Inc. | Selecting logging data for petrophysical modelling and completion optimization |
US10162081B2 (en) | 2015-08-06 | 2018-12-25 | Baker Hughes a GE Company, LLC | Downhole fluid typing |
US11119239B2 (en) | 2017-01-13 | 2021-09-14 | Baker Hughes Holdings Llc | Measuring petrophysical properties of an earth formation by regularized direct inversion of electromagnetic signals |
US10928536B2 (en) | 2017-12-07 | 2021-02-23 | Saudi Arabian Oil Company | Mapping chemostratigraphic signatures of a reservoir with rock physics and seismic inversion |
WO2020149962A1 (en) * | 2019-01-18 | 2020-07-23 | Exxonmobil Research And Engineering Company | Estimation of reservoir flow properties from seismic data |
US12084959B2 (en) * | 2020-06-18 | 2024-09-10 | Halliburton Energy Services, Inc. | Estimating borehole shape between stationary survey locations |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5487001A (en) * | 1993-05-28 | 1996-01-23 | Neff; Dennis B. | Method for determining petrophysical properties of a subterranean layer |
FR2710418B1 (fr) * | 1993-09-21 | 1995-12-15 | Inst Francais Du Petrole | Méthode d'analyse des traces sismiques utilisant une technique de calibrage statistique pour en déduire des propriétés géologiques. |
US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
US5838634A (en) | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
-
2000
- 2000-05-31 US US09/584,063 patent/US6302221B1/en not_active Expired - Lifetime
-
2001
- 2001-05-18 OA OA1200200036A patent/OA12151A/en unknown
- 2001-05-18 RU RU2002102074/28A patent/RU2271023C2/ru not_active IP Right Cessation
- 2001-05-18 AU AU64772/01A patent/AU782432B2/en not_active Ceased
- 2001-05-18 WO PCT/US2001/016443 patent/WO2001092915A1/en active IP Right Grant
- 2001-05-18 CN CNB018018734A patent/CN1213308C/zh not_active Expired - Fee Related
- 2001-05-18 UA UA2002010624A patent/UA76697C2/uk unknown
- 2001-05-18 CA CA002380463A patent/CA2380463C/en not_active Expired - Fee Related
- 2001-05-18 GB GB0203430A patent/GB2368431B/en not_active Expired - Fee Related
- 2001-05-18 DZ DZ013203A patent/DZ3203A1/xx active
- 2001-05-18 MX MXPA02001002A patent/MXPA02001002A/es active IP Right Grant
- 2001-05-29 CO CO01042596A patent/CO5221130A1/es not_active Application Discontinuation
- 2001-05-30 EC EC2001004091A patent/ECSP014091A/es unknown
- 2001-05-31 PE PE2001000502A patent/PE20011346A1/es not_active Application Discontinuation
- 2001-08-22 SA SA01220301A patent/SA01220301B1/ar unknown
-
2002
- 2002-01-30 NO NO20020494A patent/NO333339B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
CO5221130A1 (es) | 2002-11-28 |
NO20020494D0 (no) | 2002-01-30 |
CA2380463C (en) | 2006-01-10 |
ECSP014091A (es) | 2002-08-01 |
GB0203430D0 (en) | 2002-04-03 |
RU2271023C2 (ru) | 2006-02-27 |
AU6477201A (en) | 2001-12-11 |
PE20011346A1 (es) | 2002-01-31 |
CN1213308C (zh) | 2005-08-03 |
WO2001092915A1 (en) | 2001-12-06 |
OA12151A (en) | 2006-05-05 |
CN1383494A (zh) | 2002-12-04 |
UA76697C2 (uk) | 2006-09-15 |
GB2368431A (en) | 2002-05-01 |
DZ3203A1 (fr) | 2001-12-06 |
MXPA02001002A (es) | 2002-08-12 |
US6302221B1 (en) | 2001-10-16 |
NO333339B1 (no) | 2013-05-13 |
CA2380463A1 (en) | 2001-12-06 |
GB2368431B (en) | 2004-08-18 |
NO20020494L (no) | 2002-03-25 |
AU782432B2 (en) | 2005-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
SA01220301B1 (ar) | طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data | |
EP3571532B1 (en) | Systematic evaluation of shale plays | |
AU2005327937B2 (en) | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data | |
EP2846175B1 (en) | Seismic survey analysis | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
US11346833B2 (en) | Reservoir fluid characterization system | |
WO2017082897A1 (en) | Method for computing lithofacies probability using lithology proximity models | |
EP2923040A1 (en) | Energy deposit discovery system and method | |
CN103946896A (zh) | 利用基于过程的模型和动态异质性估定感兴趣的地质体积的异质性的系统及方法 | |
AU2013334868B2 (en) | Distributing petrofacies using analytical modeling | |
Masoud et al. | Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya | |
Ezebialu et al. | Field development planning based on static, dynamic and geomechanical modelling of X field, Niger delta, Nigeria | |
US20240036225A1 (en) | Thermal conductivity mapping from rock physics guided seismic inversion | |
Abe et al. | Seismic attribute analysis and 3D model-based approach to reservoir characterization of “KO” field, Niger Delta | |
CN110945384B (zh) | 流体置换 | |
Gomes et al. | Lessons learned from static reservoir modelling on complex carbonate fields, onshore UAE | |
Dubey et al. | Systematic evaluation of shale play by introducing integration of inversion geophysics, petroleum system and reservoir simulation workflows | |
Boakye et al. | A Comparative Permeability Modelling Study of the Enyenra Turonian Reservoirs, Tano Basin-Ghana. | |
Lyon et al. | Integrated study of a faulted and fractured reservoir | |
Penna et al. | 3D Seismic Flow Units, Porosity and Permeability Modelling in Brazilian Presalt Reservoirs: An Overview | |
da Silva de Souza et al. | MICP-Based Petrophysical Classification of Complex Carbonate Reservoir Rocks | |
Saidu | Integrated evaluation of wet gas reservoir: minimizing volumetric uncertainties using dynamic analysis | |
Vølstad | Permeability Prediction from Acoustic Velocity | |
Touil | Determination of Reservoir Permeability Based on Irreducible Water Saturation and Porosity from Log Data and Flow Zone Indicator From Core Data (Hassi Messaoud Region, Hassi Terfa | |
Sengupta et al. | Net-to-gross from seismic P and S impedances: estimation and uncertainty analysis using Bayesian statistics |