共查询到20条相似文献，搜索用时 453 毫秒
1.
F. W. Jones 《Pure and Applied Geophysics》1974,112(5):793-800
Summary A numerical method is used to calculate the geomagnetic fields associated with a three-dimensional conductivity anomaly. Fields associated with a two-dimensional model are also studied numerically for a range of frequencies and apparent resistivity curves for the two models are compared with that calculated for a layered earth. The apparent resistivity curves for both the three-dimensional model and the two-dimensional model differ considerably from the layered case, and it is evident that if a layered model is used for interpretation the results may be very misleading. 相似文献
2.
There exist different response characteristics in the resistivity measurements of dual laterolog (DLL) and logging while drilling (LWD) electromagnetic wave propagation logging in highly deviated and horizontal wells due to the difference in their measuring principles. In this study, we first use the integral equation method simulated the response characteristics of LWD resistivity and use the three dimensional finite element method (3D-FEM) simulated the response characteristics of DLL resistivity in horizontal wells, and then analyzed the response differences between the DLL and LWD resistivity. The comparative analysis indicated that the response differences may be caused by different factors such as differences in the angle of instrument inclination, anisotropy, formation interface, and mud intrusion. In the interface, the curves of the LWD resistivity become sharp with increases in the deviation while those of the DLL resistivity gradually become smooth. Both curves are affected by the anisotropy although the effect on DLL resistivity is lower than the LWD resistivity. These differences aid in providing a reasonable explanation in the horizontal well. However, this can also simultaneously lead to false results. At the end of the study, we explain the effects of the differences in the interpretation of the horizontal well based on the results and actual data analysis. 相似文献
3.
The difficulty to use master curves as well as classical techniques for the determination of layer distribution (ei, ρi) from a resistivity sounding arises when the presumed number of layers exceeds five or six. The principle of the method proposed here is based on the identification of the resistivity transform. This principle was recently underlined by many authors. The resistivity transform can be easily derived from the experimental data by the application of Ghosh's linear filter, and another method for deriving the filter coefficientes is suggested. For a given theoretical resistivity transform corresponding to a given distribution of layers (thicknesses and resistivities) various criteria that measure the difference between this theoretical resistivity transform and an experimental one derived by the application of Ghosh's filter are given. A discussion of these criteria from a physical as well as a mathematical point of view follows. The proposed method is then exposed; it is based on a gradient method. The type of gradient method used is defined and justified physically as well as with numerical examples of identified master curves. The practical use for the method and experimental confrontation of identified field curves with drill holes are given. The cost as well as memory occupation and time of execution of the program on CDC 7600 computer is estimated. 相似文献
4.
New Methods for Modeling Laterolog Resistivity Corrections 总被引：1，自引：1，他引：0
Jadwiga Anna Jarzyna Adam Cichy Dezső Drahos Attila Galsa Maria Joanna Bała Andrzej Ossowski 《Acta Geophysica》2016,64(2):417-442
The paper presents methods for laterolog response modeling. In Coulomb’s charges method, Laplace’s equation is solved for the electric field distribution in rock medium with internal boundaries between different resistivity layers. There, the boundary problem is reduced to Fred-holm integral equation of the second kind. The second method uses a finite element array to model apparent resistivity from laterolog. The task is treated as DC problem and the Laplace equation is solved numerically. The presented methods were applied to borehole data covering a typical stratigraphie section of the Fore-Sudetic Monocline in southwestern Poland. Apparent resistivity was calculated using the Coulomb’s charges method and alternatively modeled using a finite element method which gave similar results. Then, a series of linear corrections for borehole, shoulder bed, and filtration effects for apparent resistivity obtained by the Coulomb’s charges method demonstrated the feasibility of calculating true resistivity of virgin and invaded zones. The proposed methods provide a flexible solution in modeling which can be adapted to other logs. 相似文献
5.
本文实现了2.5维电导率正交各向异性海洋可控源电磁等参有限元数值模拟.利用傅里叶变换导出了电导率正交各向异性2.5维海洋可控源电磁法波数域电磁场耦合方程,采用伽里金加权余量法推导了相应的有限元方程;采用任意四边形单元对研究区域进行剖分,在单元中进行双二次插值,将有限元方程化为线性代数方程组;最后,求解线性方程组并进行反傅里叶变换获得空间域电磁场值.这个方法可以模拟海底起伏地形条件下地下任意形状电导率正交各向异性的复杂模型.与一维模型的数值模拟结果对比表明,电磁场数值解与解析解吻合.二维模型的计算结果与二维自适应非结构有限元模拟结果也吻合.水平海底二维地电模型考察了不同各向异性系数对海洋可控源电磁响应的影响特征.海底起伏地形地电模型的数值结果表明,电导率各向异性对海洋可控源电磁响应影响明显,有可能淹没海底地形和高阻油气藏引起的异常. 相似文献
6.
Baker Hughes INTEQ的MPR是应用较广的传播电阻率测井仪器,它有4条相位差和4条幅度比测井曲线,本文以它在垂直井中的测井响应为例开展电磁传播电阻率资料的二维反演方法研究.由于其测井响应对地层纵向边界的探测灵敏度比对其他地层参数的大得多,在其资料反演中必须重视地层纵向边界位置的提取和反演精度.本文首先应用特征识别方法由分辨率最好的测井曲线提取地层纵向边界,接着对浅、深探测的两条相位差曲线分别做纵向一维反演,利用此结果来建立二维反演的初始地层模型,采用改进的阻尼型正则化高斯 牛顿优化算法,实现了MPR资料的二维全参数反演.对含23个纵向分层的两种典型地层模型的数值模拟测井资料的二维反演得到了令人满意的结果,它可显著提高MPR的资料解释水平.发现只用它的5条测井曲线就可得到与8条测井曲线几乎相同的反演效果. 相似文献
7.
This study presents a novel mathematical model for analysis of non-axisymmetrical solute transport in a radially convergent flow field with scale-dependent dispersion. A two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates is derived based on assuming that the longitudinal and transverse dispersivities increase linearly with the distance of the solute transported from its injected source. The Laplace transform finite difference technique is applied to solve the two-dimensional, scale-dependent advection–dispersion equation with variable-dependent coefficients. Concentration contours for different times, breakthrough curves of average concentration over concentric circles with a fixed radial distance, and breakthrough curves of concentration at a fixed observation point obtained using the scale-dependent dispersivity model are compared with those from the constant dispersivity model. The salient features of scale-dependent dispersion are illustrated during the non-axisymmetrical transport from the injection well into extraction well in a convergent flow field. Numerical tests show that the scale-dependent dispersivity model predicts smaller spreading than the constant-dispersivity model near the source. The results also show that the constant dispersivity model can produce breakthrough curves of averaged concentration over concentric circles with the same shape as those from the proposed scale-dependent dispersivity model at observation point near the extraction well. Far from the extracting well, the two models predict concentration contours with significantly different shapes. The breakthrough curves at observation point near the injection well from constant dispersivity model always produce lesser overall transverse dispersion than those from scale-dependent dispersivity model. Erroneous dimensionless transverse/longitudinal dispersivity ratio may result from parametric techniques which assume a constant dispersivity if the dispersion process is characterized by a distance-dependent dispersivity relationship. A curve-fitting method with an example is proposed to evaluate longitudinal and transverse scale-proportional factors of a field with scale-dependent dispersion. 相似文献
8.
伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映. 相似文献
9.
《Advances in water resources》2005,28(8):793-806
In this paper, the numerical errors associated with the finite difference solutions of two-dimensional advection–dispersion equation with linear sorption are obtained from a Taylor analysis and are removed from numerical solution. The error expressions are based on a general form of the corresponding difference equation. The variation of these numerical truncation errors is presented as a function of Peclet and Courant numbers in X and Y direction, a Sink/Source dimensionless number and new form of Peclet and Courant numbers in X–Y plane. It is shown that the Crank–Nicolson method is the most accurate scheme based on the truncation error analysis. The effects of these truncation errors on the numerical solution of a two-dimensional advection–dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution for predicting contaminant plume distribution in uniform flow field. Considering computational efficiency, an alternating direction implicit method is used for the numerical solution of governing equation. The results show that removing these errors improves numerical result and reduces differences between numerical and analytical solution. 相似文献
10.
过套管电阻率测井通过测量套管壁电势实现测量地层视电阻率,不同的地层模型对过套管电阻率测井会有不同的测井响应. 在传输线方程系数中本文考虑了径向含多个界面地层的影响,对轴向（井轴方向）呈层状、径向阶跃变化地层模型,给出了电势分布递推公式,计算了过套管电阻率测井响应,绘出了测井响应曲线. 算例表明,该方法有较强的轴向地层界面分辨能力;对径向均匀地层计算结果非常接近地层模型电导率,能较好地反映实际地层特征;对高电导地层,低电导水泥环有明显的测井响应;当径向存在多个地层界面时计算结果则为径向各层电导率的综合反映,不是某一径向层的电导率;算例显示该方法有较快（算例运行时间在1s内）的计算速度. 相似文献
11.
最近开发了一种针对二维大地电磁野外数据进行处理解释的新反演方法。该方法以加入阻抗相位信息的一维大地电磁连续介质曲线对比法为基础 ,把一维反演得到的电阻率和相位的数据集作为二维反演的初始模型 ,使用二维有限单元法做正演模拟。在程序的后继迭代中 ,深度方向上用一维反演修改模型的电阻率和深度值 ,沿测线方向由二维有限元作修改 ,反演结束可得到一个接近真实电性分布的电阻率数据集 ,并绘制成电阻率断面图。对模型的反演实验结果显示 ,该反演方法能够较真实地反映地下电性分布 ,而且避免了偏导数矩阵的计算 ,其原理简单 ,计算速度快 ,表明该反演方案是可行的 相似文献
12.
Numerical computations using the integral equation method are presented for resistivity and IP responses due to arbitrarily shaped 3-dimensional bodies in a layered earth. The unknown surface charge density distribution is expressed as the solution of Fredholm's integral equation of the second kind. Use of moment method (with pulse basis function and point-collocation) yields the matrix equations for the unknowns. The contributions to Green's function are solved (a) analytically for the primary and (b) by convolution for the secondary contributions resulting in a fast algorithm. The further step of computing potential, apparent resistivity, chargeability etc., for any electrode system, is straightforward. Our results show a good agreement with those from finite difference methods and physical tank experiments. The CPU time is only 138 s on a super-minicomputer for an apparent resistivity pseudo-section, even with 96 elementary cells as used for discretization. A large number of models for different geological situations were studied; some are presented here. 相似文献
13.
A. M. Popov 《Pure and Applied Geophysics》1990,134(4):575-587
Experimental data show that in East Siberia resistivity curves, irrespective of their trends, are affected by galvanic (local) distortions. The preliminary step of the magnetotelluric data processing is to obtain a steady shape of resistivity curves reflecting a true deep section. For this purpose statistical averaging and different criteria of impedance rejecting were used. The available MTS curves were normalized by
level to the global magnetovariation curves. Two-dimensional modelling was performed from several sublatitudinal profiles crossing the Baikal rift zone. Three-dimensional models based on two-dimensional modelling and on induction vector distribution have been computed via programs of M. N. Yudin. Following other researchers, two conductive layers are distinguished: i) the mid- and low crustal and ii) the mantle one, with the layer surface uplifted from 100–110 km depth in the southern Baikal rift zone to 60–70 km northeastwards along the eastern Baikal coast. The top of this layer seems to correspond to the asthenospheric roof. The asthenosphere deepening in southern BRZ is likely to be related to a decrease in the asthenospheric bulge width and an increase in the rate of lithospheric thickening with mantle degasing. The origin and evolution of the Baikal rift is considered, proceeding from the model of passive rifting with regard to a long-existing lithospheric inhomogeneity between the Siberian platform and the Sayan-Baikal folded area. 相似文献
14.
An earth model with a transition layer (anisotropic inhomogeneous) is considered. The inhomogeneity in σv (vertical conductivity) of the transition layer is represented by a power law variation. Expressions for potential distribution in the upper layer, transition layer and bottom layer are obtained by solving appropriate differential equation for each layer. By utilizing the boundary conditions, expressions of apparent resistivity for Wenner and Schlumberger configurations are derived. Numerical analysis is performed for linear and quadratic variation of σv. The results are presented in the form of theoretical apparent resistivity curves for both configurations. Negative apparent resistivities are the interesting feature of this analysis. 相似文献
15.
Summary The finite element method, with triangular elements, is used to study the effect of a two-dimensional sloping contact on the surface electromagnetic fields. It is found in the case ofH-polarization and small slopes that the electric field and the apparent resistivity near the contact, on the conductive side, are higher than their asymptotic values. In the case ofE-polarization the apparent resistivity and phase values on the conductive side fall off less rapidly to their asymptotic values with decreasing slope resulting in higher apparent resistivity and phase values on the conductive side, than those expected for a vertical contact. The peak in the amplitude and phase of the normalized vertical magnetic field shifts from the resistive side for a vertical contact to the conductive side for a sloping contact. Far from the sloping contact, on the conductive side, higher values are observed for the normalized vertical magnetic field than in the case of a vertical contact. 相似文献
16.
2016年5月浚县地震台地电阻率EW向测值出现明显变化,为正确评价该变化,依据浚县地震台地质资料、岩层电性资料和EW向电测深曲线,将该台地下介质电性分布简化为二维电性结构,建立三维有限元模型。通过数值模拟,计算测区内地表铁轨干扰对地电阻率观测的影响形态和幅度。结果表明,由数值模拟所得干扰幅度、形态与实际测值的异常变化较吻合,因此从数值分析的角度确定了铁轨对浚县地震台地电阻率观测的干扰。 相似文献
17.
针对环井周的非均质地层模型,利用三维有限元方法模拟了方位梯度电极系的测井响应．结果表明,在非均质程度比较小的情况下,各个方位上的响应之间差别不大,与真实的电阻率分布差别也不大;反之,则差别很大．其特点是所测得的响应远远低于较高的真电阻率值,而接近较低的真电阻率值．因此,有必要对测井响应进行校正．考虑到传统的校正方法的弊端,本文利用非线性反演方法求取真实的电阻率分布,数值模拟的结果表明了此方法的可行性． 相似文献
18.
The advantages of the Wenner tripotential method (Carpenter 1955) for apparent resistivity profiling are described and two new data processing techniques introduced as an aid to the interpretation of apparent resistivity sections (pseudo-sections). These techniques were developed from model data computed using a two-dimensional finite difference method. Oscillatory components present in anomalies on tripotential profiles and related to electrode spacing are shown to be effectively removed by linear filtering that also simplifies their form and aids recognition. Furthermore, it is shown that the ratio of the beta- and gamma-apparent resistivities is a good indicator of resistivity variation, and is particularly sensitive to lateral change. Model data indicates that, over a wide range of conditions, enough subsurface information can be obtained by inspection of tripotential resistivity and ratio profiles, and from space sections to make possible a useful—and sometimes semi-quantitative—interpretation. A rationale for the general interpretation of tripotential data is developed. Field data are described from an area of weathered granite basement in Nigeria. A model of the subsurface is developed using parameters derived from the processed observations. The observed and calculated apparent resistivity space sections are very similar. 相似文献
19.
20.
系统地论述了用有限单元法研究复杂地形条件下三维直流电阻率的正演计算技术.首先给出了三维构造中点源电场的边值问题以及相应的变分问题;然后利用有限单元法求解变分问题,采用四面体单元对研究区域进行剖分,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后,考虑到节约计算时间,利用对称超松弛顸条件共轭梯度迭代算法求解大型线性方程组,得到了各节点的电位值,进而计算出地表的视电阻率.通过理论模型的计算检验了算法的可行性之后,给出了几种常见纯地形异常的数值模拟结果和一个组合模型的计算结果,其研究工作为研究三维直流电阻率反演奠定了基础. 相似文献