多层地电断面视电阻率的高精度计算

CALCULATION OF APPARENT RESISTIVITY FOR PROFILE OF MULTI-LAYER RESISTIVITY WITH HIGH ACCURACY

目前地电的方法在我国已应用到地震预报的研究中,为了了解与地下力学场变化相联系的地电变化,需要了解台址下勘探体积内的电阻率分层情况,从而正确地认识各种影响因素对地电阻率的作用。

In order to meet the needs of the requirement of studying underground resistivity problems related to earthquakes, we presented a formula to calculate apparent resistivity, which is deduced directly from the infinite integration for n-layer distribution of resistivity, without using the us- ual series method. According to this formula we program a Fortran procedure to compute the curve of apparent resistivity. By using the techniques both in mathematics and programming design, this procedure computes fast, stably, and has high accuracy (relative error is lower than 0.001). This has also been proved by a great deal of calculations of real problems. The method avoids the weaknesses often appeared in series method, such as lower accuracy in earthquake research and the fixed distance (usually zero) between measuring poles.