An Upward-Continuation-Based Framework for Subsurface Structural Reconstruction from Gravity Anomalies

Sintesi Tecnica: Un Framework Basato sulla Continuazione verso l'Alto per la Ricostruzione Strutturale del Sottosuolo da Anomalie di Gravità

Problematica
L'interpretazione dei dati di gravità è fondamentale per investigare l'eterogeneità della densità del sottosuolo, ma è intrinsecamente non univoca; distribuzioni di densità sostanzialmente diverse possono generare risposte di anomalia simili. Sebbene esistano tecniche di inversione e modellazione diretta, esse spesso si affidano pesantemente a vincoli geologici preesistenti, parametrizzazioni complesse o ottimizzazioni computazionalmente intensive. Metodi basati sulle trasformazioni (ad esempio, continuazione verso l'alto, segnale analitico, gradiente completo normalizzato) sono stati sviluppati per estrarre informazioni strutturali senza inversione, tuttavia la maggior parte di essi si concentra sul potenziamento dell'anomalia, sulla rilevazione dei bordi o sulla separazione delle lunghezze d'onda, piuttosto che sulla ricostruzione della geometria spaziale completa dei corpi causativi preservando la polarità dell'anomalia. Inoltre, i metodi basati sul gradiente come il Gradiente Completo Normalizzato (NFG) sono efficaci per il rilevamento dei bordi e la stima della profondità, ma non rappresentano direttamente la geometria completa o la polarità di densità della sorgente. È necessario un framework indipendente dal modello che possa ricostruire la geometria strutturale del sottosuolo direttamente dalle anomalie di gravità, fornendo sia la caratterizzazione geometrica che i vincoli di profondità senza richiedere dati geologici esterni.

Metodologia
Gli autori propongono il framework Upward Continuation Restriction to Decomposition (UCRD), integrato con il metodo Normalized Full Gradient (NFG), per affrontare tali limitazioni. Il flusso di lavoro consiste in tre fasi primarie:

1. Continuazione verso l'alto multi-livello e decomposizione spettrale:
   Il campo dell'anomalia di gravità di Bouguer viene decomposto in livelli di lunghezza d'onda ristretta (RWL) mediante continuazione verso l'alto a più livelli. Il processo inizia con il campo osservato ($U_0$) e applica successive continuazioni verso l'alto ($U_n$) con altezze esponenzialmente crescenti ($\Delta z_n$) per isolare sorgenti progressivamente più profonde. Il livello di lunghezza d'onda ristretta per ogni passaggio è calcolato come la differenza tra il campo all'altezza precedente e l'altezza corrente ($RWL_n = U_{n-1} - U_n$). Questo separa il campo di gravità in componenti associate a specifici intervalli di profondità.

2. Trasformazione del Segnale Analitico e Normalizzazione della Scala:
   Per potenziare la continuità strutturale e preservare l'interpretabilità fisica, ogni componente RWL subisce una trasformazione del segnale analitico tridimensionale. Viene quindi applicata una procedura di normalizzazione adattiva alla scala per ripristinare le unità fisiche (mGal) e preservare la polarità dell'anomalia. Ciò comporta il calcolo di un fattore di scala specifico per il livello basato sulla gamma dinamica della RWL e sull'ampiezza del segnale analitico. Il risultato è una rappresentazione Scaled Restricted Wavelength (SRWL). Le mappe SRWL catturano l'entità e la direzione dei contrasti di densità laterali ad ogni livello di profondità, fornendo una rappresentazione pseudo-3D dell'evoluzione dell'anomalia.

3. Integrazione con NFG per i Vincoli di Profondità:
   Poiché la continuazione verso l'alto non fornisce informazioni dirette sulla profondità fisica, la ricostruzione geometrica UCRD è completata dal metodo NFG. L'NFG viene applicato utilizzando una formulazione in serie di Fourier con estensione a specchio del segnale per evitare discontinuità ai bordi. Una procedura adattiva determina il numero ottimale di armoniche ($N$) tracciando le risposte massime dell'NFG attraverso intervalli armonici crescenti. Ciò fornisce vincoli di profondità quantitativi che vengono integrati con la geometria strutturale UCRD.

Contributi Chiave

• Framework UCRD: Introduzione di un metodo di decomposizione che trasforma le anomalie di gravità in immagini fisicamente interpretabili della geometria strutturale del sottosuolo senza inversione.
• Rappresentazione SRWL: Sviluppo di una procedura di scalatura che ripristina le unità fisiche e la polarità dell'anomalia ai risultati del segnale analitico, migliorando la visibilità delle caratteristiche strutturali oscurate nei dati grezzi.
• Workflow Complementare: Un approccio unificato in cui l'UCRD delinea la geometria strutturale e la polarità, mentre l'NFG fornisce vincoli di profondità indipendenti, riducendo l'ambiguità inerente alle interpretazioni basate sul gradiente.
• Indipendenza dal Modello: La metodologia opera direttamente sui dati di gravità, non richiedendo informazioni geologiche pregresse, dati sismici o ottimizzazioni computazionalmente intensive.

Risultati
Il framework è stato validato utilizzando tre modelli sintetici di crescente complessità:

1. Corpi Isolati e a Gradini: Il metodo ha delineato con successo sei corpi cubici e tre strutture a gradini con profondità e contrasti di densità variabili. L'UCRD ha ricostruito la distribuzione spaziale e la polarità, mentre l'NFG ha identificato le profondità corrispondenti. I risultati hanno mostrato una forte concordanza con il modello reale, anche per corpi profondi parzialmente mascherati da anomalie superficiali.
2. Sorgenti Sovrapposte: In un modello con prismi parzialmente sovrapposti, l'UCRD ha separato con successo le componenti strutturali dominanti delle sorgenti superficiali e profonde, risolvendo le firme composite delle anomalie. L'NFG ha fornito stime di profondità stabili per entrambe le componenti superficiale e profonda, distinguendo geometrie sovrapposte che sono difficili da interpretare tramite l'analisi della gravità standard.
3. Strutture Lateralmente Estese: Per un modello composto da venticinque blocchi distribuiti a varie profondità, l'UCRD ha recuperato la geometria complessiva della struttura su larga scala, rivelando sia le componenti superficiali che quelle profonde. L'NFG ha identificato distinti livelli di profondità e confini strutturali all'interno del sistema sorgente esteso.

In tutti i casi, le trasformazioni SRWL hanno prodotto firme di anomalia più nitide e localizzate rispetto alle componenti di lunghezza d'onda ristretta standard. Il workflow combinato UCRD–NFG ha dimostrato la capacità di recuperare sia i pattern strutturali orizzontali che quelli verticali attraverso una vasta gamma di configurazioni delle sorgenti.

Significato e Rivendicazioni
L'articolo afferma che il framework proposto UCRD–NFG offre uno strumento pratico ed efficiente dal punto di vista computazionale per l'interpretazione della gravità. Il suo significato primario risiede nella capacità di fornire una caratterizzazione preliminare della geometria del sottosuolo e della profondità senza la necessità di inversione o vincoli esterni. Gli autori pongono il metodo non come un sostituto della modellazione diretta o dell'inversione congiunta, ma come un prezioso passaggio preliminare che può generare modelli iniziali, vincoli geometrici e guida per la parametrizzazione per indagini successive più intensive. Il framework è particolarmente evidenziato per le investigazioni su scala regionale dove è richiesta una rapida caratterizzazione strutturale. Lo studio riconosce le limitazioni inerenti alla decomposizione spettrale, come la parziale sovrapposizione della lunghezza d'onda tra i livelli e gli effetti di smoothing derivanti dall'interpolazione verticale, ma asserisce che il metodo fornisce un approccio robusto e intuitivo per l'imaging strutturale.