Modellordnungsreduzierte spektrale Elemente für hochgenaue und schnelle 3‑D‑transiente elektromagnetische Vorwärtsmodellierung mit SAI‑Krylov

Blick in den Untergrund mit schnelleren virtuellen Sensoren

Das Auffinden vergrabener Erzkörper, Grundwasser oder verborgener Verwerfungen beruht oft auf feinen Störungen im elektromagnetischen Feld der Erde. Aus diesen Störungen ein Bild des Untergrunds zu erzeugen erfordert umfangreiche numerische Berechnungen an großdimensionierten 3‑D‑Modellen, die sehr zeitaufwendig sein können. Dieser Beitrag stellt eine neue Simulationsmethode für elektromagnetische Messungen vor, die die für reale Exploration nötige Genauigkeit bewahrt und gleichzeitig den Zeit‑ und Speicheraufwand am Rechner drastisch reduziert.

Warum die Simulation elektromagnetischer Messungen schwierig ist

Transiente elektromagnetische (TEM) Methoden arbeiten, indem ein Stromimpuls durch eine Drahtschleife oder ein geerdetes Kabel an der Oberfläche eingespeist wird. Beim Abschalten des Stroms werden unterirdisch Wirbelströme induziert, die zeitlich abklingen und Informationen über Leitfähigkeit und vergrabene Erzkörper tragen. Um vorherzusagen, was Empfänger messen, lösen Forschende die Maxwellschen Gleichungen auf einem 3‑D‑Gitter. Traditionelle Ansätze verwenden einfache Elemente mit nur linearer Variation in jeder Zelle. Diese niedrigordentlichen Approximationen sind leicht zu implementieren, erfassen aber feine Details nur schlecht, ohne die Anzahl der Unbekannten und die Kosten für das Durchlaufen tausender kleiner Zeitschritte extrem zu erhöhen.

Scharfere Gitter mit weniger Elementen

Die Autorinnen und Autoren nutzen eine hochordentliche spektrale Elementmethode, die jede Gitterzelle mit flexibleren mathematischen Formen anreichert, statt auf geradlinige Variationen zu setzen. Praktisch erlaubt das, das elektrische Feld glatter darzustellen bei gleicher Netzauflösung oder eine gegebene Genauigkeit mit deutlich weniger Zellen zu erreichen. Sie legen sorgfältig fest, wie Unbekannte entlang von Kanten platziert werden, wie jede physikalische Zelle auf einen einfachen Referenzwürfel abgebildet wird und wie spezielle Gauss–Lobatto‑Stützstellen und Gewichte verwendet werden. Ein wichtiger Trick, die reduzierte Integration, lockert die Exaktheit bestimmter Integrale leicht, bringt dafür aber die Leitfähigkeitsmatrix in eine strikt diagonale Form. Das erhöht die Sparsität deutlich, reduziert sowohl Speicherbedarf als auch die Kosten für das Lösen der resultierenden linearen Systeme und bewahrt gleichzeitig die hohe Genauigkeit für die interessierenden Ordnungen.

Rechnen komprimieren, Physik bewahren

Selbst mit einem präziseren Gitter bleibt das Vorwärtslaufen in der Zeit auf herkömmliche Weise teuer, weil jede Änderung des Zeitschritts eine neue Faktorisierung einer riesigen Matrix erzwingt. Stattdessen formulieren die Autorinnen und Autoren das Problem als zeitkontinuierliches System und wenden eine Modellordnungsreduktion auf Basis eines Shift‑and‑Invert‑Krylov‑Unterraumverfahrens an. Anschaulich gesagt projizieren sie das vollständige, enorme System auf eine viel kleinere Menge „repräsentativer“ Muster, die erfassen, wie sich elektromagnetische Felder tatsächlich entwickeln. Eine geschickte spektrale Transformation bündelt die wichtigen Modi, sodass der Algorithmus schnell konvergiert. Entscheidend wählt die Methode dynamisch, wie groß der reduzierte Raum sein muss, wobei eine preiswerte Residualgröße als Abbruchkriterium dient. Dadurch benötigt der Löser nur eine große Matrixfaktorisierung plus eine moderate Anzahl Rückeinsetzungen, kann aber die elektromagnetische Antwort zu beliebigen Zeiten auswerten, ohne jeden Zwischenzeitpunkt explizit durchlaufen zu müssen.

Tests der Methode

Zur Bewertung von Genauigkeit und Effizienz simulieren die Autorinnen und Autoren mehrere Benchmark‑Szenarien. Für einen einfachen Halbraum mit bekanntem analytischem Verhalten vergleichen sie verschiedene Polynomordnungen und Integrationsstrategien mit einem standardmäßigen rückwärtsgerichteten Zeitintegrationsschema. Ab der zweiten Ordnung stabilisieren sich die Fehler bei etwa einem Prozent oder weniger, und mit dem neuen reduzierten Modell lassen sich mittlere Fehler auf rund ein halbes Prozent drücken, während für höherordentliche Läufe Beschleunigungen um etwa das 16‑ bis 20‑Fache erzielt werden, ohne den Speicherbedarf in die Höhe zu treiben. In geschichteten Modellen mit starken Resistivitätskontrasten bleibt der neue spektrale‑Element‑plus‑Krylov‑Ansatz über das gesamte Zeitfenster bei etwa einem Prozent Fehler und übertrifft konventionelle Finite‑Element‑Referenzen, die auf niedrigordentliche Zellen und stärker dissipative Zeitintegration setzen. Schließlich folgt die Methode in einem vollständigen 3‑D‑Sulfiderz‑Szenario der Ausbreitung induzierter Felder von geerdeten Leitern und Schleifenquellen, ihrer Wechselwirkung mit einem komplexen Leitkörper und zeichnet so das vergrabene Ziel mit hoher räumlicher Auflösung nach.

Was das für die Erkundung des Untergrunds bedeutet

Für geophysikalische Erkundungen, die auf detaillierte 3‑D‑elektromagnetische Modellierung angewiesen sind, bietet diese Arbeit einen Weg, beide Ziele zu erreichen: hohe Genauigkeit und hohe Geschwindigkeit. Durch die Kombination hochordentlicher spektraler Elemente mit einer sorgfältig entworfenen Modellordnungsreduktion zeigen die Autorinnen und Autoren, dass sich für praktikable TEM‑Antworten Sub‑Prozent‑Fehler erzielen lassen, während die Rechenzeit gegenüber einer anerkannten Referenz um eine Größenordnung reduziert wird. Alltäglich ausgedrückt bedeutet das schnellere Durchläufe von der Rohdatenplanung bis zu interpretierbaren Bildern des Untergrunds, was die Exploration nach Mineralien, die Bewertung geohazard‑relevanter Risiken oder die Verfolgung von Grundwasser mit physikbasierten Simulationen anstelle vereinfachter Näherungen praktikabler macht.

Zitation: Fan, Y., Lu, K., Huang, Y. et al. Model-order-reduced spectral-element method for high-accuracy and fast 3-D transient electromagnetic forward modeling with SAI-Krylov. Sci Rep 16, 13356 (2026). https://doi.org/10.1038/s41598-026-44053-y

Schlüsselwörter: transiente Elektromagnetik, spektrale Elementmethode, Modellordnungsreduktion, geophysikalische Exploration, numerische Simulation