自然电场法勘探火山岩地区地下水的研究

多个火山岩地区找水的勘探实例表明,地表观测到的自然电位异常（简称SP）与地下非饱和区厚度相关。本文认为这种自然电位异常是由于雨水垂直向下流经非饱和区（渗流区）到达地下不透水界面引起,并提出产生自然电位异常的地下极化模型,进而根据静电场理论建立了水文地质模型与自电异常之间的理论关系。这一理论能够实现把一种波数域反演技术引进自然电位异常的反演,以便计算地下水的位置和深度。模型算例和勘探实例表明这一方法可以用于勘探火山岩地区地下水。     

自然电位数据阻尼最小二乘逐步解析反演研究——以桂林市寨底地下河为例

本文对利用自然电位法探测岩溶地下河做了探索性研究.选择桂林市寨底地下河发育良好的区域,通过阵列式连续时间观测方式获得了随时间连续变化的自然电位数据.提出自然电位数据阻尼最小二乘多参数逐步解析反演方法,并对获得的数据进行反演,以此提取数据中极化源的信息,包括极化源位置、极化角度和极化强度等.在所有时间道的数据反演结果中,将观测范围内极化源的空间位置等参数绘制成相应图件,与观测区地下河资料进行对比研究.研究结果表明:该反演方法能够有效提取自然电位数据中隐含的极化源信息;极化源位置可作为判断含水构造空间位置的依据,极化角度和极化强度可用来区别含水构造类型;多参数综合解释可有效定位地下河的平面位置.     

基于最小二乘正则化的自然电场场源反演成像

自然电场对多种地下污染物敏感,开展相应的场源反演成像研究可帮助有效监测污染源活动.把目标区域内的多种自然电场场源进行统一处理,利用有限单元构建二维地电模型,实现自然电位正演的高精度计算.采用最小二乘正则化反演迭代,实现自然电场场源的二维反演.在设计二维自然电场场源反演算法的基础上,利用加入高斯噪声的合成模拟数据,对反演算法进行测试.测试结果表明:基于有限单元二维模型的最小二乘正则化反演算法能有效地实现对自然电场场源的反演成像,并能准确地确定自然电位异常源的位置和深度,且算法收敛快速稳定,具有一定的抗噪声能力.     

长白山天池火山区大地电磁测深初步观测

在长白山火山区沿着两条剖面进行了１５个测点的大地电磁测深观测，对地磁感应矢量和地下电性分布维数特征进行了分析，并使用阻抗张量分解和二维自动反演等新技术对资料进行了解释。结果表明，火山区及其周围地区的地下电性结构具有二维特征，但在火山口及其附近地下存在局部三维异常体，异常体埋深约１２ｋｍ，初步推测可能是火山岩浆囊     

不同激励源井地电位成像技术三维正反演方法研究（英文）

井地电位成像是通过套管向井中供电或将电源放在井中,在地表观测电位异常的一项技术,其供电源有线源和点源两种类型。为了研究这两种电源对地下异常体产生的电位异常特征,本文针对不同激励源,采用有限差分方法进行数值模拟研究,在线性方程组求解电位时引入不完全Cholesky共轭梯度(ICCG)迭代方法,分别实现了点源和线源井地电位成像技术的三维正演。最后,基于阻尼最小二乘法实现了井地电位成像技术的电阻率三维反演。设计不同地电模型分别进行正演和反演试算,正演结果表明,供电电源的类型不同,异常体在地表的电位异常特征也不同;反演结果表明,低阻体的反演结果要好于高阻体,点源置于异常体下方时反演的电阻率对异常体边界的识别比线源更加准确。     

利用Cole-Cole模型组合得到SIP真参数的联合频谱最优化反演

利用Cole-Cole模型组合和稀释系数理论，可以描述测量得到的视复电阻率频谱与地下多个电性体SIP真参数之间的正演关系。本利用此正演关系进行反演试算，并比较了几种Cole-Cole模型反演方案的效果；提出固定极化目标体和围岩体的极化率同时反演其它各SIP真参数的反演方案。结果表明在固定极化率为真值时，利用此反演方案可以稳定、准确、快速地得出地下电性体SIP真参数。     

基于COMSOL Multiphysics的长白山天池多物理场耦合响应研究

采用COMSOL Multiphysics建立长白山天池2D多物理场模型,模拟长白山天池流场、自然电位分布及其磁场响应.研究结果表明:天池和大气降水入渗补给地下水,地下水在重力作用下向下流动时,产生电位差,形成自然电场,进而产生磁场.在火山通道外侧,电流密度方向由地表指向地下,通道周边地区磁通量密度呈现高值;在火山通道/断层破碎带内,向下运移的地下水遇到致密岩体后向上补给天池或涌出地表成泉,该过程产生向上运移的电流,电流密度方向由地下指向地表,磁通量密度呈现低值.研究长白山天池多物理场响应,有助于进一步理解长白山天池地区地下水循环过程及电磁异常出现的机理,亦可用于指导火山地区勘查、监测与开发.     

亚洲地区视磁化强度分布

利用等效偶极源方法对ＭＡＧＳＡＴ卫星地磁总强度异常数据进行分析，得到了亚洲地区地壳内部的视磁化强度分布图．该分布图基本反映了中亚和东亚地区地球深部磁性结构，表明视磁化强度异常分布与地质构造具有较好的相关性．亚洲地区各地块上（哈萨克斯坦地块、塔里木地块、扬子地块、印度地块、中朝地块、青藏地块、印度支那地块和日本海地块等）都有独立的异常显示；盆地、平原和海盆为较低磁化强度区；高原和海岭为较高磁化强度区．同时，异常边界与断裂带、一些山脉和海沟相吻合．中国中部的南北构造带两侧异常差异明显：青藏高原与印度地区交界处异常梯度大，其两边覆盖着明显不同的磁化强度异常；郯庐断裂带为正负磁化强度异常分界，这条分界线还向西南方向延伸，与三江断裂带汇合．此外，在异常图上中国东南沿海的褶皱带界线明显．     

大斜度井和水平井井地三维电阻率数值模拟和联合反演

井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向.     

三维电阻率测深数据Zohdy近似反演方法

Zohdy方法近似反演三维电阻率测深数据。正演计算采用有限单元法。反演初始模型由测量视电阻率数据给出。通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率．为了解决任意电极距测深数据的反演,采用大、小双网格剖分。大网格反映地下电性分布情况。小网格用于实际有限元正演计算．在电阻率调整公式中加入一个迭代系数,能够加快收敛速度．并对加5％随机噪声的模型理论视电阻率测深断面数据进行反演,得到的电阻率分布与模型电阻率基本一致．     

Hydrothermal system beneath Aso volcano as inferred from self-potential mapping and resistivity structure

We conducted self-potential (SP) surveys sequentially from part to part over the central cones of Aso volcano since August 1998 by December 2001. The compiled SP map revealed large SP anomalies on the central cones. The main feature of the SP map is a ‘W-shaped’ profile along the NS-transect over the central cones. It is probable that this characteristic SP profile is produced by the combination of hydrothermal upwelling in the middle and topographic effect. A positive anomaly showing a large concentric pattern has appeared after correcting the topographic effect.To evaluate this SP anomaly, we implemented a numerical code that calculates electric potential produced by arbitrarily positioned current sources and sinks in any three-dimensional resistivity structure. A layered structure obtained from a time-domain electromagnetic (TDEM) field experiment was used for the resistivity model. The estimated current source is 300 A, being located in a conductive layer around the sea level. Meanwhile, sinks were estimated to sit on a circular area corresponding to the marginal part of the conductive layer.Water and heat budget study gives a lower limit of water mass transfer from depth to the bottom of the crater lake of Nakadake. This value was used to estimate the equivalent current in either case of electro-kinetic (EK) [Mizutani, H., Ishido, T., 1976. A new interpretation of magnetic field variation associated with the Matsushiro earthquakes, J. Geomag. Geoelectr., 28, 179–188.] or rapid fluid disruption (RFD) process [Johnston, M.J.S., Byerlee, J.D., Lockner, D., 2001. Rapid fluid disruption: A source for self-potential anomalies on volcanoes, J. Geophys. Res. 106(B3), 4327-4335.]. This comparison suggests that the former process is preferable to explain the observed SP anomaly. From these results we infer a large-scale hydrothermal system beneath the central cones of Aso volcano, in which the fluid flow initiates from the surrounding area, converging to the central vent to transport the heat and materials up to the crater lake of Nakadake through a vapor-filled conduit.     

Application of Normalized Full Gradient Method to Self Potential (SP) Data

Recently, Normalized Full Gradient (NFG) method has widespread applications to natural potential fields, especially in gravity and magnetic. In this study, usage of NFG in Self-Potential (SP) data evaluation is tested. Results are compared to other SP interpretation methods. The NFG method is applied to synthetic and field SP data. As a consequence of application of the method to the anomalies of spherical, cylindrical and vertical sheet models, whose theoretical structures are explicit, the structures were found very close to their actual locations. In order to see the capability of the method in detecting the number of sources, NFG method was applied to different spherical models at different depths and locations. The least-squares inverse solution was applied to the same models and NFG method was found more powerful in detecting model structure. Sensitivity of NFG method for application to noisy data is also tested. An anomaly is generated by adding a random noise to two close sphere SP anomalies. The method seems to work for the two close spheres at high S/N ratio. Then, NFG method was applied to two field examples. The first one is the cross section taken from the SP anomaly map of the Ergani-Süleymanköy (Turkey) copper mine. The depth of the mineral deposit at that site was found about 38 m from the ground level. This result is well matched to previous studies. NFG was also applied to SP data from Seferihisar Izmir (Western Turkey) geothermal field and the location of the point source was determined. The field data of this site have already been modeled by the thermoelectric source (coupling) solution method. When these two methods are compared, they seem to support each other. It is concluded that the NFG method works perfectly when the structure model is simple. It is observed that natural potential sources close to earth’s surface are identified by the method more accurately at greater harmonics, while deep sources are identified at lesser harmonics. It produces reasonable results for noisy multi-source models than the other parameter identification methods (inverse solution, power spectrum, etc.).     

Thermal precursors in satellite images of the 1999 eruption of Shishaldin Volcano

Shishaldin Volcano, Unimak Island Alaska, began showing signs of thermal unrest in satellite images on 9 February 1999. A thermal anomaly and small steam plume were detected at the summit of the volcano in short-wave thermal infrared AVHRR (advanced very high resolution radiometer) satellite data. This was followed by over 2 months of changes in the observed thermal character of the volcano. Initially, the thermal anomaly was only visible when the satellite passed nearly directly over the volcano, suggesting a hot source deep in the central crater obscured from more oblique satellite passes. The "zenith angle" needed to see the anomaly increased with time, presumably as the thermal source rose within the conduit. Based on this change, an ascent rate of ca. 14 m per day for the thermal source was estimated, until it reached the summit on around 21 March. It is thought that Strombolian activity began around this time. The precursory activity culminated in a sub-Plinian eruption on 19 April, ejecting ash to over 45,000 ft. (13,700 m). The thermal energy output through the precursory period was calculated based on geometric constraints unique to Shishaldin. These calculations show fluctuations that can be tied to changes in the eruptive character inferred from seismic records and later geologic studies. The remote location of this volcano made satellite images a necessary observation tool for this eruption. To date, this is the longest thermal precursory activity preceding a sub-Plinian eruption recorded by satellite images in the region. This type of thermal monitoring of remote volcanoes is central in the efforts of the Alaska Volcano Observatory to provide timely warnings of volcanic eruption, and mitigate their associated hazards to air-traffic and local residents.     

Quantitative Interpretation of Self-Potential Anomalies of Some Simple Geometric Bodies

We have developed a new numerical method to determine the shape (shape factor), depth, polarization angle, and electric dipole moment of a buried structure from residual self-potential (SP) anomalies. The method is based on defining the anomaly value at the origin and four characteristic points and their corresponding distances on the anomaly profile. The problem of shape determination from residual SP anomaly has been transformed into the problem of finding a solution to a nonlinear equation of the form q = f (q). Knowing the shape, the depth, polarization angle and the electric dipole moment are determined individually using three linear equations. Formulas have been derived for spheres and cylinders. By using all possible combinations of the four characteristic points and their corresponding distances, a procedure is developed for automated determination of the best-fit-model parameters of the buried structure from SP anomalies. The method was applied to synthetic data with 5% random errors and tested on a field example from Colorado. In both cases, the model parameters obtained by the present method, particularly the shape and depth of the buried structures are found in good agreement with the actual ones. The present method has the capability of avoiding highly noisy data points and enforcing the incorporation of points of the least random errors to enhance the interpretation results.     

中国东北地区远震S波走时层析成像研究

利用中国东北流动和固定台网的234个宽频带地震仪在2009年6月-2011年5月所记录的远震波形数据,采用波形相关方法拾取了10301个有效的S震相相对走时残差数据,进一步采用两种射线走时层析成像的方法,反演获取了研究区下方深达800 km的S波速度结构,不同射线走时层析成像方法的结果对本区的S波速度异常结构起到一定的约束作用.S波成像结果与我们先前开展的P波成像研究结果整体相似：在长白山下方发现有一个高速异常结构,这可能就是俯冲到欧亚大陆板块下方的太平洋板块,由于板块的部分下沉,使得板块的形状并没有呈现出明显的板片状.长白山、阿尔山、五大连池火山下方都有低速异常体,长白山和阿尔山下的低速异常向下延伸至地幔转换带附近,可能与其上部的火山形成有关.五大连池火山下方的低速异常向下延伸至200 km左右,不同埋深的低速异常结构可能意味着五大连池与长白山和阿尔山有着不同的成因.松辽盆地呈现以高速异常为主导、高低速异常混合分布的特性,暗示松辽盆地岩石圈地幔可能遭受了改造与破坏,可能有岩石圈拆沉的过程,盆地南部下方的低速异常与长白山和阿尔山下的低速异常有连通性,可能是下地幔热物质上涌的一个通道.S波和P波相似的成像结果从另一个角度展示了中国东北地区的上地幔速度异常结构,对我们认识此区的地下结构提供了重要的约束.     

A Hybrid Model for the Summit Region of Merapi Volcano,Java, Indonesia,Derived from Gravity Changes and Deformation Measured between 2000 and 2002

Three-dimensional displacements and gravity changes were measured at 20 points between 2000 and 2002 at Merapi volcano, Java, Indonesia. Former models focused on the modeling of a single magmatic source located in the summit region of the volcano. Such models do not fit to our measurements between 2000 and 2002. A new hybrid model approach is developed consisting of an elastic-gravitational source described by a mass and energy term as well as a combined dip-slip/strike-slip fault zone in the summit region. Both nonlinear optimization problems, given by the common inversion of three-dimensional displacements and gravity changes, are solved by applying a genetic algorithm. The hybrid model fits the measurements accurately, tested by Fisher test statistics. Furthermore, our model for Merapi volcano confirms previous structural models for this region so that the new model is statistically proven as well as physically reliable.     

Estimation of depth and magnetization of magnetic sources from magnetic data: The linearized continuous inverse problem for 21/2D structures

The estimation of the depth to the top and bottom of a magnetic source from magnetic data defines a nonlinear inverse problem, while the evaluation of the distribution of magnetization determines a linear inverse problem. In this paper, these interpretation problems are resolved in the continuous case of 2_1/2D magnetized bodies with lateral magnetization variations. A formulation of the magnetic problem accounting for different directions of remanent and total magnetization vectors and including a more general definition of apparent susceptibility is presented. Differences between 2D and 2_1/2D formulations are stressed, as regards the anomaly amplitude, shape and zero-level.In order to utilize well-known continuous linear inverse methods, Fréchet derivatives of the data functionals with respect to the depth of the source top and bottom, are analytically described. Thus, using the spectral expansion inverse method (Parker, 1977) and linearizing the problem at several steps of an iterative process, the source depth is obtained within a few iterations, although the starting model is distant from the final solution. The interpretation of an anomaly in the Italian region shows the usefulness of the method.     

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

中国东北地区远震P波走时层析成像研究

利用中国东北流动和固定台网的234个宽频带地震仪记录的远震波形数据,采用波形相关方法拾取了57251个有效相对走时残差数据,进一步采用FMTT(Fast Marching Teleseismic Tomography)层析成像的方法,反演获取了研究区下方深达800 km的P波速度结构.结果显示:在长白山下方发现有一个高速异常结构,这可能就是俯冲到欧亚大陆板块下方的太平洋板块,由于板块的部分下沉,使得板块的形状并没有呈现出明显的板片状.长白山、阿尔山、五大连池火山下方都有低速异常体,长白山和阿尔山下的低速异常向下延伸至地幔转换带,可能与其上部的火山形成有关.五大连池火山下方的低速异常向下延伸至200 km左右,不同埋深的低速异常结构可能意味着五大连池与长白山和阿尔山有着不同的成因.松辽盆地呈现以高速异常为主导高低速异常混合分布的特性,暗示松辽盆地可能有岩石圈拆沉的过程,盆地南部下方的低速异常与长白山和阿尔山下的低速异常有连通性,可能是下地幔热物质上涌的一个通道.