DOI: https://doi.org/10.1007/s13146-026-01285-3
تاريخ النشر: 2026-05-12
المؤلف: Halil Bölük وآخرون
الموضوع الرئيسي: أنظمة الكارست والهيدروجيولوجيا
نظرة عامة
تبحث هذه الدراسة في تأثير السوائل الحرارية على الخصائص الفيزيائية والميكانيكية للترافرتين في حفرة كارسية، مع التركيز على التدهور المعتمد على الزمن وآثاره على إدارة المخاطر في البيئات الحرارية-الكارستية. تم إجراء تجربة تعرض في الموقع حيث تم اختبار أنوية بقطر NX في كل من ظروف المختبر وبعد وضعها في الحفرة لفترات زمنية متفاوتة (3، 6، و9 أشهر). تشير النتائج الرئيسية إلى انخفاضات منهجية في قوة الضغط غير المحصورة (UCS) وسرعة الموجة P (V_p) مع زيادة زمن التعرض، إلى جانب زيادة في الموصلية الكهربائية (EC) ودرجة حموضة قريبة من المحايدة إلى حمضية خفيفة عند المخارج الحرارية. أظهرت النمذجة العددية انخفاضًا في عامل تقليل القوة الحرجة (SRF) من 1.86 إلى 1.57 خلال فترة التعرض، مما يبرز انخفاضًا في قوة واستقرار كتلة الصخور.
تخلص الدراسة إلى أن تفاعلات السوائل الحرارية مع الصخور تسرع بشكل كبير من عملية الكارست وتساعد في توسيع التجاويف، كما يتضح من زيادة تقارب السقف واستقرار السطح حول التجاويف. يوفر دمج خرائط إزاحة حالة الخدمة ومقاييس SRF إطارًا قويًا لمراقبة الاستقرار في البيئات الحرارية-الكارستية. ومع ذلك، تعترف الدراسة بالقيود، بما في ذلك نقص نموذج تسرب متكامل تمامًا والحاجة إلى دراسات مستقبلية تشمل أخذ عينات هيدروكيميائية موسمية ونمذجة ثلاثية الأبعاد لفهم العلاقة بين الهيدروكيمياء والخصائص الميكانيكية على مر الزمن بشكل أفضل.
مقدمة
تسلط مقدمة ورقة البحث الضوء على أهمية السوائل الحرارية كمختبرات طبيعية لدراسة تفاعلات السوائل مع الصخور، خاصة في الأراضي التي تهيمن عليها الترافرتين حيث تؤثر العمليات الكارستية على المسامية والسلامة الميكانيكية. تسرع درجات الحرارة المرتفعة والقوة الأيونية العالية في هذه السوائل من التفاعلات التي يمكن أن تؤدي إلى إذابة المعادن الكربونية، مما يؤدي إلى ضعف هياكل الصخور، كما يتضح من الانخفاضات في قوة الضغط الأحادية المحور (UCS) وسرعة الموجة فوق الصوتية P ($V_p$)، إلى جانب زيادة المسامية المفتوحة. تؤكد الدراسة أن الاستجابة الجيوميكانيكية الصافية تعتمد على عوامل مختلفة، بما في ذلك تركيب السائل ودرجة الحرارة ونسيج الصخور الموجود مسبقًا.
توجد الدراسة في حقل كوزاكلي الحراري في وسط الأناضول، حيث يوفر تطور الحفر النشطة فرصة فريدة لربط النتائج المخبرية مع ظروف الميدان. تمثل حفرتان كبيرتان، تشكلتا في عامي 2007 و2018، العمليات الكارستية المستمرة وتدعمها المسوحات الجيوفيزيائية وبيانات حفر الأنوية. تهدف الدراسة إلى معالجة الفجوة في التقييمات المرتبطة بالميدان من خلال ربط التغيرات المعتمدة على الزمن في الخصائص الفيزيائية والميكانيكية مع الظروف الهيدروكيميائية والحرارية في الموقع، باستخدام النمذجة العددية (Phase2) لتعزيز التفسير الجيوميكانيكي وتقييم المخاطر. تشمل الأهداف قياس هذه التغيرات بالنسبة لخط الأساس للصخور الطازجة ودمجها في التحليلات العددية لفهم أفضل لاستقرار الطبقات تحت السطحية.
طرق البحث
تضمنت الطرق المستخدمة في هذه الدراسة سير عمل شامل شمل الاستطلاع الإقليمي ومقياس الموقع، ورسم الخرائط الهيكلية، وتجميع البيانات الجيوفيزيائية وبيانات الحفر الموجودة. ركزت الدراسة على منطقة كوزاكلي الحرارية في وسط الأناضول، التي تتميز بالنيوتيكتونيك النشط وإطار جيولوجي من الترافرتين وتشكيلات الكربونات-التبخر. تضمنت الخطوات الرئيسية تقييمات دورية لجودة المياه في الموقع لتقييم عدوانية السوائل، وجمع عينات أنوية بقطر NX للاختبار في المختبر، وبرنامج تعرض حيث تم نشر هذه الأنوية في حفرة واسترجاعها في فترات محددة مسبقًا.
شملت الاختبارات المخبرية التي أجريت على عينات الأنوية قياسات لسرعة الموجة فوق الصوتية P ($V_p$)، وقوة الضغط الأحادية المحور (UCS)، ومؤشر الحمل النقطي ($I_s(50)$)، والمسامية المفتوحة، والكثافة الكلية، ومتانة التفتت، مع الالتزام بمعايير ISRM. استخدمت الدراسة تحليلات العناصر المحدودة مع برنامج Phase2 لنمذجة السلوك الجيوتقني للمواد، مع دمج كل من الظروف السليمة والمتدهورة بناءً على النتائج المخبرية. استخدمت النمذجة نهج إجهاد مستوي لتمثيل الطبقات الجيولوجية وهندسة التجاويف، مع تطبيق نماذج تركيبية محددة على وحدات جيولوجية مختلفة، بما في ذلك نموذج موهر-كولوم المنفصل لوحدة الطين-الحطام السطحي ونموذج مميز للترافرتين. قدم دمج التصوير الجيوفيزيائي والبيانات الهيدروكيميائية إطارًا قويًا لفهم ديناميات السوائل وتكوين الحفر في المنطقة.
النتائج
تشير نتائج الدراسة إلى أن التعرض في الموقع يؤدي إلى تغييرات ملحوظة في الخصائص الميكانيكية والفيزيائية مع مرور الوقت. يكشف التحليل الدوري عن تغييرات طفيفة في الكتلة الجافة (0.06-0.11%)، بينما تزداد الكتلة المشبعة بشكل كبير (1.08-1.77%)، مما يشير إلى امتصاص محدود للمياه. تنخفض سرعة الموجة P الجافة بنسبة 9.35% إلى 17.07% بالنسبة لخط الأساس الطازج عبر ثلاث نقاط زمنية، مما يدل على زيادة في كثافة الشقوق. يظهر نسبة الفراغ زيادة تدريجية من حوالي 2.47% عند 3 أشهر إلى 9.96% عند 9 أشهر، مما يتماشى مع اتساع فتحات المسام تحت ظروف كيميائية عدوانية.
علاوة على ذلك، تظهر قوة الضغط غير المحصورة (UCS) انخفاضًا تدريجيًا من 20.34 ميغاباسكال إلى 19.11 ميغاباسكال في الظروف الجافة ومن 18.65 ميغاباسكال إلى 17.99 ميغاباسكال في الظروف المشبعة، مما يعكس اتجاهًا نحو الضعف التدريجي طوال فترة التعرض. توضح البيانات الاتجاهات المتماسكة عبر مؤشرات مختلفة، حيث تنخفض UCS من 20.90 ميغاباسكال (جاف) و15.86 ميغاباسكال (مشبع) عند خط الأساس إلى 19.11 ميغاباسكال و17.99 ميغاباسكال، على التوالي، بحلول حوالي 9 أشهر. تؤكد النتائج تأثير التعرض في الموقع على السلامة الميكانيكية للمادة، مما يبرز العلاقة بين امتصاص المياه، وتغيرات هيكل المسام، والتدهور الناتج في الخصائص الميكانيكية.
المناقشة
تسلط قسم المناقشة في ورقة البحث الضوء على العلاقة المعقدة بين العمليات الهيدروكيميائية، والتدهور الميكانيكي، والضوابط الهيكلية في البيئات الكارستية، وخاصة داخل حقل كوزاكلي الحراري. كشفت البيانات الهيدروكيميائية الميدانية التي تم جمعها خلال المواسم الجافة والرطبة عن تباينات كبيرة في درجة الحرارة ودرجة الحموضة والموصلية الكهربائية (EC)، مما يشير إلى ظروف كيميائية عدوانية تسهل إذابة الكربونات. أظهرت النمذجة العددية للدراسة، بناءً على المعلمات الجيوتقنية، أن السوائل الحرارية الغنية بـ CO₂ وH₂S لا تسرع فقط من الإذابة ولكنها تقلل أيضًا بشكل كبير من القوة الميكانيكية للترافرتين والطبقات الحجرية الأساسية على مدى فترات زمنية قصيرة. تؤكد هذه النتيجة على ضرورة دمج العوامل الهيدروكيميائية جنبًا إلى جنب مع الاعتبارات الهيكلية والطبقية في تقييم مخاطر الحفر.
كما أثبتت الدراسة أن عمليات التجوية الفيزيائية، مثل دورات التجمد-الذوبان والرطوبة-الجفاف، تتفاعل بشكل تآزري مع التدهور الكيميائي، مما يؤدي إلى تسريع تدهور الصخور. توفر النماذج العددية التي تم التحقق منها ضد الملاحظات المخبرية والميدانية إطارًا قويًا لتوقع نمو التجاويف وأنماط التشوه. ومع ذلك، يعترف المؤلفون بالقيود في نهجهم، بما في ذلك التبسيطات المتأصلة في معيار هوك-براون وإطار النمذجة ثنائية الأبعاد، الذي قد لا يعكس تمامًا تعقيدات العمليات ثلاثية الأبعاد. يُقترح أن تستخدم الأبحاث المستقبلية تقنيات نمذجة أكثر تقدمًا، مثل نمذجة العناصر المنفصلة أو المحاكاة الهيدروليكية-الميكانيكية المترابطة، لتعزيز فهم آليات الفشل التدريجي وتفاعلات السوائل المسامية في أنظمة الحرارية-الكارستية. بشكل عام، تسهم النتائج في تقديم رؤى قيمة حول الطبيعة متعددة الأوجه لتطور الحفر في البيئات الحرارية-الكارستية، مع آثار تمتد إلى ما وراء منطقة كوزاكلي إلى بيئات مماثلة عالميًا.
DOI: https://doi.org/10.1007/s13146-026-01285-3
Publication Date: 2026-05-12
Author(s): Halil Bölük et al.
Primary Topic: Karst Systems and Hydrogeology
Overview
This research investigates the impact of geothermal fluids on the physico-mechanical properties of travertine in a karst sinkhole, focusing on time-dependent degradation and its implications for risk management in geothermal-karst environments. An in situ exposure experiment was conducted where NX-diameter cores were tested in both laboratory conditions and after being emplaced in the sinkhole for varying durations (3, 6, and 9 months). Key findings indicate systematic declines in unconfined compressive strength (UCS) and P-wave velocity (V_p) with increased exposure time, alongside elevated electrical conductivity (EC) and a near-neutral to mildly acidic pH at thermal outlets. Numerical modeling demonstrated a decrease in the critical Strength Reduction Factor (SRF) from 1.86 to 1.57 over the exposure period, highlighting a reduction in rock-mass strength and stability.
The study concludes that geothermal fluid-rock interactions significantly accelerate karstification and promote cavity enlargement, evidenced by increased roof convergence and surface settlement around cavities. The integration of service-state displacement maps and SRF metrics provides a robust framework for monitoring stability in geothermal-karst settings. However, the research acknowledges limitations, including the lack of a fully coupled seepage model and the need for future studies to incorporate seasonal hydrochemical sampling and 3D modeling to better understand the relationship between hydrochemistry and mechanical properties over time.
Introduction
The introduction of the research paper highlights the significance of geothermal fluids as natural laboratories for studying fluid-rock interactions, particularly in travertine-dominated terrains where karst processes influence porosity and mechanical integrity. Elevated temperatures and high ionic strength in these fluids accelerate reactions that can lead to the dissolution of carbonate minerals, resulting in weakened rock structures, as evidenced by reductions in uniaxial compressive strength (UCS) and ultrasonic P-wave velocity ($V_p$), alongside increased open porosity. The study emphasizes that the net geomechanical response is contingent upon various factors, including fluid composition, temperature, and the pre-existing rock fabric.
The research is situated in the Kozaklı geothermal field in Central Anatolia, where active sinkhole development provides a unique opportunity to bridge laboratory findings with field conditions. Two significant sinkholes, formed in 2007 and 2018, exemplify ongoing karstification processes and are supported by geophysical surveys and core-drilling data. The study aims to address the gap in field-anchored assessments by correlating time-dependent changes in physico-mechanical properties with in situ hydrochemical and thermal conditions, utilizing numerical modeling (Phase2) to enhance geomechanical interpretation and risk assessment. The objectives include quantifying these changes relative to a fresh-rock baseline and integrating them into numerical analyses for improved understanding of subsurface stability.
Methods
The methods employed in this study involved a comprehensive workflow that included regional and site-scale reconnaissance, structural mapping, and the compilation of existing geophysical and borehole data. The research focused on the Kozaklı geothermal area in Central Anatolia, characterized by active neotectonics and a geological framework of travertine and carbonate-evaporite formations. Key steps included periodic in situ water-quality assessments to evaluate fluid aggressiveness, collection of NX-diameter core specimens for laboratory testing, and an exposure program where these cores were deployed in a sinkhole and retrieved at predetermined intervals.
Laboratory tests conducted on the core specimens included measurements of ultrasonic P-wave velocity ($V_p$), uniaxial compressive strength (UCS), point-load index ($I_s(50)$), open porosity, bulk density, and slake durability, adhering to ISRM standards. The study utilized finite-element analyses with Phase2 software to model the geotechnical behavior of the materials, incorporating both intact and degraded conditions based on laboratory findings. The modeling employed a plane-strain approach to represent the stratigraphy and cavity geometry, with specific constitutive models applied to different geological units, including drained Mohr-Coulomb for the surficial mud-rubble unit and a distinct model for travertine. The integration of geophysical imaging and hydrochemical data provided a robust framework for understanding fluid dynamics and sinkhole formation in the region.
Results
The results of the study indicate that in-situ exposure leads to notable changes in mechanical and physical properties over time. Period-wise analysis reveals minimal changes in dry mass (0.06-0.11%), while saturated mass increases significantly (1.08-1.77%), suggesting limited water uptake. The dry P-wave velocity decreases by 9.35% to 17.07% relative to the fresh baseline across three time points, indicating an increase in crack density. The void ratio exhibits a progressive increase from approximately 2.47% at 3 months to 9.96% at 9 months, which aligns with pore-throat widening under chemically aggressive conditions.
Furthermore, the unconfined compressive strength (UCS) shows a gradual decline from 20.34 MPa to 19.11 MPa in dry conditions and from 18.65 MPa to 17.99 MPa in saturated conditions, reflecting a trend of progressive weakening throughout the exposure period. The data demonstrate coherent trends across various indices, with the UCS decreasing from 20.90 MPa (dry) and 15.86 MPa (saturated) at baseline to 19.11 MPa and 17.99 MPa, respectively, by approximately 9 months. The findings underscore the impact of in-situ exposure on the mechanical integrity of the material, highlighting the relationship between water uptake, pore structure changes, and the resulting degradation of mechanical properties.
Discussion
The discussion section of the research paper highlights the intricate relationship between hydrogeochemical processes, mechanical degradation, and structural controls in karst environments, specifically within the Kozaklı geothermal field. Field hydrochemical data collected during dry and wet seasons revealed significant variations in temperature, pH, and electrical conductivity (EC), indicating chemically aggressive conditions that facilitate carbonate dissolution. The study’s numerical modeling, based on geotechnical parameters, demonstrated that geothermal fluids enriched with CO₂ and H₂S not only accelerate dissolution but also substantially reduce the mechanical strength of travertine and underlying limestones over short timescales. This finding underscores the necessity of incorporating hydrochemical factors alongside structural and lithological considerations in assessing sinkhole hazards.
The research further established that physical weathering processes, such as freeze-thaw and wetting-drying cycles, interact synergistically with chemical degradation, leading to accelerated rock deterioration. The numerical models validated against laboratory and field observations provide a robust framework for predicting cavity growth and deformation patterns. However, the authors acknowledge limitations in their approach, including the simplifications inherent in the Hoek-Brown criterion and the two-dimensional modeling framework, which may not fully capture the complexities of three-dimensional processes. Future research is suggested to employ more advanced modeling techniques, such as discrete element modeling or coupled hydro-mechanical simulations, to enhance the understanding of progressive failure mechanisms and pore-fluid interactions in geothermal-karst systems. Overall, the findings contribute valuable insights into the multi-faceted nature of sinkhole development in geothermal-karst settings, with implications extending beyond the Kozaklı region to similar environments globally.