In this article we begin The study of the problem of recovering the elastic parameters of a layered half-space from single component measurements of reflected waves, and applications to reflection seismology. We analyse the perturbational relationship between the elastic moduli (p, λ and µ and the vertical component of surface particle displacement due to a point impulsive traction. We show that this linearized problem has a unique solution and derive some stability estimates. The analysis of the linearized problem proceeds by way of the plane wave decomposition, and the uniqueness proof is constructive. To illustrate these ideas, we implement the construction numerically and recover parameter perturbations from data perturbations in a synthetic example. The solution of the linearized problem may be viewed as a step in a Newton-like scheme for solution of the full inverse problem. In considering the application of any such technique to field recorded data, however, one must confront problems such as model inadequacy, data error, source signature error, band limitation, etc. The extent to which these problems are resolvable has been explored for the simpler acoustic model using techniques compatible with those presented here. We plan to address such (practically crucial) stability considerations, and their implications for the design of useful algorithms, elsewhere. The focus of the present paper is thus the mathematical uniqueness question, which may be regarded as the most primitive stability question.
Read full abstract