基于人工源的甚低频电磁波空间传播特征统计分析

利用2005年1月至2010年11月DEMETER卫星记录的NWC发射站的VLF电场功率谱数据,采用指数拟合的方法,分析了VLF电磁波在卫星高度激发的电场空间分布和衰减特征.研究结果表明:(1)VLF电场在发射站上空及其磁共轭区有着很强的对应关系,存在南、北2个强电场中心涡;(2)相对于发射站的位置,VLF电场中心点具有经度和纬度偏移,日侧地磁经度偏移均值大于夜侧,而地磁纬度偏移均值则小于夜侧;(3)日侧VLF电场强度呈现出周期性的年变化;(4)在VLF电场中心10°范围内,电场强度随距离快速衰减,衰减常数b在长达6年的时间内保持稳定.在以上研究结果基础上初步构建的卫星高度人工源电磁波空间分布特征,将为研究地表-电离层电磁波传播机理提供基础技术支撑.     

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     

ZH-1卫星观测的VLF人工源信号特征分析与全波模拟

中国地震电磁监测试验卫星张衡一号(ZH-1)已于2018年2月2日成功发射,正在开展卫星数据在轨测试,并对卫星数据质量进行判定.本文对ZH-1卫星2018年5月至6月夜侧的VLF频段电场功率谱数据进行了分析.通过分析位于不同L值、具有不同发射频率的多个VLF人工源上空的卫星重访轨道观测数据,发现ZH-1卫星记录的人工源信号电场变化标准差与DEMETER卫星记录电场变化标准差几乎一致,说明ZH-1卫星观测数据具有较好的稳定性.通过重访轨道均值与全波模型计算结果对比,发现两者在数值上较为接近,在形态上较为一致,说明ZH-1卫星VLF频段电场功率谱数据具有一定的可靠性.此外,研究了VLF人工源上空及共轭区的电场分布特征和电波传播规律,并与DEMETER卫星的结果进行了对比,结果表明VLF人工源产生的电磁辐射穿透电离层后以导管或者非导管的哨声波模向共轭区传播,因为传播过程中的朗道阻尼,共轭区的电场能量比辐射源顶空更小.VLF人工源位于L1.5时,电磁波传播更容易发生非导管传播,VLF人工源信号导管传播模式在共轭区的电场响应相对于共轭点会发生一定程度北向偏移.     

低纬地区电离层电流的人工调制数值模拟

利用高频泵波能对低电离层进行有效的人工扰动.采用ELF/VLF调幅高频电波对电离层进行加热,电子温度会随着调制频率振荡,并引起电导率周期性变化,从而使加热区内电离层电流周期性变化,形成等效的ELF/VLF电离层虚拟天线,辐射调制频率范围内的无线电波.早期的电离层人工调制研究主要集中在高纬和极区,本文讨论低纬地区电离层人工调制的可能性.本文的理论研究和数值模拟结果表明,低纬地区低电离层电导率在周期性加热的条件下能有效地被调制,使加热区域形成ELF/VLF波的电流辐射源,并分析了不同加热参数和入射条件对调制效果的影响.     

Estimation of global lightning activity and observations of atmospheric electric field

Variations in the global atmospheric electric circuit are investigated using a wide range of globally spaced instruments observing VLF (∼10 kHz) waves, ELF (∼300 Hz) waves, Schumann resonances (4–60 Hz), and the atmospheric fair weather electric field. For the ELF/VLF observations, propagation effects are accounted for in a novel approach using established monthly averages of lightning location provided by the Lightning Image Sensor (LIS) and applying known frequency specific attenuation parameters for daytime/nighttime ELF/VLF propagation. Schumann resonances are analyzed using decomposition into propagating and standing waves in the Earth-ionosphere waveguide. Derived lightning activity is compared to existing global lightning detection networks and fair weather field observations. The results suggest that characteristics of lightning discharges vary by region and may have diverse effects upon the ionospheric potential.     

华北地区地电暴时GIC及涡旋电流响应分析

通过地电场台址近地表介质电阻率和地电场值计算大地电流,从大地电流场中分离出涡旋电流;根据平面波理论和水平导电层模型,使用地磁暴观测数据在频率域计算地磁感应电场(GIE),由GIE计算地磁暴感应电流(GIC).计算结果与实测值对比分析表明:GIE计算结果与实测地电暴具有很好的相位一致性;GIC涡旋中心相对地电暴涡旋中心存在向SE漂移约3°的现象;磁暴时地磁场Z分量的幅值分布图中极大值区域与涡旋中心重合,可能是GIC涡旋中心偏移的原因.另外,根据电磁感应原理提出的等效环电流模型,在一定程度上解释了涡旋大地电流的形成机制.本项工作可应用于地磁观测与地电观测的相互校正,同时有助于认识地电暴对大地电流分布的影响.     

NWC通信台在电离层中激发电磁响应的时变特征

本文利用DEMETER卫星VLF频段电场和磁场频谱数据对DEMETER卫星运行期间2005年至2009年澳大利亚甚低频(Very Low Frequency)通信台NWC发射的通信信号造成的电离层电磁响应的日变化、季节变化及年变化特征进行了统计分析,统计结果表明电磁响应日变化显著,夜间电场强度明显增强可达40dB,磁场变化略小也可为15dB左右,而季节变化不显著,年变化主要受太阳活动的影响,太阳活动越强,电磁响应越小.为解释数据分析结果,对地-电离层电磁波传播过程采用传递矩阵方法进行了模拟计算,模拟结果与数据分析的结果一致.我们认为这种随时间变化的特点可能由250km以下电离层电子密度分布特征导致,因此研究250km以下的电离层电子密度变化可能对寻找地震电离层电磁异常有重要意义.     

Comparison between VLF and RMT methods: A combined tool for mapping conductivity changes in the sedimentary cover

We compare the resolution power of single frequency very low frequency (VLF) electromagnetic data (real and imaginary parts of the tipper) and multi-frequency RadioMagnetoTelluric (RMT) data in delineating conductive structures typical for the sedimentary cover over crystalline basement in Scandinavia. Using VLF field data from five parallel profiles reveals that the estimated models have responses that fit the observed data well, and the models show an overall agreement with more detailed models derived from broadband RMT data. It is suggested that VLF data be used as a fast mapping tool to fill in the gaps between profiles along which more detailed RMT measurements are made. A generic model with conductive clay lenses and sandy formations over crystalline basement is used to generate synthetic data for the two cases. Using regularized inversion the corresponding estimated models clearly shows the strength and the weakness of both methods. Being inductive methods, they both have difficulties in clearly resolving the depth transition from conductive to resistive units. Especially the single frequency VLF data can be interpreted with very smooth models at depth. However, both methods resolve very well the lateral boundaries of the clay lenses and the RMT data also constrain the thickness of the clays quite well compared with the VLF models, which are less distinct at depth. Single frequency scalar VLF data emphasize those conductive structures that have dominant strikes in the direction of the transmitter. Multi-frequency VLF (tensor VLF) measurements provide the tipper vector which depends upon the underlying conductivity structure only. Real conductivity structures have significant 3D components which can be delineated easily by tensor VLF measurements. We propose that new VLF instrumentation be developed with this in mind.     

MAGNETIC WAVETILT MEASUREMENTS FOR GEOLOGICAL FRACTURE MAPPING1

The measurement of wavetilt is diagnostic for determining the electrical characteristics of the upper layers of the ground at VLF and LF frequency ranges. Theoretical and field studies have indicated that electric wavetilt using the transverse magnetic (TM) waves detects lateral inhomogeneities virtually instantly as abrupt changes in electrical properties are encountered. Theoretical studies have also indicated that magnetic wavetilt measurements using transverse electric (TE) waves are superior to electric wavetilts for such purposes. An experimental survey was conducted at two locations near Atikokan, Ontario, to verify the theoretical predictions. The survey area, forming a part of a large granitic pluton, was mapped earlier by various geophysical techniques, including the ground VLF-EM method, to detect weak conductors formed either by the presence of fractures in the bedrock filled with water and/or clay, or overburden filling bedrock depressions. A small, multi-turn, horizontal loop was used during the survey as the transmitter to generate TE waves at eleven frequencies from 10.7 to 58.5 kHz. The magnetic wavetilt measurements detected all previously known conductors at the two locations. In addition, the survey detected several weak conductors that were missed by the VLF survey. Thus, the survey indicated the usefulness of magnetic wavetilt results for detection of weak conductors at shallow depths, which may have application in engineering geophysical surveys. The multi-frequency wavetilt data also provided some indications of the depth and depth extent of such conductors.     

Airborne VLF Measurements and Variations of Ground Conductivity: A Tutorial

Airborne very low frequency (VLF) data are routinely collected by national agencies and commercial companies together with other passive geophysical measurements of the static magnetic field and radiometric data. The purpose of this paper is to demonstrate that both standard three-component VLF and tensor VLF (TVLF) data contain a lot of useful quantitative and qualitative information about the electrical conductivity distribution in the upper few hundred meters of the crystalline basement. We first give a new derivation of the fundamental transfer functions (the tipper) used in the TVLF technique. We then show that the tipper can be estimated from simultaneous measurements of the wave magnetic fields from at least two transmitters with somewhat different frequencies, and present a simple model by which the maximum error introduced by the difference in frequencies can be found. Single transmitter scalar VLF maps emphasise those conductive structures that have dominant strikes in the direction of the transmitter. Multiple transmitter transfer functions are dependent only upon the underlying conductivity structure. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic magnetotelluric (MT) soundings. In such cases three-component VLF measurements can be modelled easily upon appropriate rotation of the co-ordinate system to “strike” co-ordinates. Single frequency transfer functions (tippers) have real and imaginary parts that carry information on not only lateral contrasts in conductivity, as usually stated in text books, but, taken together, they provide a robust tool for determining the background conductivity level away from distinct conductors, and they can also be used to discriminate between deep and shallow conductors. Based upon simulations using multi-frequency data, it can be concluded that such a new development would dramatically increase the resolving power of airborne VLF measurements.     

LINEAR FILTERING OF VLF DIP-ANGLE MEASUREMENTS*

The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data. Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined. The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.     

Interpretation of ground VLF-EM data in terms of inclined sheet-like conductor models

The response of two-dimensional, inclined, sheet-like conductors with, low conductance values to plane wave electromagnetic fields in the very low frequency (VLF) range has been evaluated by using a numerical technique. The conductance values of the conductors considered are appropriate for those produced by water and/or clay-filled fracture and shear zones in the Precambrian crystalline rocks of the Canadian Shield. The surrounding host rock was assumed to be, resistive with resistivities in the 1–10 k.m range to reflect the high resistivities over the shield areas. No overburden was assumed in this analysis.The results of the computations are presented in the form of characteristic interpretation diagrams to interpret ground VLF data in the field, where facilities for direct numerical modelling may not be available. A method for interpreting ground VLF data using such characteristic diagrams has been proposed in this paper which requires a prior knowledge of the host rock resistivity and the inclination of the conductor. These two parameters may be derived from a VLF resistivity survey and from appropriate filtering of the VLF tilt angle response. The interpretation method was applied to a ground VLF anomaly obtained at a research site near Atikokan in NW Ontario, which yielded an interpretation compatible with information from geological mapping.Geological Survey of Canada Contribution No. 51888.     

地球磁层中哨声导管的物理结构

本文在理论上探讨了地球磁层中哨声导管的物理特征。利用两种能够引导甚低频电磁波的导管模式,求得哨声导管中磁场和电流的分布。当磁场偏离势场或无力场时,将会产生场向电流,该电流与导管的强度和稳定维持密切相关。根据这些结果,我们认为高纬磁层中观测到的较强的场向电流,是高纬地区地面台站频繁接收到哨声及其回波的一个重要原因。     

POTENTIAL FIELD OF A STATIONARY ELECTRIC CURRENT USING FREDHOLM'S INTEGRAL EQUATIONS OF THE SECOND KIND*

An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite.     

电离层舒曼共振的空间观测

通过功率谱分析和波阻抗函数计算,本文证实了Aureol 3卫星在电离层高度上(>600km)观测到的极低频(ELF)波场扰动是和舒曼共振相关的电磁振荡.与舒曼共振地面观测相比较,Aureol 3观测到的舒曼共振电场分量具有很好的谐振谱结构,峰值频率和各阶舒曼共振本征频率对应;磁场分量的高阶峰值频率偏离14, 20, 26Hz等舒曼共振本征频率;随着卫星高度的改变,电场与磁场谐振的一阶最大能量峰值并不会发生在同一频率,结合本文分析的数据,分别位于78Hz和10Hz;水平方向的磁场分量更接近南北方向的线极化而不是地球－电离层空腔中的椭圆极化;波阻抗随频率表现出不太规则的准正弦振荡,它会随着频率增加和飞行高度上升呈现减小的趋势.虽然舒曼共振信号和电离层密度梯度间的非线性作用可以解释舒曼共振空间观测的部分特征,但需加入其他机制,如电离层不稳定性,传播模式的耦合,进一步了解电离层高度上舒曼共振各种特征产生的原因.     

Integrated Geophysical Studies in the East-Indian Geothermal Province

Integrated geophysical surveys using vertical electrical sounding (VES), very low frequency (VLF) EM, radiation counting, total magnetic field and self-potential (SP) measurements are carried out to characterize the geothermal area around a hot spring in the Nayagarh district, Orissa, India that lies in the East Indian geothermal province. The study was performed to delineate the fracture pattern, contaminated groundwater movement and possible heating source. VES interpretations suggest a three- to four-layer structure in the area. Resistivity survey near the hot spring suggests that weathered and fractured formations constitute the main aquifer system and extend to 60 m depth. Current flow measured at various electrode separations normalized by the applied voltage suggests that fractures extend to a greater depth. Detailed VLF study shows that fractures extend beyond 70 m depth. VLF anomaly has also very good correlation with the total magnetic field measured along the same profiles. Study results suggest that a gridded pattern of VLF survey could map the underground conductive fracture zones that can identify the movement of contaminated groundwater flow. Therefore, precautionary measures can be taken to check further contamination by delineating subsurface conducting structures. Self potential (SP) measured over the hot spring does not show a large anomaly in favor of the presence of a sulphide mineral body. A small positive (5–15mV) SP anomaly is measured which may be streaming potential due to subsurface fluid flow. A high radiation is measured about four kilometers from the hot spring, suggesting possible radiogenic heating. However, the exact nature of the heating source and its depth is not known in the area. Deep resistivity followed by a magneto-telluric survey could reveal the deeper structures.     

Review of Electric and Magnetic Fields Accompanying Seismic and Volcanic Activity

New observations of magnetic, electric and electromagnetic field variations, possibly related to recent volcanic and seismic events, have been obtained on Mt. Unzen in Japan, Reunion Island in Indian Ocean, the Long Valley volcanic caldera in California, and for faults in China and Russia, California and several other locations. For volcanic events, contributions from different physical processes can be identified during the various eruption stages. Slow processes (weeks to months) include near-surface thermal demagnetization effects, piezomagnetic effects, and effects from rotation/displacement of magnetized material. Rapid processes (seconds to days) include piezomagnetic effects from instantaneous stress redistribution with explosive eruptions and electrokinetic effects from rupture of high pressure fluid compartments commonly encountered in volcanic regions. For seismic events, the observed coseismic offsets are instantaneous, provided care has been taken to ensure sensors are insensitive to seismic shaking and are in regions of low magnetic field gradient. Simple piezomagnetic dislocation models based on geodetically and seismically determined fault parameters generally match the observed signals in size and sign. Electrokinetic effects resulting from rupture of fluid filled compartments at hydrostatic to lithostatic pore pressures can generate transient signals in the frequency band 100 Hz to 0.01 Hz. However, large-scale fluid driven processes are not evident in near-field measurements in the epicentral region minutes to weeks before large earthquakes. The subset of ionospheric disturbances generated by trapped atmospheric pressure waves (also termed gravity waves and/or acoustic waves, traveling ionospheric disturbances or TID's) that are excited by earthquakes and volcanic eruptions are common and propagate to great distances. These are known and expected consequences of earthquakes, volcanic explosions (and other atmospheric disturbances), that must be identified and their effects removed from VLF/ULF electromagnetic field records before associating new observations of ionospheric disturbances with earthquake activity.     

多辐射场源地空瞬变电磁法多分量全域视电阻率定义

本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势.     

不同方位甚低频电磁波衰减的初步研究*

利用2005—2010年DEMETER卫星记录NWC发射站的19.8kHz电场功率谱数据,采用统计和线性拟合等方法,研究了NWC站发射的电磁波在顶部电离层及磁共轭区激发的电场效应及其在不同方位的衰减特性。结果发现:(1)在NWC发射站上空,卫星记录电场呈椭圆状扩散分布,电场最强中心点位置相对地面人工源位置有一定偏移;(2)研究区上空电场北部衰减梯度最小,东部衰减梯度最大;(3)在其磁共轭区,南部电场衰减梯度最小,北部衰减梯度最大。综合分析认为人工源(NWC)发射的VLF电磁波传播到电离层高度后,受地磁场影响,电磁波主要沿磁力线方向传播,空间电场最强中心点位置相对发射站的位置发生偏移,向磁赤道方向倾斜,而且偏于磁赤道方向的电场衰减梯度最小。