Modelling strategies and model assessment for wide-angle seismic traveltime data

Strategies for modelling seismic refraction/wide-angle reflection traveltim es to obtain 2-D velocity and interface structure are presented along with methods for assessing the reliability of the results. Emphasis is placed on using inverse methods, but a discussion of arrival picking and classificat ion, data uncertainties, traveltime reciprocity, crooked line geometry and the selection of a starting model is also applicable to trial-and-error for ward modelling. The most important advantages of an inverse method are the ability to derive simpler models for the appropriate level of fit to the da ta, and the ability to assess the final model in terms of resolution, param eter bounds and non-uniqueness. Given the unique characteristics of each da ta set and the local earth structure, there is no single approach to modell ing wide-angle data that is best. This paper describes the best modelling s trategies according to (1) the model parametrization, (2) the inclusion of prior information, (3) the complexity of the earth structure, (4) the chara cteristics of the data, and (5) the utilization of coincident seismic refle ction data. There are two natural end-member inversion styles: (1) a regula r, fine-grid parametrization when seeking a minimum-structure model, and (2 ) an irregular grid, minimum-parameter model when considering certain forms of prior information. The former style represents the 'pure' tomography ap proach. The latter style is closer to automated forward modelling, and can be applied best with a parameter-selective algorithm, that is, one that all ows any subset of model parameters to be selected for inversion. If there i s strong lateral heterogeneity in the near-surface only, layer stripping wo rks well. If there is complexity at all depths, all model parameters should be determined simultaneously after careful construction of a starting mode l that allows the appropriate rays to be traced to all pick locations. The lateral spacing of model nodes to use will depend on the type of inversion and whether detailed prior information is included, but a general guideline based on model resolution when seeking a minimum-parameter model is a mode l node spacing equal to the shot spacing (receiver spacing for typical mari ne data), except perhaps in the upper layers where about half this may be n ecessary; node spacing is not an issue when using smoothing constraints, pr ovided it is small enough to resolve the earth structure of interest. Trave ltimes picked from pre-stack, unmigrated or migrated coincident reflection data can be (1) used to develop the starting model, (2) inverted simultaneo usly with the wide-angle data, or (3) inverted after modelling the wide-ang le data to constrain interfaces that 'float' within the velocity model. Mod el assessment establishes the reliability of the final model. Presenting mo del statistics, traveltime fits, ray diagrams and resolution kernels is use ful, but can only indirectly address this issue. Direct model assessment te chniques that derive alternative models that satisfactorily fit the real da ta are the best means of establishing the absolute bounds on model paramete rs and whether a particular model feature is required by the data.