SA01220301B1 - طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data - Google Patents

طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data Download PDF

Info

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
Application number
SA01220301A
Other languages
English (en)
Inventor
جفري جي .هامان
دونالد اتش. كالدويل
ستيفن دي. ويلسون
Original Assignee
ماراثون أويل كومباني
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ماراثون أويل كومباني filed Critical ماراثون أويل كومباني
Publication of SA01220301B1 publication Critical patent/SA01220301B1/ar

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis 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)

  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 ‏تين قيمة زلزالية تجريبية للمعاوقة الصوتية من التوزيع المذكور لنموذج الحجم‎ . ‏الجزئي المذكور‎ ٠١ . ‏تحديد خاصية لصخر أو سائل مرتبطة بالحجم الجيولوجي المذكور‎ ١
    Yo
    ‎yy‏ تعيين قيمة تنبؤية أولى للخاصية المذكورة لصخر أو لسائل لنموذج الحجم الجزئي ‎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‏ على استخدام القيمة التنبؤية ‎٠‏ الثانية المذكورة أو قيمة تنبؤية تالية مضبوطة للخاصية المعينة المذكورة لصخر أو
    Ye
    ‎Pliny 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 ¢
SA01220301A 2000-05-31 2001-08-22 طريقة للتنبؤ بالقيم الكمية لخصائص صخر أو سائل في خزان باستخدام بيانات خاصة بالزلازل seismic data SA01220301B1 (ar)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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