Propagation of singularities and inverse problems for the viscoacoustic wave equation

1. Il Problema Studiato

L'articolo affronta un problema inverso per l'equazione delle onde viscoacustiche, un modello integro-differenziale che descrive la propagazione delle onde in mezzi viscoacustici con memoria.

• Modello Fisico: Il mezzo è caratterizzato da una velocità del suono spazialmente variabile $c(x)$ e da un nucleo di memoria $G(x, t)$ dipendente sia dallo spazio che dal tempo. L'equazione governante è:
  $$Pu := -\frac{1}{2} \left( \partial_t^2 u - \nabla \cdot (c(x)\nabla u) + \int_0^t \nabla \cdot (G(x, t-s)\nabla u(x, s)) \, ds \right) = 0$$
  Il termine integrale rappresenta il comportamento viscoso del mezzo, dove lo stress dipende dall'intera storia del gradiente della pressione.
• Configurazione Geometrica: Le onde sono generate da sorgenti $f$ supportate all'esterno della regione di interesse $\Omega$. Le misurazioni sono effettuate all'esterno e codificate tramite la mappa sorgente-soluzione ($S: f \mapsto u|_{M' \setminus M}$).
• Obiettivo: Determinare se la mappa sorgente-soluzione esterna determina univocamente la velocità del suono $c(x)$ all'interno di $\Omega$ e tutte le derivate temporali del nucleo di memoria $G(x, t)$ valutati a $t=0$ (cioè $\partial_t^j G(x, 0)$).

2. Metodologia

Gli autori adottano un approccio microlocale avanzato, superando le limitazioni dei metodi classici basati su sorgenti delta (che falliscono in presenza di caustiche).

• Costruzione di Soluzioni Approssimate (Quasimodi):
  • Vengono costruite soluzioni concentrate vicino a geodetiche fisse, note come quasimodi a fascio gaussiano (Gaussian beam quasimodes).
  • Queste soluzioni sono espresse come sovrapposizioni di pacchetti d'onda con fase complessa e ampiezze asintotiche.
  • A differenza del caso puramente elastico, la presenza del termine di memoria richiede l'introduzione di termini di ampiezza aggiuntivi ($a''$) ottenuti risolvendo equazioni di Volterra.
• Propagazione delle Singolarità:
  • Viene stabilito un teorema di propagazione delle singolarità per l'insieme frontiera d'onda semiclassico ($WF_{scl}$).
  • Un risultato cruciale è che il termine di memoria, pur essendo un operatore integro-differenziale, non altera la traiettoria delle singolarità di tipo ondulatorio (che seguono le geodetiche della metrica $c(x)^{-2}g_0$), influenzando solo le ampiezze.
• Riduzione a Problemi Inversi Geometrici:
  • I dati esterni permettono di recuperare la relazione di scattering e i tempi di viaggio (dati di lente o lens data).
  • Sotto opportune condizioni geometriche (convessità stretta, varietà semplici o foliate da ipersuperfici convesse), i dati di lente determinano univocamente la metrica (e quindi $c(x)$) tramite risultati di rigidità delle lenti e trasformate di Ray X geodetiche.
  • Le ampiezze delle soluzioni, una volta noto $c(x)$, permettono di recuperare le derivate temporali di $G$ tramite la trasformata di Ray X geodetica applicata a quantità specifiche.

3. Risultati Chiave

• Teorema 1.1 (Costruzione di Quasimodi): Esistenza di soluzioni approssimate della forma $u(x,t) = \sum (-ik)^{-l} \int e^{ik\phi} a_l d\sigma$ che soddisfano l'equazione con errore $O(k^{-N})$. Queste soluzioni concentrano l'energia lungo le bicharatteristiche nulle.
• Teorema 1.2 (Propagazione delle Singolarità): Le singolarità delle soluzioni dell'equazione viscoacustica si propagano lungo le stesse geodetiche del caso elastico (onda pura). Il termine di memoria non modifica la traiettoria delle singolarità, ma solo l'ampiezza.
• Teorema 1.3 (Unicità del Problema Inverso):
  • Se le mappe sorgente-soluzione coincidono per due coppie di parametri $(c_1, G_1)$ e $(c_2, G_2)$, allora:
    1. $c_1 = c_2$ in tutto il dominio.
    2. $\partial_t^j G_1(x, 0) = \partial_t^j G_2(x, 0)$ per ogni ordine $j \in \mathbb{N}$.
  • Questo risultato vale senza restrizioni sulla semplicità della metrica (a differenza di lavori precedenti), purché siano soddisfatte le condizioni geometriche (A1)-(A3) che evitano la presenza di geodetiche chiuse o caustiche patologiche.
• Teorema 1.4 (Applicazione al Modello di Maxwell Esteso):
  • Viene applicato il risultato generale al Modello di Maxwell Esteso (Extended Maxwell Model), dove il nucleo di memoria è una somma di esponenziali decrescenti.
  • Si dimostra che la conoscenza delle derivate temporali di $G$ fino a un certo ordine (dipendente dal numero di termini $N$ del modello) è sufficiente a recuperare univocamente tutti i parametri fisici (costanti elastiche e viscose) del modello.

4. Contributi Principali

1. Generalità Geometrica: Il lavoro estende i risultati di unicità a geometrie non semplici (dove possono esistere caustiche), un caso tipico nelle applicazioni sismiche reali, utilizzando un approccio basato su fasci gaussiani invece di sorgenti puntuali.
2. Trattamento del Termine di Memoria: Dimostrazione rigorosa che il termine di memoria di ordine superiore non influenza la propagazione delle singolarità, permettendo di separare il recupero della velocità del suono da quello dei parametri di memoria.
3. Metodologia Asintotica: Integrazione di tecniche di analisi semiclassica, equazioni di Riccati matriciali e equazioni di Volterra per costruire soluzioni adatte all'inversione in presenza di memoria.
4. Applicazione Pratica: Fornisce una base teorica solida per l'inversione di parametri in modelli di attenuazione sismica complessi (come il modello di Maxwell esteso), superando le limitazioni delle approssimazioni puramente elastiche.

5. Significato e Impatto

Questo studio è significativo per la sismologia e l'imaging geofisico. Nella pratica, i materiali terrestri (rocce, fluidi nei pori) mostrano un comportamento viscoelastico con memoria, che causa attenuazione e dispersione delle onde.

• I metodi tradizionali di inversione spesso trascurano la memoria o assumono geometrie semplici, portando a errori nella stima della velocità del suono e delle proprietà del sottosuolo.
• Questo lavoro dimostra che, teoricamente, è possibile recuperare sia la struttura elastica ($c$) che i parametri di attenuazione viscosa ($G$) dai dati esterni, anche in scenari complessi con focalizzazione delle onde (caustiche).
• L'applicazione al modello di Maxwell esteso offre un percorso verso la ricostruzione univoca dei parametri fisici microscopici (costanti di molla e smorzatori) che governano il comportamento macroscopico dei materiali geologici.

In sintesi, l'articolo fornisce un quadro teorico completo e rigoroso per risolvere problemi inversi in mezzi viscoacustici complessi, combinando analisi microlocale, geometria differenziale e teoria delle equazioni integrali.