福州八一水库-尚干断裂的高密度电法和瞬变电磁法试验探测

为了检验高密度电法和瞬变电磁法在城市活断层探测中的有效性 ,在福州市区及周围进行了试验探测① ,文中主要介绍对八一水库尚干断裂的探测结果。八一水库尚干断裂是一条穿过福州市区的主要断裂带 ,其主断裂在八一水库测区位于一条由新店通往古城的水泥马路东侧 ,在横向上表现为东西宽几十米的电阻率低异常带 ,在纵向上从盖层底部向下延伸到基底深处。试验表明 ,高密度电法和瞬变电磁法可以互相补充和验证 ,从而提高整体探测效果     

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.     

Simple relative space–time scaling of electrical and electromagnetic depth sounding arrays: implications for electrical static shift removal and joint DC-TEM data inversion with the most-squares criterion

A simple scaling relationship is shown to facilitate comparison, correlation and integration of data recorded using the common experimental configurations in electrical and electromagnetic depth sounding. Applications of the scheme to field data from typical geological and landfill environments show that it is robust and, where transient electromagnetic (TEM) data are available, enables easy identification and quantification of electrical static shift (galvanic distortion) in magnetotelluric and direct current (DC) sounding curves. TEM‐based procedures are suggested for both the direct removal of static shift in DC sounding curves and effective joint data inversion with the most‐squares criterion in the presence of static shift. A case study of aquifer characterization using sounding data from borehole sites in the Vale of York in England shows that static shift is a common problem in this glacial‐covered terrain and demonstrates the effectiveness of the proposed joint DC‐TEM inversion strategy in handling distorted soundings.     

Appraisal of equivalence and suppression problems in 1D EM and DC measurements using global optimization and joint inversion*

Global optimization with very fast simulated annealing (VFSA) in association with joint inversion is performed for 1D earth structures. The inherent problems of equivalence and suppression in electromagnetic (EM) and direct current (DC) resistivity methods are studied. Synthetic phase data from multifrequency sounding using a horizontal coplanar coil system and synthetic apparent resistivity data from Schlumberger DC resistivity measurements are inverted individually and jointly over different types of layered earth structures. Noisy data are also inverted. The study reveals that global optimization of individual data sets cannot solve inherent equivalence or suppression problems. Joint inversion of EM and DC measurements can overcome the problem of equivalence very well. However, a suppression problem cannot be solved even after combination of data sets. This study reveals that the K-type earth structure is easiest to resolve while the A-type is the most difficult. We also conclude that the equivalence associated with a thin resistive layer can be resolved better than that for a thin conducting layer.     

基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟

直流电阻率法被广泛应用在工程和环境及水文地球物理、野外采矿、地热探测等领域.地下岩石常具有层理面和裂缝等具有方向依赖性的结构,岩石电导率常常呈现各向异性特征,因此研究复杂直流电阻率各向异性问题的高精度正演算法具有迫切的理论和学术需求.本文利用面向目标的自适应有限元算法和非结构化网格相结合的方式,解决了带地形任意复杂直流电阻率各向异性问题的高精度正演这一难题.有别于前人的研究成果,本文提出了一种特别的二次虚拟场算法实现带源的任意起伏地形问题模拟;另外,本文第一次基于电流密度连续性条件构建适合直流电阻率各向异性问题的后验误差估计算法,有效地驱动面向目标有限元网格自适应加密过程.最后,通过三组电阻率各向异性模型验证本文提出算法的正确性和适应性,测试结果表明：对于任意复杂直流电阻率各向异性问题,本文提出的算法具有精度高、适应性强等特点;另外,我们还发现电流密度连续性条件可用于设计直流电阻率问题的有效后验误差估计算法.     

高效多路直流开关电源

本文介绍一种新型的直流开关电源,该电源是一种直流电压输入,三种直流电压输出,并能够达到相当高的效率。这种电源适用于数字化设备的无人职守台站。本文也论述了实现的技术思路和达到的技术要求。     

Three-dimensional distributed DC/IP method for altered tectonite-type gold ore deposits exploration: a case study of the Jiaojia gold metallogenic belt,Eastern China

To develop an effective method to identify ore-controlling faults, we studied the Jiaojia gold metallogenic belt, a most typical altered tectonite-type gold metallogenic belt in the Jiaodong Peninsula, China, and conducted experiments using the 3D distributed direct current-induced polarization (DC/IP) method. Firstly, we tested the ability of using 3D distributed DC/IP method to identify altered tectonite-type gold ore deposits by 3D synthetic modelling. We then collected real data of the Sizhuang gold deposit using the 3D distributed DC/IP method. The resistivity model obtained of this region is generally consistent with the known geological setting. Moreover, to obtain the information about the southern extension of the Jiaojia gold metallogenic belt, we conducted a 3D distributed DC/IP experiment in the Shijia area in the southern segment of the Jiaojia fault. The southern extension of the Jiaojia fault and tectonic evolution of shallow magma in this region were inferred from the 3D resistivity and chargeability models. Based on all the information above, we concluded that the 3D distributed DC/IP method has the advantages of 3D observations, high spatial resolution and great detection depth and will be one of the most effective methods for detecting altered tectonite-type gold ore deposits.     

Tectonic analysis of mine tremor mechanisms from the Upper Silesian Coal Basin

Fault network of the Upper Silesian Coal Basin (USCB) is built of sets of strike-slip, oblique-slip and dip-slip faults. It is a typical product of force couple which acts evenly with the parallel of latitude, causing horizontal and anti-clockwise movement of rock-mass. Earlier research of focal mechanisms of mine tremors, using a standard fault plane solution, has shown that some events are related to tectonic directions in main structural units of the USCB. An attempt was undertaken to analyze the records of mine tremors from the period 1992–1994 in the selected coal fields. The digital records of about 200 mine tremors with energy larger than 1×10⁴ J (M _L >1.23) were analyzed with SMT software for seismic moment tensor inversion. The decomposition of seismic moment tensor of mine tremors was segmented into isotropic (I) part, compensated linear vector dipole (CLVD) part and double-couple (DC) part. The DC part is prevalent (up to 70%) in the majority of quakes from the central region of the USCB. A group of mine tremors with large I element (up to 50%) can also be observed. The spatial orientation of the fault and auxiliary planes were obtained from the computations for the seismic moment DC part. Study of the DC part of the seismic moment tensor made it possible for us to separate the group of events which might be acknowledged to have their origin in unstable energy release on surfaces of faults forming a regional structural pattern. The possible influence of the Cainozoic tectonic history of the USCB on the recent shape of stress field is discussed.     

关于人工源极低频电磁波发射源的讨论——均匀空间交流点电流源的解

大功率人工源极低频电磁波（Control Source of Extremely Low Frequency method, CSELF）技术中,辐射天线由布设在高大地电阻率区的接地长导线源构成,通过接地点向大地注入数百安培的强电流.目前有关强电流在大地中的流动特性和分布规律,尚未完全弄清楚.本文提出一个新思路,将CSELF发射天线分解为两部分:交变的接地长导线源和接地处的交流点电流源,二者响应的叠加构成总场.本文重点讨论了均匀空间下交流点电流源的求解,对比研究了交流点电流源响应、直流点电流源响应之间的差异.结果表明,在小范围内（场源距几十、几百米以内）,均匀空间下交流点电流源场和直流源场差别较小,可相互近似,但在大范围内,交流点电流场比直流场的衰减要快得多.由于CSELF辐射天线跨度上百公里,场源距很大,其电流分布规律遵循交变场规则,不能近似为直流场.     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell’s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.     

Polar ionic conductivity profile in fair weather conditions. Terrestrial test of the Huygens/Hasi-PWA instrument aboard the Comas Sola balloon

The permittivity wave and altimetry (PWA) instrument is a part of the CASSINI/HUYGENS HASI experiment and was designed to determine the electrical parameters of the atmosphere of Titan in 2004. In December 1995, a balloon campaign was conducted in León, Spain, to test the HASI onboard hardware and software using a HUYGENS probe mock-up in an electromagnetic-disturbance-free environment (mainly from power emission lines at 50 Hz). This work is concerned with the measurements of small ion polar conductivities and DC fields using the PWA relaxation probes (RP). The two RP electrodes were periodically set to ±5 V relative to the conductive surface of the mock-up and allowed to discharge in the surrounding atmosphere. The polar components of conductivity are calculated from the discharge time, and the DC field from the floating potential differences once the electrodes reach equilibrium. In spite of some observed effects, such as mock-up charging or oscillations in the measurement of potential, the conductivity measurements are coherent and in good agreement with the obtained results in other experiments. The conductivity data were collected in ‘fair-weather’ conditions, up to 30 km during a 4-h flight, every 72 s, giving an altitude resolution better than 400 m. We also discuss the DC field data that do not lead, in presence of charging effects, to a straightforward measurement of the natural DC field. The Comas Solá balloon flight, first real test of the PWA experiment in the terrestrial atmosphere, confirmed the validity of the ionic conductivity measurements but raised the problem of a reliable interpretation of the DC field.     

Higher-order statistics correlation stacking for DC electrical data in the wavelet domain

DC (direct current) electrical and shallow seismic methods are indispensable to the near surface geophysical exploration, but the near surface areas are very difficult environments for any geophysical exploration due to the random noise caused by near surface inhomogeneities. As a new algorithm based on higher-order statistics theory, the higher-order correlation stacking algorithm for seismic data smoothing in the wavelet domain has been developed and applied efficiently to filter some correlation noise that the conventional second-order correlation stacking could not inhibit. In this paper, this higher-order statistics correlation stacking technology is presented for DC electrical data in wavelet domain. Taking into account the single section and multiple section data, we present two new formulations of correlation stacking for DC electrical data. Synthetic examples with Gaussian noise are designed to analyze the overall efficiency of the new algorithm and to determine its efficacy. Meanwhile, comparison with the traditional least-squares optimization inversion method for field examples from electrical imaging surveys and time-domain IP measurement in China shows its significant advantages. The quality of the new algorithm also has been assessed by physical simulation experiments. This new technology in DC electrical exploration measurements provides a new application in engineering and mining investigation.     

Using direct current resistivity sounding and geostatistics to aid in hydrogeological studies in the Choshuichi alluvial fan, Taiwan

Ground water reservoirs in the Choshuichi alluvial fan, central western Taiwan, were investigated using direct-current (DC) resistivity soundings at 190 locations, combined with hydrogeological measurements from 37 wells. In addition, attempts were made to calculate aquifer transmissivity from both surface DC resistivity measurements and geostatistically derived predictions of aquifer properties. DC resistivity sounding data are highly correlated to the hydraulic parameters in the Choshuichi alluvial fan. By estimating the spatial distribution of hydraulic conductivity from the kriged well data and the cokriged thickness of the correlative aquifer from both resistivity sounding data and well information, the transmissivity of the aquifer at each location can be obtained from the product of kriged hydraulic conductivity and computed thickness of the geoelectric layer. Thus, the spatial variation of the transmissivities in the study area is obtained. Our work is more comparable to Ahmed et al. (1988) than to the work of Niwas and Singhal (1981). The first "constraint" from Niwas and Singhal's work is a result of their use of linear regression. The geostatistical approach taken here (and by Ahmed et al. [1988]) is a natural improvement on the linear regression approach.     

井中震电效应的模拟实验研究(Ⅱ):孔隙饱和岩石部分

使用小功率直流电源(32V,3A),正交于声波传播方向供电,实验考察了直流电场对4块砂岩岩心的声学性质的影响,发现了一些以往文献没有报导的实验现象,总体上,有1)随通电时间增长,电流强度下降,最后趋于“饱和”;2)岩石的声速有先降后升的趋势,相对变化一般小于2％;3)纵波幅度呈现下降趋势,相对变化幅度最大可达到80％;4)声速和声幅的变化与岩石的渗透率有正相关关系,例如D7#20由于渗透率很低,因此声速和声幅的变化很小．结合岩石电阻率增大以及对含水饱和度影响岩石声学性质的一般认识,认为实验现象与电解导致的岩石孔隙内部含水饱和度的下降有关．研究结果有助于丰富对目前电法—地震勘探实验某些现象的认识。     

Evaluation of DC, FDEM and IP resistivity methods for imaging perched saltwater and a shallow channel within coastal tidal flat sediments

Presently, electrical resistivity methods are applied in a wide variety of geological and environmental site investigations. Geologically, the coastal tidal flat sediments formed shallow channel-like features at the northern part of Germany. Three geoelectrical methods are applied to image the near surface sediments including a shallow conductive zone within the tidal deposits at the North Sea coast. These methods, direct current (DC) resistivity, frequency domain electromagnetic (FDEM) and spectral induced polarization (SIP), are evaluated to show which one can provide the required spatial resolution under study area conditions. This evaluation also includes a synthetic modeling to assess the DC resistivity imaging technique.The results constitute an encouraging example using these geophysical methods in characterizing the coastal aquifers. The inversion results show that the subsurface resistivity distribution of tidal sediments can change rapidly within a short distance. A thin high conductive layer is observed above the peat and clay layers reflecting a perched saltwater. The 2D IP section shows that the perched saltwater is restricted to patched forms above an impermeable layer of clay. According to the IP images the boundaries of the clay layer are recognized with a good resolution due to the high membrane polarization of the clays. The EM and DC profiles show a shallow channel-like feature within tidal deposits. In this paper, the best FDEM field parameters and the role of EM in lithologic studies are emphasized. Two main limitations can be observed from DC synthetic modeling: (a) A smearing in the lower boundary of the perched saltwater; (b) an amplification of the lateral effect of the highly conductive layer. These limitations decrease the resolution of DC imaging for accurate defining our targets. Because the IP response depends on microgeometry, fluid chemistry and saturation, the 2D IP results demonstrate the suitability of this method to characterize the tidal deposits in the coastal area with a good resolution. In this study, the success of SIP method supports further investigations into studying the hydraulic parameters of tidal deposits in this area. The obtained results during this investigation provide an overview of the coastal aquifer and they can serve as a basis for refining the conceptual model of morphological elements and sedimentary sequences of the coastal tidal flat.     

滦县地磁场井下观测系统抗干扰效果

近距离铁磁性物质干扰(农田耕作、基建等)、直流供电干扰(高压直流输电、地电阻率观测等)是地磁场地面观测的主要干扰方式。在滦县地震台布设地磁场井下观测系统,检验井下观测的抗干扰效果,结果发现,地磁场井下观测可抑制近距离铁磁性物质干扰,但无法抑制直流供电干扰(可通过避让一定距离解决),该分析结果为未来地震监测手段的选取提供技术参考指标和依据。     

综合物探方法在铝土矿勘查中的应用

区域地层为较稳定的沉积岩层,中二叠系梁山组含铝岩系广泛分布于区域上向斜构造内,为铝土矿的找矿奠定了物质基础.区内车盘向斜构造控制着铝土矿的分布,使探矿工作成为可能.由于铝土矿无明显的磁性、密度差异等物性特征,铝土矿的地球物理勘探仍然是盲区,单一的地球物理方法难以提取出其异常特征.本文阐述了综合利用音频大地电磁法(AMT)和直流电测深(DC)两种物探方法在沉积岩中查找铝土矿,并获得良好的效果.音频大地电磁测深法可以提供大尺度的岩性分布特征.直流电测深能准确的划分出地下介质层位的电阻率分布,且对细小异常的分辨能力强.据地质资料推测,目标矿层铝土矿位于石灰岩与砂泥岩分界面中.AMT能划分出分界面,但难以识别覆存其中的目标矿体,利用DC可查出该异常.结合钻孔资料可以验证异常,同时对比进行物探数据的处理和解释,这样的组合勘探方法对于深部找矿及后期的开发是行之有效的.     

城市活断层探测中电磁噪音和环境干扰对浅层电磁方法的影响

城市噪音主要指由于城市设施和人类活动 (如电网、管线、通讯设备、建筑物、交通等 )对电磁信号的干扰。克服城市噪音、提高信噪比是保证城市活断层的电磁探测能否成功的关键。因而 ,在福州市活断层的电磁试验探测中 ,噪音观测与分析被列为试验的一项重要内容 ① 。文中将主要讨论不同类型的噪音对高密度电法 (DC)和瞬变电磁法 (TEM)的影响 ,并提出一些实用的解决办法。试验表明 ,影响高密度电法观测质量的主要因素是地下噪音源 ,如地下管道、地下电缆等 ,地面以上噪音的影响相对较小。对于瞬变电磁法 ,不管地面以下还是地面以上的噪音源都将产生影响。但是 ,只要仪器的性能好 ,对不同噪音源采取相应的措施 ,可以获得满意的结果     

Numerical Simulation of the Effect of a DC Electric Field on Seismic Wave Propagation with the Pseudospectral Time Domain Method

We have modeled the effect of a direct current (DC) electric field on the propagation of seismic waves by the pseudospectral time domain (PSTD) method, based on a set of governing equations for the poroelastic media. This study belongs to the more general term of the seismoelectric coupling effect. The set of physical equations consists of the poroelastodynamic equations for the seismic waves and the Maxwell's equations for the electromagnetic waves; the magnitude of the seismoelectric coupling effect is characterized by the charge density, the electric conductivity, the Onsager coefficient, a function of the dielectric permittivity, the fluid viscosity, and the zeta potential. The poroelastodynamic vibration of a solid matrix generates an electric oscillation with the form of streaming current via the fluctuation of pore pressure. Meanwhile, fluctuating pore pressure also causes oscillatory variation of the electric resistivity of the solid matrix. The simulated poroelastic wave propagation and electric field variation with an existing background DC electric field are compared with the results of a physical experiment carried out in an oilfield. The results show that the DC electric field can significantly affect the propagating elastic energy through the seismoelectric coupling in a wide range of the seismic frequency band.