a partition modelling approach to tomographic problems thomas bodin & malcolm sambridge research...
TRANSCRIPT
![Page 1: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/1.jpg)
A Partition Modelling Approach to
Tomographic Problems
Thomas Bodin & Malcolm Sambridge
Research School of Earth Sciences,
Australian National University
![Page 2: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/2.jpg)
Outline
Parameterization in Seismic tomography
Non-linear inversion, Bayesian Inference and Partition Modelling
An original way to solve the tomographic problem
• Method
• Synthetic experiments
• Real data
![Page 3: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/3.jpg)
2D Seismic Tomography
We want
A map of surface wave velocity
![Page 4: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/4.jpg)
2D Seismic Tomography
source
receiver
time
distv
We want
A map of surface wave velocity
We have
Average velocity along seismic rays
![Page 5: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/5.jpg)
We want
A map of surface wave velocity
We have
Average velocity along seismic rays
2D Seismic Tomography
![Page 6: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/6.jpg)
2D Seismic Tomography
We want
A map of surface wave velocity
We have
Average velocity along seismic rays
![Page 7: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/7.jpg)
Regular ParameterizationCoarse grid Fine grid
Bad GoodResolution
Constrain on the model
Good Bad
![Page 8: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/8.jpg)
Regular ParameterizationCoarse grid Fine grid
Bad GoodResolution
Constraint on the model
Good Bad
Define arbitrarily more constraints on the model
![Page 9: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/9.jpg)
9
Irregular parameterizations
Chou & Booker (1979); Tarantola & Nercessian (1984); Abers & Rocker (1991); Fukao et al. (1992); Zelt & Smith (1992); Michelini
(1995); Vesnaver (1996); Curtis & Snieder (1997); Widiyantoro & van der Hilst (1998); Bijwaard et al. (1998); Bohm et al. (2000);
Sambridge & Faletic (2003).
Nolet & Montelli (2005)
Sambridge & Rawlinson (2005)Gudmundsson & Sambridge (1998)
![Page 10: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/10.jpg)
Voronoi cells
Cells are only defined by their
centres
![Page 11: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/11.jpg)
11
QuickTime™ and a decompressor
are needed to see this picture.
Voronoi cells are everywhere
![Page 12: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/12.jpg)
12
QuickTime™ and a decompressor
are needed to see this picture.
Voronoi cells are everywhere
![Page 13: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/13.jpg)
13
QuickTime™ and a decompressor
are needed to see this picture.
QuickTime™ and a decompressor
are needed to see this picture.
Voronoi cells are everywhere
![Page 14: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/14.jpg)
Voronoi cells
Problem becomes highly nonlinear
Model is defined by:
* Velocity in each cell* Position of each cell
![Page 15: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/15.jpg)
Non Linear Inversion
X2
X1
Sampling a multi-dimensional function
X1X2
![Page 16: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/16.jpg)
Non Linear Inversion
Optimisation Bayesian Inference
Solution : Maximum Solution : statistical distribution
X2
X1
X2
X1
X2
X1
X2
X1
(e.g. Genetic Algorithms, Simulated Annealing)
(e.g. Markov chains)
![Page 17: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/17.jpg)
Partition Modelling(C.C. Holmes. D.G.T. Denison, 2002)
• Cos ? • Polynomial function?
Regression Problem
A Bayesian technique used for classification and Regression problems in Statistics
![Page 18: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/18.jpg)
n=3
The number of parameters is
variable
Dynamic irregular parameterisation
Partition Modelling
![Page 19: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/19.jpg)
n=6 n=11
n=8 n=3
Partition Modelling
![Page 20: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/20.jpg)
Mean. Takes in account all the
models
Partition Modelling
Bayesian Inference
Mean solution
![Page 21: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/21.jpg)
Adaptive parameterisation
Automatic smoothing
Able to pick up discontinuities
Partition Modelling
Can we apply these concepts to tomography ?
Mean solution
True solution
![Page 22: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/22.jpg)
Synthetic experiment
True velocity model Ray geometry
Data Noise σ = 28 sKm/s
![Page 23: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/23.jpg)
Iterative linearised tomography
Inversion step Subspace method (Matrix inversion)
Fixed Parameterisation
Regularisation procedure
Interpolation
Inversion step Subspace method (Matrix inversion)
Fixed Parameterisation
Regularisation procedure
Interpolation
Ray geometry
Ray geometry
Observed travel timesObserved
travel times
Forward calculationFast Marching Method
Forward calculationFast Marching Method
Solution Model
Solution Model
Reference Model
Reference Model
![Page 24: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/24.jpg)
Regular grid Tomographyfixed grid (20*20 nodes)
Damping
Smoothing
Km/s
20 x 20 B-splines nodes
![Page 25: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/25.jpg)
Iterative linearised tomography
Inversion step Subspace method (Matrix inversion)
Fixed Parameterisation
Regularisation procedure
Interpolation
Inversion step Subspace method (Matrix inversion)
Fixed Parameterisation
Regularisation procedure
Interpolation
Ray geometry
Ray geometry
Observed travel timesObserved
travel times
Forward calculationFast Marching Method
Forward calculationFast Marching Method
Solution Model
Solution Model
Reference Model
Reference Model
![Page 26: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/26.jpg)
Iterative linearised tomography
Inversion step
Partition Modelling
Adaptive Parameterisation No regularisation procedure No interpolation
Inversion step
Partition Modelling
Adaptive Parameterisation No regularisation procedure No interpolation
Ray geometry
Ray geometry
Observed travel timesObserved
travel times
Forward calculationFast Marching Method
Forward calculationFast Marching Method
Ensemble of ModelsEnsemble of Models
Reference Model
Reference Model
Point wise spatial average
![Page 27: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/27.jpg)
Description of the method
I. Pick randomly one cell
II. Change either its value or its position
III. Compute the estimated travel time
IV. Compare this proposed model to the current one
)(
)(,1min)(
current
proposed
mP
mPacceptP
Each stepKm/s
![Page 28: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/28.jpg)
Description of the method
Step 150 Step 300 Step 1000
![Page 29: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/29.jpg)
Solution
Maxima Mean
Best model sampled Average of all the models sampled
Km/s
![Page 30: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/30.jpg)
Regular Grid vs Partition Modelling
200 fixed cells 45 mobile cells
Km/sKm/s
![Page 31: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/31.jpg)
Model Uncertainty
Standard deviation
1
0
Average Cross Section
True modelAvg. model
![Page 32: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/32.jpg)
Computational Cost Issues
Monte Carlo Method cannot deal with high dimensional problems, but …
Resolution is good with small number of cells.
Possibility to parallelise.
No need to solve the whole forward problem at each iteration.
![Page 33: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/33.jpg)
Computational Cost Issues
When we change the value of one cell …
![Page 34: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/34.jpg)
Computational Cost Issues
When we change the position of one cell …
![Page 35: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/35.jpg)
Computational Cost Issues
When we change the position of one cell …
![Page 36: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/36.jpg)
When we change the position of one cell …
Computational Cost Issues
![Page 37: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/37.jpg)
When we change the position of one cell …
Computational Cost Issues
![Page 38: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/38.jpg)
Real Data
(Erdinc Saygin ,2007)
Cross correlation of seismic
ambient noise
![Page 39: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/39.jpg)
Real Data
Maps of Rayleigh
waves group velocity at
5s.
Damping
Smoothing
Km/s
![Page 40: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/40.jpg)
40
Changing the number of Voronoi cells
The birth step
Generate randomly the location of a new cell nucleus
![Page 41: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/41.jpg)
Real Data
Variable number of Voronoi cells
Average model (Km/s) Error estimation (Km/s)
![Page 42: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/42.jpg)
Real Data
Variable number of Voronoi cells
Average model (Km/s)
![Page 43: A Partition Modelling Approach to Tomographic Problems Thomas Bodin & Malcolm Sambridge Research School of Earth Sciences, Australian National University](https://reader038.vdocument.in/reader038/viewer/2022103112/551af59b550346f70d8b4ff0/html5/thumbnails/43.jpg)
Conclusion
Adaptive Parameterization
Automatic smoothing and regularization
Good estimation of model uncertainty