Numerical Model of a Multiple-Input-Multiple-Output Distributed Acoustic Sensing System with Joint Phase and Birefringence Estimation

Sintesi Tecnica: Modello Numerico di un MIMO-DAS con Stima Congiunta di Fase e Birifrangenza

Definizione del Problema
La sensoristica distribuita in fibra ottica (DAS) sfrutta l'esistente infrastruttura delle telecomunicazioni per applicazioni di monitoraggio come la sorveglianza del traffico, il monitoraggio della salute strutturale e la rilevazione dei terremoti. Gli approcci convenzionali si concentrano tipicamente su uno dei due parametri seguenti:

1. OTDR di fase differenziale ($\Delta\phi$-OTDR): stima la fase ottica comune della luce retrodiffusa, che è altamente sensibile alla deformazione longitudinale, ma spesso trascura gli aspetti di polarizzazione, rendendola incapace di rilevare perturbazioni puramente trasversali o anisotrope.
2. Polarization-OTDR (P-OTDR): sfrutta le proprietà di polarizzazione per rilevare deformazioni trasversali (flessioni, torsioni, stress trasversale). Tuttavia, la P-OTDR tradizionale richiede spesso molteplici Stati di Polarizzazione (SOP) in ingresso, il che porta a tempi di misurazione più lunghi che limitano il monitoraggio dinamico. Inoltre, tipicamente non utilizza la fase comune, risultando cieca rispetto alle perturbazioni assialsimmetriche.

L'articolo affronta la necessità di un approccio di sensing dinamico multi-parametrico che possa stimare simultaneamente sia la fase ottica comune che le quantità legate alla polarizzazione (specificamente la birifrangenza) per discriminare tra diversi tipi di eventi di deformazione (longitudinale vs. trasversale/anisotropa) con alta sensibilità e risoluzione spaziale.

Metodologia
Gli autori hanno sviluppato e validato sperimentalmente un modello numerico per un sistema di Distributed Acoustic Sensing MIMO (MIMO-DAS). La metodologia comprende tre componenti principali:

1. Modello di Propagazione della Fibra e di Scattering:

   • La fibra è modellata come una serie di segmenti, ciascuno trattato come una lamina d'onda (waveplate) con birifrangenza lineare.
   • Condizioni Statiche: La birifrangenza intrinseca è generata utilizzando un modello in cui le componenti del vettore di birifrangenza di Stokes seguono processi di Langevin indipendenti, risultando in una distribuzione di Rayleigh della magnitudo della birifrangenza e in un orientamento casuale lungo la fibra.
   • Perturbazioni Dinamiche: Il modello incorpora l'effetto fotoelastico e l'allungamento della fibra. Distingue tra due tipi di eventi:
     • Deformazione Longitudinale: Modellata come allungamento della fibra ($\epsilon_z$) con contrazione di Poisson ($\epsilon_x = \epsilon_y = -\nu\epsilon_z$). Ciò altera la fase comune tramite l'effetto fotoelastico e la variazione di lunghezza, ma non induce birifrangenza aggiuntiva (poiché $\epsilon_x = \epsilon_y$).
     • Deformazione Trasversale/Anisotropa: Modellata con deformazioni trasversali disuguali ($\epsilon_x \neq \epsilon_y$), che induce birifrangenza aggiuntiva. Il modello calcola la somma vettoriale della birifrangenza intrinseca e di quella indotta dall'evento nello spazio di Stokes, notando che la magnitudo risultante è non lineare e dipende dall'orientamento relativo degli assi rapidi.

2. Architettura del Sistema MIMO-DAS:

   • Il sistema utilizza sequenze di Golay ortogonali mappate BPSK con divisione di polarizzazione (PDM) per sondare la fibra. Questo schema di codifica riduce la potenza ottica di picco, consentendo la coesistenza con la trasmissione dati WDM pur mantenendo il Rapporto Segnale-Rumore (SNR).
   • Il sistema impiega un ricevitore coerente a doppia polarizzazione per catturare il segnale retrodiffuso.
   • Elaborazione del Segnale: Il sistema stima le matrici di Jones round-trip per ogni segmento. Da queste matrici, estrae congiuntamente:
     • Fase Comune ($\phi$): Derivata dal determinante della matrice di Jones.
     • Birifrangenza Effettiva ($\Delta\beta_{eff}$): Derivata dagli autovalori del prodotto di matrici di Jones consecutive, che rappresenta il ritardo accumulato su una lunghezza di gauge.

3. Strategia di Validazione:

   • Simulazione Numerica: Il modello è stato testato su una fibra di 3 km con eventi di deformazione longitudinale e trasversale simulati.
   • Validazione Sperimentale: Esperimenti di laboratorio sono stati condotti utilizzando due bobine di fibra (235 m e 1956 m) e una configurazione che applicava deformazioni longitudinali controllate (tramite uno stadio piezoelettrico) e trasversali (tramite un compressore di fibra).

Contributi Chiave

• Framework di Stima Congiunta: L'articolo propone un metodo per stimare congiuntamente la fase ottica comune e la magnitudo della birifrangenza lineare effettiva dalle matrici di Jones distribuite in un sistema MIMO-DAS dinamico.
• Discriminazione degli Eventi: Lo studio dimostra che il sistema può discriminare tra i tipi di evento basandosi sulle loro firme distinte:
  • Deformazione Longitudinale: Influenza significativamente la fase comune ma lascia la magnitudo della birifrangenza sostanzialmente invariata.
  • Deformazione Trasversale/Anisotropa: Influenza sia la fase comune che la magnitudo della birifrangenza.
• Caratterizzazione della Birifrangenza Statica: Il modello stima con successo la distribuzione statica della birifrangenza effettiva, mostrando una distribuzione di Rayleigh coerente con le aspettative teoriche per le fibre standard monomodali.
• Verifica Sperimentale: Il modello numerico è stato validato rispetto agli esperimenti di laboratorio, confermando la capacità del sistema di rilevare, localizzare e discriminare tra deformazioni longitudinali e trasversali con una risoluzione spaziale media di circa 1,3 metri.

Risultati

• Condizioni Statiche: Il sistema ha stimato la magnitudo della birifrangenza effettiva lungo la fibra. La distribuzione dei valori stimati per una bobina di 1956 m si è adattata a una distribuzione di Rayleigh con un valore medio di $\approx 0,57$ rad/m, coerente con le previsioni teoriche. Il modello ha anche distinto con successo tra bobine concatenate con diversi valori medi di birifrangenza.
• Deformazione Longitudinale Dinamica: Simulazioni ed esperimenti hanno mostato che una deformazione longitudinale sinusoidale (800 Hz in simulazione, 100 Hz in esperimento) ha prodotto una chiara firma nella stima della fase differenziale, ma era non rilevabile nelle tracce di ritardo (birifrangenza).
• Deformazione Trasversale Dinamica: La deformazione trasversale anisotropa (simulata e applicata tramite un compressore di fibra) ha prodotto firme rilevabili sia nella fase differenziale che nelle tracce di ritardo.
• Risoluzione Spaziale: Utilizzando un fattore di bassa risoluzione (selezionando i segmenti con la riflessione più alta tra gruppi di 5), il sistema ha raggiunto una risoluzione spaziale media di 1,28 m (circa 1,3 m). Gli autori hanno osservato un compromesso tra sensibilità e risoluzione, dove fattori di risoluzione più elevati (es. 10) migliorano il contrasto ma riducono la risoluzione spaziale.

Significatività e Rivendicazioni
L'articolo sostiene che, sfruttando sia la stima della fase che quella della polarizzazione in un'architettura MIMO-DAS codificata, è possibile migliorare la discriminazione degli eventi e la sensibilità alle perturbazioni trasversali rispetto ai sistemi a parametro singolo allo stato dell'arte. Gli autori sottolineano che questo approccio consente il rilevamento sia di disturbi isotropi (longitudinali) che anisotropi (trasversali).

Tuttavia, l'articolo mantiene una posizione modesta riguardo alle capacità quantitative del sistema. Gli autori affermano esplicitamente che la birifrangenza effettiva stimata non varia linearmente con l'ampiezza della perturbazione; la sua risposta dipende in modo non lineare dall'orientamento relativo tra la birifrangenza intrinseca della fibra e la birifrangenza indotta dall'evento. Di conseguenza, sebbene il sistema possa discriminare la natura della deformazione (longitudinale vs. trasversale), non fornisce attualmente una quantificazione precisa della magnitudo della deformazione senza una conoscenza ulteriore dell'orientamento. Gli autori notano inoltre limitazioni pratiche, come il basso trasferimento delle deformazioni trasversali dal cavo alla fibra nelle installazioni reali, che potrebbe limitare il rilevamento di piccoli eventi trasversali. Il lavoro futuro suggerito si concentra sulla comprensione degli effetti del rumore e sulla conduzione di prove sul campo.