The neural network approximation method for solving multidimensional nonlinear inverse problems of geophysics

The iterative approximation neural network method for solving conditionally well-posed nonlinear inverse problems of geophysics is presented. The method is based on the neural network approximation of the inverse operator. The inverse problem is solved in the class of grid (block) models of the medium on a regularized parameterization grid. The construction principle of this grid relies on using the calculated values of the continuity modulus of the inverse operator and its modifications determining the degree of ambiguity of the solutions. The method provides approximate solutions of inverse problems with the maximal degree of detail given the specified degree of ambiguity with the total number of the sought parameters ~n × 10³ of the medium. The a priori and a posteriori estimates of the degree of ambiguity of the approximated solutions are calculated. The work of the method is illustrated by the example of the three-dimensional (3D) inversion of the synthesized 2D areal geoelectrical (audio magnetotelluric sounding, AMTS) data corresponding to the schematic model of a kimberlite pipe.     

Direct and inverse problems in electromagnetic sounding of three-dimensional heterogeneous medium

The methods for solving three-dimensional (3D) direct and inverse problems of electromagnetic sounding are considered. It is shown that the method of integral equations is an efficient instrument for mathematical modeling of electromagnetic fields in 3D heterogeneous media. Adaptation of the integral equation technique to the solution of inverse problems is demonstrated.     

Magnetotellurics and radio-wave interference sounding

The plane harmonic electromagnetic fields are considered in the theory of magnetotelluric methods in the range of frequencies from 0.0001 Hz to 20 kHz. These fields are natural by their origin and contain information on the depths from tens of meters up to 100 km and more. The magnetotelluric soundings, which use the fields of radio stations, expand the frequency band almost up to 1 MHz and make it possible to study the depths from the first few meters. The method of radio-wave interference sounding supplements geoelectric prospecting on plane waves into the range of even higher frequencies (up to 100 MHz). In this case, the conduction and displacement currents become comparable, which makes it possible to distinguish objects both by their electrical conductivity and by their dielectric permittivity. For the two-layered model of a medium, there exist simple kinematic methods of data interpretation of a radio-interferometry sounding. Within multilayer, and especially horizontally heterogeneous, media, methods for solving equations of electrodynamics and inverse problems of geophysics are required. In the present paper, the foundations of the theory of radio-interferometry sounding, the methodology, its role in geoelectric prospecting, and the opportunities for the solution of geological problems are discussed.     

Trefftz method for the solution of three-dimensional forward and inverse problems of geoelectrics

The Trefftz method is applied to the numerical solution of the three-dimensional (3D) forward problem for an electromagnetic field harmonically oscillating with time in a 3D environment, which ensures the possibility to solve a 3D inverse problem. Here, the known idea of simultaneous joint interpretation of the tangential components of electric and magnetic fields measured on the Earth's surface, which was suggested by A.N. Tikhonov for the development of the electromagnetic sounding method, is employed in the context of the numerical-analytical representation of the forward problem solution in accordance with the Trefftz method. Application of this method allows enables one to generalize the basic plane-stratified vertically 1D model of the medium by allowing for lateral variations in electric conductivity.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

Geoelectric section of the Central Tien Shan: Sequential inversion of the magnetovariational and magnetotelluric data along the Naryn Line

The paper presents the results of 2D inversion of deep magnetotelluric (MT) and magnetovariational (MV) soundings along the Naryn Line. The method of partial (sequential) inversions is used. According to this method, at the first stage, magnetovariation responses are used for the localization of deep anomalies of electrical conductivity, and then the magnetotelluric sounding data are invoked to refine the structure of the host medium and the structural details in the upper part of the section. It is shown that this approach enables one to estimate the informativeness of separate components of the electromagnetic field, to reduce the distorting influence of the near-surface geoelectric inhomogeneities, and to increase the stability of the final solution of the inverse problem.     

大地电磁有限差分数值解对比

基于网格精度、形成系统方程的方式、边界条件以及预条件线性算子等,文中对大地电磁（MT）有限差分数值解作了对比.对不同网格剖分方式下的三个均匀半空间模型的一维MT响应对比显示,在降低首层厚度的同时保持层间厚度变化在合理范围可以同时提高主场和辅助场的精度.在利用正常中心网格法（主场和辅助场都定义在单元顶面的中心）计算二维（2-D）TM模式响应时,应该从Maxwell一次差分方程开始组建二次差分方程,这样可以更充分考虑模型电阻率的变化.在对边界值如何影响数值解的测试表明,仅仅提高一维（1-D）边界值的精度对提高2-D MT有限差分数值解的精度是有限的.线性算子对提高MT解的效率十分重要,简单的对比进一步表明合适的预条件再配合好的线性算子（如文中求解2-D MT时所采用的DILU-BICGSTAB方法）不仅可以加速收敛,而且可以降低迭代次数.     

Neural network inversion of data in classes of parametrized geoelectric sections

A method of approximate magnetotelluric sounding (MTS) data inversion is developed on the basis of the representation of the inverse operator by an artificial neural network in classes of geoelectric structures. A methodology of the neural network inversion of magnetotelluric data is proposed for a family of classes of geoelectric structures and the uncertainty of the inferred results is estimated. A neural network algorithm of MTS data inversion is tested using synthetic 2-D data.     

二维水平电偶极变频测深阻抗视电阻率的有限元正演计算

在二维地电断面情况下,电偶源变频测深阻抗视电阻率的有限元计算,存在着两个困难:一个是对源的处理;另一个是边界条件的确定。本文提出了解决这两个问题的方法,构成了含源波动方程的边值问题,用有限元法求出了相应的电磁场分量,从而计算出二维地电断面的阻抗视电阻率。对各种地电断面计算的结果表明,计算方法是正确的,可对实际地电断面进行有限元模拟,求得其电磁响应。     

Computer modeling of geoelectromagnetic fields in three-dimensional media by the finite element method

This work describes the basic approaches to the solution of three-dimensional (3D) problems of geoelectromagnetism with the use of the finite element method and the possibilities of the GeoEM program complex for its implementation. The methods of modeling geoelectromagnetic fields for the most widely known types of controllable sources and the mathematical apparatus for the solution of problems of magnetotelluric soundings are considered. Examples of the calculations of 3D fields for the horizontal and vertical electric lines in planning electric exploration works on the shelf are presented, as well as an example of the 3D interpretation of array transient electromagnetic sounding data during the search for deep-seated target objects in conditions of a mostly heterogeneous upper part of the section.     

一种基于亥姆霍兹分解的大地电磁测深有限元正演预条件解法

本研究针对大地电磁测深法有限元数值模拟中,迭代法求解线性方程组效率较低的问题,利用亥姆霍兹分解原理,将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题,并在四面体非结构化棱边元离散的情况下,借助节点元辅助网格离散上述预条件问题,进一步利用代数多重网格方法（AMG）实施求解,最终实现预条件算法.利用经典的COMMEMI理论模型进行试算并与前人的积分方程解进行对比,验证了本文数值模拟程序与预条件方法的正确性和可靠性.此外,利用不同自由度规模的实验模型对这一预条件算法的效率进行了测试.结果表明,这一算法可以有效地提升大地电磁测深法棱边有限元数值模拟迭代法的收敛性,计算效率较通用的不完全LU分解预条件算法明显更高;在较大自由度网格（>1000万）数值模拟计算中,其算法效率及内存占用相对直接解法有较大优势,也使小型工作站上利用较大自由度的有限元网格进行大地电磁测深数值模拟计算成为可能.     

实数编码遗传算法在大地电磁测深二维反演中的应用

作为全局非线性优化的新方法之一的遗传算法,近年来已从生物工程流行到大地电磁测深资料解释中.然而,大地电磁反演问题具有不适定性,解的非唯一性.通过结合求解不适定问题的Tikhonov正则化方法,本文采用实数编码遗传算法求解大地电磁二维反演问题.此算法在构建目标函数时引入正则化的思想,利用遗传算法求解最优化问题.常规的基于局部线性化的最优化反演方法易使解陷入局部极小值,而且严重的依赖初始模型的选择.与传统线性化的迭代反演方法相比,实数编码遗传算法能够克服传统方法的不足且能获得更好的反演结果.通过对大地电磁测深理论模型进行计算,结果表明:该算法具有收敛速度快、解的精度高和避免出现早熟等优点,可用于大地电磁资料解释.     

Combined method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-approximations in solving the problems of geophysics and geomorphology

The interrelation between different modifications of the method of linear integral representation is studied. Combined approximations of the topography and geopotential fields enable more refined tuning of the method in solving inverse problems of geophysics and geomorphology and provide a more complete allowance for the a priori information about the surface elevation data and elements of anomalous fields. A technique for finding the numerical solution for the inverse problem for determining the mass distributions equivalent in terms of the external field is presented. The results of the mathematical experiment are discussed.     

地震偏移反演成像的迭代正则化方法研究

利用伴随算子L^*,直接的偏移方法通常导致一个低分辨率或模糊的地震成像.线性化偏移反演方法需求解一个最小二乘问题.但直接的最小二乘方法的数值不稳定,为目视解译带来困难.本文建立约束正则化数学模型,研究了地震偏移反演成像问题的迭代正则化求解方法.首先对最小二乘问题施加正则化约束,接着利用梯度迭代法求解反演成像问题,特别是提出了共轭梯度方法的混合实现技巧.为了表征该方法的可实际利用性,分别对一维,二维和三维地震模型进行了数值模拟.结果表明该正则偏移反演成像方法是有效的,对于实际的地震成像问题有着良好的应用前景.     

Inverse problems of gravity prospecting as a decision-making problem under uncertainty and risk

A new class of algorithms for solving the inverse problems of gravity prospecting is considered. The best interpretation is selected from the set Q of the admissible versions by the optimality criteria that are borrowed from the solution-making theory and adapted for the geophysical problems. The concept of retrieving the information about the sources of gravity anomalies, which treats the result of the interpretation as a set of locally optimal solutions of the inverse problem but not as a single globally optimal solution is discussed. The locally optimal solutions of the inverse problem are sort of singularity points of set Q. They are preferable to the other admissible solutions by a certain criterion formulated in terms of the geologically important information about the anomalous bodies. The admissible versions of the interpretation of the gravimetry data that meet the criteria of the decision-making theory are the primary candidates for the singularity points. The results of the numerical calculations are presented. The set of the admissible solutions from which the locally optimal versions of interpretation are selected is formed by the modifications of the assembly method developed by V.N. Strakhov.     

Forward computation of the magnetic field of a 3D body with arbitrary boundary and continually varying magnetization

The forward computation of the gravitational and magnetic fields due to a 3D body with an arbitrary boundary and continually varying density or magnetization is an important problem in gravitational and magnetic prospecting. In order to solve the inverse problem for the arbitrary components of the gravitational and magnetic anomalies due to an arbitrary 3D body under complex conditions, including an uneven observation surface, the existence of background anomalies and very little or no a priori information, we used a spherical coordinate system to systematically investigate forward methods for such anomalies and developed a series of universal spherical harmonic expansions of gravitational and magnetic fields. For the case of a 3D body with an arbitrary boundary and continually varying magnetization, we have also given the surface integral expressions for the common spherical harmonic coefficients in the expansion of the magnetic field due to the body, and a very precise numerical integral algorithm to calculate them. Thus a simple and effective method of solving the forward problem for magnetic fields due to 3D bodies of this kind has been found, and in this way a foundation is laid for solving the inverse problem of these magnetic fields. In addition, by replacing the parameters and unit vectors in the spherical harmonic expansion of a magnetic field by gravitational parameters and a downward unit vector, we have also derived a forward method for the gravitational field (similar to that for the magnetic case) of a 3D body with an arbitrary boundary and continually varying density.     

Delineation of layer boundaries from smooth models obtained from the geoelectrical sounding data inversion

Two techniques have been presented for the delineation of boundaries from smooth models obtained by smooth inversion techniques of geoelectrical sounding data, such as straightforward inversion scheme, Occam’s and Zohdy’s methods. The smooth model consists of a large number of equally spaced layers, wherein the real geological boundaries are missing. The techniques proposed here suppress the geologically irrelevant boundaries and support the real structural boundaries present in the geoelectrical data. In the first technique, solution of linear inverse problem is improved iteratively through weighted minimum norm inverse, the weight being taken from the current solution. The technique is referred as Iterative Straightforward Inversion Scheme. The second method is analytical, based on the application of smoothing filter, referred in the literature as edge-preserving smoothing. A few examples of theoretical magnetotelluric, dc resistivity and field sounding data have been presented to demonstrate the capabilities of the techniques. The methodologies also reduce the conspicuous oscillations in the smooth solutions caused by the conversion of sharp boundaries to the smooth ones.     

Solution of the one-dimensional inverse magnetotelluric problem

Summary Methods of solving the inverse magnetotelluric problem are compared. Basic relations for Newton's method, the least-squares method and Marquardt's method are presented and the convergence properties of these methods are studied. The high effectiveness of Marquardt's method is demonstrated and its application to practical magnetotelluric data is discussed.