الانعكاس الزلزالي

📅 آخر تحديث: ٢٥ يوليو ٢٠٢٦
✍️ نُشر في: ٢٥ يوليو ٢٠٢٦
الموسوعة الحرة

في الجيوفيزياء (وخاصةً في استكشاف وتطوير النفط والغاز )، يُعرف الانعكاس الزلزالي بأنه عملية تحويل بيانات الانعكاس الزلزالي إلى وصف كمي لخصائص الصخور في المكمن . قد يكون الانعكاس الزلزالي قبل أو بعد التراكم ، حتميًا أو عشوائيًا أو إحصائيًا جغرافيًا ؛ ويتضمن عادةً قياسات أخرى للمكمن مثل سجلات الآبار والعينات اللبية. [ 1 ]

مقدمة

يُجري الجيوفيزيائيون بشكل روتيني مسوحات زلزالية لجمع معلومات حول جيولوجيا حقول النفط والغاز . تسجل هذه المسوحات الموجات الصوتية التي تنتقل عبر طبقات الصخور والسوائل في باطن الأرض. ويمكن تقدير سعة وتردد هذه الموجات لإزالة أي تأثيرات جانبية أو تأثيرات ضبط التردد [ 2 ] الناتجة عن الموجة الصغيرة.

يمكن فحص البيانات الزلزالية وتفسيرها بشكل مستقل دون الحاجة إلى عملية عكسية، إلا أن ذلك لا يوفر رؤية شاملة لطبقات ما تحت السطح، وقد يكون مضللاً في بعض الحالات. ونظرًا لكفاءتها وجودتها، تستخدم معظم شركات النفط والغاز حاليًا تقنية العكس الزلزالي لزيادة دقة البيانات وموثوقيتها، ولتحسين تقدير خصائص الصخور، بما في ذلك المسامية وصافي الإنتاج. [ 3 ]

توجد العديد من التقنيات المختلفة المستخدمة في الانعكاس الزلزالي. [ 4 ] ويمكن تصنيفها تقريبًا إلى فئتين:

  1. قبل التكديس أو بعد التكديس
  2. الدقة الزلزالية أو دقة تسجيل الآبار

يُنتج الجمع بين هذه الفئات أربعة مناهج تقنية لحل مشكلة الانعكاس، ويعتمد اختيار التقنية المناسبة على الهدف المنشود وخصائص الصخور الجوفية. ورغم أن الترتيب المعروض يعكس التطورات التي شهدتها تقنيات الانعكاس خلال العشرين عامًا الماضية، إلا أن لكل مجموعة استخدامات فعّالة في مشاريع محددة أو كجزء من سير عمل أوسع.

تقدير الموجات الصغيرة

تتطلب جميع طرق الانعكاس الزلزالي الحديثة بيانات زلزالية وموجة صغيرة مُقدَّرة من هذه البيانات. عادةً، تُستخدم سلسلة معاملات الانعكاس من بئر ضمن حدود المسح الزلزالي لتقدير طور الموجة الصغيرة وترددها. يُعدّ التقدير الدقيق للموجة الصغيرة أمرًا بالغ الأهمية لنجاح أي عملية انعكاس زلزالي. قد يؤثر الشكل المُستنتج للموجة الزلزالية الصغيرة تأثيرًا كبيرًا على نتائج الانعكاس الزلزالي، وبالتالي على التقييمات اللاحقة لجودة المكمن.

يتم تقدير طيف سعة وطور الموجات الصغيرة إحصائيًا إما من البيانات الزلزالية وحدها أو من مزيج من البيانات الزلزالية وبيانات التحكم في الآبار باستخدام آبار ذات منحنيات صوتية وكثافية متاحة . بعد تقدير الموجة الزلزالية الصغيرة، تُستخدم لتقدير معاملات الانعكاس الزلزالي في عملية الانعكاس الزلزالي.

عندما تتوافق المرحلة المُقدَّرة (الثابتة) للموجة الإحصائية مع النتيجة النهائية، يتقارب تقدير الموجة بشكل أسرع مما لو بدأ بافتراض مرحلة صفرية . يمكن تطبيق تعديلات طفيفة وتقنية "التمديد والضغط" على البئر لتحسين محاذاة الأحداث. يتطلب التقدير الدقيق للموجة ربطًا دقيقًا لسجل المعاوقة بالبيانات الزلزالية. قد تؤدي أخطاء الربط بين البئر والبيانات الزلزالية إلى ظهور تشوهات في المرحلة أو التردد في تقدير الموجة. بمجرد تحديد الموجة، يحسب الانعكاس الزلزالي سجلًا اصطناعيًا لكل مسار زلزالي. ولضمان الجودة، تُجرى عملية التفاف بين نتيجة الانعكاس والموجة لإنتاج مسارات زلزالية اصطناعية تُقارن بالبيانات الزلزالية الأصلية. [ 4 ]

مكونات الانعكاس

تشمل عملية الانعكاس بيانات الحقل الزلزالي وبيانات الآبار، حيث تُستخدم بيانات الآبار لإضافة الترددات العالية أسفل النطاق الزلزالي ولتحديد نطاق الانعكاس. تُعالج سجلات الآبار أولاً لضمان وجود علاقة مناسبة بين سجلات المعاوقة والخصائص المطلوبة. ثم تُحوّل السجلات إلى بيانات زمنية، وتُصفّى لتقريب عرض النطاق الزلزالي، وتُعالج لتأثيرات البئر، وتُوازن وتُصنّف حسب الجودة.

تُعدّ البيانات الزلزالية محدودة النطاق الترددي، مما يقلل من دقتها وجودتها. ولتوسيع نطاق التردد المتاح، تُستخلص بيانات التردد المنخفض من بيانات التسجيلات، أو سرعات الموجات الزلزالية المُرحّلة زمنيًا أو عمقيًا قبل التراكم، و/أو من التدرج الإقليمي. [ 5 ] ويمكن استخلاص بيانات التردد العالي من بيانات التحكم في الآبار أو التحليل الجيولوجي الإحصائي.

تُجرى عمليات الانعكاس الأولية عادةً بقيود مُخففة، بدءًا من البيانات الزلزالية ثم إضافة بيانات الاتجاهات المحدودة من الآبار. يوفر هذا نظرة عامة تقريبية على المكمن بطريقة غير متحيزة. من الأهمية بمكان في هذه المرحلة تقييم دقة الربط بين نتائج الانعكاس وبيانات الآبار، وبين البيانات الزلزالية الأصلية والبيانات المُصنّعة. من المهم أيضًا التأكد من تطابق الموجة مع طور وتردد البيانات الزلزالية.

بدون استخدام الموجات الصغيرة، لا يكون الحل فريدًا. تعالج عمليات الانعكاس الحتمية هذه المشكلة بتقييد الإجابة بطريقة ما، عادةً باستخدام بيانات سجلات الآبار. أما عمليات الانعكاس العشوائية فتعالج هذه المشكلة بتوليد مجموعة من الحلول المحتملة، والتي يمكن تضييق نطاقها من خلال اختبار أفضل تطابق مع قياسات مختلفة (بما في ذلك بيانات الإنتاج).

انعكاس دقة البيانات الزلزالية بعد التراكم

من الأمثلة على تقنيات عكس البيانات الزلزالية بعد التراكم تقنية عكس النبضات المتباعدة المقيدة (CSSI). تفترض هذه التقنية عددًا محدودًا من معاملات الانعكاس ذات سعة أكبر. ينتج عن عملية العكس المعاوقة الصوتية (AI)، وهي حاصل ضرب كثافة الصخور وسرعة الموجة الأولية (p) . على عكس بيانات الانعكاس الزلزالي (التي تُعد خاصية سطحية)، فإن المعاوقة الصوتية هي خاصية صخرية. يتميز النموذج الناتج بجودة أعلى، ولا يعاني من التشويش والتداخل الناتج عن الموجة.

يحوّل برنامج CSSI البيانات الزلزالية إلى سجل معاوقة صوتية افتراضية عند كل مسار. تُستخدم المعاوقة الصوتية لإنتاج تفسيرات هيكلية وطبقية أكثر دقة وتفصيلاً مما يمكن الحصول عليه من التفسير الزلزالي (أو تفسير السمات الزلزالية ). في العديد من البيئات الجيولوجية، ترتبط المعاوقة الصوتية ارتباطًا وثيقًا بالخصائص الفيزيائية للصخور ، مثل المسامية، والتركيب الصخري ، وتشبع السوائل.

تُنتج خوارزمية CSSI الجيدة أربعة مجسمات عالية الجودة للممانعة الصوتية من بيانات الزلازل الكاملة أو ما بعد التراكم: ممانعة كاملة النطاق، وممانعة محدودة النطاق ، ونموذج الانعكاسية، ومكون التردد المنخفض. يمكن فحص كل مكون من هذه المكونات لمعرفة مساهمته في الحل وللتحقق من جودة النتائج. ولمزيد من ملاءمة رياضيات الخوارزمية لسلوك الصخور الحقيقية في باطن الأرض، تستخدم بعض خوارزميات CSSI نهجًا مختلطًا للمعيار، وتسمح بمعامل ترجيح بين تقليل تباين الحل وتقليل عدم تطابق الآثار المتبقية.

عكس دقة البيانات الزلزالية قبل التراكم

يُستخدم الانعكاس قبل التراكم غالبًا عندما يفشل الانعكاس بعد التراكم في التمييز بشكل كافٍ بين السمات الجيولوجية ذات بصمات المعاوقة الأولية المتشابهة. [ 6 ] يحل الانعكاس المتزامن المعاوقة الثانوية والكثافة، بالإضافة إلى المعاوقة الأولية. في حين أن العديد من السمات الجيولوجية قد تُظهر خصائص معاوقة أولية متشابهة، إلا أن القليل منها يشترك في سمات معاوقة أولية وثانية مُجتمعة (مما يسمح بتحسين الفصل والوضوح). غالبًا ما تُشير دراسة الجدوى باستخدام سجلات الآبار إلى ما إذا كان من الممكن تحقيق فصل النمط الصخري المطلوب باستخدام المعاوقة الأولية وحدها أو ما إذا كانت المعاوقة الثانوية مطلوبة أيضًا. وهذا ما يُحدد ما إذا كان الانعكاس قبل التراكم أو بعده ضروريًا.

الانعكاس المتزامن (SI) هو أسلوب ما قبل التراكم يستخدم عدة مجموعات فرعية زلزالية بإزاحة أو زاوية، بالإضافة إلى الموجات الصغيرة المرتبطة بها، كمدخلات. ينتج هذا الأسلوب معاوقة الموجة P، ومعاوقة الموجة S، والكثافة كمخرجات (مع العلم أن دقة مخرجات الكثافة نادراً ما تكون بنفس دقة المعاوقات). يساعد هذا في تحسين التمييز بين أنواع الصخور، والمسامية، وتأثيرات السوائل. يتم تقدير موجة صغيرة فريدة لكل مجموعة جزئية مُدخلة. تُدخل جميع النماذج، والمجموعات الجزئية، والموجات الصغيرة إلى خوارزمية انعكاس واحدة، مما يُمكّن الانعكاس من التعويض بفعالية عن تأثيرات الطور، وعرض النطاق، والضبط، وتمدد NMO، التي تعتمد على الإزاحة. [ 7 ]

تعمل خوارزمية الانعكاس عن طريق تقدير انعكاسات الموجة P المعتمدة على الزاوية لمجموعات الإدخال الجزئية. بعد ذلك، تُستخدم هذه الانعكاسات مع معادلات زويبريتز الكاملة (أو تقريبات، مثل معادلة آكي-ريتشاردز، في بعض الخوارزميات) لإيجاد انعكاسات مرنة محدودة النطاق. تُدمج هذه الانعكاسات بدورها مع نظيراتها منخفضة التردد من النموذج، وتُكامل للحصول على الخصائص المرنة . ثم تُحسّن هذه النتيجة التقريبية في عملية انعكاس نهائية لممانعة الموجة P، وممانعة الموجة S، والكثافة، مع مراعاة قيود صارمة ومرنة مختلفة. يمكن لأحد هذه القيود التحكم في العلاقة بين الكثافة وسرعة الموجة الانضغاطية؛ وهذا ضروري عندما لا يكون نطاق الزوايا واسعًا بما يكفي لتشخيص الكثافة.

يُعدّ تقدير الموجات الزلزالية جزءًا أساسيًا من عملية الانعكاس. ويتم ذلك بحساب مُرشِّح يُحسِّن شكل معاملات انعكاس سجلات الآبار المعتمدة على الزاوية في المنطقة محل الاهتمام، وذلك بمقارنتها مع مجموعة الإزاحة المقابلة عند مواقع الآبار. تُحسب معاملات الانعكاس من سجلات الموجات الصوتية الأولية (P-sonic) والثانوية (S-sonic) وسجلات الكثافة باستخدام معادلات زويبريتز . تُدخَل الموجات، بسعات تُمثِّل كل مجموعة إزاحة، مباشرةً إلى خوارزمية الانعكاس. ونظرًا لحساب موجة مختلفة لكل حجم إزاحة، يتم التعويض تلقائيًا عن تأثيرات عرض النطاق الترددي والقياس والضبط المعتمدة على الإزاحة. ويمكن استخدام موجة قريبة من مجموعة الإزاحة كنقطة بداية لتقدير موجة الزاوية البعيدة (أو الإزاحة).

لا تتوفر أي معلومات مسبقة عن المعاملات المرنة والكثافة خارج نطاق الحلول المحددة بأي قيود صارمة عند مواقع الآبار. وهذا يجعل مقارنة سجلات الآبار المُفلترة ومخرجات الانعكاس عند هذه المواقع بمثابة آلية طبيعية لمراقبة الجودة. يتم استبدال الترددات المنخفضة من الانعكاس بمعلومات من النموذج الجيولوجي نظرًا لضعف تحديدها بواسطة البيانات الزلزالية. عند تطبيقها في الوضع العالمي، يُضاف حد تحكم مكاني إلى دالة الهدف، ويتم عكس مجموعات فرعية كبيرة من المسارات في وقت واحد. تأخذ خوارزمية الانعكاس المتزامن مجموعات بيانات زلزالية متعددة مكدسة بزوايا مختلفة، وتُنتج ثلاثة مجلدات للمعاملات المرنة كمخرجات.

تُمثل المعاملات المرنة الناتجة خصائص الصخور الحقيقية التي يمكن ربطها مباشرةً بخصائص الخزان. تستخدم الخوارزميات الأكثر تطورًا معادلات نوت-زويبريتز الكاملة، مع مراعاة كاملة لتغيرات السعة والطور مع الإزاحة. ويتم ذلك من خلال استخلاص موجات صغيرة فريدة لكل مجموعة جزئية من المدخلات. يمكن تقييد المعاملات المرنة نفسها مباشرةً أثناء عملية الانعكاس الزلزالي، كما يمكن تطبيق علاقات فيزياء الصخور لتقييد أزواج المعاملات المرنة ببعضها البعض. تُعيد نماذج المعاملات المرنة النهائية إنتاج البيانات الزلزالية المدخلة على النحو الأمثل، حيث يُعد ذلك جزءًا من تحسين عملية الانعكاس الزلزالي.

الانعكاس الإحصائي الجغرافي بعد التراكم

يدمج التحليل الجيولوجي الإحصائي بيانات الآبار عالية الدقة مع بيانات الزلازل ثلاثية الأبعاد منخفضة الدقة، ويُنتج نموذجًا بتفاصيل رأسية عالية بالقرب من نقطة التحكم في البئر وبعيدًا عنها. يُولّد هذا نماذج للمكامن ذات أشكال جيولوجية معقولة، ويُقدّم تحديدًا كميًا واضحًا لعدم اليقين لتقييم المخاطر. كما تُنتج نماذج بتروفيزيائية عالية التفصيل، جاهزة للإدخال في محاكاة تدفق المكمن.

يختلف علم الإحصاء الجيولوجي عن الإحصاء العادي في أنه يُقرّ بأن بعض النتائج فقط هي الممكنة جيولوجيًا. يدمج التحليل الجيولوجي الإحصائي العكسي البيانات من مصادر متعددة، ويُنشئ نماذج ذات دقة أعلى من البيانات الزلزالية الأصلية، وتتوافق مع الأنماط الجيولوجية المعروفة، ويمكن استخدامها لتقييم المخاطر والحد منها.

تُمثَّل البيانات الزلزالية وسجلات الآبار وغيرها من البيانات المدخلة بدالة كثافة احتمالية (PDF)، والتي توفر وصفًا إحصائيًا جغرافيًا قائمًا على المدرجات التكرارية والمخططات التباينية . وتحدد هذه البيانات مجتمعةً احتمالات قيمة معينة في موقع معين، والنطاق الجيولوجي المتوقع والتركيب الجيولوجي في جميع أنحاء المنطقة التي تم نمذجتها.

بخلاف خوارزميات الانعكاس والنمذجة الجيولوجية التقليدية، يعتمد الانعكاس الجيولوجي الإحصائي على نهج أحادي الخطوة، حيث يقوم بحل المعاوقة وأنواع الخصائص المنفصلة أو الليثوفاسيس في آن واحد. ويؤدي اتباع هذا النهج إلى تسريع العملية وتحسين دقتها.

تُدمج دوال كثافة الاحتمال الفردية باستخدام تقنيات الاستدلال البايزي ، مما ينتج عنه دالة كثافة احتمال لاحقة مشروطة بمجموعة البيانات الكاملة. تحدد الخوارزمية وزن كل مصدر بيانات، مما يزيل أي تحيز محتمل. ثم تُدخل دالة كثافة الاحتمال اللاحقة إلى خوارزمية مونت كارلو لسلسلة ماركوف لإنشاء نماذج واقعية للممانعة الكهربائية والوحدات الصخرية، والتي تُستخدم بدورها لمحاكاة خصائص الصخور مثل المسامية. تُكرر هذه العمليات عادةً حتى يظهر نموذج يطابق جميع المعلومات. حتى مع أفضل نموذج، يبقى بعض عدم اليقين. يمكن تقدير عدم اليقين باستخدام بذور عشوائية لإنشاء نطاق من النتائج. يُعد هذا مفيدًا بشكل خاص عند التعامل مع المعلمات الحساسة للتغير؛ إذ يُمكّن هذا النوع من التحليل من فهم أعمق لمخاطر التطوير.

عكس تفاصيل السجل قبل التراكم

يدمج الانعكاس الجيولوجي الإحصائي لسعة الإشارة مقابل الإزاحة (AVO) عملية الانعكاس المتزامن لسعة الإشارة مقابل الإزاحة (AVA) في خوارزمية الانعكاس الجيولوجي الإحصائي، مما يتيح الحصول على دقة عالية، وإحصاءات جيولوجية، وبيانات AVO في طريقة واحدة. تتوافق نماذج الإخراج ( النتائج ) مع معلومات سجلات الآبار، وبيانات AVO الزلزالية، وتراعي علاقات خصائص الصخور الموجودة في الآبار. كما تُنتج الخوارزمية في آنٍ واحد الخصائص المرنة (معاوقة الموجة P، ومعاوقة الموجة S، والكثافة) وأحجام الليثولوجيا، بدلاً من حساب الليثولوجيا أولاً ثم ملء الخلية بقيم المعاوقة والكثافة. ولأن جميع نماذج الإخراج تتطابق مع جميع بيانات الإدخال، يمكن تقييم عدم اليقين كمياً لتحديد نطاق احتمالات الخزان ضمن البيانات المتاحة.

يستخدم برنامج AVA للتحليل الجيولوجي الإحصائي تقنيات متطورة، بما في ذلك أخذ عينات مونت كارلو لسلسلة ماركوف (MCMC) ونمذجة الليثولوجيا متعددة التوزيعات الغاوسية. وبذلك، يُمكن استغلال "التآزر المعلوماتي" لاستخلاص تفاصيل تُخفيها أو تُغفلها تقنيات التحليل الحتمية. ونتيجةً لذلك، يُحقق علماء الجيولوجيا نجاحًا أكبر في إعادة بناء كلٍ من البنية العامة والتفاصيل الدقيقة للمكمن. كما يُتيح استخدام مجسمات الزلازل متعددة الزوايا في برنامج AVA للتحليل الجيولوجي الإحصائي إجراء تقييم أدق لخصائص الصخور المرنة والليثولوجيا المحتملة أو السحنات الزلزالية وتوزيعات السوائل.

تبدأ العملية بتحليل بتروفيزيائي مفصل ومعايرة سجلات الآبار. تستبدل عملية المعايرة قياسات الصوت والكثافة غير الموثوقة والمفقودة بقيم مُصنّعة من نماذج بتروفيزيائية وفيزيائية صخرية مُعايرة. تُستخدم معلومات سجلات الآبار في عملية الانعكاس لاستخلاص الموجات الصغيرة، وتوفير مُكوّن التردد المنخفض غير الموجود في البيانات الزلزالية، وللتحقق من النتائج النهائية وتحليلها. بعد ذلك، تُستخدم بيانات الأفق وسجلات الآبار لبناء الإطار الطبقي للمعلومات الإحصائية اللازمة لبناء النماذج. وبهذه الطريقة، تُستخدم بيانات سجلات الآبار فقط لتوليد إحصائيات ضمن أنواع الصخور المتشابهة في الطبقات الطبقية للأرض.

يُجرى تحليل الموجات الصغيرة باستخراج مرشح من كل حجم من الأحجام الزلزالية باستخدام معاوقة البئر المرنة (الزاوية أو الإزاحة) كمخرج مطلوب. وتعتمد جودة نتيجة الانعكاس على الموجات الزلزالية المستخرجة. ويتطلب ذلك سجلات دقيقة للموجات الصوتية الأولية والثانوية والكثافة، مرتبطة بالأحداث المناسبة في البيانات الزلزالية. تُستخرج الموجات الصغيرة بشكل فردي لكل بئر. ثم تُستخرج موجة صغيرة نهائية "متعددة الآبار" لكل حجم باستخدام أفضل روابط الآبار الفردية، وتُستخدم كمدخل لعملية الانعكاس.

يتم إنشاء المدرجات التكرارية والتباينات لكل طبقة جيولوجية ونوع صخري، وتُجرى محاكاة أولية على مناطق صغيرة. ثم يُشغَّل برنامج AVA للتحليل الجيولوجي الإحصائي العكسي لإنشاء العدد المطلوب من النماذج، والتي تتطابق مع جميع بيانات الإدخال. تُضبط جودة النتائج من خلال مقارنة مباشرة بين أحجام خصائص الصخور المُستنتجة وسجلات الآبار. تشمل عملية ضبط الجودة الإضافية مراجعة جميع معلمات الإدخال ونتائج المحاكاة من قِبل فريق متعدد التخصصات. ينتج عن تحليل النماذج المتعددة مكعبات أو خرائط للخصائص المتوسطة (P50). في أغلب الأحيان، تكون هذه المكعبات خاصة بالنوع الصخري أو السحنات الزلزالية، بالإضافة إلى احتمالات مُتوقعة للنوع الصخري أو السحنات ، ولكن من الممكن أيضًا الحصول على مخرجات أخرى. كما يتم إنشاء مكعبات مختارة للنوع الصخري والسحنات لاحتمالات P15 وP85 (على سبيل المثال). يتم تمثيل الأجسام ثلاثية الأبعاد للخزانات الحاملة للهيدروكربونات مع خصائصها الصخرية المقابلة، ويتم تحديد عدم اليقين في حجم الخزان وخصائصه.

انظر أيضاً

مراجع

  1. ↑ تشين، يانغكانغ؛ تشين، هانمينغ؛ شيانغ، كوي؛ تشين، شياوهونغ (2017). "استيفاء سجلات الآبار الموجه بالبنية الجيولوجية لعكس الموجة الكاملة عالي الدقة" . المجلة الجيوفيزيائية الدولية . 209 (1): 21-31 . doi : 10.1093/gji/ggw343 .
  2. ↑ معجم مصطلحات حقول النفط ، تم الاطلاع عليه بتاريخ 2011-06-03.
  3. ↑ بيندريل، ج.، "الانعكاس الزلزالي - أداة حاسمة في توصيف الخزانات"، مجلة النفط والغاز الإسكندنافية، العدد 5/6، 2006، ص 19-22.
  4. 1 2 سين، إم كيه، "الانعكاس الزلزالي"، جمعية مهندسي البترول، 2006.
  5. ↑ Latimer, R., Davison, R., van Riel, P., "دليل المفسر لفهم والعمل مع بيانات المعاوقة الصوتية المستمدة من الزلازل"، The Leading Edge، مارس 2000، ص 242-256.
  6. ↑ بيندريل، ج.، "الانعكاس الزلزالي - أفضل أداة لتوصيف الخزانات"، مسجل CSEG.
  7. ↑ Pendrel, J., Dickson, T., "Simultaneous AVO Inversion to P Impedance and Vp/Vs", SEG.

للمزيد من القراءة

  • Caulfield, C., Feroci, M., Yakiwchuk, K. "الانعكاس الزلزالي لتخطيط الآبار الأفقية في غرب ساسكاتشوان"، الجيوفيزياء المتطورة من خلال الابتكار، ص  213-214.
  • Chakrabarty, C., Fossey, J., Renard, G., Gadelle, C. "عملية SAGD في حقل شرق سينلاك: من توصيف الخزان إلى التطبيق الميداني"، رقم 1998.192.
  • كونتريراس، أ.، توريس-فيردين، س.، تشيسترز، و.، كفين، ك.، غلوب، م.، "الانعكاس العشوائي المشترك للسجلات البتروفيزيائية وبيانات الزلازل ثلاثية الأبعاد قبل التراكم لتقييم الاستمرارية المكانية لوحدات السوائل بعيدًا عن الآبار: تطبيق على خزان هيدروكربوني في المياه العميقة في خليج المكسيك"، ندوة SPWLA السنوية السادسة والأربعون للتسجيل، 26-29 يونيو 2005.
  • دي باروس، ديتريش، م.، "انعكاس الموجة الكاملة لتجمعات الطلقات من حيث المعلمات المرنة المسامية"، EAGE، لندن، يونيو 2007.
  • دويتش، سي.، نمذجة الخزانات الجيولوجية الإحصائية، نيويورك: مطبعة جامعة أكسفورد، 2002، 376 صفحة.
  • فرانسيس، أ.، "محدودية ومزايا الانعكاس الزلزالي الحتمي والعشوائي"، مجلة CSEG Recorder، فبراير 2005، ص  5-11.
  • حسنوسي، د.، أدهيتياوان، إ.، باسير، أ.، ليسابالي، ل.، فان إيكينهوف، ر.، "الانعكاس الزلزالي كأداة مثيرة لتحديد توزيع السحنات في خزانات تياكا الكربوناتية، سولاويزي - إندونيسيا"، وقائع المؤتمر السنوي الحادي والثلاثين لجمعية البترول الإندونيسية، مايو 2007.
  • راسل، ب.، هامبسون، د.، "القديم والجديد في الانعكاس الزلزالي"، مسجل CSEG، ديسمبر 2006، ص  5-11.
  • ستيفن، ك.، ماكبيث، سي.، "تقليل عدم اليقين في التنبؤ بالخزانات عن طريق تحديث نموذج عشوائي باستخدام مطابقة التاريخ الزلزالي"، تقييم هندسة الخزانات SPE، ديسمبر 2008.
  • فارغاس-ميليزا، ل.، ميغتشون، ج.، فاسكيز، ج.، "تقدير الخصائص الفيزيائية للصخور من خلال دمج AVO، والانعكاس الزلزالي، وتحليل السمات المتعددة في حجم ثلاثي الأبعاد لبلاويلا، فيراكروز"، المؤتمر الدولي AAPG: 24-27 أكتوبر 2004، كانكون، المكسيك.
  • وانغ، إكس.، وو، إس.، شو، إن.، تشانغ، جي.، "تقدير تشبع هيدرات الغاز باستخدام عكس النبضات المتفرقة المقيدة: دراسة حالة من شمال بحر الصين الجنوبي"، علوم الأرض والغلاف الجوي والمحيطات، المجلد 17، العدد 4، 799-813، ديسمبر 2006.
  • واتسون، آي.، لاينز، إل.، "الانعكاس الزلزالي في قمة بايكس، ساسكاتشوان"، تقرير بحثي من CREWES، المجلد 12، 2000.
  • ويتفيلد، ج.، "العلاقة بين صافي الأجر والسعة مقابل تدرجات الإزاحة: دراسة حالة خليج المكسيك"، رسالة ماجستير من جامعة هيوستن، 1993.
  • Zou, Y., Bentley, L., Lines, L., "دمج محاكاة الخزان مع نمذجة الزلازل ذات الفاصل الزمني"، المؤتمر الوطني لجمعية الجيوفيزيائيين الجيولوجيين الكندية لعام 2004.