ON THE BREMMER SERIES DECOMPOSITION: EQUIVALENCE BETWEEN TWO DIFFERENT APPROACHES*

A Bremmer Series decomposition of the solution y(t) to the lossless wave equation in layered media is where the y_j(t) are physically meaningful constituents (i.e., y₁(t) are primaries, y₂(t) are secondaries, etc.). This paper reviews Mendel's state space models for generating the constituents; reviews Bremmer's integral equation models for generating the constituents; and demonstrates how Mendel's state space models can be obtained by a careful decomposition of Bremmer's integral equation models. It shows that Mendel's equations can be viewed as approximate numerical solutions of Bremmer's integral equations. In a lossless homogeneous medium, the approximations become exact.     

APPLICATION OF ELASTIC MODELLING IN PROCESSING AND INTERPRETATION OF VSP DATA: A CASE HISTORY1

Data from offshore Norway is used to study applications of elastic VSP modelling in detecting shear waves and observing the effects of successive mode conversion in field-recorded VSP data. The shear-wave velocities and densities from log data are used in conjunction with compressional wave velocities determined from surface seismic and log data in the VSP modelling. The time domain non-normal incidence elastic VSP modelling technique of Aminzadeh and Mendel is used as the modelling algorithm. Two surface seismograms are computed first. One is the vertical component and the other is the horizontal component for plane waves that have specified incident angles. A downward continuation method is then applied to generate seismograms at different depth points. The collection of these seismograms constitutes non-normal incidence VSPs. Both vertical and horizontal components of VSP data can be obtained by this procedure. In this paper non-normal incidence VSPs are generated for a 12.5° incident plane wave. The modelling results of layered earth systems of thin layers and thick layers are both compared with field data, and the effect of mode conversions in thin layers is observed. Several events in the field data can be explained by this elastic VSP modelling. Comparison of the model data and field data enabled a probable tube wave or out-of-plane event to be identified, the removal of which significantly improved the final VSP section. This study also shows how the VSP data helped the interpretation of the surface 3D data.     

DERIVATION OF LAYER PARAMETERS OF AN ELASTIC MEDIUM FROM REFLECTION COEFFICIENT MATRICES*

In this paper we develop a recursive algorithm to obtain the layer parameters of an elastic medium (density, P-wave velocity, S-wave velocity) from reflection coefficient matrices in terms of energy flux ratios for a non-normal incidence case. We define a layer impedance matrix, analogous to the impedance of an acoustic medium. Next we derive a matrix relationship between the layer impedance matrix of the n+ 1st layer and the reflection coefficient and parameter matrices of the nth layer. This relationship leads to recursively computing the parameters of the subsurface. We show that the elastic case—unlike the acoustic case—allows one to recover the layer parameters from the impedance matrix for non-normal incidence. The results of this work play a key role in the solution of the inverse problem with non-normal-incidence plane-wave seismic data when using a downward continuation technique.     

NON-NORMAL INCIDENCE STATE-SPACE MODEL AND LINE SOURCE REFLECTION SYNTHETIC SEISMOGRAM*

This paper is directed at modeling layered media. We extend the plane-wave normal-incidence state-space model developed by Mendel, Nahi and Chan in 1979, to the non-normal incidence case. To do this we introduce a shifting principle, a zero-offset wavefront, and zero-offset travel times for different layers. We also develop an algorithm for obtaining a synthetic line source reflection seismogram. In this algorithm non-normal incidence plane-wave seismograms are summed over a range of incident angles. The algorithm is based on a modified version of Sommerfield's (1896) theorem. Simulations of acoustic and elastic media are included which illustrate the applicability of our plane-wave and line source seismograms for both elastic and acoustic cases.     

FILTER DESIGN FOR SUPPRESSION OF SURFACE MULTIPLES IN A NON-NORMAL INCIDENCE SEISMOGRAM*

In this paper we design a non-linear filter to suppress surface multiples in a non-normal incidence plane wave seismogram of a horizontally-layered elastic (or acoustic) earth model. Our filter is optimal in a least-squares sense and is very efficient in suppressing surface multiples, especially for small incidence angles. The design is based on an extension of earlier work by Mendel on normal incidence Bremmer series decomposition to the case of non-normal incidence, and relies heavily on work by Aminzadeh on a non-normal incidence state space model.     

Variation of Rhenium Contents in Molybdenites from the Sar Cheshmeh Porphyry Cu‐Mo Deposit in Iran

This paper presents the results of a study on the variation of Re in the Sar Cheshmeh porphyry copper‐molybdenum deposit in Iran. The Re content was determined in monomineralic molybdenite separates from veinlets and from molybdenum milled floatation concentrates. Mineralogical data and chemical analyses indicate two generations of veinlet molybdenites in the Sar Cheshmeh deposit. The first generation (molybdenite I), which contains low Re, is found as intergrowth with quartz in silicified rocks with low total sulfide content, and the second generation (molybdenite II), which contains high Re, is found in veinlets with low quartz content. There is a positive correlation between the elevation of molybdenite deposition and Re content. Decrease in the temperature and pH of the hydrothermal fluids are effective for the accumulation Re in molybdenite. X‐ray diffraction study of 14 molybdenite samples indicate that the molybdenites are of the 2H polytype.