Interpretation of total intensity magnetic anomalies when anomalous field is arbitrarily large

Summary Not infrequently, in mining geophysics, the anomalous field of the magnetized body is appreciably large and it varies from the direction of the earth's normal field within the vicinity of the magnetized body. Total magnetic intensity data collected on the ground over shallow magnetized bodies cannot be interpreted quantitatively, since all the available methods of interpretation assume that the resultant field lies in the direction of earth's normal field as the anomalous field is small. For geological bodies that can be magnetically represented by infinite line pole, point pole, infinite line dipole, and point dipole, this paper gives quantitative interpretation procedures for total magnetic intensity anomalies without the above assumption using the characteristic points from the theoretical curves (not reproduced). These characteristic curves can be used for the direct determination of depths and offsets by using information from suitably chosen field profiles.Contribution No. N.G.R.I.-71-233.     

成层粘弹性土中桩土耦合纵向振动时域响应研究

从三维轴对称角度出发,采用粘性阻尼粘弹性连续土介质模型,考虑桩土相互作用效应,对成层土中桩土纵向耦合振动时的桩顶时域响应进行了解析研究。求解时,首先建立定解问题,然后利用拉氏变换先对底部土层进行求解得到其振动位移形式解,然后利用桩土接触界面连续条件来考虑桩土耦合作用,分析底层土中桩段的动力反应,然后利用桩段阻抗函数的传递性,进行逐层递推求解,最终得到桩顶时域和频域响应的半解析解。通过参数影响分析和与工程实测曲线的对比,讨论分析了成层土中桩土耦合振动的响应特性,验证了本文解。基于本文研究可为桩基抗震、防震设计、桩基动力检测提供新的理论支持。     

Tracer measurements in groyne fields for the quantification of mean hydraulic residence times and of the exchange with the stream

Tracer experiments were carried out in groyne fields (GF) of the river Elbe near Havelberg (Germany) in order to estimate the hydraulic connectivity with the river channel. The characteristic times of the five groyne fields, which were estimated from the exponentially declining tracer curves in 43 runs, ranged between 15 min and 69 min and did not correlate with the water level. Methodological investigations show that single point injection and two measurements (in the outflowing water and in the dominant region) are sufficient to provide robust in‐situ tracer curves. Using simplified mathematical simulations with connected stirred tanks, the conditions are investigated for the development of breaks in tracer curves and for the occurrence of significant errors in the estimation of intrinsic residence times. It was shown that an initial uniform dye distribution is not mandatory for the estimation during steady states. In special cases, point injections are more advisable. Moreover, the mean hydraulic residence time was found to be not equivalent to the estimated characteristic time. In fact, it is mostly overestimated by tracer experiments. The degree of overestimation depends on mixing and the volumetric proportions between the different parts of the GF and can be calculated from measured dye concentration differences. For example, an overestimation of 32% was calculated for a groyne field with a commonly found circulation flow pattern.     

区域地电观测网地电场日变化模型研究

应用泰勒多项式和曲面样条函数法建立了我国地电场观测台网中华东、 华北两个区域台网的地电场日变化(K_p≤5)模型,并对拟合结果进行了对比分析. 结果表明: ① 尽管地电场变化复杂,但用这两种方法建立的日变化模型相当吻合,模型曲线随时间的变化符合样本曲线随时间的变化; ② 用模型曲线可以描述其它年份同农历日期、 同等磁情时实测的地电场日变化; ③ 用模型曲线描述实测地电场日变化时,低磁情时得出的结果要优于高磁情时的结果; ④ 用多日样本曲线得到的模型曲线更能清晰地展示出地电场日变化的两次起伏波形、 幅度以及极值集中在午前午后等变化特点. 今后建立全国各区域台网不同磁情、 不同农历日期的地电场日变化模型可定量评价区域台网的地电场观测资料的质量,同时根据该模型能够分离出地电场的非正常变化,进而研判真实的异常,以服务于地震等灾害事件的预测以及地球电磁环境变化的监测等方面.      

A simple formula for shape and depth determination from residual gravity anomalies

This paper presents a simple method for shape and depth determination of a buried structure from residual gravity anomalies along profile. The method utilizes the anomaly values of the origin and characteristic points of the profile to construct a relationship between the shape factor and depth of the causative source. For fixed points, the depth is determined for each shape factor. The computed depths are then plotted against the shape factor representing a continuous monotonically increasing curve. The solution for the shape and depth of the buried structure is then read at the common intersection point of the depth curves. This method is applied to synthetic data with and without random errors. Finally, the validity of the method is tested on two field examples from the USA.     

Aquifer parameter estimation using an incremental area method

Theoretical well functions have been derived over the years to predict ground water level behaviour in aquifer systems under stress owing to groundwater extraction. The drawdown data collected during pump tests are typically analysed using graphical curve‐matching procedures to estimate aquifer parameters based on these well functions. Difficulty in aquifer characteristic identification and parameter estimation may arise when the field data do not perfectly match the drawdown curves obtained from the well functions. The present study provides a new method for the interpretation of aquifer pump tests which supplements the existing curve‐matching procedures in case ideal conditions do not exist; the proposed method provides a greater degree of flexibility in the data analysis for diagnostic tool purposes. The method, referred to as the Incremental Area Method (IAM) is based on integrating the logarithmic‐based drawdown curves within a discrete time and matching the results with a corresponding time integral of the Theis ( 1935 ) Well Function which governs ideal confined aquifers. The application of the proposed method to synthetically generated data and field data showed that IAM represents a viable method which yields information on potential non‐idealness of the aquifer and provides aquifer parameter estimates thus potentially overcoming drawdown data curve‐matching difficulties. Copyright © 2011 John Wiley & Sons, Ltd.     

Electromagnetic induction in a half-space with a cylindrical inhomogeneity

Summary The problem of electromagnetic induction in a half-space with a cylindrical inhomogeneity is treated. Solutions for TE and TM polarization of the exciting electromagnetic field are given, which can be used for computing sounding and profiling curves. The anisotropy of the surface impedance is pointed out, as well as a whole series of other interesting properties of the magnetotelluric field in this model; some of them are also demonstrated on the computed theoretical curves.     

磁性层状介质的CSAMT电场分量响应特征

在传统的勘探电磁学理论中,往往把岩层磁导率近似看成空气中的磁导率,从而得到简化后的电磁场响应式,主要用于指导非磁性介质情况下的电磁探测.本文开展了磁性介质的电磁响应理论、数值模拟和特征分析等研究,以便实现对磁性矿体的电磁精细探测.首先推导出磁性层情况下的电磁响应表达式;然后,计算了几种典型磁性介质地质断面的响应曲线,并对其响应特征进行了分析;研究结果表明:当磁性层位于地表时,当表层相对磁导率μr1.2,磁性层对电磁测深曲线具体明显的影响;当岩层为中间层时,只要μr1.2,对H型、A型曲线中段会有影响.对H型曲线,主要表现在使极小值处的曲线抬高、变宽,以致使曲线形态发生变化;当磁性层处于最后一层时,不论岩层的磁导率是多少,对各类二、三层曲线来讲,不会影响曲线的形态,只是使渐近线前第一个极值点的位置右移,极值点至渐近线之间的线段变得较为陡峭.     

Determination of dispersive phase velocities for SASW method using harmonic wavelet transform

The spectral analysis of surface waves (SASW) method is an in situ, seismic method for determining the shear wave velocity (or maximum shear modulus) profile of a site. The SASW test consists of three steps: field testing, evaluation of dispersion curve by phase unwrapping method, and determination of shear modulus profile by inversion process. In general, field testing and dispersion curve evaluation are regarded as simple work. However, because of characteristic of Fourier transform used in the conventional phase unwrapping method, dispersion curve is sensitive to background noise and body waves in the low frequency range. Furthermore, under some field conditions such as pavement site, the usual phase unwrapping method can lead to erroneous dispersion curve. To overcome problem of the usual phase unwrapping method, in this paper, a new method of determining dispersion curve for SASW method was applied using time–frequency analysis based on harmonic wavelet transform as an alternative method of a current phase unwrapping method. To estimate the applicability of proposed method to SASW method, numerical simulations at various layered soil and pavement profiles were performed and the dispersion curves by proposed method are more reliable than those by the usual phase unwrapping method.     

基于泰勒多项式及曲面样条函数的地电场变化分析

选取江苏及邻区地电台网13个台站地电场观测分钟值,采用泰勒多项式和曲面样条函数模型建立地电日变化模型,使用模型曲线描述不同年份、同农历日期、同等磁情天气的实测地电日变化,结果表明:使用分钟值曲线建立的地电日变化模型,更能清晰反映日变化特点,并可反映一天当中台站观测环境、测量系统等变化引起的地电变化;使用同一区域地电台网中若干台站同方向测道的样本曲线,建立各测道电场变化模型,可反映其自身特征及共性特征。     

Wavelet analysis of characteristic length scales and orientation of two-dimensional heterogeneous porous media

The global wavelet energy spectrum (GWES) has been used on one-dimensional transects of porous media to identify the characteristic length scale of a material property. The characteristic length scale is a measure of the distance over which property values are correlated. We extend the wavelet analysis technique to two-dimensional stationary porous media for identifying both the characteristic length scales and orientation of heterogeneity. We theoretically develop and evaluate the GWES for isotropic and anisotropic random fields using the two-dimensional Mexican hat wavelet. The relationship between the wavelet scale and the characteristic length scale of a two-dimensional random field is investigated for different covariance structures. The ability of the GWES to identify multiple characteristic length scales and orientations in a random field with nested covariance structures is also investigated. We use the technique to identify the characteristic length scales of a log hydraulic conductivity field used in a laboratory experiment.     

THE INTERPRETATION OF INDUCED-POLARIZATION SOUNDING CURVES IN THE FREQUENCY DOMAIN *

Two-layer type curves of apparent frequency effect for the Wenner configuration are presented. The formulation is based on the normal definition of frequency effect in terms of resistivities measured at different frequencies plus the definition of apparent resistivity over two horizontal layers as a function of first and second layer resistivities. The use of these type curves in the interpretation of multilayer apparent frequency-effect curves is described and some field examples are given.     

Magnetic-doublet method of interpretation of arbitrarily large total intensity anomalies

Summary In mining geophysics, analysis of total intensity magnetic anomalies in terms of equivalent magnetic-doublets has been extensively used for approximate determination of depth, offset and the length of the disturbing bodies. Total intensity data collected on the ground over shallow magnetized bodies cannot be interpreted quantitatively by the available methods of interpretation which assume that the resultant field lies in the direction of Earth's normal field as the anomalous field vector is small when compared to the normal field vector. In the present investigation, the above assumption has been eliminated and total intensity magnetic anomalies due to disturbing bodies are analyzed in terms of their magnetic-doublet equivalents. From sets of theoretical total magnetic anomaly curves (some are reproduced) Tables have been prepared using the characteristic points and are shown in Tables 1 to 7. These Tables can be used for the direct determination of depth, offset, etc., by using information from suitably chosen field profiles.Contribution No. NGRI 286.     

A Numerical Approach to Depth Determination from Residual Gravity Anomaly Due to Two Structures

—The residual anomaly expression due to two simple geologic structures can be represented by the sum of fields produced by two known source bodies whose depths are to be determined. For a few fixed characteristic points and distances along the residual anomaly profile, the depth to the upper structure is determined iteratively using a simple formula for each depth of the lower structure. The computed depths to the lower structure are plotted against the depths to the upper structure representing a depth curve. The solution for the depths to both structures is read at the common intersection of the depth curves. Procedures are also formulated to estimate the radii of the two buried structures. The method is a powerful technique for separating the composite residual anomaly into two residual components due to two structures. The method is applied to theoretical data with and without random errors and tested on a field example from U.S.A.     

CM4模型数据与台站实测数据的对比研究

本文利用第四代地磁场综合模型(Comprehensive Model 4,CM4),计算了1982—2002年隆尧地磁台站的磁层源磁场及其感应场、电离层源磁场及其感应场的地磁北向分量X、东向分量Y、垂直分量Z,分析了各场源磁场随时间的变化特征。磁层源磁场及其感应场呈现出11年和27天的周期性变化,有些年的27天周期性变化显著,有些年则不太显著;电离层源磁场及其感应场具有明显的季节变化,不同年相同季节的变化形态一致但幅度不同;日变化分析显示,磁静日和磁扰日期间的模型数据与台站实测数据变化一致性较好,相关性较高。     

Direction and relative paleointensity of the geomagnetic field in the middle and later pleistocene from paleomagnetic data of the Roxolany section (Ukraine)

The lower part of the Roxolany section (Ukraine) is studied, and paleomagnetic and petromagnetic characteristics of rock samples are obtained for the time interval 300-75 ka. Detailed curves of variations in angular parameters of the geomagnetic field are constructed in the entire time interval, and curves of variations in the relative paleointensity are obtained for the interval 300-180 ka. Using the values of the angular parameters and the deviations of the virtual geomagnetic pole from the position of the stationary field, anomalous directions compatible with the Jamaica excursion are identified in the intervals 250-249 and 221-220 ka. The geomagnetic field evolution is studied by methods of wavelet analysis, and the field generation process is shown to vary in the interval 300-180 ka.     

高光谱遥感技术在建 (构) 筑物震害识别中的应用

高光谱遥感作为20世纪空间对地观测技术重大进步的产物,通过其较高的光谱分辨率,为人们提供了丰富的地球表面信息,在各个研究领域得到了快速发展和广泛应用,并取得了卓越的成果。尽管高光谱遥感具有独特的优势,但是针对其在震害评估领域中应用的相关研究较少。本文在总结高光谱遥感的特征、优势及不同领域应用现状的基础上,开展了其在震害评估领域的应用研究。基于ASD地物波谱仪获取的建（构）筑物光谱曲线构建可用于震害分析所需的光谱特征库,对比光谱库中地物曲线之间的差异后,发现高光谱遥感在震害评估领域中的应用是可行的,因不同震害地物之间的光谱特征曲线存在差异,依据这种差异可区分不同的震害信息,从而对震害遥感图像进行精细分类。     

Investigating peat hydrological properties using field and laboratory methods: application to the Lanoraie peatland complex (southern Quebec,Canada)

Flow dynamics within a peatland are governed by hydraulic parameters such as hydraulic conductivity, dispersivity and specific yield, as well as by anisotropy and heterogeneity. The aim of this study is to investigate hydraulic parameters variability in peat through the use of different field and laboratory methods. An experimental site located in the Lanoraie peatland complex (southern Quebec, Canada) was used to test the different approaches. Slug and bail tests were performed in piezometer standpipes to investigate catotelm hydraulic conductivity. Combined Darcy tests and tracer experiments were conducted on cubic samples using the modified cube method (MCM) to assess catotelm hydraulic conductivity, anisotropy and dispersivity. A new laboratory method is proposed for assessing acrotelm hydraulic conductivity and gravity drainage using a laboratory experimental tank. Most of slug tests' recovery curves were characteristic of compressible media, and important variability was observed depending on the initial head difference. The Darcy experiments on cubic samples provided reproducible results, and anisotropy (K_h > K_v) was observed for most of samples. All tracer experiments displayed asymmetrical breakthrough curves, suggesting the presence of retardation and/or dual porosity. Hydraulic conductivity estimates performed using the experimental tank showed K variations over a factor of 44 within the upper 40 cm of the acrotelm. The results demonstrate that the intrinsic variability associated with the different field and laboratory methods is small compared with the spatial variability of hydraulic parameters. It is suggested that a comprehensive assessment of peat hydrological properties can be obtained through the combined use of complementary field and laboratory investigations. Copyright © 2007 John Wiley & Sons, Ltd.     

黏弹介质中勒夫波频散问题的统一解及其动态特征分析

目前完全弹性介质中面波频散特征的研究已较为完善,多道面波分析技术(MASW)在近地表勘探领域也取得了较好的效果,但黏弹介质中面波的频散特征研究依然较少.本文基于解析函数零点求解技术,给出了完全弹性、常Q黏弹和Kelvin-Voigt黏弹层状介质中勒夫波频散特征方程的统一求解方法.对于每个待计算频率,首先根据传递矩阵理论得到勒夫波复频散函数及其偏导的解析递推式,然后在复相速度平面上利用矩形围道积分和牛顿恒等式将勒夫波频散特征复数方程的求根问题转化为等价的连带多项式求解问题,最后通过求解该连带多项式的零点得到多模式勒夫波频散曲线与衰减系数曲线.总结了地层速度随深度递增和夹低速层条件下勒夫波频散特征根在复相速度平面上的运动规律和差异.证明了频散曲线交叉现象在复相速度平面上表现为:随频率增加,某个模式特征根的移动轨迹跨越了另一个模式特征根所在的圆,并给出了这个圆的解析表达式.研究还表明,常Q黏弹地层中的基阶模式勒夫波衰减程度随频率近似线性增加,而Kelvin-Voigt黏弹地层中的基阶模式勒夫波衰减程度随频率近似指数增加,且所有模式总体衰减程度强于常Q黏弹地层中的情况.