שיטת אינטרפולציה POCS מונעת-אילוץ לקווי הסתגלי

תמונות תת-קרקעיות חדים יותר מתוך נתונים חסרים

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

מדוע מחסור בעקבות חשוב

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

דרך חכמה יותר להשתמש במידע השכן

המחברים מבססים את עבודתם על מסגרת מתמטית ידועה בשם projection onto convex sets (הטלה על קבוצות קמורות), שמלטפת שוב ושוב את הנתונים לעבר תנאים "מותרתים" שונים עד שכולם מתיישבים. הגרסה הקלאסית אוכפת שני תנאים עיקריים: שהמדידות הידועות לא ישתנו, ושהנתונים ייראו קומפקטיים ומסודרים בטרנספורמציה ספאורסית כגון תחום תדירות–גל-רוחב. השיטה החדשה מוסיפה מרכיב שלישי: אילוץ לטרלי מפורש המעודד עקבות סמוכות להשתנות בצורה חלקה במקום שבו הגאולוגיה רציפה, אך עדיין מאפשר שינויים חדים כאשר התת-קרקע באמת נשברת, כמו בשברים או בגבולות סלע חדים.

הסתגלות לגאולוגיה מורכבת

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

איזון בין דיוק ליעילות

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

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

הצוות בוחן את הגישה שלהם על שני סוגי נתונים. ראשית, הם משתמשים ברשומות סייסמיות סינתטיות שנוצרו ממודל גאולוגי מורכב של Marmousi, עם עקבות שהוסרו באקראי הכוללים 30%, 50% ואפילו 70% אובדן נתונים. בכל המקרים, השיטה החדשה משחזרת אירועים שהם רציפים יותר באופן לטרלי, עם רעש נראה פחות ותוכן תדירות קרוב יותר לנתוני המקור המלאים בהשוואה לאלגוריתם הישן. שנית, הם מיישמים את השיטה על נתוני פוסט-סטאק אמיתיים ממאגר נפט וגז במזרח סין. כאן גם הגישה החדשה מניבה חתכים נקיים יותר, המשכיות טובה יותר באזורים משתנים בעדינות ושימור משופר של פרטים גאולוגיים חשובים, כפי שאושש באמצעות מדדים כמותיים של שגיאה, יחס אות/רעש ודמיון מבני.

מה משמעות הדבר לראיית התת-קרקע

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

ציטוט: Qin, Z., Pan, S., Chen, J. et al. Adaptive lateral constraint-driven POCS interpolation method. Sci Rep 16, 11518 (2026). https://doi.org/10.1038/s41598-026-39281-1

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