Skip navigation
Skip navigation

Geophysical inversion with a neighbourhood algorithm - I. Searching a parameter space

Sambridge, Malcolm

Description

This paper presents a new derivative-free search method for finding models of acceptable data fit in a multidimensional parameter space. It falls into the same class of method as simulated annealing and genetic algorithms, which are commonly used for global optimization problems. The objective here is to find an ensemble of models that preferentially sample the good data-fitting regions of parameter space, rather than seeking a single optimal model. (A related paper deals with the quantitative...[Show more]

dc.contributor.authorSambridge, Malcolm
dc.date.accessioned2015-12-13T23:35:29Z
dc.identifier.issn0956-540X
dc.identifier.urihttp://hdl.handle.net/1885/93927
dc.description.abstractThis paper presents a new derivative-free search method for finding models of acceptable data fit in a multidimensional parameter space. It falls into the same class of method as simulated annealing and genetic algorithms, which are commonly used for global optimization problems. The objective here is to find an ensemble of models that preferentially sample the good data-fitting regions of parameter space, rather than seeking a single optimal model. (A related paper deals with the quantitative appraisal of the ensemble.) The new search algorithm makes use of the geometrical constructs known as Voronoi cells to derive the search in parameter space. These are nearest neighbour regions defined under a suitable distance norm. The algorithm is conceptually simple, requires just two 'tuning parameters', and makes use of only the rank of a data fit criterion rather than the numerical value. In this way all difficulties associated with the scaling of a data misfit function are avoided, and any combination of data fit criteria can be used. It is also shown how Voronoi cells can be used to enhance any existing direct search algorithm, by intermittently replacing the forward modelling calculations with nearest neighbour calculations. The new direct search algorithm is illustrated with an application to a synthetic problem involving the inversion of receiver functions for crustal seismic structure. This is known to be a non-linear problem, where linearized inversion techniques suffer from a strong dependence on the starting solution. It is shown that the new algorithm produces a sophisticated type of 'self-adaptive' search behaviour, which to our knowledge has not been demonstrated in any previous technique of this kind.
dc.publisherBlackwell Publishing Ltd
dc.sourceGeophysical Journal International
dc.subjectKeywords: geophysics; inverse analysis; mathematical method Numerical techniques; Receiver functions; Waveform inversion
dc.titleGeophysical inversion with a neighbourhood algorithm - I. Searching a parameter space
dc.typeJournal article
local.description.notesImported from ARIES
local.description.refereedYes
local.identifier.citationvolume138
dc.date.issued1999
local.identifier.absfor040312 - Structural Geology
local.identifier.ariespublicationMigratedxPub25365
local.type.statusPublished Version
local.contributor.affiliationSambridge, Malcolm, College of Physical and Mathematical Sciences, ANU
local.description.embargo2037-12-31
local.bibliographicCitation.issue2
local.bibliographicCitation.startpage479
local.bibliographicCitation.lastpage494
local.identifier.doi10.1046/j.1365-246X.1999.00876.x
dc.date.updated2015-12-12T09:40:06Z
local.identifier.scopusID2-s2.0-0032753561
CollectionsANU Research Publications

Download

File Description SizeFormat Image
01_Sambridge_Geophysical_inversion_with_a_1999.pdf1.29 MBAdobe PDF    Request a copy


Items in Open Research are protected by copyright, with all rights reserved, unless otherwise indicated.

Updated:  22 January 2019/ Responsible Officer:  University Librarian/ Page Contact:  Library Systems & Web Coordinator