瞬变电磁法对存在山谷地形时的多个异常体的探测能力研究

以层状地层中的2个三维导电薄板为例,考虑了山谷地形存在时,瞬变电磁法(TEM)分辨地下多个三维异常体的能力。模型计算的结果表明,当有山谷地形存在时,总的异常场由薄板的异常和地形的异常叠加在一起,表明异常场在空间上具有可加性质。在地形平坦时,TEM方法对相同尺寸的近源异常体的分辨能力要优于远源,而且离源越远的异常体异常越小,最后直到无法分辨。但是当存在山谷地形并且把源放在其中时,它们的总异常场叠加后会“放大”远源异常体的异常场,在一定意义上说,这对于分辨远处的异常体更加有利。当异常体在接收器下方埋深不是太大时,不论山谷地形在源处还是接收器处,TEM法对这些异常体的探测能力都是比较好的。2个异常体相距太近时,其异常会叠加在一起,给TEM的分辨增加一些困难。总体表明,在有山谷地形存在时,TEM方法对多个异常体仍具有较好的探测能力     

不规则回线瞬变电磁法一维烟圈反演研究

瞬变电磁法中最常用矩形或圆形回线,但有时由于地形限制,只能使用不规则回线发射和接收.如果仍用现有的常规回线理论进行处理解释,得到的结果会存在很大的偏差.本文首先进行了不规则回线瞬变电磁法一维正演理论研究,基于电偶极子源的频率域响应公式,通过沿回线积分和时频转换,推导出不规则回线源在水平层状介质中的时间域响应公式.采用欧拉算法、高斯积分和快速汉克尔变换,计算不规则回线内任意一点处的磁场响应.利用改进的二分搜索法计算全区视电阻率,并在此基础上利用烟圈法反演电阻率和深度.通过四个典型地电模型的正反演计算,表明烟圈法反演能够有效反映地电模型的大致形态,可以用于不规则回线瞬变电磁数据的快速反演解释.     

ATTEM系统中电流关断期间瞬变电磁场响应求解的研究

在瞬变电磁法中,由于发射电流关断时间不为零、接收线圈的谐振频率有限,早期瞬变电磁信号发生畸变,只能舍弃,因此存在着探测盲区. 针对这一问题,研究了瞬变电磁方法中发射电流关断期间总磁场的形成过程,论证了一次场、二次场和总瞬变场的关系,分析了接收线圈的频率特性和关断时间对瞬变电磁场的影响,提出从总磁场中剔除一次磁场影响的方法,从而获得电流关断期间和电流关断后的早期瞬变电磁场. 采用吉林大学自主研制的瞬变电磁测量系统(ATTEM)在长春市伊通河活断层进行勘探,进一步验证了算法的有效性,缩短了瞬变电磁法的勘探盲区,实现了近地表4 m以下的勘探,可以清晰地分辨近地表的低阻异常,提高了浅层探测精度和分辨率.     

考虑关断时间的地面瞬变电磁三维带地形反演

起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性.     

多通道瞬变电磁法2D有限元模拟

多通道瞬变电磁法的反演工作是在大地脉冲响应的基础上进行的,因此本文用2D有限元正演来求取大地脉冲响应,实现多道瞬变电磁法的正演模拟研究.研究工作首先在有源大地电磁法2D有限元模拟的基础上进行,原有2D有限元法的频带相对于多通道瞬变电磁法的频带是窄频的,为此采用校正法将其扩展到所需的宽度,实现了频率域宽频带有源电磁勘探方法的正演模拟;然后通过频时变换变换到时间域,得到瞬变电磁法的阶跃响应;最后通过求取阶跃响应的时间导数,得到大地脉冲响应.针对多通道瞬变电磁方法的装置特点,用有一定埋深的顺层成矿模型分析了大地脉冲响应特征,发现对有一定埋深的顺层产出的矿脉模型来说,频率域电磁场结果以及时间域大地脉冲响应结果均需移除背景场才能突出异常体的存在,然而,在实际工作中,未知的背景场限制了其实用性.为此,本文利用多通道瞬变电磁方法数据量大的特点,提出用不同偏移距的等时曲线与等时断面来展示脉冲响应模拟结果,事实表明这两种展示方式可以清晰地分辨出矿体引起的异常,且能准确定位矿体的中心在地面的投影位置,从而说明多通道瞬变电磁方法相对于频率域有源电磁方法,对有一定埋深的顺层矿床有更高的分辨能力.     

井中磁源瞬变电磁响应特征研究

井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率...     

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

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

多孔径瞬变电磁场物理模拟

为了提高瞬变电磁法对地探测精度,提出了一种新颖的瞬变电磁装置一多孔径发射装置,并进行了多孔径瞬变电磁场物理模拟观测.分别使用四孔径发射装置和常规回线发射装置,进行不同收一发距情况下的一次场观测,并对铜板在不同埋藏深度情况下的二次响应进行观测,通过比较两种不同装置下的响应曲线和计算观测数据误差,最终发现使用四孔径发射装置可使一次场强度提高34%,可使二次场强度提高31%.结论表明,多孔径发射源可以改变电磁场的空间分布,认为瞬变电磁场存在类相干现象.研究成果对提高瞬变电磁的探测精度具有很重要的意义.     

瞬变电磁法的探测深度问题

用解析分析、时域有限差分、时-频分析的方法,以地面中心回线装置和阶跃脉冲激励源为例,分析讨论了瞬变电磁测深法的勘探深度问题,以便为野外勘探工作设计提供依据,达到预期的探测目的.解析计算证实了瞬变场在地下以有限速度传播,数值模拟表示出了准静态条件下瞬变场的反射.研究结果表明,由于时间域电磁场遵循因果律,瞬变电磁法的探测深度主要由观测时间决定. 瞬变电磁场的初始传播速度与大地电阻率无关,继后在大地色散作用下,阶跃脉冲前沿逐渐变得平缓,各频率分量的传播速度与电阻率有关,在低阻地层中探测同样的深度需要较长的观测时间. 最大探测深度是在给定时间内电磁波往返地下某一深度的单程距离,最小探测深度受仪器性能的限制,但是埋藏较浅的异常体也有可能在晚时段被观测到.从时-频密度谱中可得到瞬变电磁场信号时间与频率的关系.     

海底油气藏地质模型的冲激响应

海洋可控源电磁法(mCSEM)的时间域冲激响应特征可以反映海底油气高阻薄层.本文计算了水平电偶极子源均匀大地半空间,海洋均匀双半空间和海洋四层模型的阶跃响应和冲激响应,提出了瞬变冲激时刻的概念.分析了水平电偶源瞬变冲激时刻与介质电导率的指示关系.对于海底油气高阻薄层宜采用多偏移距同时测量方式,由于在低电导率介质中电磁能量传播得要快,在适当的收发距瞬变冲激时刻会提前到达,提出的瞬变冲激时刻道间变化量可以明确指示高阻薄层的存在及埋深.文中还分析了海水深度对瞬变冲激时刻的影响.由于“天波”干扰,瞬变冲激响应受到一定收发距观测的限制.消除 “天波”影响是时间域和频率域mCSEM数据处理的研究热点.     

三维起伏地形条件下航空瞬变电磁响应特征研究

航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论.     

A first application of a marine inductive source electromagnetic configuration with remote electric dipole receivers: Palinuro Seamount,Tyrrhenian Sea

We study a new marine electromagnetic configuration that consists of a ship‐towed inductive source transmitter and a series of remote electric dipole receivers placed on the seafloor. The approach was tested at the Palinuro Seamount in the southern Tyrrhenian Sea, at a site where massive sulphide mineralization has been previously identified by shallow drilling. A 3D model of the Palinuro study area was created using bathymetry data, and forward modelling of the electric field diffusion was carried out using a finite volume method. These numerical results suggest that the remote receivers can theoretically detect a block of shallowly buried conductive material at up to ～100 m away when the transmitter is located directly above the target. We also compared the sensitivity of the method using either a horizontal loop transmitter or a vertical loop transmitter and found that when either transmitter is located directly above the mineralized zone, the vertical loop transmitter has sensitivity to the target at a farther distance than the horizontal loop transmitter in the broadside direction by a few tens of metres. Furthermore, the vertical loop transmitter is more effective at distinguishing the seafloor conductivity structure when the vertical separation between transmitter and receiver is large due to the bathymetry. As a horizontal transmitter is logistically easier to deploy, we conducted a first test of the method with a horizontal transmitter. Apparent conductivities are calculated from the electric field transients recorded at the remote receivers. The analysis indicates higher apparent seafloor conductivities when the transmitter is located near the mineralized zone. Forward modelling suggests that the best match to the apparent conductivity data is obtained when the mineralized zone is extended southward by 40 m beyond the zone of previous drilling. Our results demonstrate that the method adds value to the exploration and characterization of seafloor massive sulphide deposits.     

An improved interpolation scheme at receiver positions for 2.5D frequency-domain marine controlled-source electromagnetic forward modelling

An improved interpolation scheme is presented for 2.5-dimensional marine controlled-source electromagnetic forward modelling. For the marine controlled-source electromagnetic method, due to the resistivity contrast between the seawater and seafloor sedimentary layers, it is difficult to compute the electromagnetic fields accurately at receivers, which are usually located at the seafloor. In this study, the 2.5-dimensional controlled-source electromagnetic responses are simulated by the staggered finite-difference method. The secondary-field approach is used to avoid the source singularities, and the one-dimensional layered background model is used for calculating the primary fields excited by the source quasi-analytically. The interpolation of electromagnetic fields at the cell nodes for the whole computational domain to the receiver locations is discussed in detail. Numerical tests indicate that the improved interpolation developed is more accurate for simulating the electromagnetic responses at receivers located at the seafloor, compared with the linear or rigorous interpolation.     

海底油气藏及天然气水合物的时频电磁辨识

提出了伪随机多频海洋电磁法观测方案.采用伪随机多频信号作为激励场源,多偏移距同线偶极-偶极同时观测,相关辨识海底地电系统的频率特性和冲激响应,可以在时间域和频率域同时辨识海底高阻薄层.在时间域,瞬变冲激时刻可以直接指示海底地层电导率的变化;在频率域,利用多个频率的电场响应计算的频散率及其道闻变化量,相对相位道间变化量对高阻薄层有很好的反映.从而实现对海底油气及天然气水合物的多参数辨识.     

Transient electromagnetic surveys with unimodal transverse magnetic field: ideas and results

The theory behind transient electromagnetic surveys can be well described in terms of transverse magnetic and transverse electric modes. Soundings using transverse magnetic and transverse electric modes require different source configurations. In this study, we consider an alternating transverse magnetic field excitation by a circular electric dipole. The circular electric dipole transmitter is a horizontal analogue of the vertical electric dipole. Offshore surveys using circular electric dipole might represent an alternative to the conventional marine controlled‐source electromagnetic method at shallow sea and/or for exploring relatively small targets. Field acquisition is carried out by recording either electric or magnetic responses. Electric responses bear information on the 1D structure of a layered earth and successfully resolve high‐resistivity targets in marine surveys. Land‐based circular electric dipole soundings are affected by induced polarisation. On the contrary, magnetic responses are absent on the surface of a 1D earth, and as a result, they are very sensitive to any and even very small 3D conductivity perturbations. In addition, they are sensitive to induced polarisation or some other polarisation effects in the subsurface. At present, circular electric dipole transmitters and magnetic receivers are successfully used in on‐land mineral and petroleum exploration.     

Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data

Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.     

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

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

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

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