Implications of low-pass filtering on power spectra and autocorrelation functions of turbulent velocity signals

Filtering either through the electronics of an instrument or through digital procedure is performed routinely on geophysical data. When velocity fluctuations are measured in turbulent flows using electromagnetic current meters (ECMs), a builtin lowpass Butterworth filter of order n usually attenuates fluctuations at high frequencies. However, the effects of this filter may not be acknowledged in turbulence studies, thus impeding comparisons between data collected with different ECMs. This paper explores the implications of the filters on the characteristics of velocity signals, mainly on variance, power spectra, and correlation analyses. Variance losses resulting from filtering can be important but will vary with the order n of the Butterworth filter, decreasing as n increases. Knowing the filter response, it is possible to reconstruct the original signal spectrum to evaluate the effect of filtering on variance and to allow comparisons between data collected with different instruments. The autocorrelation function also is affected by filtering which increases the value of the coefficients in the first lags, resulting in an overestimation of the integral length scale of coherent structures. These important effects add to those related to size and shape differences in ECM sensors and must be taken into account in comparative studies.     

Electromagnetic signals associated with stick–slip of quartz-free rocks

We recorded clear transients in the electric and magnetic fields upon sudden slip in stick–slip experiments on dry, cylindrically shaped, quartz-free rock specimens of basalt and peridotite with a 30° saw-cut (representing a fault) at confining pressures of up to 120 MPa. The amplitudes of the measured electric field signals were always higher at the electrode pair oriented parallel to the strike of the fault than at the pair oriented perpendicular. This anisotropy suggests a preferred electric polarization normal to the slip surface. The transients in the electric and magnetic fields were observed only when the fault slip occurred by stick–slip mode, not by a stable mode of the sliding, and the amplitudes of the electric field signals increased with increasing stress drop. It is suggested that the generation process of the electromagnetic signals is closely related to the characteristic behavior of the fault at the time of the initiation of slip during stick–slip events, probably with respect to the intensity of the signals. We propose that one or both of the following two processes characteristic of the fault at the time of the initiation of slip during stick–slip events are essential for the generation of detectable electromagnetic signals: rapid slip along the simulated fault and separation of the rock masses across the fault.     

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Acoustic and electromagnetic properties of soils saturated with salt water and NAPL

Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

云闪放电通道发展及其辐射特征

利用闪电宽带干涉仪系统,对中国南方(广东)地区云闪时空演变特征、辐射及其相应电场变化特征进行分析研究。根据云闪电场变化波形,云闪放电过程可划分成活跃阶段和最后阶段,辐射源定位结果表明,云闪放电起始于向上发展的负击穿过程,通道向上发展的速度约为3～3.3×105m·s-1。云闪放电的主通道在活跃阶段形成,该期间辐射源随时间演变和相应电场变化表明,云内电荷结构具有上正下负的偶极性电荷结构。云闪的最后阶段辐射源主要在早期形成的通道内出现,其辐射源活动特征与地闪的回击过程比较相似;云闪辐射能量主要集中在2～3MHz以下的低频段,且辐射强度随频率增加迅速减弱。     

对波达波夫和Pride震电波方程组的对比分析

用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义，解释了第一类和第二类震电效应的意义，在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较，说明了二者在描述第二类震电效应方面的异同点. 同时指出 ：波达波夫方程组忽略了流体与固体的耦合质量；方程中的黏性耗散项丢掉了一个孔隙度因 子，依据该方程组计算出的弹性波和转换电场的幅度将偏大；边界条件之一存在错误，会影 响对波在界面上的反射透射规律的描述.     

多相流动电磁波成像测井测量敏感场计算

在多相流动电磁波成像测井研究中，测量敏感场是设计测量仪器和研究图像重建算法的前提工作. 本文用数值方法计算多相流动电磁波成像测井测量敏感场. 根据测量探头电极阵列的对称性，将三维敏感场问题划分为轴向和横向两个平面问题，然后应用有限元算法求解两个平面电势分布问题，对于不同介质模型，计算各种测量条件下测量区域内电势分布，进而根据电势分布计算测量敏感场. 结果表明测量敏感区域为发射电极到测量电极的弧形区域，在发射电极到测量电极的电势梯度线上，测量敏感函数近似于一个指数函数.     

MDCB型地震前兆监测仪映震效果分析

对ＭＤＣＢ－Ⅱ型地震前兆监测仪观测获得的类地应力值，异常频度和８方位异常值等参数进行了比较深入的分析及其与地震关系的研究。资料反映的事实表明，在一定距离，一定强度的地震的短临阶段，确定存在着电磁波异常现象，其动态趋势，变化次数，异常突变和异常方向性反映等参数，可为地震预测预报研究提供科学判据。