Improved NMO correction with a specific application to shallow-seismic data

The application of traditional NMO‐correction techniques in the processing of seismic data may result in severe distortion. This distortion is observed as a decrease in frequency content (NMO stretch) or even a total loss of data for steep velocity gradients. It is shown that the specific nature of the observed distortion is introduced by the particular implementation of the traditional NMO‐correction technique. It is also shown that other techniques that have been proposed in the literature suffer from similar artefacts. An alternative approach is suggested, based on the correction of tapered blocks of seismic data, followed by a coherence filter to compensate for the specific artefacts thus introduced. Correction of seismic amplitudes (e.g. geometrical spreading and attenuation) is implemented as an integral part of the NMO‐correction method thus allowing for both dynamic and kinematic reconstruction of interfering events. From the application of this method to synthetic and field data, it is concluded that the proposed technique may be particularly useful in the processing of shallow seismic (shear‐wave) data. Although not illustrated in this paper, it is emphasized that application of the proposed method is by no means restricted to shear waves; application to P‐wave data will prove equally useful, although the stretch effects involved may not be as severe as in the examples given in this paper.     

Analysis of the resolution function in seismic prestack depth imaging

We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions.     

Iterative resolution estimation in least‐squares Kirchhoff migration

We apply iterative resolution estimation to least‐squares Kirchhoff migration. Reviewing the theory of iterative optimization uncovers the common origin of different optimization methods. This allows us to reformulate the pseudo‐inverse, model resolution and data resolution operators in terms of effective iterative estimates. When applied to Kirchhoff migration, plots of the diagonal of the model resolution matrix reveal low illumination areas on seismic images and provide information about image uncertainties. Synthetic and real data examples illustrate the proposed technique and confirm the theoretical expectations.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

Fractal properties of medium and seismoelectric phenomena

Electrokinetic phenomena in a water-porous medium with a fractal structure above percolation threshold are theoretically investigated. Fracture zone with space-variable porosity is considered as a model of an earthquake hypocenter zone in which the electrokinetic current results from fluid filtration in a fractal pore network. A critical exponent of the streaming potential coefficient is found to depend on both the transport critical exponent and correlation length critical exponent. In this model, logarithmic dependence of electric field amplitude E on the earthquake magnitude M is derived which is compatible with the one observed by the VAN group. Without fractal properties, this form of dependence contradicts the empirical data. The electromagnetic field far from the hypocenter is calculated, which leads to the prediction of weak magnetic field variations. To explain the observed amplitude of VAN's Seismic Electric Signals (SES), the electric source must be at a distance of about 10 km from the registration point if the medium is homogeneous. Therefore, some conductive channel(s) are needed to explain the long distance selective SES transmission.     

Paleolatitudes of the Kerguelen hotspot: new paleomagnetic results and dynamic modeling

The Kerguelen Plateau, a Large Igneous Province in the southern Indian Ocean, was formed as a product of the Kerguelen hotspot in several eruptive phases during the last 120 Myr. We obtained new paleolatitudes for the central and northern Kerguelen Plateau from paleomagnetic investigations on basalts, which were drilled during ODP Leg 183 to the Kerguelen Plateau-Broken Ridge. The paleolatitudes coincide with paleolatitudes from previous investigations at the Kerguelen Plateau and Ninetyeast Ridge (the track of the Kerguelen hotspot) and indicate a difference between paleolatitudes and present position at 49°S of the Kerguelen hotspot. We show that true polar wander, the global motion between the mantle and the rotation axis, cannot explain this difference in latitudes. We present numerical model results of plume conduit motion in a large-scale mantle flow and the resulting surface hotspot motion. A large number of models all predict southward motion between 3° and 10° for the Kerguelen hotspot during the last 100 Myr, which is consistent with our paleomagnetic results.     

Development and testing of numerical methods for two-way nested air pollution modelling

A new 3-D model REGINA (REGIonal high resolutioN Air pollution model) is under development at the National Environmental Research Institute (NERI). The model is based on pieces from several models developed over the last decades at NERI. The aim of the work is to develop a nested model which can operate with very high resolution in both space and time. To fulfill this aim the choice and implementation of accurate numerical methods is crucial. The model will be applied for studying air pollution phenomena (both monitoring, forecasting and scenarios) over Denmark. The present paper is focussed on the modification, implementation and testing of a numerical method for treating the horizontal advection in the model as well as the implementation of two-way nesting techniques. The horizontal transport in the model is solved using an accurate space derivative algorithm. This method is traditionally implemented with periodic boundary conditions, however, this is not an option for nested modelling. A new method for calculating non-periodic boundary conditions has been developed in order to overcome this problem. Extensive testing of the numerical solution of the advection and the coupling of the solution of advection and chemistry in the model using Molenkamp–Crowley rotation tests have been carried out. The results show that the model with the current implementation of numerical methods is suitable for calculating air pollution levels with high resolution.     

Use of Chemotaxonomy To Study the Influence of Benzalkonium Chloride on Bacterial Populations in Biodegradation Testing

Contradictory results are reported for the behaviour of quaternary ammonium compounds (QACs) in sewage treatment plants (STPs). QACs may sorb onto activated sludge. Only little information is available with respect to effects of QACs against bacteria in STPs. Only 5 to 15 % of bacteria present in sewage sludge can be detected by means of culture dependent microbiological methods. The shift of the bacterial populations due to effects of test compounds have not been studied up to now with culture independent methods. The microbial populations shift was studied in situ using culture independent chemotaxonomy profiling ubiquinones and polyamines. Additionally, toxic effects of QACs against bacteria present in the test vessels of the Zahn‐Wellens test (OECD 302 B) were assessed with a toxicity control in the test. The ubiquinone profiles representing changes in Gram‐negative populations mainly showed that the activated sludge was affected only in test vessels containing benzalkonium chloride. According to chemotaxonomy Acinetobacter or/and some members of Pseudomonas spp. have been selected by benzalkonium chloride after some adaptation period (8 to 12 days).