تشخيص الزلازل لإنشاء النفاذية من مقاييس السنتيمتر إلى ما دون الكيلومتر في الصخور البلورية خلال التحفيز القصي
Seismicity diagnostic of permeability creation from centimeter to subkilometer scales in crystalline rock during shear stimulation

المجلة: Science Advances، المجلد: 12، العدد: 2
DOI: https://doi.org/10.1126/sciadv.ady5201
PMID: https://pubmed.ncbi.nlm.nih.gov/41512060
تاريخ النشر: 2026-01-09
المؤلف: Pengliang Yu وآخرون
الموضوع الرئيسي: دراسات الزلازل والتكتونيات

نظرة عامة

تبحث الدراسة في التطور الديناميكي للاختراق في قشرة الأرض، وخاصة كيف يتأثر بالزلازل الصغيرة المتقطعة التي تولد المسامية والاختراق من خلال إنشاء أو إعادة تنشيط الشقوق. تؤسس الدراسة علاقة مقاييس بين اللحظة الزلزالية ($M_0$) وتوليد الاختراق الإضافي ($\Delta k$) باستخدام تجارب مختبرية توفر قيودًا مطلقة على اللحظة الزلزالية. تكشف النتائج عن علاقة تناسبية بين $\Delta k$ و $M_0$ على مقياس السنتيمتر، وقوة تنبؤية لوصف المقاييس بموجب المعادلة $\Delta k = \lambda M_0^{2/3}$ التي تنطبق من السنتيمتر إلى مقاييس دون الكيلومتر.

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

مقدمة

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

يشير المؤلفون إلى أن الزلازل الصغيرة (MEQs) التي ت triggered بواسطة حقن السوائل توفر رؤى حول تطور الاختراق وفعالية التحفيز. يستكشفون كيف يمكن ربط MEQs بالتغيرات في شكل الشقوق وميزانيات الطاقة، مما يشير إلى أن الطاقة المنبعثة خلال MEQs يمكن أن تعمل كمؤشر على تعزيز الاختراق. لقد أقامت الدراسات السابقة روابط بين بيانات MEQ والاختراق، بما في ذلك استخدام كثافة الزلازل وتقديرات الانتشار الهيدروليكي. على الرغم من هذه التقدمات، لا تزال العلاقة الحتمية التي تربط خصائص MEQ بتغيرات الاختراق بعيدة المنال. تهدف هذه الدراسة إلى معالجة هذه الفجوة من خلال التحقيق في العلاقة المقياسية بين تغيير الاختراق ($\Delta k$) واللحظة الزلزالية ($M_0$) من خلال تجارب مختبرية على الشقوق في الصخور البلورية، مع فحص ما إذا كانت العلاقة يمكن التعبير عنها كـ $\Delta k \propto M^\alpha_0$ مع $\alpha = 1$.

طرق

في هذه الدراسة، تم إجراء تجارب مختبرية على شق مفرد مائل مقطوع بالمنشار داخل نوى أسطوانية من جرانيت ويستلي وجرانيت يوتا FORGE، موجهة بزاوية 30° إلى المحور الطولي. تم تجعيد سطح الشق باستخدام مسحوق جلخ بحجم حبيبات #60 لضمان إعادة الإنتاج عبر التجارب. تم عزل النوى هيدروليكيًا وتعرضت لحمل هيدروستاتيكي قدره 10 ميجا باسكال، مع تطبيق ضغط قصي مسبق على حوالي 80% من قوتها. تم استخدام حالتين حدوديتين: ضغط قصي ثابت (CSS)، يحاكي التحميل التكتوني المستمر دون تخفيف الضغط القصي، وعدم تشوه (ZD)، يمثل استرخاء الضغط أثناء الانزلاق. تم تحفيز إعادة تنشيط الشق من خلال زيادة تدريجية في ضغوط المسام upstream مع الحفاظ على ضغوط downstream ثابتة.

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

نتائج

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

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

مناقشة

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

يستكشف المؤلفون أيضًا العلاقة المقياسية بين اللحظة الزلزالية ($M_0$) وتغيير الاختراق ($\Delta k$)، ويجدون ارتباطًا خطيًا عبر تجارب متعددة. يقترحون قانون مقياسي تشخيصي، $\Delta k \propto M_0^{2/3}$، الذي يأخذ في الاعتبار الاعتماد المتبادل بين معامل القص، خشونة الشق، وانخفاض الضغط، مما يشير إلى أن هذه العلاقة تنطبق عبر مقاييس مختلفة من المختبر إلى البيئات الميدانية. تؤكد الدراسة على أهمية فهم هذه الديناميات للتنبؤ بتغييرات الاختراق في الصخور القشرية ذات الاختراق المنخفض، وربطها بالزلازل وتقديم إطار للبحوث المستقبلية في الجيوميكانيكا وتحفيز الخزانات.

Journal: Science Advances, Volume: 12, Issue: 2
DOI: https://doi.org/10.1126/sciadv.ady5201
PMID: https://pubmed.ncbi.nlm.nih.gov/41512060
Publication Date: 2026-01-09
Author(s): Pengliang Yu et al.
Primary Topic: earthquake and tectonic studies

Overview

The research investigates the dynamic evolution of permeability in the Earth’s crust, particularly how it is influenced by episodic microearthquakes that generate porosity and permeability through the creation or reactivation of fractures. The study establishes a scaling relationship between seismic moment ($M_0$) and incremental permeability generation ($\Delta k$) using laboratory experiments that provide absolute constraints on seismic moment. The findings reveal a proportional relationship between $\Delta k$ and $M_0$ at the centimeter scale, and a predictive power-law scaling described by the equation $\Delta k = \lambda M_0^{2/3}$ that holds from centimeter to subkilometer scales.

The prefactor $\lambda$ is influenced by factors such as stress drop, fault roughness, and deformation modulus, exhibiting a wide range of variability across ten decades for natural parameter ranges. However, the observed permeabilities are more tightly constrained, with $\lambda$ varying only over two orders of magnitude across the specified length scales. This suggests interdependencies among the controlling variables, highlighting the relationship’s potential as a diagnostic tool for predicting fluid flow in crustal reservoirs.

Introduction

The introduction of the research paper discusses the relationship between permeability in the shallow crust and various geological processes. It highlights that permeability decreases with depth due to increasing confining pressure and the healing of tectonic strain-induced damage. Although crustal temperatures rise with depth, the recovery of geothermal heat through forced convection is limited by low permeabilities. The study emphasizes the role of shear stimulation, which enhances permeability through the reactivation of faults and fractures, particularly during tectonic earthquakes. This process is crucial for the sustainable recovery of deep geothermal resources.

The authors note that microearthquakes (MEQs) triggered by fluid injection provide insights into permeability evolution and the effectiveness of stimulation. They explore how MEQs can be linked to changes in fracture morphology and energy budgets, suggesting that the energy released during MEQs can serve as an indicator of permeability enhancement. Previous studies have established connections between MEQ data and permeability, including the use of seismicity density and hydraulic diffusivity estimates. Despite these advancements, a deterministic relationship linking MEQ characteristics to permeability changes remains elusive. This study aims to address this gap by investigating the scaling relationship between permeability change ($\Delta k$) and seismic moment ($M_0$) through laboratory experiments on faults in crystalline rocks, specifically examining whether the relationship can be expressed as $\Delta k \propto M^\alpha_0$ with $\alpha = 1$.

Methods

In this study, laboratory experiments were conducted on a saw-cut single inclined fault within cylindrical cores of Westerly granite and Utah FORGE granitoid, oriented at 30° to the longitudinal axis. The fault surface was roughened using #60 grit abrasive powder to ensure reproducibility across experiments. The cores were hydraulically isolated and subjected to a hydrostatic load of 10 MPa, with shear prestress applied to approximately 80% of their strength. Two boundary conditions were employed: Constant Shear Stress (CSS), simulating continuous tectonic loading without shear stress relief, and Zero Deformation (ZD), representing stress relaxation during slip. Fault reactivation was triggered by incrementally increasing upstream pore pressures while maintaining constant downstream pressures.

The flow rate along the fault was analyzed using the cubic law to derive the overall equivalent bulk permeability, \( k \), which accounts for the influence of fracture spacing on fluid flow. The study acknowledges limitations in its assumptions of idealized flow conditions and homogeneous aperture distributions, which may not accurately reflect the complexities of real fault mechanics. Acoustic Emission (AE) events were recorded using calibrated PZT sensors, providing insights into fault slip dynamics and the relationship between fluid pressure increments and shear displacement. The results indicate that fault reactivation occurs when the balance between frictional resistance and driving shear forces is achieved, allowing for a steady flow regime that facilitates permeability assessment. Detailed calculations of effective normal and shear stress, along with seismic moment evaluations, are provided in supplementary texts.

Results

The “Results” section of the research paper presents key findings derived from the conducted experiments or analyses. It details the outcomes of various tests, highlighting significant trends and patterns observed in the data. The results are often accompanied by statistical analyses, including p-values and confidence intervals, to validate the findings.

Additionally, the section may include graphical representations, such as charts or tables, to illustrate the relationships between variables or the effectiveness of interventions. These visual aids serve to enhance the clarity of the results and support the conclusions drawn by the researchers. Overall, the findings contribute to the broader understanding of the research question and may suggest implications for future studies or practical applications.

Discussion

In this section, the authors discuss laboratory experiments designed to investigate the relationship between fluid injection-induced earthquakes and changes in permeability within fault systems. Using calibrated Piezoelectric transducers (PZTs), they establish a mechanistic link between seismic moment and permeability enhancement. The experiments demonstrate that shear slip is the primary factor influencing permeability changes, significantly more so than changes in effective normal stress. Specifically, permeability changes during shear slip are approximately an order of magnitude greater than those caused by normal stress variations alone. The findings indicate that shear-induced dilation enhances permeability by creating or enlarging fluid pathways, while initial slip may lead to temporary reductions in permeability due to wear product generation.

The authors further explore the scaling relationship between seismic moment ($M_0$) and permeability change ($\Delta k$), finding a linear correlation across multiple experiments. They propose a diagnostic scaling law, $\Delta k \propto M_0^{2/3}$, which accounts for the interdependencies of shear modulus, fault roughness, and stress drop, suggesting that this relationship holds across various scales from laboratory to field settings. The study emphasizes the importance of understanding these dynamics for predicting permeability changes in low-permeability crustal rocks, linking them to seismicity and offering a framework for future research in geomechanics and reservoir stimulation.