深海热液硫化物矿体3D瞬变电磁正演

深海热液硫化物矿体瞬变电磁的正演是考虑深海环境的全空间条件下三维体的涡流电磁响应.采用全空间矢量有限元法模拟计算深海热液硫化物矿的三维瞬变电磁响应,对硫化物矿体采用矩形单元模型剖分,应用Galerkin法推导有限元方程,先计算频率域响应,再通过Fourier反变换将其转换至时间域,得出深海热液硫化物矿矿体的瞬变电磁响应.并用双半空间模型的解析解检验了全空间矢量有限元法模拟计算算法和程序的正确性,最后按照等比例缩小电磁物理实验原则,比对数值计算和物理实验结果论证了全空间3D模型数值的正确性.结果表明:对于海水、矿体以及围岩复杂电磁边界,应用全空间矢量有限元法模拟计算深海热液硫化物矿瞬变电磁响应异常与物理模拟结果一致,而且计算方法简单精确,异常幅值明显,边界清晰.     

华北煤田瞬变电磁勘探深度研究

为研究瞬变电磁勘探在华北煤田的勘探深度,利用瞬变电磁感应电压晚期段表达式,通过指定仪器最小可分辨电压,推导出瞬变电磁勘探最晚可检测出有效信号时间及最大探测深度表达式,计算了特定情况下的极限探测深度和最晚有效信号检出时间:在发射磁矩为10002×20 A.m2,地层电阻率分别为10和100Ω.m时,瞬变电磁勘探最晚有效信号检出时间为209.6 ms和52.6 ms,最大探测深度为1823 m和2889 m;以实测典型华北煤系地层瞬变电磁曲线为例,其计算机反演深度达1400 m;结果表明:采用合适的仪器和正确的工作装置,瞬变电磁勘探可以满足华北煤田目前开采深度的勘探需要.     

井中瞬变电磁旁测范围初探

为查明井中瞬变电磁旁测能力,从室内数值模拟和野外勘查实例两个角度开展了井中瞬变电磁旁测范围研究.在室内,以单边长度100 m、厚度10 m、电导率为100 S/m的水平板体为研究对象,通过建立相同工作参数、不同空间位置的地电模型,分析其理论响应认为:当目标体埋深小于1500 m时,其与钻孔距离在200 m以内时能够被井中瞬变电磁所识别;且随着深度变浅,井中瞬变电磁旁测范围扩大.在野外,利用大官厂地区的井中瞬变电磁实测,准确识别钻孔南西方位约71 m、埋深约652 m的隐伏火山口一处.室内数值模拟及野外现场实测表明,井中瞬变电磁具备在深部发现井旁低阻地质异常体的能力,其旁测范围甚至可达井旁200 m,可推广应用于我国的深部找矿工程.     

线源瞬变电磁法在熊耳山区多金属找矿中的应用研究

通过理论研究和野外试验,总结了线源瞬变电磁法在熊耳山区深部勘查时使用的装置参数,并把该装置应用到龙门店矿K9矿脉的深部勘查中.线源瞬变电磁法装置最主要的参数是偏离距和线源的长度.理论上偏离距的大小只影响到瞬变衰减曲线首支幅值的高低和衰减的快慢.但在野外试验时发现当偏离距为500 m时,该方法产生了强激发极化效应,当偏离距为1000 m时,激发极化效应减弱.说明强激发极化效应仅在收发距较小时产生,实测表明偏离距大于840 m可以有效压制激发极化干扰.线源AB长度一般取1000~2000 m,当信噪比大时取小值.该装置在龙门店矿区的应用表明,线源瞬变电磁法装置轻便灵活,探测深度大,分辨率高,可多台接收机同时工作,经济快速,适合山区深部有色金属地质勘查工作.     

超慢速扩张西南印度洋中脊硫化物成矿模型

洋中脊多金属硫化物已经成为人类重要的战略资源,科学的成矿模型是对其调查研究和勘探的重要依据.相较快速、慢速扩张洋中脊,超慢速扩张洋中脊在岩浆供给、构造和围岩等特征均存在明显差异,但目前对其热液循环及硫化物成矿模型缺乏系统梳理,制约了其资源勘探评价与研究的有效进程.本文系统总结了超慢速扩张西南印度洋中脊热液活动分布以及典型热液区的构造、热源、热液通道、围岩类型、流体性质和硫化物等特征,根据其成矿地质背景的差异性特点将该洋中脊赋存的热液系统分为局部强岩浆控制型、单向拆离/高角度大偏移距断层控制型以及双向拆离控制型三类,根据岩浆供给率(M值)的大小进一步将其划分为五种类型,从而建立了超慢速扩张西南印度洋中脊的局部强热供给-深大断裂控制硫化物成矿模型.超慢速扩张西南印度洋中脊扩张速率整体变化不大(14～18mm/a),岩浆供给呈分段不均匀性.通过近20年的调查研究,发现其发育类型多样的热液系统和硫化物.在岩浆供给充足的洋脊段,发育局部强岩浆供给条件下的深部岩浆房(4～9km).而在岩浆供给贫瘠的洋脊段,发育长期持续活动的深大拆离断层(可达13km),并沿拆离断层形成成矿带.因而超慢速扩张洋脊具...     

浅水区的瞬变电磁法：一维数值模拟结果分析

在浅水环境中,利用频率域海洋可控源电磁法探测油气藏薄层面临着挑战,这是因为空气波支配着电磁响应而它只含有少量海底电阻率结构的信息.本文研究了空气波对时间域可控源电磁(CSEM)响应的影 响. 空气波到达时间取决于海水深度.在浅水区空气波到达早,而在深海水区空气波到达晚,空气波与来自深部电阻体的电磁波出现在不同时间段.虽然在中等深度的海水区空气波和来自深部电阻体的电磁信号几乎同时到达,但与不含有油气藏的背景模型相比,仍可以看到明显的异常.浅水区勘探的好处在于可以使用表面拖曳系统.     

浅海瞬变电磁全波形响应特征及探测能力分析

瞬变电磁法在浅海工程勘探等领域受到了越来越多的关注.目前浅海瞬变电磁仍处于应用初期,相关研究少且未有成熟装备问世,有必要研究其探测能力并为观测系统选取最佳观测参数.本文以几种典型发射波形为例,采用褶积算法细致分析了不同发射波形条件下浅海瞬变电磁on-/off-time响应受海底介质电导率、磁导率及发射波形脉宽等参数的影响特征与规律;通过三维正演并设定极限探测深度阈值,进一步分析不同发射波形on-/offtime期间浅海瞬变电磁探测能力及对典型三维目标体的极限探测深度.基于本文研究成果,可为浅海瞬变电磁探测装置设计、观测系统的参数选取及试验参数的选取等提供了一些有价值的理论借鉴.     

地空瞬变电磁法研究进展

随着瞬变电磁法的研究逐渐转向高海拔山区等地形地质条件复杂地区进行深部勘探,地空瞬变电磁法得到了更多的关注.地空瞬变电磁法结合了地面瞬变电磁法和航空瞬变电磁法各自的特点,相比地面瞬变电磁法,能有效提高野外工作效率;相比航空瞬变电磁法,勘探深度更大、采集信号信噪比更高.近几年来,地空瞬变电磁法得到了较快的发展.本文回顾了地空瞬变电磁法的历史发展与现状,总结了地空瞬变电磁法研究的最新进展,并对地空瞬变电磁法的各种装置类型进行了相关评述和总结,指出目前地空瞬变电磁法的问题主要在于解释水平跟不上,大都仍采用基于地面瞬变电磁的数据解释方法.为了进一步丰富地空瞬变电磁勘查体系,从仪器装备、数据预处理和处理技术、正反演及数据解释等方面对地空瞬变电磁法的发展方向进行了总结.     

海洋拖曳式水平电偶源数值模拟与电场接收机研制

海洋天然气水合物资源调查中,电磁法已成为地震勘探最有效的补充手段之一,而我国海洋水合物电磁探测中的拖曳式电场接收机此前是空白.本文从方法原理、仪器设备和海洋试验等方面,对海洋拖曳式水平电偶极-偶极方法进行研讨.建立三维电阻率模型,进行数值模拟,讨论收发距、发射频率、拖曳装置离海底高度等参数与电场异常响应的关系,针对特定模型给出了最佳观测参数以指导仪器硬件设计与海洋测试.在此基础上,采用高级嵌入式计算机的智能化控制、前端低噪声斩波放大、高精度同步采集等技术,研制了由拖体、承压舱、采集电路、低噪声电场传感器等组成的海洋拖曳式水平轴向电场接收机,其系统噪声达到2 nV/√HZ @(1~100Hz).将自主研制的电场接收机用于国内首次拖曳式电偶极-偶极方法的海洋试验与评估,并对海试资料进行了处理与分析.研究结果表明,所研发的仪器达到了设计指标,能够用于海洋天然气水合物调查.     

瞬变电磁法正反演问题研究进展

对瞬变电磁法的方法发展概况和仪器研制状况做出了综述性评价.对瞬变电磁法正反演问题的研究成果进行了系统总结.目前的数值模拟正演方法主要有一维滤波系数法,三维积分方程法,二维,三维有限差分法,2.5维有限元法等,主要的反演方法有:一维浮动薄板解释法,人机对话自动反演法,烟圈理论解释法,神经网络反演法,成像类反演等,论述了瞬变电磁法各种计算方法的特点.瞬变电磁法的正反演发展趋势主要是研究三维正反演的计算方法和目标体成像系统.     

The behavior of deep convective clouds over the warm pool and connection to the Walker circulation

As the deep convective clouds (DCCs) over the western Pacific and Indian Ocean warm pool may play different roles in the climate system, variations in DCC properties over these two sectors are investigated and compared. The DCC intensity and area varies more significantly in the Indian Ocean than the western Pacific sector, while the DCC frequency is comparable in both sectors at the seasonal scale. Although the Indian Ocean sector is strongly dominated by the seasonal evolution, the interannual variations in the two sectors are comparable for all three DCC properties (frequency, intensity, and area). Besides, Walker circulation is closely correlated with the interannual variability of DCCs in both sectors. The Walker circulation strengthens (weakens) as the DCCs shift eastward (westward) over the Indian Ocean sector and westward (eastward) over the western Pacific sector. When more or stronger DCCs occur over the Indian Ocean sector (western Pacific sector), the Walker circulation becomes stronger (weaker) and shifts westward (eastward). Interestingly, the response of the Walker circulation to DCC variability over the warm pool is asymmetry. The asymmetry response of the Walker circulation to the negative and positive DCC anomaly may be related to the non-linearity internal variability of the atmosphere. DCCs over the Indian Ocean sector have a much weaker nonlinear correlation with the Walker circulation than DCCs over the western Pacific sector.     

Effects of a variety of Indian Ocean surface temperature anomaly patterns on the summer monsoon circulation: Experiments with the NCAR general circulation model

The time mean response of the summer monsoon circulation, as simulated by the 2.5° latitude-longitude resolution, July version of the National Center for Atmospheric Research (NCAR) General Circulation Model (GCM), to a variety of Indian Ocean surface temperature anomaly patterns is examined. In separate experiments, prescribed changes in surface temperature are imposed in the Western Arabian Sea, the Eastern Arbian Sea or the Central Indian Ocean. The influence of these anomaly patterns on the simulated summer monsoon circulation is evaluated in terms of the geographical distribution of the prescribed change response for any field of interest. This response is defined as the grid point difference between a 30-day mean from a prescribed change experiment and the ensemble average of the 30-day means from the control population for which the same set of climatological ocean surface temperatures are used in each simulation. The statistical significance of such a prescribed change response is estimated by relating the normalized response (defined as the ratio of the prescribed change response to the standard deviation of 30-day means as estimated from the finite sample of control cases) to the classical Student'st-statistic. Using this methodology, the most prominent and statistically significant features of the model's response are increased vertical velocity and precipitation over warm anomalies and typically decreased vertical velocity and precipitation in some preferred region adjacent to the prescribed change region. In the case of cold anomalies, these changes are of opposite sign. However, none of the imposed anomaly patterns produces substantial or statistically significant precipitation changes over large areas of the Indian sub-continent. The only evidence of a major nonlocal effect is found in the experiment with a large positive anomaly (+3°C) in the Central Indian Ocean. In this instance, vertical velocity and precipitation are reduced over Malaysia and a large area of the Equatorial Western Pacific Ocean. Thus, while these anomaly experiments produce only a local response (for the most part), it is hoped, as one of the purposes of the planned Monsoon Experiment (MONEX), that the necessary data will be provided to produce detailed empirical evidence on the extent to which Indian Ocean surface temperature anomalies correlate with precipitation anomalies over the Indian subcontinent—a correlation which generally does not appear in these GCM results.The National Center for Atmospheric Research is sponsored by the National Science Foundation     

On the interpretation of near-bottom water temperature anomalies

A positive water temperature anomaly of 0.11°C and an inverse gradient of potential temperature of 1.5 × 10^?2°C/m has been measured at the TAG hydrothermal field in the rift valley of the Mid-Atlantic Ridge at latitude 26°N by means of a thermistor array towed between 2 and 20 m above the seafloor. This anomaly appears to be associated with hydrothermal discharge from the oceanic crust. The temperature data are interpreted in terms of (1) a steady, turbulent thermal plume rising in a homogeneous, neutrally buoyant medium, and (2) turbulent diffusion in the ocean-bottom boundary layer. The calculations indicate that the thermal output of the TAG anomaly area is of the order of several megawatts, which is of the same order of magnitude as some continental geothermal systems. The thermal output from the TAG anomaly area represents a significant fraction of the total heat loss resulting from the generation of new lithosphere at the Mid-Atlantic Ridge at 26°N.     

电磁辐射"EMOLS"仪观测结果原理及震例

“九五”期间研制的“EMAOS”电磁辐射仪,其接收原理是按大地电磁测深的原理改进而研制的。选择接收频段为0．01～10Hz,该频段探测深度在0～20km,是浅源地震发生的层位。记录的电磁扰动是在应力场作用下,动电效应产生的电磁场。电场和磁场同时记录增强了信息可信度的判断能力,记录震例表明,异常规律与已有模拟结果基本相同。因此,已总结的模拟记录预报指标等可在数字记录结果中使用。     

A finite-element time-domain forward-modelling algorithm for transient electromagnetics excited by grounded-wire sources

A three-dimensional finite-element time-domain forward-modelling algorithm is developed to simulate transient electromagnetics excited by grounded-wire sources. The main advantage of this finite-element time-domain algorithm is that full transmitting-current waveforms and complex-shaped sources resulting from topography can be directly dealt with in this algorithm. The models used to test this algorithm include a homogeneous half-space model, a stratified-medium model, the model of a complex conductor at a vertical contact and the Ovoid Zone massive sulfide deposit at Voisey's Bay, Canada. The homogeneous half-space model is used to determine the truncation boundary for a computational domain, and to compare with the electromagnetic responses excited by step-off, step-on and direct current waveforms. For the stratified-medium model, results demonstrate that full transmitting waveforms have strong effects on the observed electromagnetic responses. The model of a complex conductor at a vertical contact is designed for the grounded electrical source airborne transient electromagnetic method and is also used to examine the effectiveness of the broadside and inline configurations for such a vertical, thin plate embedded in the subsurface. The area of the Ovoid Zone massive sulfide deposit possesses non-negligible topography, the effects of which on the shapes of the grounded-wire sources must be taken into account when implementing the finite-element time-domain solver. The results show that both the broadside and inline electromagnetic responses are strongly affected by the massive conductive ore body.     

Equipment for carrying out quantitative geoelectromagnetic model experiments

Summary The paper describes equipment designed to study the effect of a conducting hostrock/overburden on the electromagnetic anomaly of sulphide ore bodies embedded therein. The model was constructed strictly according to the theory of electromagnetic similitude so as to constitute a direct reading reproduction of the field vectors. The experiments were conducted at a fixed, crystal-controlled frequency of 100 kHz using mainly a graphite sheet immersed in a dilute hydrochloric acid solution of predetermined conductivity to simulate vein and manto type of ore deposits surrounded by a partially conducting zone. Both the inphase and quadrature components of the anomaly were measured in terms of the primary field after elimination of the regional anomaly by means of a measuring bridge and a compensator. A comparison of the anomaly profiles obtained over the ore model when it is held in air with those obtained for corresponding situations inside a conducting solution shows a general enhancement of response [1–3].     

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Réunion)

 A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Réunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60–100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos, the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise. Received: 3 November 1999 / Accepted: 15 September 1999     

Three‐dimensional inversion of towed streamer electromagnetic data

A towed streamer electromagnetic system capable of simultaneous seismic and electromagnetic data acquisition has recently been developed and tested in the North Sea. We introduce a 3D inversion methodology for towed streamer electromagnetic data that includes a moving sensitivity domain. Our implementation is based on the 3D integral equation method for computing responses and Fréchet derivatives and uses the re‐weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present two model studies relevant to hydrocarbon exploration in the North Sea. First, we demonstrate the ability of a towed electromagnetic system to detect and characterize the Harding field, a medium‐sized North Sea hydrocarbon target. We compare our 3D inversion of towed streamer electromagnetic data with 3D inversion of conventional marine controlled‐source electromagnetic data and observe few differences between the recovered models. Second, we demonstrate the ability of a towed streamer electromagnetic system to detect and characterize the Peon discovery, which is representative of an infrastructure‐led shallow gas play in the North Sea. We also present an actual case study for the 3D inversion of towed streamer electromagnetic data from the Troll field in the North Sea and demonstrate our ability to image all the Troll West Oil and Gas Provinces and the Troll East Gas Province. We conclude that 3D inversion of data from the current generation of towed streamer electromagnetic systems can adequately recover hydrocarbon‐bearing formations to depths of approximately 2 km. We note that by obviating the need for ocean‐bottom receivers, the towed streamer electromagnetic system enables electromagnetic data to be acquired over very large areas in frontier and mature basins for higher acquisition rates and relatively lower cost than conventional marine controlled‐source electromagnetic methods.