利用积分方程法计算双侧向测井仪的井眼校正曲线

对电测井积分方程中的各积分项进行了理论推导，利用任意闭合曲面对曲面外一点所张立体角为零的特点将电势积分项由对柱体侧面的积分转化为对柱体底面的积分，从而使在纵向采用线性插值条件下给出了电势积分项的解析解．对电流积分项进行了化简，使之由曲面积分化为一维数值积分．利用所得结果计算了双侧向测井仪的井眼校正曲线．     

三维地形频率域人工源电磁场的边界元模拟方法

提出了一种用边界元法计算频率域人工源三维地形电磁场的数值模拟方法.首先用矢量积分理论和电磁场边界条件，将上半空间（空气）和下半空间（地下介质）两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中，假设单元中电磁场为无限大气空间电磁场与地形影响的叠加,并假设地形影响为常项,这样既保证了计算精度又使得计算方法简便.通过分解和计算，每一个矢量面积分方程分解为对应三个坐标方向的三个常量线性方程，这些线性方程组成了对角占优的线性方程组，可用SSOR方法求解.文中给出了垂直磁偶源的垂直磁场地形影响的例子.     

时间二阶积分波场的全波形反演

通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异.     

New digital linear filters for Hankel J0 and J1 transforms

The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J₀ transform (61-point and 120-point operators), and two for the Hankel J₁ transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics.     

Analytical approximations of Arctic basin depths and their spectral analysis

Using the method of linear integral representations, analytical approximations and spectral analysis of Arctic depths are performed on the basis of a large data set. Results of calculations performed with the use of computer technologies developed by the authors are presented.     

Analytical Models of the Physical Fields of the Earth in Regional Version with Ellipticity

Izvestiya, Physics of the Solid Earth - Abstract—The regional version of the method of linear integral representations is used to interpret the data on the anomalous physical fields of the...     

The method of linear integral representations for solving 3D structural inverse problems of potential theory

The algorithm for the solution of a three-dimensional (3D) structural inverse problem in potential theory is described. The algorithm is based on the method of linear integral representations and on the method of extending compacts. The proposed approach is verified on model examples for an anomalous gravity field.     

适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.     

非线性土木结构振动控制的逆系统方法

采用多变量非线性控制的逆系统理论,将多变量、非线性的复杂被控对象精确线性化成伪线性系统,对于该伪线性系统很容易基于线性系统理论设计综合最优控制器。克服了传统的非线性结构近似线性化控制方法不够精确的缺点,实现了非线性结构的动态精确线性化控制。给出了一个仿真实例,将其控制效果与一般的线性和近似线性化最优控制作了比较,达到同样控制效果时本方法所需控制力显著减小。     

A BRIEF INTRODUCTION TO THE NEW METHOD FOR RIVER PROFILE ANALYSIS: Integral Approach

The topography and geomorphology of active orogens result from the interaction of tectonics and climate. In most orogens, a fluvial channel is most sensitive to the coupling between tectonics, lithology, and climate. Meanwhile, the related signals have been recorded by both the drainage geometry and channel longitudinal profile. Thus, how to extract tectonic information from fluvial channels has been a focused issue in geologic and geomorphologic studies. The well known stream-power river incision model bridges the gap between tectonic uplift, river incision and channel profile change, making it possible to retrieve rock uplift pattern from river profiles. In this model, the river incision rate depends on the rock erodibility, contributing drainage area and river gradient. The steady-state form of the river incision model predicts a power-law scaling between the drainage area and channel gradient. Via a linear regression to the log-transformed slope-area data, the slope and intercept are channel concavity and steepness indices, respectively. The concavity relates to lithology, climatic setting and incision process while the channel steepness can be used to map the spatial pattern of rock uplift. For its simple calculation process, the slope-area analysis has been widely used in the study of tectonic geomorphology during past decades. However, to calculate river slope, the coarse channel elevation data must be smoothed, re-sampled, and differentiated without any reasonable smooth window or rigid mathematical fundamentals. One may lose important information and derive stream-power parameters with high uncertainties. In this paper, we introduce the integral approach, a procedure that has been widely used in the latest four years and demonstrated to be a better method for river profile analysis than the traditional slope-area analysis. Via the integration to the steady-state form of the stream-power river incision equation, the river longitudinal profile can be converted into a straight line of which the independent variable is the integral quantity χ with the unit of distance and the dependent variable is the relative channel elevation. We can calculate the linear correlation coefficient between elevation and χ based on a series of concavity values and find the best linear fit to be the reasonable channel concavity index. The slope of the linear fit to the χ value and elevation is simply related to the ratio of the uplift rate to the erodibility. Without calculating channel slope, the integral approach makes up for the drawback of the slope-area analysis. Meanwhile, via the integral approach, a steady-state river profile can be expressed as a continuous function, which can provide theoretical principle for some geomorphic parameters (e.g., slope-length index, hypsometric integral). In addition, we can determine the drainage network migration direction using this method. Therefore, the integral approach can be used as a better method for tectonogeomorphic research.     

Forward plane-wave electromagnetic model in three dimensions using hybrid finite volume–integral equation scheme

We present a concept of the hybrid finite volume–integral equation technique for solving Maxwell's equation in a quasi-static form. The divergence correction was incorporated to improve the convergence and stability of the governing linear system equations which pose a challenge on the discretization of the curl–curl Helmholtz equation. A staggered finite volume approach is applied for discretizing the system of equations on a structured mesh and solved in a secondary field technique. The bi-conjugate gradient stabilizer was utilized with block incomplete lower-upper factorization preconditioner to solve the system of equation. To obtain the electric and magnetic fields at the receivers, we use the integral Green tensor scheme. We verify the strength of our hybrid technique with benchmark models relative to other numerical algorithms. Importantly, from the tested models, our scheme was in close agreement with the semi-analytical solution. It also revealed that the use of a quasi-analytical boundary condition helps to minimize the runtime for the linear system equation. Furthermore, the integral Green tensor approach to compute at the receivers demonstrates better accuracy compared with the conventional interpolation method. This adopted technique can be applied efficiently to the inversion procedure.     

On the problem of determining the vertical gradients of gravity anomalies

Summary Green's theorem on harmonic functions makes it possible to determine the integral relationship between the harmonic function and its derivative with respect to the normal on a closed Lyapunov surface. The conditions of solvability are given by Fredholm's theory of integral equations. The solution for a sphere was presented by Molodenskii[3] and the general solution with the help of Molodenskii's parameter k by Ostach[4]. The present paper indicates a possibility of solving this problem with the help of a system of linear algebraic equations, a simplified modification of the Ostach-Molodenskii solution and, finally, a method, based on Eremeev's solution of the fundamental integral equation[5].     

Combined Method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-Approximations of Increased Dimensionality in Solving the Problems of Geophysics and Geomorphology

The interrelation between different variants of the method of linear integral representations in the spaces of an arbitrary dimension is considered. The combined approximations of the topography and geopotential fields allows the selection of the optimal parameters of the method in solving a wide range of inverse problems in geophysics and geomorphology, as well as a most thorough use of the a priori information about the elevations and elements of anomalous fields. A method for numerically solving an inverse problem on finding the equivalent, in terms of the external field, mass distributions in the ordinary three-dimensional (3D) space and in the four-dimensional (4D) space is described.     

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.     

A METHOD FOR CALCULATING IP ANOMALIES FOR MODELS WITH SURFACE POLARIZATION*

A numerical method is given for calculating resistivity and induced polarization anomalies produced by a surface polarization model. Surface polarization is generated when a purely electronic conductor is located in an electrolyte environment. The system that develops on the boundary between the conductor and the electrolyte is described macroscopically by a net surface charge distribution and an electric double layer. An integral equation is derived for the potential by assuming that the electronic conductor forms an equipotential system and that the polarization impedance across the boundary is linear. The integral equation is solved by means of the method of subsections. As an application some numerical modeling results are presented. The surface impedance values used in calculations are based on laboratory measurements that are briefly described. Implications of the results for scale modeling are discussed.     

COMPUTATION OF TYPE CURVES FOR ELECTROMAGNETIC DEPTH SOUNDING WITH A HORIZONTAL TRANSMITTING COIL BY MEANS OF A DIGITAL LINEAR FILTER*

The computation method described in this paper is based on the existence of a linear relationship between the mutual coupling ratio and the kernel function in the integral expression for it. Accordingly, the mutual coupling ratio can be determined by first computing sample values of the kernel function and then subjecting these to a digital linear filter. In the present paper the appropriate sampling distance is determined and the values of the digital filter coefficients are computed, both for electromagnetic sounding with horizontal coils and for electromagnetic sounding with perpendicular coils.     

一种新的三维大地电磁积分方程正演方法

采用规则六面体单元和并矢Green函数奇异积分等效积分技术,已有的大地电磁积分正演方法具有不能有效模拟地下复杂地质体和计算精度偏低的缺点.本文提出了一种新的三维大地电磁积分方程正演技术,即采用四面体单元、解析的并矢Green函数奇异积分表达式,达到既能模拟地下复杂异常体,又能有效提高已有积分方程法计算精度的目的.首先,采用四面体网格技术离散地下复杂异常体,获得四面体单元上的大地电磁积分方程.然后,利用针对四面体单元开发的新的奇异值积分的解析表达式,准确计算线性方程中的并矢Green函数的奇异积分,从而获得精确的线性方程.借助于PARDISO高性能并行直接求解器,实现了三维大地电磁问题的高精度求解.最后,基于国际标准3D-1模型和六棱柱模型,通过与其他方法结果的对比分析,验证了本文方法的正确性、处理高电导率对比度的能力(1000:1)和处理复杂模型的能力.     

POTENTIAL FIELD OF A STATIONARY ELECTRIC CURRENT USING FREDHOLM'S INTEGRAL EQUATIONS OF THE SECOND KIND*

An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite.     

基于谱烈度改进算法的研究与应用

本文基于Housner谱烈度公式中存在的不足和问题提出了一种改进算法,首先将相对速度反应谱积分公式上下限范围扩大为0.1～10.1 s,通过统计分析集集地震、汶川地震和芦山地震反应谱积分值和宏观烈度之间的关系,给出8组分档线性公式来计算改进谱烈度值I_s,并采用普通克里金插值法绘制谱烈度分布图。本文提出的改进算法在地震动加速度时程积分与谱烈度之间建立起一组长周期、连续性和精细化的函数关系,经鲁甸地震对比验证,该改进算法得出的谱烈度分布图与宏观烈度图整体变化趋势一致,谱烈度值与宏观烈度吻合率(偏差±1度以内)约为90%,具有一定的科学性和准确性,可为灾后应急救援以及决策部署等提供科学的数据支撑。     

Precise integration methods based on the Chebyshev polynomial of the first kind

This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method （IFM） and the homogenized initial system method （HISM）. In both methods, nonlinear variable loadings within time intervals are simulated using Chebyshev polynomials of the first kind before a direct integration is performed. Developed on the basis of the integral formula, the recurrence relationship of the integral computation suggested in this paper is combined with the Crout decomposed method to solve linear algebraic equations. In this way, the IFM based on Chebyshev polynomial of the first kind is constructed. Transforming the non-homogenous initial system to the homogeneous dynamic system, and developing a special scheme without dimensional expansion, the HISM based on Chebyshev polynomial of the first kind is able to avoid the matrix inversion operation. The accuracy of the time integration schemes is examined and compared with other commonly used schemes, and it is shown that a greater accuracy as well as less time consuming can be achieved. Two numerical examples are presented to demonstrate the applicability of these new methods.