تحكم خشونة العيوب في الزلازل الناتجة عن الحقن
Fault roughness controls injection-induced seismicity

المجلة: Proceedings of the National Academy of Sciences، المجلد: 121، العدد: 3
DOI: https://doi.org/10.1073/pnas.2310039121
PMID: https://pubmed.ncbi.nlm.nih.gov/38215182
تاريخ النشر: 2024-01-12
المؤلف: Lei Wang وآخرون
الموضوع الرئيسي: دراسات الزلازل والتكتونيات

نظرة عامة

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

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

طرق

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

أشارت النتائج إلى أن خشونة العيوب تؤثر بشكل كبير على ديناميات الانزلاق. بالنسبة للعيب السلس، بدأ الانزلاق بعد تأخير قدره حوالي 73 ثانية، مما أدى إلى سرعات انزلاق قصوى تبلغ حوالي 3.8 ميكرومتر/ثانية، بينما أظهرت العيوب الخشنة سرعات انزلاق قصوى أقل تبلغ حوالي 1.0 ميكرومتر/ثانية وتأخيرات بدء أقصر تتراوح بين 40 إلى 50 ثانية. قام المؤلفون بتحديد معدلات ضعف الانزلاق، ووجدوا أن العيب السلس كان له معدل قدره $K_{\text{smooth}} = 61 \, \text{MPa/mm}$، مقارنة بـ $K_{\text{rough}} = 48 \, \text{MPa/mm}$ للعيوب الخشنة. تشير هذه النتائج إلى أن الخشونة الهندسية تعزز استقرار الانزلاق الناتج عن الحقن، مع آثار لفهم سلوك العيوب تحت ظروف حقن السوائل.

نتائج

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

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

نقاش

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

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

Journal: Proceedings of the National Academy of Sciences, Volume: 121, Issue: 3
DOI: https://doi.org/10.1073/pnas.2310039121
PMID: https://pubmed.ncbi.nlm.nih.gov/38215182
Publication Date: 2024-01-12
Author(s): Lei Wang et al.
Primary Topic: earthquake and tectonic studies

Overview

This study investigates the influence of fault structure and heterogeneity on fluid-induced seismicity, particularly in the context of industrial activities such as wastewater disposal and geothermal energy production. Through controlled fluid injection experiments on laboratory faults and corresponding numerical modeling, the research reveals that the geometrical roughness of faults significantly impacts the dynamics of injection-induced fault slip. Specifically, increased roughness is shown to slow down fault slip and alter the distribution of seismic events, including their frequency-magnitude characteristics and source mechanisms.

The findings underscore the critical roles of fault roughness and stress heterogeneity in facilitating the transition from aseismic to seismic slip under elevated fluid pressure conditions. This transition is associated with the occurrence of large induced seismic events, which tend to be localized around areas of high stress, referred to as asperities. Overall, the study contributes valuable insights into the mechanisms governing induced seismicity, emphasizing the need for careful consideration of fault characteristics in the management of industrial fluid injection practices.

Methods

In this section, the authors describe their experimental methodology for investigating injection-induced slip behavior on faults with varying roughness. They conducted fluid injection experiments on three cylindrical samples of Bentheim sandstone, which included one smooth fault and two rough faults. The samples were subjected to a constant confining pressure of 35 MPa and an initial fluid pressure of 5 MPa. The experiments involved a stepwise increase in fluid pressure from 5 MPa to 29 MPa, with a pressurization rate of 2 MPa/min, while monitoring shear stress and slip behavior.

The results indicated that fault roughness significantly influences slip dynamics. For the smooth fault, slip initiated after a delay of approximately 73 seconds, leading to peak slip velocities of about 3.8 μm/s, while the rough faults exhibited lower peak slip velocities of around 1.0 μm/s and shorter initiation delays of 40 to 50 seconds. The authors quantified the slip-weakening rates, finding that the smooth fault had a rate of $K_{\text{smooth}} = 61 \, \text{MPa/mm}$, compared to $K_{\text{rough}} = 48 \, \text{MPa/mm}$ for the rough faults. These findings suggest that geometrical roughness enhances the stability of injection-induced slip, with implications for understanding fault behavior under fluid injection conditions.

Results

In this study, cylindrical samples featuring an inclined fault interface with varying surface roughness were subjected to triaxial compression to investigate the impact of geometrical roughness on fault slip. The experiments involved a smooth saw-cut fault and two rough faults, with the latter exhibiting significant surface irregularities characterized by a maximum height difference of up to 7 mm. The results demonstrated that the smooth fault exhibited a uniform and flat surface topography, while the rough faults displayed pronounced surface roughness.

To analyze the roughness quantitatively, fast Fourier transform was employed to compute the power spectral density (PSD) of the elevation profiles both parallel and perpendicular to the fault dip. The findings indicated that the rough faults exhibited broad-band self-affine roughness with a Hurst exponent of approximately 0.65, suggesting a complex fractal nature. In contrast, the PSD of the smooth saw-cut fault was nearly flat across the same wavelength range, indicating a lack of self-affine characteristics. These results underscore the significant influence of surface roughness on the mechanics of fault slip under controlled stress conditions.

Discussion

The discussion section of the research paper highlights the influence of fault roughness on injection-induced acoustic emissions (AEs) and fault slip behavior. Utilizing a network of 16 piezoelectric transducers, the study recorded a total of 17,000 AEs across different fault surface roughness conditions. The spatial distribution of AEs revealed that smooth faults exhibited a narrow band of hypocenters, while rough faults displayed broader zones with distinct clustering patterns. The fractal dimension of AEs was lower for rough faults, indicating more localized activity correlated with surface topography, which suggests that fault roughness contributes to stress concentration and the occurrence of larger magnitude events.

The findings also indicate that fault roughness affects the kinematics of AEs, with a notable increase in non-double-couple source types in rough faults, suggesting complex rupture mechanisms. The study employed a numerical model to simulate the effects of fluid injection on slip behavior, revealing that rough faults led to high spatial variability in slip rates and stress distributions. The results imply that roughness stabilizes slow stick-slip events, with larger slip-weakening rates observed in smooth faults. Overall, the research underscores the critical role of fault roughness in modulating induced seismicity, with implications for understanding natural earthquake processes and improving monitoring strategies for fluid injection activities.