基于遗传神经网络的大地电磁反演

为进一步提高大地电磁非线性反演的稳定性、运算效率及准确度,将遗传神经网络算法引入大地电磁反演.首先针对大地电磁二维地电模型建立BP（Back Propagation）神经网络基本框架进行学习训练,网络输入为已知地电模型的视电阻率参数,输出为该地电模型参数;再利用遗传算法对神经网络学习训练过程进行优化,计算出多种地电模型网络连接权值和阈值的最优解;最后将最优连接权值和阈值对未知模型进行反演测试,网络输入为未知地电模型的视电阻率参数,输出为该地电模型参数.模型实验表明：遗传神经网络算法充分结合了遗传算法的全局寻优性和神经网络的局部寻优性,相比单一神经网络算法,在网络学习训练中提高了解的收敛成功率和计算速度,在反演测试中能更准确地逼近真实模型.将遗传神经网络算法与最小二乘正则化反演进行对比,理论模型和实测数据都验证了遗传神经网络算法在大地电磁反演中的可行性和有效性.     

带相位信息的一维大地电磁曲线对比反演法

在大地电磁连续介质曲线对比法反演的基础上 ,引入了相位信息 ,即在反演的模型视电阻率的修改量中加入了相位改正量 ,从而提高了收敛速度 ,稳定了反演过程。与其他方法的模型研究比较 ,加入相位信息的曲线对比法反演的结果 ,能够较真实地反映模型的电性分布 ,并可为多维反演提供初始     

大地电磁相位资料的一维反演

近年来,由于仪器研制水平的提高,特别是高精度数字仪器的出现,使相位资料达到了与视电阻率相同的精度,从而可用相位资料进行反演解释,获得地球的层状电性模型。相位曲线变化幅差大,求得的地电模型参数将更可靠。 本文采用马奎特法反演相位资料。选择的目标函数为     

前兆地电阻率资料处理的剥皮法

在台址下电阻率分层均匀的前提下，针对我国绝大多数前兆地电台单极距观测的事实与现状，提出一种试图从这种观测所获得的视电阻率时间序列中识别电阻率表层变化和深部变化的等效模型方法．分析表明，在保持布极方式和观测结果不变的前提下，任意多层的地电结构都可以视为一个二层地电结构——等效模型．本文中，等效模型上层的地电参数（电阻率和厚度）与实际模型第一层相同，下层的地电参数（电阻率）则根据等效要求，由实际模型第二层及其以下各层的地电参数经计算得到，因此，等效模型的下层包含了实际模型第二层及其以下各层的综合影响．初步研究表明，这种方法是可行的．      

CSAMT单分量数据解释方法

可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.     

带相位信息的一维大地电磁曲线对比反演法

在大地电磁连续介质曲线对比法反演的基础上，引入了相位信息，即在反演的模型视电阻率的修改量中加入了相位改正量，从而提高了收敛速度，稳定了反演过程。与其他方法的模型研究比较，加入相位信息的曲线对比法反演的结果，能够较真实地反映模型的电性分布，并可为多维反演提供初始模型。     

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

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

方位梯度电极系的正反演

针对环井周的非均质地层模型,利用三维有限元方法模拟了方位梯度电极系的测井响应．结果表明,在非均质程度比较小的情况下,各个方位上的响应之间差别不大,与真实的电阻率分布差别也不大;反之,则差别很大．其特点是所测得的响应远远低于较高的真电阻率值,而接近较低的真电阻率值．因此,有必要对测井响应进行校正．考虑到传统的校正方法的弊端,本文利用非线性反演方法求取真实的电阻率分布,数值模拟的结果表明了此方法的可行性．     

线源三维井地电法异常边界增强及地形校正技术

基于非结构网格有限元方法开展了三维复杂地电模型的线源井地电法的高效正演模拟研究,探讨了通过求取电场响应导数来刻画目标体边界范围、采用差异场地形校正技术来消除地形影响等措施对井地电法成像的效果和精度的影响。并通过对比与解析解,验证了本文数值解算法的有效性。模型计算结果表明：积水巷道的空间位置和走向均会引起视电阻率的显著变化,视电阻率变化率的极值准确且清晰地指示了巷道边界的位置;电位的归一化总水平导数极大地提高了井地电法对目标体复杂边界位置的识别能力;地形对井地电场分布的影响也很大,其视电阻率响应与地形形状近似呈对称关系,利用差异场技术能有效地削弱地形对井地电法高精度成像的影响。     

双感应测井资料的快速近似迭代反演

本文给出一种双感应资料快速近似迭代反演技术. 首先建立Fréchet导数的快速算法,保证在反演过程中能够同时获得测井响应相对于地层电阻率和层界面的偏导数,并给出用规范化处理与奇异值分解技术进行迭代反演的具体过程. 为了对理论模拟和井场实际资料进行反演,利用综合分层技术从双感应曲线中提取层界面初始位置和地层电阻率初值,通过单独迭代反演中感应资料,修改层界面和地层电阻率实现中感应资料的最佳拟合,得到探测深度相对较浅的地层电阻率,然后固定层界面位置,再迭代反演深感应资料,得到另一组探测深度相对较深的地层电阻率. 理论和实际资料处理结果证明,两个不同探测深度的电阻率反演结果的相对大小能够准确地反映地层真实的侵入特征.此外,由于深感应仪器具有较深的探测深度,不论在高侵或低侵地层上,深感应反演结果与地层原状电阻率的差异大大小于视电阻率的差异,所以利用反演结果也能得到更好的地层原状电阻率的估计值.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.     

Assessment of the reliability of 2D inversion of apparent resistivity data

The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity.     

A new procedure for the interpretation of VES data: 1.5-D simultaneous inversion method

In the paper the principles of a new inversion method as well as the results of its application are described by using synthetic and field VES data. The basic idea of the method says that the horizontal changes in the layer thicknesses and the resistivities of the 2-D geological structure can be described by (expanding in series) functions of one variable. The coefficients of the functions are determined from the VES data by simultaneous inversion method using a least-squares technique. The local thicknesses and the resistivities of the geological structure are calculated from the coefficients from point to point along the profile. Using this method, the equivalence can considerably be reduced, and by using the functions chosen on the basis of a priori knowledge, the whole geological model can be determined by this inversion method. In the inversion the local 1-D forward modelling is used and this approximation is often applied in the VES practice. To qualify the results of the inversion, the correlation matrix for the coefficients are calculated and the error values for the local model parameters are introduced.     

THE AUTOMATIC FITTING OF A RESISTIVITY SOUNDING BY A GEOMETRICAL PROGRESSION OF DEPTHS *

The problem of interpretation of geoelectrical resistivity soundings has been studied and a compromise method of interpretation developed. A simple depth sequence of increasing bed thickness is taken and resistivities assigned to the various layers so that a good fit is obtained to the observed apparent resistivity curve. An initial estimate of the resistivities is made by a direct method but the model is then altered to give an improved fit by the indirect method of curve matching. The whole method is computerized and input data to the computer program consists of spacings and apparent resistivities. The method is extremely rapid, about one third of the computation time is devoted to obtaining the resistivity transform by the method of convolution. A normal type of solution, incorporating additional information if necessary, can also be obtained.     

大地电磁视电阻率和阻抗相位资料的联合解释

本文叙述了大地电磁视电阻率和阻抗相位资料联合解释的原理与计算方法。通过具体算例的分析,表明该方法是正确有效的,而且在缩小多解性方面,要比单纯的视电阻率资料反演法或阻抗相位资料反演法优越     

THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA

On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained.     

CONTRIBUTION A L'ETUDE DE L'EQUIVALENCE EN PROSPECTION ELECTRIQUE (COURANT CONTINU ET MAGNETOTELLURIQUE)*

In all electrical methods (EM, MT, DC), a theoretical one-to-one relation links the apparent resistivity diagram at a station with the local distribution of the resistivities in a tabular earth. Unfortunately, the actual field diagrams derive from a few measurements of limited accuracy and the one-to-one relation vanishes. Calculating the distribution of the earth resistivities from the field diagrams is thus an undetermined problem, which experience has shown to be partly controlled by approximate laws of “equivalence”. The use of these rough rules can produce only a very restricted number of the resistivity distributions that correspond to calculated diagrams which do not differ from a given field curve by more than a tolerable discrepancy. On the other hand, the full set of such “equivalent” solutions is given automatically by a computer program. Tests made on field examples prove that a few geological data generally reduce the wide range of undetermination.     

大地电磁测深资料的联合反演研究

本文从大地电磁扬的色散关系出发,阐述了进行大地电磁资料一维联合反演的物理意义。论证了对于以色散关系相联系的一对大地电磁响应函数,特别是阻抗实部和虚部视电阻率进行联合反演的必要性和可能的应用前景。对比研究了大地电磁阻抗各单参量及其联合反演的效果。结果表明,在加快反演迭代速度,限定反演解的非唯一性、保证反演迭代的稳定进行等方面,联合反演比单参量反演具有明显的优越性;在大多数情况下,阻抗实部视电阻率和阻抗虚部视电阻率联合反演的效果要比通常所采用的常规视电阻率和阻抗相位的联合反演效果好。实际大地电磁资料的试验表明,各种不同平均阻抗的利用和综合反演解释将有助于获得更为可靠的结果。