长短极距比值干扰源定位方法及其在地电场观测中的实际应用

提出长短极距比值干扰源定位方法用来定位地电观测中出现的地表干扰源。该方法基于均匀半空间模型且地表干扰为点电源,利用地表不同位置的干扰源所引起的同测向长短极距的电场变化差异,联合多测向比值来定位干扰源位置,并应用到徐庄子台漏电干扰一例和高邮台偶极源试验中。结果表明:利用长短极距比值法,能有效定位干扰点源,且能定位偶极干扰源的其中一个电极。     

水平层状大地中任意方向偶极源激发的电磁场

目前已存在多种物理模型用于解释地震前兆电磁异常,不同机制下产生电磁波场的源可等效为不同类型偶极源(磁型/电型)的组合.针对不同偶极源激发的电磁响应特征的研究对地震孕育过程中观测的电磁信号的提取与解释具有重要意义.基于此,本文推导了一套正演算法,可用于研究任意方向(方位角和倾角)、大小和位置的一个或多个偶极源在地壳或岩石圈-大气圈-电离层(Lithosphere-Atmosphere-Ionosphere, LAI)模型中产生的电磁波场的传播特征.基于多源组合,该算法可进一步拓展到线电流源和面电流源的电磁响应的计算.本文给出了水平、垂直和倾斜电偶极源响应计算的实例,并针对不同类型的偶极源激发下的地中、地面以及空间电磁场响应特征进行了初步分析.模拟结果显示,相比于单纯水平或垂直的电偶极源,地中倾斜电偶极源激发的电磁响应在地表的各场量辐射花样虽然仍可反映偶极子的方位角和空间位置,但不再具有对称性,其模式更复杂.本文提出的正演算法可用于地中低频电磁辐射源的正演模拟和反演定位.     

点源干扰源定位方法及其在地电场观测中的实际应用

我国地震地电场观测时的漏电干扰问题会影响到观测质量和地震监测效能。针对这一实际问题,本文提出了点源干扰源定位方法。该方法首先建立地表点源对地电场观测影响模型,根据模型计算点源对地电场观测的干扰幅度,然后求解干扰源相对于中心电极的位置。将定位方法应用于高邮地震台和汉王地震台实际干扰源排查工作中,定位结果与实际干扰源位置基本一致,测区外干扰源定位误差小于2%,测区内干扰源定位误差小于20 m,并对影响计算误差的因素进行了分析。两个台站的实际应用表明：本文提出的定位方法为地电观测实践中查找点源漏电干扰源提供解决方案,能够有效地定位点源干扰源位置,提高干扰排查效率。     

基于牛顿迭代法的地电场观测人为漏电干扰源反演及其应用

地电场异常变化是地震前兆观测中的一个重要研究方向.我国地震台站地电场观测资料经常受到各类干扰的影响,人为漏电亦是一种常见干扰.本工作通过构建漏电偶极源二维地电场模型,正演偶极源对观测电极产生的电位差比例系数,采用牛顿迭代算法,实现漏电位置的精确求解.在模型基础上,进行模拟位置漏电及台站实际漏电实验,分别对反演算法进行测试.测试结果表明:基于牛顿迭代的地电场人为漏电反演算法,无论是在模拟实验还是台站的实测检验中,均能实现对漏电干扰源的有效定位,且算法收敛快速稳定.地电场人为漏电干扰源定位可执行程序可通过邮件向作者获取.     

水平两层均匀介质中井下电阻率观测信噪比的理论计算

在点电源水平两层均匀介质模型下计算了在不同地电断面中观测时地表干扰电流源对观测的影响。得到:地表干扰电流源对地电阻率观测电势差的影响取决于地电断面类型和参数、供电电极和测量电极的埋深以及避开干扰源的距离。本文研究结果对实施井下电阻率观测中台址电性结构选择、电极埋深、干扰源避让有参考意义。     

井下地电阻率观测中地表电流干扰影响计算

本文把地电阻率台站地下介质简化为水平三层均匀介质模型,以点电流模拟地表干扰电流,针对对称四极观测装置,计算了在不同电性结构中的不同深度观测时,地表电流干扰源对对称四极装置地电阻率观测的影响,计算得到:地表干扰源对电阻率观测的影响取决于电性结构的类型和层参数、供电电极和测量电极的埋深以及避开干扰源的距离.本文研究结果对实施井下电阻率观测中台址电性结构的选择、电极埋深、干扰源避让距离等有参考意义.     

地表点电流对地电阻率观测的影响

将地电阻率台站地下介质简化为水平层状均匀介质模型,以点电流模拟地表干扰电流源,针对对称四极观测装置,计算干扰源电流增减及位置变化对地电阻率观测的影响。结果表明,地电阻率的受干扰程度取决于地表干扰点电流的大小、在测区的方位及距中心点的距离。     

地磁场长期变化特征及机理分析

将地磁场的总变化分为三部分：偶极场自身变化，非偶极场自身变化及非偶极场磁斑区通过对核幔边界（CMB）层环形电流的调制来影响偶极场的变化. 本文利用国际地磁参考场模型IGRF）1900～2000计算分析了地球不同深度地磁场分布及长期变化特征，且讨论了变化的可能机制. 可以推论，地磁场西漂和倒转不仅是非偶极场引起，同时与偶极场有密切关系.     

大功率长偶极与环状电流源电磁波响应特征对比

为了在深部第二成矿带找矿,需要进一步提高电磁法可控源的强度,因此,一种新的大功率固定发射源的电磁法正在孕育之中.该方法采用长偶极源或者环状电流源,为了提高发射源的强度,偶极长度L或者环状源半径R需要非常大.然而,电磁场信号的强度不仅仅依赖于发射源,当L或者R为数十公里时,接收点的位置远离发射源,电离层对电磁场的影响不可以忽略.电离层的影响既可能是增强信号,也可能是削弱信号,这依赖于不同类型源产生的电磁波.本文利用三维积分方程法对包含电离层、空气层、固体地球层(简称地-电离层模型)的线状和环状大功率电流源电磁波响应进行了三维数值模拟,对两种源的电磁波的响应特征进行了对比研究.研究表明虽然磁性源的发射效率高于电性源,但是由于发射源在波导区产生的直达波经电离层的反射波要产生相长相消干涉,波导区的场特征既和发射源的特征有关,也和电离层有关.对比研究发现对于x方向的水平长偶极源电离层对电磁场 E_x 和 H_y 的影响是使场的幅值增强,即波导场的 E_x 和 H_y 的衰减是减缓的.而对于水平环状电流源,电离层对水平方向的电磁场的影响是使场的幅值减弱,即波导场的 E_x 和 H_y 的衰减是增大的.为此,认为采用大功率固定源在波导场工作时应该采用线状电流源.     

基于Weimer模型的高纬地磁扰动预测

地磁扰动是空间天气中的重要现象,对地基技术系统具有重要的影响.准确预报地磁扰动可以有效避免重大灾害发生.本文基于Weimer电势和磁势模型发展了高纬地区地磁扰动的模拟方法,并与地面台站观测数据进行了比较.地表磁场扰动主要受电离层电流系统的影响,利用Weimer模式计算出电离层等效电流分布后,基于毕奥-萨伐尔定律推导了地磁扰动三分量与电流的关系,最终计算出地磁扰动量.模型的输入参数为太阳风速度、太阳风密度、行星际磁场和磁偶极倾角.模型计算结果与不同纬度和经度的地磁台站观测结果对比表明本文的计算方法能有效地模拟地磁暴期间地磁扰动特征.本文结果对今后发展高纬地区地磁场预报模型奠定了重要基础.     

EFFICIENT COMPUTATION OF TRANSIENT SOUNDING CURVES FOR WIRE SEGMENTS OF FINITE LENGTH USING AN EQUIVALENT DIPOLE APPROXIMATION1

Calculations of emf induced in a vertical axis coil by a current step in a straight wire segment of finite length were carried out using an equivalent dipole placed along the source wire. An approximation is valid for homogeneous and layered earth models. The location of the equivalent dipole is calculated by forcing the early- and late-time asymptotes of the transient sounding curve of the equivalent dipole to match those of the finite wire. This approximation works because the early-time asymptote of the emf depends on the component of the receiver position which is parallel to the wire while the late-time asymptote does not. Analytical integration of the early-time asymptote yields an expression for the equivalent dipole location in closed form. The coil can be placed anywhere in the vicinity of the finite wire. Square or rectangular loop sources can be simulated by one or more finite wire segments depending on the source-receiver geometry. The equivalent radius calculated for central loop soundings agrees well with the value derived using a circular loop with the same area as the square loop. Results show that acceptable sounding curves can be generated by the equivalent dipole for coils placed as close as 0.5 source lengths from the finite wire segment. Higher accuracy can be obtained by splitting the finite wire into two or more subsegments. Results for layered models are slightly better than homogeneous earth models when the resistivity increases with depth and slightly worse for models with resistivity decreasing with depth. Approximate calculations are about 10 times less expensive than exact calculations depending on the method used for the numerical integration.     

Global anomalies of the main geomagnetic field and the dynamic model of their sources

The model of the main geomagnetic field sources, including 12 more dipoles in addition to the main dipole, is presented. The source parameters have been determined for the 100-year-long period at an interval of 5 years. The problem is solved using the iteration method proposed and developed by the authors. Six more sources have been included in the iteration process at its current step. This made it possible to substantially smooth the time variations in the dipole parameters obtained previously and, on the other hand, to describe all global anomalies of the main geomagnetic field using this model. A change in the parameter values can be approximated by rather smooth functions describing specific features of development of the sources operating for 20–30 years. However, the short duration of the obtained series (21 points) makes the problem of prediction substantially more difficult. The parameters of 13 dipoles, extrapolated for 2000 using the linear prediction filter, are in good agreement only with the general development of the sources.     

复杂场源形态的海洋可控源电磁三维正演

在使用电偶极发射源的可控源电磁法(CSEM)勘探中,发射源的方位、长度、形状等对观测数据有重要的影响,然而现有的大部分三维数值模拟方法没有全面地将这些因素考虑进来,很多都只能应对非常简单的场源形态,例如单一方位的点电偶极子,这有可能显著降低模拟结果的准确性.本文实现了基于交错网格有限体积(FV)离散的海洋CSEM三维正演算法,能够模拟形态相对复杂的场源,包括任意方位的有限长直导线和弯曲导线发射源.该算法使用一次场/二次场方法,只需对二次场使用FV法求解,避免了场源的奇异性问题;一次场的计算为一维正演问题,使用准解析法求解,并且只要在计算一次场时考虑复杂的场源形态便可以实现同样场源的三维正演.通过与一维理论模型的解析解对比验证了三维程序的准确性,并针对三维理论模型进行了一系列正演测试,初步考察了场源形态对三维正演结果的影响.     

多极源随钻声波测井实验分析

本文针对随钻声波测井中钻铤波干扰地层声波测量的问题,设计了小尺寸随钻声波测井探头,在实验室内开展了多极源(单极源、偶极源和四极源)随钻声波测井实验研究.先在水池中对裸露的探头进行了随钻测量,记录到了沿钻铤传播的直达钻铤波,并获得了不同声源激发的钻铤波速度.之后在砂岩和有机玻璃模型中进行了随钻声波测井实验,观测到多极源随钻测井的全波波形,并通过与水池中实验结果的对比,分析了井中钻铤波的传播特性.特别地,在偶极随钻测井实验中不仅记录到了偶极一阶钻铤波,还观测到对弯曲波干扰较大的偶极二阶钻铤波,进而研究了不同声源频率下该波群的响应特性及其对弯曲波测量的影响.此外,基于单极源和偶极源随钻声波测井实验数据,本文发现:随着声源频率的增加,单极和偶极钻铤波的传播特性不同,但它们在测井全波中的相对幅度均降低,进而可从测井全波中较好地提取地层的声波速度.本文实验结果对随钻声波测井仪器设计及测井数据解释具有重要意义.     

一种地电阻率场地干扰的快速检测方法

对2018年洛阳地震台地电阻率观测区域出现的1次典型漏电干扰进行分析。经初步分析、干扰源模拟定位、干扰源现场调查等,确定干扰源的位置。通过对此次干扰的分析与排查,归纳总结出一种便利、快速的干扰检测定位方法。     

层状垂直各向异性介质海洋CSEM发射源及接收站姿态和位置对电磁响应影响研究

众所周知,电阻率各向异性对海洋可控源电磁(CSEM)响应有显著影响.在模拟和解释海洋可控源电磁资料时,通常假设接收站布放于平稳的海底,发射源以理想的水平电偶极源形式按预设路径拖曳前进.然而,在实际海洋作业中,自由下沉的接收站可能沉降在局部倾斜的海底面上,在海底的方位可能是任意取向的.同时,由于洋流等因素的影响,发射源偶极子会发生旋转、倾斜,其位置也会偏离预设的拖曳路径,所以发射源和接收站姿态的变化都将对电磁资料产生影响.本文计算了电阻率各向异性介质中,发射源任意姿态任意位置激发、接收站任意姿态任意位置接收的海洋CSEM响应,分析观测系统引起数据误差的机制,分别在电阻率各向同性、围岩电阻率各向异性和高阻储层电阻率各向异性三种模型中,讨论发射源和接收站的姿态及位置变化对响应的影响.计算结果表明,电磁场振幅误差主要来源于发射源和接收站的倾斜角,而相位的变化为发射源位置和接收站倾斜角、位置共同影响的结果.     

一种改进的基于网格搜索的微地震震源定位方法

震源定位是微地震监测技术要解决的主要问题.目前,井下微地震监测多采用走时拟合法计算震源位置.常规方法受到环境噪声、初至拾取误差、速度模型误差等因素的影响,定位结果存在一定误差.为了提高定位精度,本文提出了一种改进的基于网格搜索的微地震震源定位方法.本文方法根据P波的偏振特征参数计算概率密度函数求取震源方位角,并采用改进的目标函数和搜索算法计算震源的径向距离和深度.模型数据和实际资料的处理结果表明,本文方法具有较强的抗噪性,计算得到的震源方位角更加接近真实值;与常规目标函数相比,本文方法采用的目标函数具有更好的收敛性,其定位结果受初至拾取误差和速度模型误差的影响更小;本文提出的搜索算法能够消除由于错误拾取造成的观测到时中的异常值对定位结果的影响.     

Locating near-surface scatterers using non-physical scattered waves resulting from seismic interferometry

We use controlled-source seismic interferometry (SI) and inversion in a unique way to estimate the location of near-surface scatterers and a corner diffractor by using non-physical (ghost) scattered surface and body waves. The ghosts are arrivals obtained by SI due to insufficient destructive interference in the summation process of correlated responses from a boundary of enclosing sources. Only one source at the surface is sufficient to obtain the ghost scattered wavefield. We obtain ghost scattered waves for several virtual-source locations. To determine the location of the scatterer, we invert the obtained ghost traveltimes by solving the inverse problem. We demonstrate the method using scattered surface waves. We perform finite-difference numerical simulations of a near-surface scatterer starting with a very simple model and increase the complexity by including lateral inhomogeneity. Especially for the model with lateral variations, we show the effectiveness of the method and demonstrate the estimation of the subsurface location of a corner diffractor using S-waves. In all models we obtain very good estimations of the location of the scatterer.     

Comparison of the time-domain electromagnetic field from an infinitesimal point charge and dipole source

An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the dipole source.     

Acoustic source characterization of impulsive Strombolian eruptions from the Mount Erebus lava lake

We invert for acoustic source volume outflux and momentum imparted to the atmosphere using an infrasonic network distributed about the erupting lava lake at Mount Erebus, Ross Island, Antarctica. By modeling these relatively simple eruptions as monopole point sources we estimate explosively ejected gas volumes that range from 1,000 m³ to 24,000 m³ for 312 lava lake eruptions recorded between January 6 and April 13, 2006. Though these volumes are compatible with bubble volumes at rupture (as estimated from explosion video records), departures from isotropic radiation are evident in the recorded acoustic wavefield for many eruptions. A point-source acoustic dipole component with arbitrary axis orientation and strength provides precise fit to the recorded infrasound. This dipole source axis, corresponding to the axis of inferred short-duration material jetting, varies significantly between events. Physical interpretation of dipole orientation as being indicative of eruptive directivity is corroborated by directional emissions of ejecta observed in Erebus eruption video footage. Although three azimuthally distributed stations are insufficient to fully characterize the eruptive acoustic source we speculate that a monopole with a minor amount of oriented dipole radiation may reasonably model the primary features of the recorded infrasound for these eruptions.