Adaptive regridding in 3-D reflection tomography

3D reflection tomography allows the macro-model of complex geological structures to be reconstructed. In the usual approach, the spatial distribution of the velocity field is discretized by regular grids. This choice simplifies the development of the related software, but introduces two serious drawbacks: various domains of the model may be poorly covered, and a relevant mismatch between the grid and a complex velocity field may occur. So the tomographic inversion becomes unstable, unreliable and necessarily blurred. In this paper, we introduce an algorithm to adapt the grid to the available ray paths and to the velocity field in sequence: so we get irregular grids with a locally variable resolution. We can guide the grid fitting procedure interactively, if we are going to introduce some geological a priori information; otherwise, we define a fully automatic approach, which exploits the Delauny triangles and Voronoi polygons.

La tomografia a riflessione in 3D permette di ricostruire un macro modello di strutture geologiche complesse. Nell’approccio comune, la distribuzione spaziale del campo di velocità è discretizzata da grigie regolari. Questa scelta semplifica lo sviluppo del software relativo, ma introduce due inconvenienti seri: diversi domini del modello possono essere scarsamente coperti, e può esserci una notevole discrepanza tra la griglia ed un campo di velocità complesso. Così l’inversione tomografica diviene instabile, non attendibile, e necessariamente poco definita. In quest’articolo presentiamo un algoritmo per adattare la griglia in sequenza ai percorsi dei raggi a disposizione ed al campo di velocità.: così otteniamo griglie irregolari con una risoluzione localmente variabile. Possiamo guidare l’adattamento della griglia interattivamente, se vogliamo introdurre informazioni geologiche note a priori, o al contrario definiamo un approccio totalmente automatico , basato sui triangoli di Delauny ed i poligoni di Voronoi.