Fisica

Blog

CasaCasa / Blog / Fisica

Nov 29, 2023

Fisica

Scientific Reports volume 13, numero articolo: 13133 (2023) Cita questo articolo 6798 Accessi a 3 dettagli delle metriche alternative Previsione a breve termine della magnitudo massima stimata (\({\widehat{M}}_{max}\))

Rapporti scientifici volume 13, numero articolo: 13133 (2023) Citare questo articolo

6798 accessi

3 Altmetrico

Dettagli sulle metriche

La previsione a breve termine della magnitudo massima stimata (\({\widehat{M}}_{max}\)) è fondamentale per mitigare i rischi di sismicità indotta durante la stimolazione dei fluidi. La maggior parte dei metodi precedenti richiede dati di iniezione in tempo reale, che non sono sempre disponibili. Questo studio propone due approcci di deep learning (DL), insieme a due metodi di partizionamento dei dati, che si basano esclusivamente su precedenti modelli di sismicità. Il primo approccio prevede \({\widehat{M}}_{max}\) utilizzando direttamente DL; il secondo incorpora vincoli fisici utilizzando DL per prevedere il tasso di sismicità, che viene poi utilizzato per stimare \({\widehat{M}}_{max}\). Questi approcci vengono testati utilizzando un set di dati di monitoraggio delle fratture idrauliche del Canada occidentale. Troviamo che il DL diretto apprende dai precedenti modelli di sismicità per fornire una previsione accurata, anche se con un ritardo che ne limita l’utilità pratica. L'approccio basato sulla fisica prevede accuratamente i cambiamenti nel tasso di sismicità, ma a volte sottostima (o sovrastima) \({\widehat{M}}_{max}\). Proponiamo che un superamento significativo di \({\widehat{M}}_{max}\) possa annunciare l'inizio di una rottura fuori controllo della faglia.

La fratturazione idraulica (HF), un metodo di stimolazione dei fluidi per migliorare la permeabilità producendo fratture nelle rocce serbatoio a bassa permeabilità1, produce tipicamente microterremoti (MEQ) con magnitudo del momento MW < 0. Tuttavia, l'HF può anche indurre terremoti moderati (MW > 4)2 ,3,4,5,6, che sono associati all'attivazione di guasti preesistenti7. Ottenere una stima probabilistica della massima magnitudo attesa dell'evento (\({\widehat{M}}_{max}\)) per una data operazione HF è importante per la valutazione del pericolo8 e potrebbe informare strategie proattive di mitigazione in tempo reale per la sismicità indotta che sono richiesti in alcuni sistemi di monitoraggio avanzati9,10.

Sono stati sviluppati vari approcci per stimare la \({\widehat{M}}_{max}\) della sismicità indotta dai fluidi. Ad esempio, la distribuzione attesa delle magnitudo dei terremoti può essere espressa in termini di volume netto di fluido iniettato (∆V) e di indice sismogenico (∑), un parametro sismotettonico proposto specifico per l'area che caratterizza il livello di attività sismica atteso in risposta al fluido. iniezione11. Questa espressione è stata utilizzata per sviluppare una stima probabilistica per la magnitudo massima12, che scala linearmente con log10 ∆V. La stessa relazione di scala volumetrica è stata derivata utilizzando un approccio teorico diverso, basato sul criterio di equilibrio della fessura di Griffith13. In questo caso, la stima della magnitudo massima si applica al caso di rottura arrestata, un concetto in cui la zona di rottura della faglia è confinata in una regione sotterranea in cui la pressione è perturbata dall'iniezione di fluido. Questo concetto è stato utilizzato anche per sviluppare un vincolo geometrico per la magnitudo massima, basato sulla distribuzione spaziale dei MEQ14. In un'altra formulazione, il momento sismico massimo atteso per un terremoto indotto dall'iniezione è espresso come il prodotto del modulo di taglio del mezzo e del volume netto del fluido iniettato15. Con l'eccezione dell'approccio geometricamente vincolato14, che richiede la determinazione delle posizioni dell'ipocentro MEQ, tutti questi metodi utilizzano il volume netto iniettato ∆V come parametro per la stima di \({\widehat{M}}_{max}\).

Durante le operazioni HF, le osservazioni sismiche possono essere utilizzate per identificare MEQ operativi1 nonché eventi sismici indotti che si verificano su faglie vicine16,17,18,19. I MEQ operativi si verificano tipicamente in cluster che si estendono lontano dal pozzo, solitamente perpendicolarmente alla direzione dello stress orizzontale minimo20,21. In alcuni casi, una faglia riattivata è caratterizzata da un evento ritardato rispetto al momento di inizio di una fase di iniezione, accoppiato con un orientamento obliquo dei cluster di sismicità rispetto alle principali direzioni di stress16,17,18,19. La riattivazione della faglia è spesso contrassegnata da un aumento del tasso di sismicità, accompagnato da una diminuzione del valore b di Gutenberg-Richter22,23. Sebbene tali cambiamenti nel modello spaziotemporale della sismicità possano essere impercettibili, il loro rilevamento utilizzando metodi di deep learning (DL) potrebbe fornire una strada per una migliore previsione a breve termine.

 1000 MEQs)30. The fixed-window method with a floating estimate of b appears to track temporal fluctuations for small seismic magnitudes (MW < 2), but it fails to forecast larger events (Fig. 4a). This can be ameliorated by fixing b to unity, which leads to a forecast that approximates the upper limit for most seismic events but still fails to provide a forecast envelope for the largest observed events. For the cumulative data partition method, \({\widehat{M}}_{max}\) increases monotonically over time (Fig. 4c), as expected. In all cases the forecast has a low R2 value (Fig. 4b,d), indicating that for this approach the calculated value is not suitable for a direct forecast, although it could provide a forecast of the envelope of \({\widehat{M}}_{max}\)./p> 10,000 events, with a maximum magnitude of 3.1 MW. Overall, the observed seismicity is characterized by b >  > 1, as expected for operational MEQs1; however, individual event clusters associated with fault activation show a marked drop in b-value45. Based on the b-value stability method28,29 and the maximum likelihood method27, we determine Mc using the first 1000 MEQs in the catalogue finding Mc = – 0.15 (Fig. S1). Since the sensors used in the study are fixed and the event depths remain approximately the same throughout the HF program46, we assume that Mc is fixed at this value (– 0.15) for the duration of the experiment./p>