שיטת היפוך תלת‑ממדית ל‑TEM המשלבת אופטימיזציה PSO‑NLCG ויסות אדפטיבי

לראות תצורות חבויות מתחת לרגלינו

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

להקשיב לאותות הדועכים באדמה

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

איזון בין התאמת נתונים לתמונה תת‑קרקעית סבירה

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

שילוב חיפוש רחב וכיוונון עדין

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

בדיקה על פונקציות תקן ומודלים סינתטיים של הקרקע

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

הוכחה בשטח באתר בניה אמיתי

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

מפות תת‑קרקעיות ברורות יותר להחלטות בטוחות יותר

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

ציטוט: Jianqiang, C., Feng, Z., Xuhai, D. et al. A 3D inversion method of TEM combining PSO-NLCG optimization and adaptive regularization. Sci Rep 16, 16115 (2026). https://doi.org/10.1038/s41598-026-48117-x

מילות מפתח: אלקטרומגנטיות חולפת, היפוך תלת‑ממדי, הדמיה גיאופיזיקלית, אלגוריתם אופטימיזציה, תצורות תת‑קרקעיות