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.     

Trefftz method for the solution of three-dimensional forward and inverse problems of geoelectrics

The Trefftz method is applied to the numerical solution of the three-dimensional (3D) forward problem for an electromagnetic field harmonically oscillating with time in a 3D environment, which ensures the possibility to solve a 3D inverse problem. Here, the known idea of simultaneous joint interpretation of the tangential components of electric and magnetic fields measured on the Earth's surface, which was suggested by A.N. Tikhonov for the development of the electromagnetic sounding method, is employed in the context of the numerical-analytical representation of the forward problem solution in accordance with the Trefftz method. Application of this method allows enables one to generalize the basic plane-stratified vertically 1D model of the medium by allowing for lateral variations in electric conductivity.     

A numerical study of the topographic effect on electromagnetic fields in a three-dimensional conductivity model

Summary The topographic effect on electromagnetic fields in a three-dimensional conductivity model is investigated using a finite-difference numerical method. The model is that of a conducting block buried in the host Earth beneath a conducting hill. The altitude effect as well as the effect of conductivity of the hill has been studied by comparing the results for the field components at the surface of the hill with those at the surface of a flat half-space. The results indicate that both altitude and the hill conductivity influence the behaviour of the electromagnetic fields. For certain traverses beyond the edge of the block, it is seen that the variations in some field components are mainly due to the topographic effect and that the contribution from the block itself is negligible.     

含起伏地形地层中点源激发的震电响应研究

地震诱导电磁现象是国内外地学领域十分关注的前沿问题,前人对地震波和电磁场耦合波场的认识主要是基于规则模型获得的.为研究含起伏地表和地下界面的地层中震电波场激发、传播特性,本文采用有限元软件COMSOL Multiphysics模拟点震源激发的电磁场.首先给出频率域二维SHTE模式震电耦合方程组,然后利用COMSOL软件建立计算模型,并求解出点力源激发震电波场的频率域响应,最后利用FFT变换得到地震波场和电磁场的时间域波形.模拟结果表明,震电波场中存在三种类型的电磁信号,第一种是震源直接激发的电磁波;第二种是地震波在分界面处激发的电磁波(包括自由表面、地下不同介质分界面);第三种是伴随地震波的同震信号,前两种电磁波比地震波更早到达远处观测台站,对地震预警有重要意义.此外,研究还发现:当地震波传播至地表并沿着地表传播时,在地表附近空气层中同样记录到了伴随地震波传播的电磁扰动信号,该信号与相同水平源距条件下、地下观测点接收到的电磁信号相同,这与前人的一些观测结果相符.本文研究结果为今后地震电磁信号的解释提供了理论证据.     

Electromagnetic induction of a three-dimensional conductivity inhomogeneity in a two-layer earth

Summary The results of numerical computations of the electromagnetic field induced in a two-layer model of the Earth with a three-dimensional inhomogeneity (a block) in the subsurface layer are given. Several recommendations are given which have enabled the solution of a complicated system of integral equations and the computation of the field at the Earth's surface in an effective way. The analysis of the obtained solution has proved that, in an anomalous electromagnetic field, the field of the horizontal electric dipole, oriented in the direction of the exciting field, is predominant. A number of practically usable diagrams and approximative formulae is given.     

Controlled‐source electromagnetic monitoring of reservoir oil saturation using a novel borehole‐to‐surface configuration

To advance and optimize secondary and tertiary oil recovery techniques, it is essential to know the areal propagation and distribution of the injected fluids in the subsurface. We investigate the applicability of controlled‐source electromagnetic methods to monitor fluid movements in a German oilfield (Bockstedt, onshore Northwest Germany) as injected brines (highly saline formation water) have much lower electrical resistivity than the oil within the reservoir. The main focus of this study is on controlled‐source electromagnetic simulations to test the sensitivity of various source–receiver configurations. The background model for the simulations is based on two‐dimensional inversion of magnetotelluric data gathered across the oil field and calibrated with resistivity logs. Three‐dimensional modelling results suggest that controlled‐source electromagnetic methods are sensitive to resistivity changes at reservoir depths, but the effect is difficult to resolve with surface measurements only. Resolution increases significantly if sensors or transmitters can be placed in observation wells closer to the reservoir. In particular, observation of the vertical electric field component in shallow boreholes and/or use of source configurations consisting of combinations of vertical and horizontal dipoles are promising. Preliminary results from a borehole‐to‐surface controlled‐source electromagnetic field survey carried out in spring 2014 are in good agreement with the modelling studies.     

Electromagnetic response of a buried cylindrical conductor for sheet and line current sources

A major problem in electromagnetic induction studies in regions of localized source fields, such as the auroral and equatorial electrojet regions, is the source effect. Using an analytical model, the electromagnetic response of a buried conducting cylinder to sheet current and line current excitations has been studied for the period rangeT=5 s to 24 h. The validity of the numerical results obtained from the analytical model are compared with the numerical results obtained from a finite difference model. The results show that for periods less than 30 min, there is no significant difference in the response of the cylinder to both source fields. However, significant differences are observed at longer periods. It was also observed that the equivalent height at which a uniform sheet current at 100 km above the earth's surface can be approximated by a line current varies as a function of the source period.     

Asthenospheric parameters of dipping plate regions by magnetotelluric investigations

The multifrequency resolving capability of the electromagnetic surface impedance parameter, employed in magnetotelluric investigations, has been examined for an upper mantle model incorporating (1) a deep-seated lithosphere-asthenosphere boundary and (2) preferred orientations of olivine crystals in the lithosphere, usually associated with the subducting plate boundary zones. Numerical results display the quantitative errors in the resolution of various anisotropic and dipping plate conditions. It is found that (1) conductivity anisotropy, particularly with higher dips, causes a significant dispersion in the impedance values and higher resolving power, and (2) variations in the geometrical/physical parameters are reflected predominantly in the phase component of the surface impedance. The study has relevance to the understanding of the electrical nature of materials, its relationship with physical properties and associated geothermal and seismic activities in the dipping plate regions.     

南极普里兹湾海冰厚度的电磁感应探测方法研究

本文介绍一种由EM31-ICE型电磁感应仪和激光测距仪组合而成的船载电磁感应海冰厚度探测系统. 针对海冰和海水的电学特征,运用电磁感应技术提取探测系统至海冰下底面的距离,运用激光测距仪测量冰面粗糙度和探测系统至海冰上表面的距离,两组数据结合,实现了海冰厚度的探测. 通过南极现场探测数据分析,并与钻孔实测冰厚数据对比研究,定量分析了探测系统距离冰面的高度效应,建立了该系统冰厚测定值随高度变化的修正关系式,并对船载航行数据进行了系统校正. 与SCAR ASPeCt的冰厚数据对比分析,表明该系统能够获得可靠的海冰厚度并具有较高的精度,且能满足对极区大范围海冰厚度观测的需求.     

可控源电磁三维频率域有限元模拟

本文采用电磁场的磁矢量位和电标量势,将Maxwell方程组化为位势的类似于Helmholtz型方程,并引入罚项及稳定化方法克服了电磁三维有限元算法中的伪解及数值不稳定性;采用人工边界把计算区域局域化,将均匀半空间中水平电偶极子源产生的位势值作为人工边界上的第一类边界条件以表示源的作用,减少了实际的计算区域.理论模型和复杂模型的计算结果均表明,可控源电磁三维有限元数值模拟给出了稳定、可靠的电磁场分布.     

Ulf-wave diagnostics of a mid geomagnetic latitude model of the ionosphere

Summary Wave diagnostics of several alternatives of a constructed model of the daytime and night ionosphere of the mid geomagnetic latitudes were carried out in the ULF range of electromagnetic waves. The altitude profiles of local wave parameters of the ionosphere, the wavelengths and wave attenuation, were analysed at a fixed frequency. The frequency dependence of the amplitude and polarization characteristics of the total wave at the Earth's surface were also analysed in dependence on the frequency-amplitude spectrum of geomagnetic pulsations. The effect of the low regions of the ionosphere (h<100 km) was tested on the damping of the propagating wave and on the frequency-amplitude and polarization characteristics of the wave at the Earth's surface.     

植被电磁散射的半空间模型研究

本文给出了地面植被电磁散射的半空间模型研究方法.在以往的相关文献中,均采用自由空间格林函数求解单个叶片散射体的散射场,本文利用半空间并矢格林函数求出了单个散射体的散射场,然后运用Monte Carlo方法模拟生成地面植被层,得到了在半空间下植被层的电磁散射特性,计算出了单、双站雷达散射截面,并与自由空间下的散射场做了相互对比.结果表明,在半空间格林函数下得到的散射场更为适用于描述地面植被的电磁散射特性.     

电离层-均匀地球模型中地表水平电偶极子激发的电磁场

本文在电离层-均匀地球模型中使用Debye势方法,在电磁场一般性边界条件下导出了地表水平电偶极子源在地球介质中激发的电磁场的表达式.对球Bessel函数使用指数函数近似,并采用加速算法,计算了地表及地下的电磁场.其正确性得到理论计算和实际观测数据的多方面支持.结果表明,在收发距离大于1000 km,电磁场强度很稳定,几乎没有衰减的迹象;电磁场受电离层的影响很大,由于电离层高度升高或电导率减小引起的电磁场幅度的减少可以达到50%,但是由电场和磁场得到的视电阻率变化却很小.     

时变海面与其上方导体平板的复合电磁散射研究

应用当前的数值方法求解动态海面与目标的复合电磁散射,由于海面的变化,在不同时刻需要对阻抗矩阵各元素及海面表面电流重新求解,因而要耗费大量内存和运算量.为了克服这一问题,本文应用物理光学(PO)近似求解了导体海面表面电流及导体平板的一阶散射场,应用基尔霍夫近似给出了海面的后向散射场,同时借助互易性定理降低了求解平板和海面之间二次耦合散射场的难度,讨论了平板尺寸、风速等对后向复合散射场的影响.另外,本文还推导出了耦合散射场多普勒谱频移的理论公式,详细分析了复合后向散射场的Doppler 频谱特性.     

Three‐dimensional controlled‐source electromagnetic modelling with a well casing as a grounded source: a hybrid method of moments and finite element scheme

Steel well casings in or near a hydrocarbon reservoir can be used as source electrodes in time‐lapse monitoring using grounded line electromagnetic methods. A requisite component of carrying out such monitoring is the capability to numerically model the electromagnetic response of a set of source electrodes of finite length. We present a modelling algorithm using the finite‐element method for calculating the electromagnetic response of a three‐dimensional conductivity model excited using a vertical steel‐cased borehole as a source. The method is based on a combination of the method of moments and the Coulomb‐gauged primary–secondary potential formulation. Using the method of moments, we obtain the primary field in a half‐space due to an energized vertical steel casing by dividing the casing into a set of segments, each assumed to carry a piecewise constant alternating current density. The primary field is then substituted into the primary–secondary potential finite‐element formulation of the three‐dimensional problem to obtain the secondary field. To validate the algorithm, we compare our numerical results with: (i) the analytical solution for an infinite length casing in a whole space, excited by a line source, and (ii) a three‐layered Earth model without a casing. The agreement between the numerical and analytical solutions demonstrates the effectiveness of our algorithm. As an illustration, we also present the time‐lapse electromagnetic response of a synthetic model representing a gas reservoir undergoing water flooding.     

层状介质高频电磁场计算方法及结果分析

本文在研究了高密度采样线性滤波算法计算层状介质高频电磁场的电磁响应的基础上,系统比较了Anderson等数值计算方法的速度和精度,给出了不同的地电参数和接收距对该计算方法精度的影响.通过比较表明,高密度采样算法不仅速度快,而且精度高.本文还对水平回线源层状介质当发射线圈位于地面时的表达式的第二项积分进行了公式推导,为电磁场的计算提供了清晰的思路.     

基于共形辛Euler算法的非金属地下管道精细化高效探地雷达正演模型

我国城市地下管线种类繁多、规模巨大.近年来,因城市管网家底不清、位置不明导致的管道挖断事故频发.探地雷达(GPR)是目前最为有效的地下设施探测定位设备,通过对探地雷达回波信号进行反演分析,可获取管道埋深、位置等信息.然而,地下管道等地下圆形目标,传统正演模型中阶梯近似方法会产生一定误差.本文结合辛Euler算法与曲面共形技术,建立探地雷达电磁波在含非金属管道地下结构中传播的精细化数值模型,并通过图形处理器(GPU)加速技术进一步提升模型计算效率.数值模拟结果显示,共形网格技术有效降低了由于阶梯近似造成的虚假绕射波,并通过结合GPU并行计算,在不增加网格密度的前提下,实现了对地下管道电磁响应的高效精细化计算.     

地-井瞬变电磁响应特征数值模拟分析

井中瞬变电磁法(Bore-hole transient electromagnetic method-BHTEM)是指接收线圈在钻井中观测瞬变场响应用以勘查深部矿产资源的勘探方法,其中以地-井(地面激发井中接收)组合方式研究最多、应用最广.本文应用时域有限差分法(FDTD),建立包含薄板导体的均匀半空间二维数学模型,采用线源为激发源,选用Mur吸收边界条件,对矩形回线源在半空间中产生的瞬变电磁场进行数值模拟,计算了低阻板状导体在均质半空间和有低阻覆盖层影响情况下的地-井瞬变电磁异常响应,并对响应的特征及规律进行研究分析,为研究地-井TEM提供参考.     

2.5维起伏地表条件下频率域地空电磁正演模拟

频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特性进行了对比.     

坑(井)-地多参数电磁接收系统

坑（井）-地多参数电磁接收系统在地面、坑道、井中三维空间观测天然场源及人工源电磁信号,观测装置接近或穿过矿体,借助动源、大功率发射可获得更加明显的异常值,旨在加大探测深度的基础上提高深部分辨率,同时获取电阻率、极化率、复阻抗等参数.为坑（井）-地电磁成像方法现场数据采集提供硬件设备支持,针对坑道、井中电磁观测的特殊需求,开发了由不极化电极、非接触电极、小型感应式线圈、微型三分量音频磁传感器、地面-坑道电磁接收机、井中电磁接收机等组成的坑（井）-地多参数电磁接收系统,解决了传统不极化电极坑道硬岩表面接地困难和传统感应式磁传感器难以适应坑道狭窄空间的问题;通过多层次的室内、弱干扰条件下野外测试验证了接收系统的各项技术指标,接收机通道全频段噪声水平接近10 nV/rt（Hz）,微型三分量音频磁传感器体积压缩至32 cm×32 cm×32 cm;强干扰条件下的矿山试验验证了接收系统的可靠性、适用性及先进性.