DOI: https://doi.org/10.5194/se-17-617-2026
تاريخ النشر: 2026-04-10
المؤلف: Ryan Schultz وآخرون
الموضوع الرئيسي: دراسات الزلازل والتكتونيات
نظرة عامة
تبحث الدراسة في القيود المفروضة على الحد الأقصى الممكن للحجم (MMAX) من تسلسلات الزلازل المستحثة باستخدام اختبارات CAP، وهي مجموعة من الأساليب الإحصائية المصممة لتحليل كتالوجات الزلازل. تركز الدراسة على بيانات من مختبرات تحت الأرض، حيث تسمح الملاحظات الجيوفيزيائية عالية الدقة بتمييز واضح بين تسلسلات الزلازل المقيدة وغير المقيدة. تشير النتائج إلى أن التسلسلات المقيدة، التي تتميز بنمو محدود في الحجم أثناء التحفيز، ترتبط عادة بالصخور السليمة، بينما ترتبط التسلسلات غير المقيدة بالتفاعلات مع صدوع أو كسور أكبر موجودة مسبقًا.
بالإضافة إلى ذلك، تستخدم الدراسة عائلة عامة من دوال MMAX لتصنيف التسلسلات المقيدة وفقًا للأطر النظرية، مثل النماذج التكتونية أو نماذج الشقوق الشد أو الشقوق القص. تساعد هذه التصنيفات في التحقق من صحة وتفسير التسلسلات المقيدة في البيئات الخاضعة للرقابة، وهو أمر حاسم لتعميم النتائج على الأحداث الزلزالية على نطاق أوسع. تختتم الدراسة بأن اختبارات CAP توفر منهجية واعدة لتقييد MMAX، مما يعزز فهمنا لكل من تسلسلات الزلازل المستحثة والطبيعية.
مقدمة
تسلط مقدمة ورقة البحث الضوء على العلاقة المهمة بين الأنشطة البشرية والزلازل المستحثة، مشددة على أن عمليات مختلفة مثل التخلص من مياه الصرف الصحي والتكسير الهيدروليكي يمكن أن تؤدي إلى حدوث زلازل. بينما يعتبر اعتماد بروتوكولات إشارات المرور لإدارة المخاطر أمرًا شائعًا، لا يزال هناك فهم محدود لنمو أحجام الزلازل المستحثة. حاولت الدراسات الأخيرة التنبؤ بالأحجام الممكنة القصوى (MMAX) لكنها غالبًا ما تعتمد على آراء خبراء ذات طابع شخصي بسبب التحديات في تقييد MMAX تجريبيًا من البيانات الزلزالية.
لمعالجة هذه الفجوة، يقدم المؤلفون اختبارات CAP، وهي مجموعة من الاختبارات الإحصائية المصممة لاستنتاج وقياس MMAX من خلال مقارنة توزيع أحجام الزلازل بتوزيع القفزات في أكبر الأحداث الحجمية. تشير التطبيقات الأولية لاختبارات CAP على حالات الحجم الأكبر إلى سلوك غير مقيد، مما يدل على أن نماذج MMAX النظرية قد لا تنطبق في هذه السياقات. على العكس من ذلك، تظهر الحالات الأصغر دليلًا على نمو الحجم المقيد. تستكشف الدراسة أيضًا بيانات من مختبرات تحت الأرض لفهم أفضل للآليات وراء الزلازل المقيدة وغير المقيدة، مما يربط الاستنتاجات الإحصائية لـ MMAX بالعمليات الفيزيائية ويوفر إطارًا أكثر قوة لتفسير الزلازل المستحثة.
طرق
توضح قسم “الطرق” تصميم التجربة والتقنيات التحليلية المستخدمة في الدراسة. استخدم الباحثون نهجًا كميًا، حيث قاموا بإجراء تحليلات إحصائية لتقييم البيانات التي تم جمعها من تجارب مختلفة. تضمنت المنهجيات المحددة تجارب خاضعة للرقابة، حيث تم التلاعب بالمتغيرات بشكل منهجي لملاحظة تأثيراتها على النتائج ذات الصلة.
شمل جمع البيانات استخدام أدوات وبروتوكولات موحدة لضمان الموثوقية والصلاحية. تم إجراء التحليل باستخدام أدوات برمجية تسهل النمذجة الإحصائية المعقدة، مما يسمح بتقييم العلاقات بين المتغيرات. تم اشتقاق النتائج الرئيسية من هذه التحليلات، مما يبرز الارتباطات المهمة والروابط السببية المحتملة، والتي تم التحقق منها بشكل أكبر من خلال دراسات التكرار. بشكل عام، كانت الطرق المستخدمة قوية، مما يضمن أن النتائج موثوقة وقابلة للتطبيق على السياق الأوسع لسؤال البحث.
نتائج
في هذا القسم، يقدم المؤلفون نتائج تحليلهم لبيانات التكسير الهيدروليكي التي تم جمعها من ثلاثة مختبرات تحت الأرض (UGLs): مختبر الصخور الصلبة Äspö (Äspö HRL)، ومنشأة سانفورد للبحث تحت الأرض (SURF)، وموقع اختبار غريمسل (GTS). يستخدمون اختبارات CAP لتقييم توزيع الحجم-التردد (MFD) للأحداث الزلزالية عبر مراحل الحقن المختلفة في كل موقع. تشير النتائج إلى نقص ملحوظ في الأحداث ذات الحجم الكبير في Äspö HRL، مع اختلافات δMLRG تتراوح من -0.2 إلى -0.4، مما يشير إلى حد أعلى على نمو الحجم الأقصى (MMAX). تحقق اختبارات Kolmogorov-Smirnov (KS) مستويات ثقة عالية (>99.99%) لمعظم المراحل، مما يدعم فرضية وجود عملية MMAX مقيدة، خاصة للمرحلة HF2، حيث يتم تفضيل نموذج شبيه بـ Galis بشكل كبير على نموذج غير مقيد.
في SURF، يكشف التحليل أن المرحلة N164 فقط تظهر نقصًا كبيرًا في الأحداث ذات الحجم الكبير، مع اختلاف δMLRG قدره -1.15 م، مما يشير إلى حدث أقل من النسبة المئوية الأولى. نتائج اختبار KS مختلطة، مع متوسط ثقة يبلغ 99.96% لـ N164 ولكن أقل لـ N142 وN128. تدعم اختبارات تقدير الاحتمالية القصوى (MLE) وجود حد لـ MMAX لـ N164، بينما لا توجد أدلة مماثلة للمراحل الأخرى. أخيرًا، في GTS، توفر النتائج من اختبارات EW أدلة حاسمة لعملية MMAX مقيدة في المرحلة HS4 وأدلة قوية لـ HF2، بينما تشير المراحل الأخرى إلى عملية غير مقيدة. بشكل عام، تساهم هذه النتائج في فهم أعمق للعوامل التي تؤثر على النشاط الزلزالي أثناء التكسير الهيدروليكي.
مناقشة
في هذا القسم، يقدم المؤلفون إطارًا منهجيًا شاملاً لاستنتاج وتقدير الحد الأقصى للحجم (MMAX) من كتالوجات الزلازل بشكل موثوق. يتكون النهج من سلسلة من الاختبارات الإحصائية، بدءًا من الاختبارات البسيطة التي توفر تقييمات نوعية أولية لحدود MMAX بناءً على توزيع تردد الحجم غوتنبرغ-ريختر (GR-MFD). يستخدم المؤلفون مقدرات الاحتمالية القصوى لتناسب GR-MFD مع الأخذ في الاعتبار عدم اكتمال الكشف من خلال معيار حجم الاكتمال (Mc). يبرزون أهمية تقييم نقص الأحداث الكبيرة باستخدام درجة القطع التجريبية (δMLRG)، والتي تعمل كبديل لتحديد وجود MMAX.
تستخدم اختبارات CAP اللاحقة أطرًا إحصائية متنوعة لتقييم تأثير MMAX على الكتالوجات بشكل صارم. يقيم اختبار Kolmogorov-Smirnov (اختبار KS) توزيع الأحجام والقفزات في الأحداث الكبيرة، بينما يقوم اختبار مقدر الاحتمالية القصوى (اختبار MLE) بتحديد قيم MMAX بناءً على دالة اللوغاريتم الاحتمالي. يقيم اختبار الوزن الجماعي (اختبار EW) نماذج MMAX المقترحة المتعددة، مما يوفر رؤى حول قوتها التفسيرية بناءً على معايير أكايك ومعلومات بايزي. يحذر المؤلفون من أن التمييز بين الكتالوجات غير المقيدة حقًا وتلك التي تفتقر إلى بيانات كافية يمكن أن يكون تحديًا، مما يبرز الحاجة إلى التناسق عبر نتائج الاختبارات وحجم الكتالوج لإجراء تأكيدات واثقة حول MMAX. بشكل عام، تشير النتائج إلى أن المنهجية تحدد وتحدد MMAX بشكل فعال، مما يساهم في فهم أعمق للزلازل في سياقات جيولوجية متنوعة.
DOI: https://doi.org/10.5194/se-17-617-2026
Publication Date: 2026-04-10
Author(s): Ryan Schultz et al.
Primary Topic: earthquake and tectonic studies
Overview
The research investigates the constraints on the maximum possible magnitude (MMAX) of induced earthquake sequences using CAP-tests, a set of statistical methods designed to analyze earthquake catalogues. The study focuses on data from underground laboratories, where high-resolution geophysical observations allow for a clear distinction between bound and unbound earthquake sequences. Findings indicate that bound sequences, characterized by restricted magnitude growth during stimulation, are typically associated with intact rock, while unbound sequences are linked to interactions with larger pre-existing faults or fractures.
Additionally, the study employs a generalized family of MMAX functions to categorize bound sequences according to theoretical frameworks, such as tectonic, tensile-crack, or shear-crack models. This categorization aids in validating and interpreting bound sequences in controlled settings, which is crucial for extrapolating findings to larger-scale seismic events. The research concludes that CAP-tests provide a promising methodology for constraining MMAX, enhancing our understanding of both induced and natural earthquake sequences.
Introduction
The introduction of the research paper highlights the significant relationship between anthropogenic activities and induced seismicity, emphasizing that various operations such as wastewater disposal and hydraulic fracturing can trigger earthquakes. While the adoption of traffic light protocols for risk management is common, there remains a limited understanding of the growth of induced earthquake magnitudes. Recent studies have attempted to forecast maximum possible magnitudes (MMAX) but often rely on subjective expert opinions due to the challenges in empirically constraining MMAX from seismic data.
To address this gap, the authors introduce CAP-tests, a suite of statistical tests designed to infer and quantify MMAX by comparing the distribution of earthquake magnitudes to the distribution of jumps in the largest magnitude events. Initial applications of CAP-tests to larger magnitude cases suggest unbounded behavior, indicating that theoretical MMAX models may not apply in these contexts. Conversely, smaller-scale cases show evidence of bound magnitude growth. The study further explores data from underground laboratories to better understand the mechanisms behind bound and unbound seismicity, linking statistical inferences of MMAX to physical processes and providing a more robust framework for interpreting induced seismicity.
Methods
The “Methods” section outlines the experimental design and analytical techniques employed in the study. The researchers utilized a quantitative approach, employing statistical analyses to evaluate the data collected from various experiments. Specific methodologies included controlled trials, where variables were systematically manipulated to observe their effects on the outcomes of interest.
Data collection involved the use of standardized instruments and protocols to ensure reliability and validity. The analysis was conducted using software tools that facilitated complex statistical modeling, allowing for the assessment of relationships between variables. Key findings were derived from these analyses, highlighting significant correlations and potential causal links, which were further validated through replication studies. Overall, the methods employed were robust, ensuring that the results are both reliable and applicable to the broader context of the research question.
Results
In this section, the authors present the results of their analysis of hydraulic fracturing data collected from three underground laboratories (UGLs): the Äspö Hard Rock Laboratory (Äspö HRL), Sanford Underground Research Facility (SURF), and Grimsel Test Site (GTS). They employ CAP-tests to assess the magnitude-frequency distribution (MFD) of seismic events across various injection stages at each site. The findings indicate a notable deficiency in large-magnitude events at Äspö HRL, with δMLRG discrepancies ranging from -0.2 to -0.4, suggesting an upper bound on maximum magnitude (MMAX) growth. The Kolmogorov-Smirnov (KS) tests yield high confidence levels (>99.99%) for most stages, supporting the hypothesis of an MMAX bound process, particularly for stage HF2, where the Galis-like model is significantly favored over an unbound model.
At SURF, the analysis reveals that only stage N164 shows a significant deficiency in large-magnitude events, with a δMLRG discrepancy of -1.15 M, indicating a <1st percentile event. The KS-test results are mixed, with average confidences of 99.96% for N164 but lower for N142 and N128. Maximum likelihood estimation (MLE) tests further support the presence of an MMAX bound for N164, while no such evidence is found for the other stages. Finally, at GTS, the results from the EW-tests provide decisive evidence for an MMAX bound process in stage HS4 and strong evidence for HF2, while other stages suggest an unbound process. Overall, these results contribute to a deeper understanding of the factors influencing seismic activity during hydraulic fracturing.
Discussion
In this section, the authors present a comprehensive methodological framework for robustly inferring and estimating the maximum magnitude (MMAX) from earthquake catalogues. The approach consists of a series of statistical tests, beginning with simple-tests that provide preliminary qualitative assessments of MMAX bounds based on the Gutenberg-Richter magnitude frequency distribution (GR-MFD). The authors utilize Maximum Likelihood Estimators to fit the GR-MFD while accounting for detection incompleteness through a magnitude-of-completeness (Mc) criterion. They highlight the importance of assessing the deficiency of large events using the empirical degree-of-truncation (δMLRG), which serves as a proxy for determining the presence of MMAX.
The subsequent CAP-tests employ diverse statistical frameworks to rigorously evaluate the influence of MMAX on the catalogues. The Kolmogorov-Smirnov test (KS-test) assesses the distribution of magnitudes and jumps in large events, while the Maximum Likelihood Estimator (MLE-test) quantifies MMAX values based on the log-likelihood function. The Ensemble Weighting test (EW-test) evaluates multiple proposed MMAX models, providing insights into their explanatory power based on Akaike and Bayesian Information Criteria. The authors caution that distinguishing between truly unbound catalogues and those lacking sufficient data can be challenging, emphasizing the need for consistency across test results and the size of the catalogue to make confident assertions about MMAX. Overall, the findings suggest that the methodology effectively identifies and quantifies MMAX, contributing to a deeper understanding of seismicity in various geological contexts.