共查询到20条相似文献，搜索用时 15 毫秒
1.
Analytical solutions of vertical electrical soundings (VES) have mostly been applied to groundwater exploration and monitoring groundwater quality on terrains of fairly simple geology and geomorphology on which the electrode arrays are symmetrical (e.g. Schlumberger or Wenner configurations). The sounding interpretation assumes flat topography and horizontally stratified layers. Any deviations from these simple situations may be impossible to interpret analytically. The recently developed GEA-58 geoelectrical instrument can make continuous soundings along a profile with any colinear electrode configuration. This paper describes the use of finite-difference and finite-element methods to model complex earth resistivity distributions in 2D, in order to calculate apparent resistivity responses to any colinear current electrode distribution in terrains in which the earth resistivities do not vary along the strike. The numerical model results for simple situations are compared with the analytical solutions. In addition, a pseudo-depth section of apparent resistivities measured in the field with the GEA-58 is compared with the numerical solution of a real complex resistivity distribution along a cross-section. The model results show excellent agreement with the corresponding analytical and experimental data. 相似文献
2.
3.
PRABHAKAR S. NAIDU 《Geophysical Prospecting》1970,18(Z1):774-785
A point source or a dipole source of electric current, placed on a randomly layered semi-infinite medium, produces an inhomogeneous random potential field on the surface. The variance of the random potential, normalized with reference to the normal field (that is, a field due to a point source or a dipole source on a homogeneous medium), falls off inversely as the distance from the source. The conductivity of the layers is assumed to vary randomly about a mean value (unity) such that the variations can be represented by a zero mean homogeneous random function. Further the variations are assumed to be small compared to the mean so that the first order perturbation is adequate. The analysis shows that the dipole field is more sensitive than the single pole field to the conductivity variations. 相似文献
4.
井地电阻率成像法利用井套管作电流源向井下供入大功率直流电流,在地表测量由地下介质的电性变化形成的电位分布,通过反演可得到地下介质的电阻率分布.针对大斜度井和水平井开展井地三维电阻率数值模拟和反演研究,对油田注水及压裂效果监测具有重要意义.基于井地电阻率成像法原理,采用有限差分法和不完全切勒斯基共轭梯度法进行了三维正演模拟,研究了大斜度井和水平井的井地电位响应特征.提出了采用层状约束阻尼最小二乘法由浅到深地进行大斜度井和水平井的多层联合反演,并对实际水平井井地电位各个层段数据进行了三维反演.模拟结果表明,倾斜线源和水平线源会对地面电位响应产生明显影响,在反演中需要考虑线源形态.实际水平井井地电位反演成像表明,考虑倾斜线源或者水平线源的联合反演得到了准确的水平井三维注水层成像图,得出注水层的真电阻率分布,能够判断注水运移方向. 相似文献
5.
针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性. 相似文献
6.
7.
8.
地震波在地下含流体孔隙介质中传播时,会引起中观尺度的"局域流",进而产生地震波震电效应.基于Biot(1941)固结理论的准静态方程,在频率域中采用空间有限差分方法,正演模拟虚岩石物理岩样的地震波衰减和震电效应.与时间域虚岩石物理方法相比,该方法既可以直接求取任一频率下的地震波衰减和电势,便于应用于实际岩样的预测分析,也避免了讨论岩样外表面施加的力源函数表达式及时间剖分稳定性条件等问题.首先利用周期性层状介质模型验证了本文所描述方法的有效性,并进一步求取分析了周期性层状介质两种不同特征单元的渗流电流密度及电势,数值模拟结果表明由中观尺度"局域流"引起的震电效应电势振幅数量级在实验室测量范围之内,随后,分析研究了四种不同高渗介质占比值的地震衰减及震电效应特征.最后,将本文提出的震电效应数值计算方法推广至二维,并求取了二维斑块饱和模型的地震波衰减、速度频散、电势的振幅和相位角数值结果. 相似文献
9.
While finite-difference methods have been used extensively for many years to model wave propagation in elastic media, some of the more subtle effects observable in such models are very inadequately documented in the geophysical literature, especially in regard to their practical numerical consequences. In addition to the intended travelling waves, and the undesirable exponential instability revealed by the von Neumann test, typical second-order-time finite-difference equations also support drifting linear solutions, as can be verified, both theoretically and by numerical experiment. The necessity of these solutions, and their relationship to the incompleteness of the set of travelling-wave eigenfunctions of the finite-difference operator, can be exposed by a matrix-based analysis, and exact expressions for them can be obtained by using standard algebraic techniques. A further peculiarity of the finite-difference formulation is numerical anisotropy, which emerges in a grid of more than one spatial dimension, even when the modelled medium is intended to be isotropic. This anisotropy can be explained and quantified in terms of the exact eigenfunction solutions to the finite-difference equation, which, it is found, can be obtained in a simple, closed form, for a typical modern 3D staggered scheme. 相似文献
10.
R.D. Hibbs F.W. Jones V. Ramaswamy H.W. Dosso 《Physics of the Earth and Planetary Interiors》1978,16(4):327-340
This work compares experimental analogue model measurements and finite-difference numerical calculations of the electric and magnetic fields for a highly conducting slab embedded in a poorly conducting host earth for three different source field configurations. Measurements and calculations were carried out for a uniform source, a sheet current source with a y exp(?ay) current intensity distribution, and a horizontal magnetic dipole source. The results indicate reasonable agreement with some exceptions between the analogue and numerical methods. The source field is found to have an important effect on the field anomalies at the interface of the highly conducting slab and the poorly conducting host medium. 相似文献
11.
本文从一阶方程组形式的波动方程出发,发展了一种计算二维不均匀介质中点源P-SV波响应的近似方法。该方法通过引入线分布的应力作为震源,利用二维有限差分方法计算出线源响应,然后再经过波形校正和几何扩散校正得出相应的近似点源响应。通过把波形和振幅与精确解比较表明,该方法具有较好的精度。由于有限差分方法对于介质中速度和密度的分布没有特殊要求,另一方面,本文所给出的震源可以适用于位错点源、爆炸源或集中力源,因此上述方法十分适合于研究横向不均匀介质中的近场强地运动、爆炸振动或地震勘探等问题。 相似文献
12.
D. PATELLA 《Geophysical Prospecting》1977,25(4):699-729
The electrical potential generated by a point source of current on the ground surface is studied for a multi-layered earth formed by layers alternatively characterized by a constant conductivity value and by conductivity varying linearly with depth. The problem is accounted for by solving a Laplace's differential equation for the uniform layers and a Poisson's differential equation for the transitional layers. Then, by a simple algorithm and by the introduction of a suitable kernel function, the general expression of the apparent resistivity for a Schlumberger array placed on the surface is obtained. Moreover some details are given for the solution of particular cases as 1) the presence of a infinitely resistive basement, 2) the absence of any one or more uniform layers, and 3) the absence of any one or more transitional layers. The new theory proves to be rather general, as it includes that for uniform layers with sharp boundaries as a particular case. Some mathematical properties of the kernel function are studied in view of the application of a direct system of quantitative interpretation. Two steps are considered for the solution of the direct problem: (i) The determination of the kernel function from the field measurements of the apparent resistivity. Owing to the identical mathematical formalism of the old with this new resistivity theory, the procedures there developed for the execution of the first step are here as well applicable without any change. Thus, some graphical and numerical procedures, already published, are recalled. (ii) The determination of the layer distribution from the kernel function. A recurrent procedure is proposed and studied in detail. This recurrent procedure follows the principle of the reduction to a lower boundary plane, as originally suggested by Koefoed for the old geoelectrical theory. Here the method differs mainly for the presence of reduction coefficients, which must be calculated each time when passing to a reduced earth section. 相似文献
13.
A numerical experiment carried out to investigate the structural model of the Domenico lignite site is discussed. The model is a 2D structure containing several lignite layers at different depths, and a low-velocity layer at the top of the model. The experiment consists in simulating a measured CDP section by two independent techniques, based on completely different concepts: the finite-difference method and the ray method. Due to the incompleteness of the ray synthetic wave field, as well as to numerical problems of the finite differences at higher frequencies, the agreement between the synthetic seismogram sections for the individual shot points is poor. However, the CDP stacked sections modelled by the ray and finite-difference methods agree rather well. This is because the main differences between the wave fields computed by the two methods are due to the presence of the low-velocity layer (ground roll, head waves, etc.), and just these parts of the wave field can be suppressed by routine data processing such as f–k filtration. Synthetic ray and finite-difference CDP stacks agree relatively well with the observed data. They confirm three lignite seams and a fault in the shallower one. The synthetic data also indicate that many apparent horizons of the measured section may be due to the multiple reflections within the subsurface low-velocity layer. 相似文献
14.
To be able to make use of all existing probe interpretation techniques, it is desirable that field resistivity observations be conducted in such a way that it is possible to construct an equivalent curve for the surface variation of potential about a point source. Further, the usual parallel layer interpretation will only be justified in as much as the observed potential curve is compatible with such a subsurface geometry. Thus, whilst a potential curve may be constructed from suitable finite potential differences obtained using a Wenner configuration of electrodes, it can be shown, using the tri-potential technique of measurements, that these potential differences may arise partly from lateral resistivity variations. In this paper, a ladder network technique is employed to display these lateral effects and an adjustment method proposed to reduce them. The adjusted potential differences obtained form a consecutive series suitable for summation to give the potential function. These values are subject to further minimum adjustments required for them to comply with slope and curvature conditions for a layered medium. After forming the potential curve by summation, a final numerical smoothing process is carried out. From this smoothed potential curve, corresponding Wenner and Schlumberger curves can readily be derived for interpretation. The method proposed thus attempts to extract from the observed data the maximum part which can be reconciled with a purely depth variation of resistivity. A method is also proposed for compounding the ‘errors of closure’ of the network to provide a Lateral Inhomogeneity Index which gives a measure of the departure of the observed data from the basic interpretation requirements. The method involves no great labour, but can readily be programmed for a computer if desired. Examples are given of the application of the method to field observations. In the interest of objectivity, the final smoothing has been confined to a single stage but in certain very extreme examples a further stage may prove desirable. 相似文献
15.
Borehole-to-surface electrical imaging (BSEI) uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient (ICCG) method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low-resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source. 相似文献
16.
E.H. ELORANTA 《Geophysical Prospecting》1986,34(6):856-872
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. 相似文献
17.
针对地下工程领域隧道超前预报地震波波场传播与成像中存在的问题,通过数值模拟,构建二维含低速异常的隧道介质模型,研究隧道弹性波场传播规律和异常体边界成像准确性.首先,利用一阶速度-应力波动方程和高阶交错网格有限差分计算方法,导出隧道超前预报数值模拟的稳定性条件和边界条件,对上述隧道模型进行数值模拟,识别波场特征;其次,利用叠前逆时偏移成像方法,对压制噪音干扰后的波场在互相关成像条件下,对隧道模型中的异常体边界进行逆时偏移成像.研究结果表明：采用高阶交错网格有限差分正演获得异常体边界清晰的反射波和角点产生的散射波;逆时偏移算法获得隧道内异常体准确成像结果,从而大大提高隧道超前预报的分辨率与准确性;靠近掌子面单一震源、多道接收观测系统对异常体成像效果最佳,为隧道内高效数据采集提供理论依据. 相似文献
18.
通过对近似解析离散化(NAD)方法的分析,给出了一种求解声波和弹性波方程的带权重的近似解析离散化(WNAD)方法,并用WNAD方法、Lax-Wendroff 修正格式(LWC)和二阶中心差分方法计算了二维波动方程初值问题的应力场数值误差.结果表明WNAD方法具有更高的数值精度.用WNAD方法、LWC和四阶交错网格法对二维非均匀介质中弹性波传播的应力场进行了数值模拟.应力场快照和地表地震记录表明,即使是在粗网格条件下WNAD方法的模拟结果仍无可见的数值频散和源噪声.另一方面,由于WNAD方法同时计算了地震位移和梯度场,使得应力的计算更为便捷和精确,而且WNAD方法中波位移梯度局部连接关系的使用使得应力在间断处能够自动近似地满足应力连续性. 相似文献
19.
本文研究了三层介质中膨胀球模拟震源的应力、应变场以及由这个应力、应变场产生的电阻率变化量和地表视电阻率变化的关系。结果表明,对于一定强度的“震源”,在地电装置探测范围内的应力场随具有不同弹性参数的层位有显著的变化;而应变场却看不出有明显的突变;由地表视电阻率变化量的计算结果揭示了地电异常受台站介质电性结构与力学结构的综合影响。因此作为预报地震的地电阻率法一方面要寻找具有高电阻率——应变灵敏度岩层的台址,另一方面台址的视电阻率变化对各层真电阻率变化要具有好的响应。 相似文献
20.
高精度频率域弹性波方程有限差分方法及波场模拟 总被引：18，自引：4，他引：14
有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,本文采用25点优化差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入最佳匹配层,构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作（如成像、反演）提供了研究基础. 相似文献