本文提出了一种用于模拟非轴对称套管井中扇区水泥胶结测井100kHz3维瞬态声场的并行算法.采用时域有限差分方法,以及非裂化的完全匹配层吸收边界条件,基于信息传递接口(MPI)并行环境,在计算机集群上实现并行计算.对内边界上的差分方程进行参数平均处理,提高了计算稳定性.通过在轴对称情况下将有限差分模拟结果与实轴积分法计算结果对比,验证了方法的正确性.使用并行有限差分算法,模拟了水泥环有不同大小和在不同方位缺失时非轴对称套管井中8扇区的扇区水泥胶结测井仪的测井波形,并分析了套管波幅度随水泥缺失一换能器夹角的变化.根据模拟和分析结果,得到了非轴对称套管井中新的声场传播规律,并建议了一种用于反演水泥缺失大小和方位的方法.
Using finite difference time domain (FDTD) method with nonsplitting perfectly matched layer absorbing boundary condition, a parallel computation scheme based on Message Passing Interface on computer-based cluster, is proposed to simulate high frequency (100 kHz) 3D transient acoustic field of Segmented Bond Tool logging in nonaxisymmetric cased hole. Average parameters of difference equations on the internal interfaces are proposed to improve the computational stability. This method is verified by the good agreement of the simulated waveforms with those obtained from the real-axis integration method for axisymmetric cased hole. The waveforms of 8-sector SBT logging in nonaxisymmetric cased hole with cement groove of various size and azimuth are simulated using the parallel FD algorithm, and the amplitude characteristics of casing wave versus included angle between cement groove and transducer are studied. New wave phenomena in nonaxisymmetric cased hole are presented. An idea to determine the size and azimuth of the cement groove is proposed.