Signal detectability of marine electromagnetic methods in the exploration of resistive targets

We compare selected marine electromagnetic methods for sensitivity to the presence of relatively thin resistive targets (e.g., hydrocarbons, gas hydrates, fresh groundwater, etc.). The study includes the conventional controlled‐source electromagnetic method, the recently introduced transient electromagnetic prospecting with vertical electric lines method, and the novel marine circular electric dipole method, which is still in the stage of theoretical development. The comparison is based on general physical considerations, analytical (mainly asymptotic) analysis, and rigorous one‐dimensional and multidimensional forward modelling. It is shown that transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods represent an alternative to the conventional controlled‐source electromagnetic method at shallow sea, where the latter becomes less efficient due to the air‐wave phenomenon. Since both former methods are essentially short‐offset time‐domain techniques, they exhibit a much better lateral resolution than the controlled‐source electromagnetic method in both shallow sea and deep sea. The greatest shortcoming of the transient electromagnetic prospecting with vertical electric lines and marine circular electric dipole methods comes from the difficulties in accurately assembling the transmitter antenna within the marine environment. This makes these methods significantly less practical than the controlled‐source electromagnetic method. Consequently, the controlled‐source electromagnetic method remains the leading marine electromagnetic technique in the exploration of large resistive targets in deep sea. However, exploring laterally small targets in deep sea and both small and large targets in shallow sea might require the use of the less practical transient electromagnetic prospecting with vertical electric lines and/or marine circular electric dipole method as a desirable alternative to the controlled‐source electromagnetic method.     

EQUIVALENCE IN ELECTROMAGNETIC (FREQUENCY) SOUNDING*

An analysis of the equivalence phenomenon, analogous to that encountered in resistivity sounding, has been made for electromagnetic (frequency) sounding with the following systems: horizontal coplanar coils, perpendicular coils, vertical coplanar coils, and vertical coaxial coils. Cases analyzed are three-layer H and K type resistivity distributions. The theoretical responses for the EM sounding systems have been computed by the digital linear filter method using short filters. An analysis has been made of the resolution by the four EM systems considered f the equivalence effect. It is concluded that from the equivalence point of view electromagnetic (frequency) sounding has relatively better resolution compared to resistivity sounding for the H type three-layer cases. For the K type three-layer cases electromagnetic (frequency) sounding provides a serious problem in resolution similar to that encountered in resistivity sounding. Empirical relationships have been established for the H and K type cases respectively in order to numerically quantify the equivalence involved in electromagnetic (frequency) sounding.     

LEAST-SQUARES FILTERS WITHOUT TRANSIENT ERRORS: AN EXAMINATION OF THE ERRORS IN LEAST-SQUARES FILTER DESIGN*

A new approach has been developed for the design of cross-equalization filters by the least-squares method. The filters estimated by this new exact method are subject to only two types of error: bias and random error. Cross-equalization filters estimated by a more conventional least-squares method are further subject to “transient error”. This type of error becomes important when designing filters from a data gate of a length comparable with the length of the filter, i.e., less than four times the length of the filter. The effect of altering various design parameters has been investigated for the new method. It has been found that the proportion of bias in the filter decreases as the effective filter length increases, whereas the random error in the filter decreases with increase in either the signal-to-noise ratio of the data or the ratio of the data duration to the filter length. The level of whitening applied to the auto-correlation matrix before inversion was not found to be a critical design parameter. Also, two techniques have been tested for reducing any anomalous d.c. component in the calculated filter.     

Hankel transform filters for dipole antenna radiation in a conductive medium

We discuss the Hankel transforms related to a particular application, i.e. the dipole antenna radiation in conductive media, such as the antenna radiation in sea-bed electromagnetic applications. In this application, the electromagnetic wavefields decay very rapidly with distance. A good filter means that it can be used to evaluate weak fields. Exponential sampling transforms a Hankel transform into a convolution equation, which must be solved to obtain the filter coefficients. Here, we use a direct matrix inversion method to solve the convolution equation in the sample domain, instead of the Fourier transform method and the Wiener–Hopf method, previously used to solve the convolution equation. This direct method is conceptually simple and is suitable for our optimization process: by using the Sommerfeld identity, we search for the optimum sampling interval, which corresponds to the minimum wavefield, evaluated for a given length filter. The performances of the new filters obtained are compared with some well-known filters. We find that our filters perform better for our application; that is, for the same length filters, our filters are able to calculate weaker fields. For users working in similar applications, three sets of filters with lengths 61, 121, 241 are available from the author.     

New digital linear filters for Hankel J0 and J1 transforms

The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J₀ transform (61-point and 120-point operators), and two for the Hankel J₁ transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics.     

A NOTE ON THE LINEAR FILTER METHOD OF COMPUTING ELECTROMAGNETIC SOUNDING CURVES *

The values of the filter coefficients used for the computation of electromagnetic sounding curves are studied in conjunction with the values of the input function to the filter, or the range of values which the input function may assume, and the filters are broken off at such a place that the error in the sum of the products of filter coefficient and input function does not exceed a prescribed value. This analysis is carried out for the horizontal coils system, the perpendicular coils system, and the vertical coplanar coils system. The lengths of the filters so derived depend on the layer parameters, the frequency and the coil spacing. Even in the most unfavourable cases the filters are shorter than the filters published by Koefoed, Ghosh, and Polman (1972).     

VLF RESISTIVITY MAPPING AND VERTICALIZATION OF THE ELECTRIC FIELD1

For over 20 years, powerful VLF transmitters have been used as electromagnetic sources for subsurface investigations in mining exploration. Measurements initially concerned the vertical component of the magnetic field or the inclination of the field and were later extended to measurement of the horizontal electric field in the direction of the transmitter, to determine the resistivity of the terrain. Measurement of the electric field is usually performed with electric lines, grounded or not, with lengths of at least 5 m. This paper presents the concept of a VLF resistivity meter with a very short electric sensor (1 m) and the results obtained with it. This technique improves the measurement of the electric field, which is in principle a point value. It also permits a higher spatial sampling rate and, by closely linking the electric sensor with the magnetic sensor on a lightweight mount, makes it possible for the instrument to be used by a single operator. In addition, transformation of the electric field data, analogous to reduction to the pole in magnetism, is proposed to correct the horizontal deformation of the anomalies created by polarization of the primary field. Comparison with direct current electrical measurements shows highly satisfactory correlations. This transformation, validated for VLF, can be extended to any electrical or electromagnetic method using a uniform primary field, i.e. gradient array in direct current or low-frequency magnetotellurics. We call this verticalization of the electric field. Resistivity measurements and mapping using the VLF frequency range can be applied not only to mining but also to a wide range of shallow geophysical studies (hydrology, civil engineering, etc.) and are not limited to problems concerning the location of conductive targets     

Electric and electromagnetic processes in the near-Earth atmosphere before earthquakes on Kamchatka

Simultaneous records of the diurnal variations in the quasistatic electric field in the near-Earth atmosphere, fluxes of discrete electromagnetic pulses in the VLF band, source azimuths, narrowband filter output emission intensity at frequencies of 4.65 and 5.3 kHz, and time forms and spectra of VLF pulses have been analyzed. The anomalous behavior of these parameters in October 2002 and August 2004 with different time delay was accompanied by earthquakes near the southeastern coast of Kamchatka at distances of 250–400 km from the registration points. Based on the results of a fine frequency-time analysis of the broadband records of VLF signals, it has been indicated for the first time that discrete electromagnetic pulses observed in anomalous fluxes before earthquakes were signals of local thunderstorm processes.     

Comparison of deterministic ensemble Kalman filters for assimilating hydrogeological data

Groundwater models are critical decision support tools for water resources management and environmental remediation. However, limitations in site characterization data and conceptual models can adversely affect the reliability of groundwater models. Therefore, there is a strong need for continuous model uncertainty reduction. Ensemble filters have recently emerged as promising high-dimensional data assimilation techniques. Two general categories of ensemble filters exist in the literature: perturbation-based and deterministic. Deterministic ensemble filters have been extensively studied for their better performance and robustness in assimilating oceanographic and atmospheric data. In hydrogeology, while a number of previous studies demonstrated the usefulness of the perturbation-based ensemble Kalman filter (EnKF) for joint parameter and state estimation, there have been few systematic studies investigating the performance of deterministic ensemble filters. This paper presents a comparative study of four commonly used deterministic ensemble filters for sequentially estimating the hydraulic conductivity parameter in low- and moderately high-dimensional groundwater models. The performance of the filters is assessed on the basis of twin experiments in which the true hydraulic conductivity field is assumed known. The test results indicate that the deterministic ensemble Kalman filter (DEnKF) is the most robust filter and achieves the best performance at relatively small ensemble sizes. Deterministic ensemble filters often make use of covariance inflation and localization to stabilize filter performance. Sensitivity studies demonstrate the effects of covariance inflation, localization, observation density, and conditioning on filter performance.     

Defining parameters for Eckhardt's digital baseflow filter

Digital filters are useful tools for assessing the contribution of groundwater to total river flow. Several of those filters have been proposed in the last decades. One of the last contributions on this subject was given by Ekchardt (2005) who proposed a more general form of a digital baseflow filter and showed that some of the most used filters are special cases of this general form. This new filter has the inconvenience of having two parameters, one of them may be obtained directly from recession analysis, but the other (maximum baseflow index (BFImax)) is routinely estimated by a priori defined values according to the predominant geological characteristics of the drainage basin. In this short communication, we propose a method to estimate BFImax by a backwards filtering operation. The method was applied using data from 15 gauging stations in Brazil, with a varied range of groundwater contribution to streamflow. Results of the new estimation method for the BFImax parameter are coherent with the values which would be adopted by considering geological characteristics of the river basins. Copyright © 2012 John Wiley & Sons, Ltd.     

A Multilevel Ground Water Monitoring System for Casing Advance and Bedrock Drilling Methods

A single-hole multilevel sampling piezometer system (MLSPS) has been designed by the Geological Survey of Canada (GSC) to be installed using drilling systems that continuously core (e.g., Rotosonic) or continuously sample (e.g., hollow-stem auger, Becker hammer) overburden and that have the flexibility of allowing additional coring (diamond drilling) or sampling (hammer drilling) of bedrock. The GSC-MLSPS (under license to Solinst Canada Ltd.) uses a patented GSC dry injection system for accurate emplacement of filter packs and seals. This system permits (a) the use of variable screen lengths; (b) the complete evacuation of piezometers before introduction of new ground water (no bailing); (c) the use of a number of types of hydraulic tests (e.g., slug, withdrawal/recovery, vacuum, pressure-pulse); (d) ground water sampling under a nitrogen atmosphere; (e) dissolved gas sampling; (f) a great deal of flexibility in the use of design materials; and (g) the elimination of bridging and collapse of filter packs and seals.     

Factors affecting sediment trapping in vegetated filter strips: simulation study using VFSMOD

Soil and water conservation practices have been promoted for a long time, in order to sustain agricultural activities and prevent environmental pollution. Vegetated filter strips (VFS) have been used to reduce sediment pollution into water bodies at or near the pollutant source. However, factors effecting VFS performance under natural conditions have not been well understood owing to the physical, time and financial limitations of field experiments. The use of well‐validated simulation models to understand the performance of VFS and factors affecting sediment deposition is highly justified. The objective of this research is to investigate sediment trapping in VFS and to study various factors affecting VFS performance using the simulation model VFSMOD, which was developed by researchers at University of North Carolina. Recently, VFSMOD has been validated successfully by using 21 filters with varying length, slope and vegetated cover. A wide range of five parameters was selected for the simulations, namely filter length, filter slope, manning roughness coefficient, soil type and characteristics of incoming sediment from adjacent fields. Computer simulations revealed that the length of filter is the most significant factor affecting sediment trapping in VFS. The relative increase in trapping efficiencies was not linearly related to an increase in filter length. Inflow sediment class also has a major influence on sediment trapping in VFS. The trapping efficiency of clay sediments in a 15 m length VFS was 47% compared with 92% for silt from incoming sediment. Manning roughness coefficient had a moderate effect on sediment trapping and was more significant in short filters. Land slope and soil type of VFS had a minor influence on the performance of VFS. Copyright © 2001 John Wiley & Sons, Ltd.     

A biological oil adsorption filter

A new oil adsorption method called adsorption filtration (AF) has been developed. It is a technology where by oil residues can be cleaned from water by running it through a simple filter made from freeze treated, dried, milled and then fragmented plant material. By choosing suitable plants and fragmentation sizes it is possible to produce filters, which pass water but adsorb oil. The aim of this study was to investigate the possibilities of manufacturing oil adsorbing filter materials from reed canary grass (Phalaris arundinacea), flax (Linum usitatissimum L.) or hemp fibre (Cannabis sativa L.). The oil (80 ml) was mixed with de-ionised water (200 ml) and this mixture was filtered through 10 or 20 g adsorption filters. Fine spring harvested hemp fibre (diameter less than 1 mm) and reed canary grass fragments adsorb 2-4 g of oil per gram of adsorption material compared to 1-3 g of water. Adsorption filtration is thus a novel way of gathering spilled oil in shallow coastal waters before the oil reaches the shore.     

On the Use of Crushed Recycled Glass Instead of Silica Sand in Dual‐Media Filters

Crushed recycled glass was evaluated as an alternative to silica sand in dual‐media filters. Pilot scale inline filtration experiments were carried out using raw waters from three different water sources with turbidities between 6.0 and 14.0 NTU. Two physically identical filter columns were operated in parallel in the experiments. One filter consisted of 62.5 cm silica sand and 41.5 cm anthracite coal, whereas the other filter contained 62.5 cm crushed recycled glass plus 41.5 cm anthracite coal. The total bed depth was 104 cm for both filters. The properties of the media were as follows: Glass effective size = 0.77 mm, uniformity coefficient = 1.41. Sand effective size = 0.79 mm, uniformity coefficient = 1.33. Coal effective size = 1.45 mm, uniformity coefficient = 1.39. Experiments were repeated five times as follows: (i) Without the use of a coagulant, (ii–iii) with 5 and 10 mg/L of alum, and (iv–v) with 5 and 10 mg/L of ferric chloride. The filtration rate used was 11.5 m/h. Turbidity, particle counts, and head losses were measured and compared as functions of time. The following were observed: (i) Effluent turbidities and particle counts of the two filters were very close, i.e., essentially the same effluent quality was obtained when crushed glass was used instead of silica sand. (ii) In the majority of the tests, the filter with crushed glass generated both a smaller clean‐bed head loss and smaller clogging head losses than those of the filter containing sand. It is concluded that crushed glass may be a good alternative to silica sand in dual‐media filtration.     

Bias aware Kalman filters: Comparison and improvements

This paper reviews two different approaches that have been proposed to tackle the problems of model bias with the Kalman filter: the use of a colored noise model and the implementation of a separate bias filter. Both filters are implemented with and without feedback of the bias into the model state. The colored noise filter formulation is extended to correct both time correlated and uncorrelated model error components. A more stable version of the separate filter without feedback is presented. The filters are implemented in an ensemble framework using Latin hypercube sampling. The techniques are illustrated on a simple one-dimensional groundwater problem. The results show that the presented filters outperform the standard Kalman filter and that the implementations with bias feedback work in more general conditions than the implementations without feedback.     

Experimental investigation of runoff reduction and sediment removal by vegetated filter strips

The impact of vegetated filter strips (VFS) on sediment removal from runoff has been studied extensively in recent years. Vegetation is believed to increase water infiltration and decrease water turbulence thus enhancing sediment deposition within filter media. In the study reported here, field experiments have been conducted to examine the efficiency of vegetated filter strips for sediment removal from cropland runoff. Twenty filters with varying length, slope and vegetated cover were used under simulated runoff conditions with an average sediment concentration of 2700 mg/L. The filters were 2, 5, 10 and 15 m long with a slope of 2·3 and 5% and three types of vegetation. Three other strips with bare soil were used as a control. The experimental results showed that the average sediment trapping efficiency of all filters was 84% and ranging from 68% in a 2‐m filter to as high as 98% in a 15‐m long filter compared with only 25% for the control. The length of filter has been found to be the predominant factor affecting sediment deposition in VFS up to 10 m. Increasing filter length to 15 m did not improve sediment trapping efficiency under the present experimental conditions. The rate of incoming flow and vegetation cover percentage has a secondary effect on sediment deposition in VFS. Copyright © 2004 John Wiley & Sons, Ltd.     

Anomalous propagation of VHF radiowaves behind the horizon in the seismic region

A theory of electromagnetic radiation generation by random electric discharges in the troposphere and VHF radiowave scattering by these discharges has been developed. The discharge model, which makes it possible to calculate the spatiotemporal distribution of the discharge channel conductivity depending on the electric current value in this discharge, has been obtained. The electromagnetic radiation spectrum in the troposphere occupied by random discharges has been calculated. VHF electromagnetic wave scattering by random electric discharges in the troposphere has been considered. Equations have been derived, and the method for calculating the average value of the electromagnetic wave field scattered by random discharges has been developed. The calculations indicated that the scattered wave field amplitude is much larger than the diffraction wave field amplitude behind the horizon. The theoretical results agree with the observations of the electromagnetic radiation and VHF transmitter signals behind the horizon relative to the earthquake epicenter during the earthquake preparation.     

Three-dimensional EM computer simulation on sprite initiation above a horizontal lightning discharge

The purpose of this paper is to study the three-dimensional (3D) effects of the source current and the electromagnetic (EM) pulse on the distribution of upper atmospheric electric field and ionization created. A new lightning model has been employed; i.e., we have included a horizontal channel in addition to the conventional vertical channel. The effects of the horizontal lightning channel are summarized as follows: (1) the effect of a vertical channel plays a fundamental role in the whole view of sprite initiation; (2) the position of a sprite is shifted from the position of its parent vertical channel in response to the length of a horizontal channel; and (3) we observe very fine structures with local maxima and minima in the “reduced” electric field (the electric field divided by the neutral gas density). These theoretical inferences are discussed and compared with the observational facts so far reported (such as lateral shift of sprites, morphological difference of sprites, etc.).     

Distributed electric field sensing using fibre optics in borehole environments

In the past decade, rapid advances in distributed optical fibre sensing technologies have made it possible to record various geophysical data (e.g. strain, temperature and pressure) continuously in both time and space along the fibre, providing an unprecedented quantity and spatial density of data compared to traditional geophysical measurements as well as reducing data acquisition cost. To date, no distributed fibre-based electromagnetic field sensing system has been implemented although electromagnetic sensing could have a broad range of applications to geophysical imaging and monitoring in borehole environments. The goal of this paper is to provide a theoretical feasibility study regarding the design and use of an electromagnetic sensing optical fibre for geophysical applications. First, we present the sensitivity analysis of a ‘hypothetical’ optical fibre coated with polyvinylidene fluoride, a polymer that provides relatively high piezoelectric properties, yet unlike ceramics, is flexible. Using a two-dimensional electromagnetic modelling algorithm, we simulate the earth electric-field-to-fibre-strain transfer function and estimate the theoretical sensitivity of the optical fibre to electric fields. Given the state-of-the-art distributed acoustic sensing strain sensitivities in the picometres strain range, our numerical modelling analysis suggests that a perfectly coupled polyvinylidene fluoride–coated optical fibre can measure electric field values in the mV/m to V/m amplitude range. We then apply a cylindrically symmetric modelling algorithm to simulate numerical models demonstrating the applicability of such a fibre in an oilfield environment. Scenarios investigated employ an electric field source and suggest that the measurements can be used to distinguish the oil versus water ratio with a fibre mounted inside a producing steel cased oil well as well as distinguishing between brine and hydrocarbon filled reservoir zones with a fibre located outside of the casing.     

分离重磁区域场与局部场的维纳滤波器

本文从最佳线性滤波理论出发,对目前重磁资料数据处理中分离区域场与局部场的两种滤波器--匹配滤波和维纳滤波的频率响应特性作了分析比较,指出了匹配滤波只是一般维纳滤波的一个特例。将该两种滤波器与一般情况的维纳滤波器的误差作了对比,并通过简单的理论试例,说明它们的局限性和应用范围。