מודל מכני של חלחול לא-ליניארי וניתוח סינון זוגי עבורCoupling נוזל-מוצק במסת סלע שסועה

מדוע המים הנסתרים בסלעים סדוקים חשובים

הרחק מתחת לרגלינו זורמים מים דרך שכבות סדוקות של סלע במכרות, מנהרות ומאגרים נפט וגז. האופן שבו המים מפעילים לחץ על הסלע — וכיצד הסלע מגיב חזרה — יכול להניח את היסוד להצפות פתאומיות, קריסות או רעידות אדמה. המאמר הזה מפתח שיטה חדשה לתיאור ולחיזוי האינטראקציה ההדדית בין זרימת מים והתעוותות סלע בסלעים שסועים, וחושף מדוע מערכות כאלה יכולות להיראות קבועות במשך זמן רב ואז להכשל בפתאומיות.

מים החוצים את האבן השבורה

ברבות מהפרויקטים התת-קרקעיים, מסת הסלע אינה מוצקה כבלוק אלא דומה לספוג סדוק: אין-ספור שברים ורווחים יוצרים נתיבי מעבר למים. לחץ חיצוני מהסלע העילי ולחץ פנימי של המים בנקבים מעצבים כל העת את אותם נתיבים. כשהסלע נלחץ, הנקבים מצטמצמים או נסגרים; כשהלחץ המימי עולה, הם יכולים להיפתח מחדש או להתרחב. העיצוב המתמשך הזה משנה את היכולת של המים לנוע, מה שמפעיל משוב על הסלע. המחברים טוענים שכדי להבין אסונות כגון הצפת מי מכרה או דליפה במאגרים, יש להתייחס למערכת זו כאל מערכת דינמית וצמודה ולא כאל תמונת מצב סטטית.

בניית תמונה מצומדת של סלע ומים

המחקר מתחיל בהרחבת מושג קלאסי במכניקת קרקע הנקרא "מתח יעיל", שמתאר כמה מתוך הלחץ הכולל נישא בפועל על ידי השלד המוצק. המחברים מעצבים מחדש רעיון זה כך שהוא יכלול במפורש את הנפּיחות — החלק מנפח הסלע שהוכרע על ידי חללים — כך ששינויים במרחב הנקבים יקושרו ישירות לחלוקת הלחץ בין סלע למים. הם משלבים זאת עם משוואות המתארות כיצד סלע מתכופף מעט מגיב באופן אלסטי ללחץ ועם תיאור מציאותי יותר ולא-ליניארי של זרימת מים דרך שברים, שעובר את חוק דארסי הפשוט המשמש בדגמים הנדסיים רבים.

מזרימה חלקה לשינוי פתאומי

עם המסגרת הזו, המחברים מתמקדים במקרה חד-ממדי של מים שחודרים אנכית דרך שכבת סלע שבורה. הם גוזרים זוג משוואות לא-ליניאריות שמנטרות כיצד לחץ המים וקצב הזרימה מתפתחים לאורך זמן ומרחב, כאשר הנפיחות מסתגלת בזמן שהסלע מתכווץ. פתרון המשוואות מראה כי בתנאים מסוימים המערכת אינה בעלת התנהגות יציבה יחידה: במקום זאת היא מציגה מה שמדענים מתמטיקאים קוראים לו סינון מסוג saddle-node. במילים פשוטות, כאשר פרמטר זרימה מרכזי משתנה, מצב שהיה יציב עלול להתפצל לענף בטוח וענף מסוכן, או להיעלם לחלוטין, ולגרום לקפיצה פתאומית מזרימה עדינה לזרימה בלתי נשלטת.

לחיצה אטית ויציבות מעוכבת

המחברים בוחנים כיצד המצב משתנה כאשר הלחץ בגבולות — למשל כתוצאה מהעמסה הדרגתית של הכרייה שמעל — משתנה לאורך תקופות ארוכות. סימולציות נומריות מראות שכאשר הלחץ החיצוני משתנה באיטיות, המערכת המשותפת של סלע ומים גם היא לוקחת הרבה יותר זמן להגיע למצב יציב. לחץ המים, קצב הזרימה ועיוות נפח הסלע זוחלים לעבר יציבות במקום להתאזן במהירות. עיכוב זה קורה מכיוון ששלד הסלע חייב להמשיך להתאים את מבנה הנקבים שלו בזמן שאנרגיה מוזרמת למערכת על ידי העומס המשתנה, ובכך מתארך מסלול ההגעה לשיווי משקל.

אזהרות לפני הצפה

כדי לקשר את התיאוריה למציאות, המחקר משווה את תחזיותיו למקרה אמיתי של פרץ מים מפגיעה בפיצול בפחמיית. כאשר הכרייה התקרבה לפיצול, פרמטר שהשתקף בכמה חזק הזרימה חורגת מהתנהגות דארסית פשוטה, הזיז את עצמו לתחום קריטי שבו יכלו להתקיים שתי מצבי זרימה: אחד יציב ואחד לא יציב. מדידות שדה הראו שהמהירות של המים החלה להתנדנד בין שתי רמות מובחנות לפני שפרצה בסופו של דבר לעלייה מהירה וקטסטרופלית, בדיוק כפי שדיאגרמת הסינון של המודל ניבאה. המחברים טוענים שהתנודות האלה הן סימן אזהרה ברור ומוקדם יותר להצפת מים מאשר מדדי הבטיחות המסורתיים שמתייחסים למערכת כאל ליניארית וסטטית.

מה המשמעות בשביל בטיחות תת-קרקעית

בסך הכל, המאמר מראה שסלע שסוע רווי מים מתנהג יותר כמו מערכת מורכבת ולא-ליניארית מאשר כמו צנרת פשוטה. שינויים קטנים במתח או בתנאי הזרימה יכולים לדחוף אותו מעבר לספי מפתח שבהם התנהגותו משתנה באופן איכותי, ולא רק בכמותה. על ידי קישור מפורש בין עיוות הסלע, מבנה הנקבים וזרימה לא-ליניארית, המודל החדש יכול ללכוד מצבים יציבים מרובים, מעברים פתאומיים ביניהם ורגישות חזקה לתנאי ההתחלה. למהנדסים המתכננים מכרות, מנהרות ומאגרים, משמעות הדבר היא שמעקב אחרי התפתחות הזרימה והעיוות לאורך זמן — וצפייה בתנודות דו-יציבות אופייניות — עשוי לספק אזהרות מוקדמות ומהימנות יותר על חוסר יציבות חבוי לפני שהוא מתפרץ לאסון מלא.

ציטוט: Zhengzheng, C., Mengqi, X., Tao, R. et al. Nonlinear seepage mechanical model and bifurcation analysis for fluid-solid coupling in fractured rock mass. Sci Rep 16, 9578 (2026). https://doi.org/10.1038/s41598-025-25823-6

מילות מפתח: סלע שסוע, חלחול מי תהום, צימוד נוזל–מוצק, דינמיקה לא-ליניארית, הצפה פתאומית במכרה