张军, 刘百祥, 吴燕清, 袁永榜, 许洋铖. 基于最小二乘多项式分段拟合的全区时间域电磁视电阻率直接计算[J]. 煤炭学报, 2019, 44(S1): 214-219. DOI: 10.13225/j.cnki.jccs.2019.0289
引用本文: 张军, 刘百祥, 吴燕清, 袁永榜, 许洋铖. 基于最小二乘多项式分段拟合的全区时间域电磁视电阻率直接计算[J]. 煤炭学报, 2019, 44(S1): 214-219. DOI: 10.13225/j.cnki.jccs.2019.0289
ZHANG Jun, LIU Baixiang, WU Yanqing, YUAN Yongbang, XU Yangcheng. Direct calculation of electromagnetic apparent resistivity in whole time domain based on piecewise fitting of least squares polynomials[J]. Journal of China Coal Society, 2019, 44(S1): 214-219. DOI: 10.13225/j.cnki.jccs.2019.0289
Citation: ZHANG Jun, LIU Baixiang, WU Yanqing, YUAN Yongbang, XU Yangcheng. Direct calculation of electromagnetic apparent resistivity in whole time domain based on piecewise fitting of least squares polynomials[J]. Journal of China Coal Society, 2019, 44(S1): 214-219. DOI: 10.13225/j.cnki.jccs.2019.0289

基于最小二乘多项式分段拟合的全区时间域电磁视电阻率直接计算

Direct calculation of electromagnetic apparent resistivity in whole time domain based on piecewise fitting of least squares polynomials

  • 摘要: 时间域电磁视电阻率计算中,由于感应电压与视电阻率的关系是非单调函数,无法通过感应电压直接计算。现有方法主要是晚期视电阻率近似计算法、核函数逐步逼近法和核函数平移法。晚期视电阻近似计算法在早期和中期产生计算误差; 核函数逼近法需要多次逐步计算逼近目标电阻率,计算时间较长; 核函数平移法需要通过插值方式计算各时间道所对应的实测截距在理论截距曲线上对应的时间,再根据感应电压衰减曲线的伸缩特性计算出视电阻率,技术较为复杂。为了提高计算精度,缩短计算时间,根据核函数表达式在视电阻率计算中保持不变的特点,提出基于最小二乘多项式拟合的全区时间域电磁视电阻率计算方法。采用高阶最小二乘多项式拟合核函数左右2支的反函数,推导出视电阻率多项式拟合直接计算式。为了降低拟合算法的跨度,调和数量级所导致的大数吃小数问题,避免因动态范围太大影响无法直接拟合,进行对数变换来降低动态范围,达到增加拟合效果和提高计算精度。建立均匀半空间、层状大地模型,采用12阶多项式进行拟合计算并通过MATLAB语言编程实现,对比分析核函数二分逐步逼近法和对数变换最小二乘多项式拟合计算的视电阻率和2种方法各自的计算时间,最后结合现场已知采空积水区试验,通过数值模型和实测数据验算了方法的正确性。研究表明:基于最小二乘多项式分段拟合的全区时间域电磁视电阻率计算方法能够正确计算视电阻率,与模型和实际情况基本一致; 与核函数二分逐步逼近法相比,本方法利用对数变换,计算速度提高近5倍,缩短了计算时间。

     

    Abstract: Real-time downhole time-domain electromagnetic imaging technology requires a fast calculation of apparent resistivity in downhole.As the relationship between induction voltage and apparent resistivity is not a monotonic function, the apparent resistivity cannot be calculated directly from induction voltage.The existing methods mainly consist of late apparent resistivity approximation method, kernel function step-by-step approxi-mation method and kernel function translation method.The late apparent resistance approximation method pro-duces some calculation errors in the early and middle stages.The kernel function approximation method needs to calculate the approximation resistivity step by step for many times, which takes a long time for calculation.The Kernel function translation method needs interpolation to calculate the time corresponding to the measured intercept of each time channel on the theoretical intercept curve, then the apparent resistivity is calculated according to the stretching characteristic of the induced voltage attenuation curve, it is more complex.In order to shorten the calculation time, according to the feature that the kernel function expression remains unchanged in the calculation of apparent resistivity, the authors propose a direct calculation method of electromagnetic apparent resistivity in the whole time domain based on least square polynomial fitting.The inverse function of the left and right branches of the kernel function is fitted by the higher order least square polynomial, and the direct calculation formula of apparent resistivity polynomial fitting is derived.In order to reduce the span of the fitting algorithm and reconcile the decimal problem caused by the multiple orders of magnitude, avoid too large a dynamic range to fit directly, a double logarithmic transformation is used to reduce the dynamic range, achieve the purpose of increasing fitting accuracy and speed.A homogeneous half-space and layered earth model is established, and the results of the kernel function dichotomy step-by-step approximation method and the logarithmic transformation least squares pol-ynomial fitting method are compared and analyzed by using the 12-order polynomial fitting method.Finally, the correctness of the method is verified by the field test of the known water-filled goaf.The results show that the method based on the least squares polynomial piecewise fitting can calculate the apparent resistivity correctly, which is basically consistent with the model and the actual situation.In comparison with the kernel function dichotomy step-by-step approximation method, the method proposed shortens the calculation time and increases the calculation speed for nearly five times by using logarithmic transformation.

     

/

返回文章
返回