A method of obtaining a velocity-depth envelope from wide-angle seismic data

Summary. Due to the non-uniqueness of traveltime inversion of seismic data, it is more appropriate to determine a velocity-depth ( v-z ) envelope, rather than just a v-z function. Several methods of obtaining a v-z envelope by extremal inversion have been proposed, all of which invert the data primarily from either x-p , or T-p , or both domains. These extremal inversion methods may be divided into two groups: linear extremal and non-linear extremal. There is some debate whether the linearized perturbation techniques should be applied to the inherently non-linear problem of traveltime inversion. We have obtained a v-z envelope by extremal inversion in T-p with the constraint that the inversion paths also satisfy x-p observations. Thus we use data jointly in r-p and x-p , and yet avoid the linearity assumptions.
This joint, non-linear extremal inversion method has been applied to obtain a v-z envelope down to a depth of about 30 km in the Baltimore Canyon trough using x-t data from an Expanding Spread Profile acquired during the LASE project. We have found that the area enclosed by the v-z envelope is reduced by about 15 per cent using x-p control on the T-p inversion paths, compared to the inversion without x-p control.     

Structural evolution and the partitioning of deformation during basin growth and inversion: A case study from the Mizen Basin Celtic Sea,offshore Ireland

The Celtic Sea basins lie on the continental shelf between Ireland and northwest France and consist of a series of ENE–WSW trending elongate basins that extend from St George’s Channel Basin in the east to the Fastnet Basin in the west. The basins, which contain Triassic to Neogene stratigraphic sequences, evolved through a complex geological history that includes multiple Mesozoic rift stages and later Cenozoic inversion. The Mizen Basin represents the NW termination of the Celtic Sea basins and consists of two NE–SW-trending half-grabens developed as a result of the reactivation of Palaeozoic (Caledonian, Lower Carboniferous and Variscan) faults. The faults bounding the Mizen Basin were active as normal faults from Early Triassic to Late Cretaceous times. Most of the fault displacement took place during Berriasian to Hauterivian (Early Cretaceous) times, with a NW–SE direction of extension. A later phase of Aptian to Cenomanian (Early to Late Cretaceous) N–S-oriented extension gave rise to E–W-striking minor normal faults and reactivation of the pre-existing basin bounding faults that propagated upwards as left-stepping arrays of segmented normal faults. In common with most of the Celtic Sea basins, the Mizen Basin experienced a period of major erosion, attributed to tectonic uplift, during the Paleocene. Approximately N–S Alpine regional compression-causing basin inversion is dated as Middle Eocene to Miocene by a well-preserved syn-inversion stratigraphy. Reverse reactivation of the basin bounding faults was broadly synchronous with the formation of a set of near-orthogonal NW–SE dextral strike-slip faults so that compression was partitioned onto two fault sets, the geometrical configuration of which is partly inherited from Palaeozoic basement structure. The segmented character of the fault forming the southern boundary of the Mizen Basin was preserved during Alpine inversion so that Cenozoic reverse displacement distribution on syn-inversion horizons mirrors the earlier extensional displacements. Segmentation of normal faults therefore controls the geometry and location of inversion structures, including inversion anticlines and the back rotation of earlier relay ramps.     

乌鲁木齐市低空大气逆温特征分析

为了分析乌鲁木齐低层大气逆温变化特征,通过利用乌鲁木齐市气象站2000-2004年的探空观测资料,从逐日探空原始记录中提取得到贴地逆温、脱地逆温的底部、顶部高度、温度,分析它们的变化情况。乌鲁木齐市低空大气温度层结全年以弱稳定为主,一年四季都存在逆温层。冬季是逆温特征最为显著的季节,强度最强,平均为1.06℃/100 m。冬季逆温发生频率最高,平均为89%,贴地逆温厚度最厚,平均为860 m,脱地逆温底高平均为534 m,顶高平均为1 187 m。由于蒙古反气旋的冷空气在准噶尔盆地堆积,使得乌鲁木齐冬季近地面大气稳定。     

A practical solution for linear inference computations

Summary. In this paper, we present a matrix form of Backus' theory of linear inference with multiple predictions. The Bayesian approach used by Backus allows the treatment of erroneous data and the imposition of the essential a priori bound on the model norm. We introduce a modification which involves translating the problem in accordance with an estimated model. Such a model may be known a priori or it may be constructed from the data. We are effectively able to bound the norm of all acceptable models from above and below and this results in more confining estimates of the predictions than provided by just an upper bound. In addition, the model construction approach allows us to make maximum use of the data in the inference computation. Our algorithm is robust and efficient, and estimates comparable to to those obtained from linear programming techniques have been achieved.     

Three-dimensional massively parallel electromagnetic inversion—I. Theory

An iterative solution to the non-linear 3-D electromagnetic inverse problem is obtained by successive linearized model updates using the method of conjugate gradients. Full wave equation modelling for controlled sources is employed to compute model sensitivities and predicted data in the frequency domain with an efficient 3-D finite-difference algorithm. Necessity dictates that the inverse be underdetermined, since realistic reconstructions require the solution for tens of thousands of parameters. In addition, large-scale 3-D forward modelling is required and this can easily involve the solution of over several million electric field unknowns per solve. A massively parallel computing platform has therefore been utilized to obtain reasonable execution times, and results are given for the 1840-node Intel Paragon. The solution is demonstrated with a synthetic example with added Gaussian noise, where the data were produced from an integral equation forward-modelling code, and is different from the finite difference code embedded in the inversion algorithm     

Inversion of sonic waveforms with unknown source and receiver functions

A seismogram that is several times the length of the source-receiver wavelet is windowed into two parts—these may overlap—to obtain two seismograms with approximately the same source function but different Green's functions. A similarly windowed synthetic seismogram gives two corresponding synthetic seismograms. The spectral product of the window 1 data with the window 2 synthetic is equal to the spectral product of the window 1 synthetic with the window 2 data only if the correct earth model is used to compute the synthetic. This partition principle is applied to well-log sonic waveform data from Ocean Drilling Project hole 806B, a slow formation, and used there to estimate Poisson's ratio from a single seismogram whose transmitter and receiver functions are unknown. A multichannel extension of the algorithm gives even better results. The effective borehole radius R _b, was included in the inversion procedure, because of waveform sensitivity to R _b. Inversion results for R _b agreed with the sonic caliper, but not the mechanical caliper; thus if R _b is not included in the inversion its value should be taken from the sonic caliper.     

Lithotype discrimination within a gridded model by inverting gravity data

In this paper we present a method which allows delineation of geologic structures in a bi-modal lithotype setting. We propose to use gravity data in combination with a priori information about the density contrast between the two lithotypes. The iterative method uses an objective function with five tunable parameters which need to be set. Using an efficient parameter search, suitable ranges of these are investigated to determine their optimal values, respectively, which in turn, ensures good inversion results.
The approach produces structural images of the subsurface, without the need of an a priori density model; the depth to the top of the inhomogeneity is also retrieved.
Besides synthetic simulations, the methodology has also been applied to a small gravity data set, acquired by the industry over a basinal structure. A consistent, bi-modal image of the bedrock depression is obtained from the data, which, in this case, was the goal. Other potential areas of application include delineation of salt structures and ore deposits.