מחקר על גובה החזוי של אזורי סדיקה מוליכי מי-קרקע בהתבסס על מודל BO-RFR וניתוח SHAP

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

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

נתיבים חבויים למים בקרקע

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

מכלל אצבע לחיזוי מבוסס-נתונים

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

לימוד יער של אלגוריתמים לקריאת הסלעים

כדי לפענח את הנתונים המעורבים והלא מושלמים האלה, הצוות פנה לשיטת למידת מכונה שנקראת רגרסיית יער אקראי (random forest regression), שיטה שמשלבת עצי החלטה רבים למנבא אחד ויציב. לאחר מכן השתמשו באופטימיזציה בייסיאנית — אסטרטגיית חיפוש יעילה — כדי לכוונן באופן אוטומטי את הפרמטרים הפנימיים של המודל כך שיתפקד היטב אפילו עם מספר דגימות יחסית מצומצם. מודל ה-BO-RFR הזה אומן בנפרד עבור שדות הפחם המזרחיים והמערביים, ונבדק בקפדנות ובאמצעי חצייה על נתונים שלא נראו קודם, כולל מדגמים "אילמים" ממכרות נוספים. בכל הבדיקות המודל המכוון חזה את גובהי הסדיקה בדיוק גבוה בהרבה מנוסחאות מסורתיות וממספר אלגוריתמים מתקדמים אחרים, ותפס את האופי המורכב והלא־ליניארי של האינטראקציה בין גאולוגיה לתכנון הכרייה.

פירוק מה שבאמת חשוב

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

להפוך תובנות למכרות בטוחים יותר

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

ציטוט: Qiu, M., Wen, Y., Teng, C. et al. Research on the predicted height of water-conducting fracture zones based on the BO-RFR model and SHAP analysis. Sci Rep 16, 7230 (2026). https://doi.org/10.1038/s41598-026-38043-3

מילות מפתח: בטיחות בכריית פחם, סכנות מי-תהום, סדיקות בסלע, מודלים של למידת מכונה, שדות פחם בסין