Simulated various characteristic waves in acoustic full waveform relating to cement bond on the secondary interface

Bad cement bond with water channeling or low strength cement is usually generated on the first and/or secondary interfaces which lie between the cement and casing, and the cement and formation, respectively. It is an inherently complicated problem to evaluate cement bond on the secondary interface (SICB) in a cased borehole. So we need to find a useful wave for evaluating SICB and provide its propagation features. To achieve this, we simulated the acoustic full waveform, 2-D spectrum and dispersion curves in cased boreholes based on the real axis integration method, and a propagator matrix. Simulation results indicated that “casing-cement mode” wave can be utilized to evaluate SICB, as it is most sensitive to water channeling on the secondary interface. Velocity of this mode wave is in between that of the casing wave and the formation wave, so that it is easy to extract it in time. The dispersion curves showed that the casing-cement wave is a dispersion wave, but not a pseudo-Rayleigh wave. Its velocity is obviously larger than v_s of cement. Simulation results further indicated that the amplitude of casing-cement wave is decreasing with increasing cement sheath thickness, but not water channeling thickness, and its travel time is influenced by casing and cement together. Moreover, we investigated factors influencing the casing–cement wave and found that casing diameter and peak frequency are major factors influencing the amplitude and the travel time, while other factors, such as P-wave velocities of cement and wall thickness of casing, are minor factors. But none of these factors greatly influence bond index which is computed from the amplitude. So the bond index has an advantage over the amplitude during evaluating SICB. Formation wave is a helpful signal to evaluate low strength cement without water channeling, and its amplitude is mainly related to the velocities (or impedances) of both cement and formation except for formation attenuation and disturbance wave. Experimental data in calibration pits and field data detecting channeling demonstrated the simulation results.     

套管井贴壁声源激发的非轴对称声场的有限差分模拟及结果分析

套管井贴壁声波测井仪器在固井质量评价中的应用越来越广泛.本文利用三维交错网格高阶有限差分方法首次研究了贴壁声源在套管井中激发的非轴对称声场特征,给出了测井声源在贴壁固体弹性介质中激发时的数值模拟方法,模拟的全波波形与测井仪器在实际井中记录的波形吻合较好,这为利用数值模拟方式开发贴壁声波仪器的潜在应用奠定了基础.通过对套管井不同胶结状况下的数值模拟,得到了贴壁声源在套管井中激发的声场特征和传播规律,模拟结果显示加载声源的贴壁滑板的存在大大压制了套管中弯曲型Lamb波的幅度;套管中传播的拉伸波的幅度相对较强,其泄漏到水泥中的声波传播到水泥和地层界面时发生反射,反射波的幅度和相位特征携带了水泥环第二界面的胶结信息.数值计算结果为利用贴壁声源激发接收的波形实现水泥环第二界面井周不同角度上的胶结评价奠定了理论基础.     

双层套管井中单极子声波测井的响应特征

目前,在深井中采用多层套管实现层间封隔的方法已非常普遍,对多层套管井固井质量评价的需求越来越迫切,目前还未有针对外层套管与地层之间水泥胶结评价的有效手段.本文利用三维有限差分仿真了单极声波测井在双层套管井中的响应特征,探讨了利用套管波幅度和斯通利波频散特征识别外层套管与地层之间环空是否胶结水泥的可行性.数值计算结果显示在两层套管之间胶结水泥时,套管波的幅度或衰减对外层套管和地层之间是否胶结水泥非常敏感,外层环空若充填水泥则全波中的套管波幅度降低明显;在两层套管之间充填泥浆时,全波波形中首先到达的套管波幅度较大,其对外层套管的胶结状况不再敏感,通过对接收的阵列波形做频率-慢度相关分析发现,在外层套管和地层之间胶结水泥后,井孔内接收到的低于10 kHz的外斯通利波慢度明显增大,仿真波形的处理结果与解析计算的斯通利波的频散曲线吻合较好,进一步验证了可尝试利用内层环空激发的外斯通利波判断外层套管与地层之间是否胶结水泥.     

An acoustic-wave theory for casing bond evaluation using viscoelastic slip boundary modeling

Slip boundary condition is commonly utilized to model elastic wave propagation through layered earth media. The same approach is used here to characterize acoustic wave propagation along a cased borehole with various cement bond conditions. By modeling the cement layer between casing and formation as a viscoelastic slip interface with complex coupling rigidity parameters, one can not only reduce the complexity in the classical elastic wave modeling of the problem, but also efficiently model various complicated wave phenomena that are difficult for the existing modeling. More specifically, the new theory can well describe the effect of the cement bond condition change and the location of the change (i.e., whether it is in the first interface between casing and cement, or the second interface between cement and formation) on the acoustic waves, demonstrating the good modeling capability and predicting power. Application of the theory to field data shows that the theory can correctly model the acoustic wave characteristics and interpret the cement bond condition, thus providing a useful fundament theory for casing bond evaluation using acoustic logging.     

偶极子声波测井在双层套管井中激发的声场特征

多层套管井中的套管-水泥以及水泥-地层在径向上的分布较单层套管井界面多且径向分布范围更深.目前还未有针对多层套管井固井质量评价的有效方法,偶极子阵列声波测井具有工作频率低径向探测深的优势,本文通过实轴积分和三维有限差分计算了偶极子声源在不同胶结状况的双层套管井中激发的声场,旨在探讨偶极子阵列声波在双层套管井中评价外层套...     

套后介质声学性质对套管中三种模式波影响的对比分析

由于套管的厚度小于目前声波测井仪器工作频率下的波长,声波仪器工作时可在套管中激发Lamb波或SH模式波.目前评价水泥胶结质量的声波测井仪器主要是利用套管中的对称Lamb波(拉伸波,S0模式)或反对称Lamb波(泄漏弯曲波,A0模式)的相对幅度或衰减评价套管与水泥之间的胶结程度.Lamb波可看作是纵波和横波在套管边界的耦合形成的,在套管中还存在另一种非耦合的SH横波,偏振方向与套管轴向平行沿着周向传播.本文基于定向声源辐射技术,实现了对称Lamb波、反对称Lamb波的分波模拟技术,数值研究了套管后耦合不同声阻抗介质时,对称Lamb波、反对称Lamb波以及SH横波的测井响应,通过频散衰减曲线以及模式波的波形响应幅度的直观显示表明,由于套管中的SH波仅向固体水泥中辐射剪切波,这一传播特征与套管中的Lamb波相比可以极大提高套管-水泥界面胶结状况的测量灵敏度.Lamb波的固有缺陷是套管与两侧介质的较大声阻抗反差使得在井中测量套管后面水泥界面的灵敏度受到相当大的限制.     

转换波方位各向异性裂缝检测技术研究及应用

HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果.     

套管井应力集中诱导井周波速各向异性的研究

针对在套管井中应用正交偶极声测井检测地层各向异性和确定异常构造应力的需要, 本文通过求解水平构造应力作用下套管井的应力场分布,应用声弹性理论说明应力集中对套管井井周弹性波速空间分布的影响,分析了波速空间分布对加套管后的正交偶极各向异性声测量的影响.井孔、套管、水泥环和地层系统的复杂性改变了原裸眼井情况地层中的应力场分布,使应力集中的范围明显缩小到井眼附近.对应不同方位径向偏振横波的交叉点向井孔附近移动.横波波速空间分布的方位极值规律发生变化,但纵波波速空间分布的方位极值规律不变,仍能提供最大远场构造应力施加的方向.     

随钻声波测井固井质量评价理论与数值模拟研究

大斜度井、水平井的固井质量评价需求,推动了基于随钻测井仪器的固井质量评价技术的发展.本文基于弹性波传播理论,建立柱状分层介质的套管井随钻测井模型,模拟了不同水泥胶结情况的随钻声场响应,特别是其波形特征.结果表明,套管井中随钻与电缆波形特征最大的区别在于随钻波列中存在隔声衰减后的钻铤波.特别是当水泥胶结较好时,钻铤波与套管波幅度大小相近,甚至钻铤波幅度更大,而时间域上两者波形重叠,该特点使得常规电缆声波(声幅/变密度)固井质量评价方法不适用于随钻测井环境.本文考察了不同水泥胶结程度下随钻波列中套管波幅度和衰减与胶结指数的关系,形成一种同时利用幅度和衰减的随钻声波胶结指数计算方法,为随钻声波固井质量评价提供理论指导.     

非轴对称套管井中声场的并行有限差分模拟

本文提出了一种用于模拟非轴对称套管井中扇区水泥胶结测井100 kHz 3维瞬态声场的并行算法.采用时域有限差分方法,以及非裂化的完全匹配层吸收边界条件,基于信息传递接口(MPI)并行环境,在计算机集群上实现并行计算.对内边界上的差分方程进行参数平均处理,提高了计算稳定性.通过在轴对称情况下将有限差分模拟结果与实轴积分法计算结果对比,验证了方法的正确性.使用并行有限差分算法,模拟了水泥环有不同大小和在不同方位缺失时非轴对称套管井中8扇区的扇区水泥胶结测井仪的测井波形,并分析了套管波幅度随水泥缺失-换能器夹角的变化.根据模拟和分析结果,得到了非轴对称套管井中新的声场传播规律,并建议了一种用于反演水泥缺失大小和方位的方法.     

套管井声波全波列中地层波的幅度特性

套管井声全渡中的地层波可以检测固井第二界面的水泥胶结质量状况.地层波不仅受地层衰减、套管波干扰等因素影响,更主要受第二界面两侧的声阻抗差值的影响.为此,基于固体介质的声势函数和界面的边界条件,推导了固-固界面纵横波的反射系数和透射系数,以及第一临界角入射情况下反射系数和透射系数的近似公式,并且分析了透射系数与声阻抗之间的关系.数值计算结果发现,地层波幅度随第二界面两侧声阻抗差增大而减小,且这种关系是非线性的.基于临界入射的透射系数,对经过衰减校正和干扰校正的地层波进行声阻抗校正,消除地层声阻抗对地层波幅度影响.使用地层波幅度评价第二界面的水泥胶结状况,并在实际资料处理中见到了好的效果.     

基于黏弹滑移界面理论揭示套管中模式波的传播特征

目前固井质量评价主要是利用套管中模式波的幅度或衰减变化判断套管与水泥间的胶结状况.在超声测量频段（70~500 kHz）,现有声波测井仪器工作时在套管中激发的模式波类似薄板中的对称兰姆波（拉伸波）或反对称兰姆波（弯曲型兰姆波）.本文将黏弹滑移界面理论应用到套管和水泥的耦合界面,探讨了套管中的准兰姆波（包括拉伸波和弯曲型兰姆波）与套管后物质的耦合方式,拉伸波与套后物质主要以剪切刚度耦合,即便是0.01 mm的微环也会使得拉伸波的衰减明显降低,接近套后是流体时的响应特征;而弯曲型兰姆波与套管后物质主要以法向刚度耦合,这使得在套后胶结轻质水泥时弯曲型兰姆波对微环不敏感,在套后耦合常规水泥时微环的存在又会使得弯曲型兰姆波的衰减明显增强,衰减的增加说明其与套后物质的声耦合更紧密,也表明微环与常规水泥的等效声阻抗适当降低,更有利于弯曲型兰姆波传播时与套后介质实现法向刚度的匹配.拉伸波和弯曲型兰姆波对同一胶结状况的不同响应特征也为固井质量综合评价时区分不同的水泥胶结状况提供了可能.     

柱状多层准弹性介质中声压波形的数值计算与分析——声法测井理论研究（Ⅱ）

为逼近油井实际情况,本文处理圆柱状（任意）多层介质中弹性波的传播问题。采用与文[1]相同的物理模型和二维FFT算法,给出了任意层数柱状准弹性介质中全波声压波形的计算方法和程序。 对多种参数进行了声压波形的数值计算,考察了串槽对全波波形的影响,着重分析了两个胶结面串槽引起初至波的变化。结果表明,串槽导致套管滑行纵波幅度明显降低,固、流界面的反射波是影响串槽时初至波特征的一个重要机制。     

随钻电磁波传播方位电阻率仪地质导向关键技术

随钻电磁波传播方位电阻率仪器在钻井过程中可以提供地层边界的方位及距离信息,因此在地质导向应用中发挥着重要作用,是提高油气资源开采率的重要手段.该技术是国外近十年来发展起来的前沿技术,在国内尚属空白.本文首次构思并实现了一种应用"交联天线"的随钻电磁波传播方位电阻率仪器,并从天线结构,测量原理,数据处理以及地层电阻率成像及解释等几个方面详细阐述.交联天线由极化方向正交的线圈串绕组成,同时具有方位探测和地层背景电阻率测量的功能.通过分析交联天线的电压在仪器旋转过程中的变化规律,可以计算出仪器所在地层的电阻率,判断地层边界的方位以及估算地层边界相距仪器的距离.该随钻电磁波传播方位电阻率仪器还包括用于常规电磁波传播电阻率测试的多频多测距补偿天线结构.结合地层电阻率测量及其在方位上的相对变化,可以实现对地层电阻率的全方位成像.该仪器在国内油田进行了多次实井测试,测试结果证明仪器能够在水平井地质导向中准确提供仪器所在地层的电阻率,以及地层边界方位和距离信息.已有实井测试结果表明仪器在油层(100Ωm)和泥岩(5Ωm)中的地层边界探测深度分别为2.2m和1.6m,平均误差在0.2m以内.     

Numerical simulation of seismic wave in elastic and viscoelastic TTI media*

Based on the two-dimensional (2D) three-component first-order velocity-stress equation, the high order staggered mesh finite difference numerical simulation method was used to simulate the elastic and viscoelastic tilted transversely isotropic (TTI) media. The perfect matched layer (PML) absorption boundary condition was selected to eliminate the boundary effect. The results show that: ① Under the condition of fixed elastic parameters of elastic TTI medium, when the polarization angle and azimuth are 60° and 45° respectively, the degree of shear wave splitting is significantly greater than the angle of 0°; ② The influence of viscoelasticity on TTI medium is mainly reflected in the amplitude. If the quality factor decreases, the attenuation of the seismic wave amplitude increases, causing the waveform to become wider and distorted. If the quality factor increases, the viscoelastic medium becomes closer to elastic medium; ③ For TTI medium with different polarization angle and azimuth angle in the upper and lower layers, the shear wave can multiple splits at the interface of medium. The symmetry of seismograms is affected by the polarization angle and azimuth angle of TTI medium; ④ Viscoelasticity has a great influence on reflected wave, transmitted wave and converted wave in the low-velocity model. When the viscoelasticity is strong, the weaker waves may not be shown.     

Reflections from a randomly grooved interface: ultrasonic modelling and finite-diff erence calculation1

This study utilizes ultrasonic water tank modelling to examine the complexities which can result when seismic energy is diffracted by even a simple 3D interface structure. Ultrasonic laboratory modelling is carried out in our in-house water tank, using an aluminum block to simulate very simple ocean-bottom topography. The topography consists of parallel depressions which simulate valleys running parallel to a mid-ocean spreading ridge. This study shows that the phase and amplitude of primary ocean-bottom reflections can vary strongly with the azimuth of propagation. The trailing energy coda also varies dramatically as a function of propagation direction, showing the importance of developing a seismic model which includes the full 3D azimuthal distribution of interface structure.     

Observation of shear‐wave splitting in the multicomponent node data from Atlantis field,Gulf of Mexico

In 2005, a multicomponent ocean bottom node data set was collected by BP and BHP Billiton in the Atlantis field in the Gulf of Mexico. Our results are based on data from a few sparse nodes with millions of shots that were analysed as common receiver azimuthal gathers. A first‐order look at P‐wave arrivals on a common receiver gather at a constant offset reveals variation of P‐wave arrival time as a function of azimuth indicating the presence of azimuthal anisotropy at the top few layers. This prompted us to investigate shear arrivals on the horizontal component data. After preliminary processing, including a static correction, the data were optimally rotated to radial (R) and transverse (T) components. The R component shows azimuthal variation of traveltime indicating variation of velocity with azimuth; the corresponding T component shows azimuthal variation of amplitude and phase (polarity reversal). The observed shear‐wave (S‐wave) splitting, previously observed azimuthal P‐wave velocity variation and azimuthal P‐wave amplitude variation, all indicate the occurrence of anisotropy in the shallow (just below the seafloor) subsea sediment in the area. From the radial component azimuthal gather, we analysed the PP‐ and PS‐wave amplitude variation for the first few layers and determined corresponding anisotropy parameter and V_P/V_S values. Since fracture at this depth is not likely to occur, we attribute the observed azimuthal anisotropy to the presence of microcracks and grain boundary orientation due to stress. The evidence of anisotropy is ubiquitous in this data set and thus it argues strongly in favour of considering anisotropy in depth imaging for obtaining realistic subsurface images, at the least.     

基于二维三分量伪谱法模拟数据的EDA介质中横波分裂研究

为了检测定向裂隙介质中横波分裂的方位属性特征,分析地震属性随裂隙密度和方位变化,采用人工吸收边界和反周期扩展边界,用伪谱法获得不同裂隙密度和不同方位地质模型三分量地面记录;应用时频分析和剪切波偏振分析研究由于裂隙方位和密度引起的横波分裂.结果显示,裂隙密度和方位决定着横波分裂的时差和偏振.快慢横波的延迟时间随裂隙密度增大而增加,不同方位相同裂隙密度的横波分裂时差有微小的变化.在45°方位检测时间延迟时间最大.通过时频分析,可以看到不同方位的瞬时主频有显著的变化,在横波分裂处瞬时主频有明显变化.因此,瞬时主频和快横波的偏振以及延迟时间可以作为裂隙方位和密度的指示.     

Experimental study on the performance of an azimuthal acoustic receiver sonde for a downhole tool

A new azimuthal acoustic receiver sonde with a body and corresponding circuits was designed for a downhole tool. The 64‐sensor receiver sonde holds eight receiver stations that can be combined into at least 64 three‐sensor receiver subarrays. As a result, the receiver sonde can use different sensor combinations instead of different transducer types to produce multiple modes, including a phased azimuthal reception mode and conventional monopole, dipole, and quadruple modes. Laboratory measurements were conducted to study the performance of the azimuthal acoustic receiver sonde for a downhole tool, and the experimental results indicate that the receiver sonde provides a consistent reception performance. Individual sensors receive similar time‐domain waveforms, and their corresponding frequency bands and sensitivities are consistent within the measurement errors of around 5%. The direction of the reception main lobe is approximately parallel to its exterior normal direction. In addition, a receiver subarray with three sensors receives waveforms that have higher energy and narrower beamwidths. For individual sensors, the angular width of the dominant reception lobe is 191.3^° on average, whereas that of the individual receiver subarrays is approximately 52.1^° on average. The amplitude of the first arrival received by the receiver subarray centred at the primary sensor directly pointing to the source is approximately 2.2 times the average amplitude of the first arrivals received by the other receiver subarrays in the same receiver station. Thus, the maximum amplitude of the waveforms received by the receiver subarrays can be used to determine the direction of the incident waves. This approach represents a promising method for determining the reflector azimuth for acoustic reflection logging and three‐dimensional acoustic logging.     

含偏心钻铤充液井孔中正交偶极子声波测井的数值模拟研究

采用三维有限差分方法模拟了正交偶极子声源在含偏心钻铤的充液井孔中激发的声场,研究了钻铤偏心对模式波的种类、激发幅度、以及频散特征的影响.研究结果表明,钻铤偏心导致偶极子声源激发的声场的模式不唯一,观察到了除偶极模式外的单极模式波和四极模式波;钻铤偏心导致偶极模式波出现分裂现象,尤其是快速地层F2模式和慢速地层弯曲波的稍高频率的部分,且快、慢波所对应的两个方位为偏心的方位和与偏心方位垂直的方向;井孔折射横波以及快速地层F2模式的低频部分的慢度基本未受到钻铤偏心的影响,仍然能够正确反映地层的横波慢度及各向异性;对于本文研究的慢速地层井孔模型,当偏心距离l小于等于0.01 m时,弯曲波的慢度和激发幅度受钻铤偏心的影响很小,从快、慢弯曲波中提取的快、慢横波慢度基本能够反映地层的各向异性特征.