一种多根钻杆钻柱中的纵波传输特性分析.docxVIP

一种多根钻杆钻柱中的纵波传输特性分析 为了确保钻头安全，提高工作效率，降低钻孔成本，必须实时测量和传输钻头设计参数和地质参数，从而促进钻孔和井罐评估领域越来越多的应用需求。尽管国外利用随钻测井技术几乎实现了所有电缆测井的功能,但遥测技术仍然是随钻测井技术发展的瓶颈。利用钻井液脉冲波的方法实现信息从钻头到地面的传输速率很低,远远小于人们的期望值。电磁遥测钻杆系统的传输速率理论上可以高达2 Mbit/s,但这种遥测方式需要对现有钻杆和接箍进行特殊设计和改造,花费巨大。通过由多节钻杆和接箍构成的空心金属钻柱中传播的弹性波来传递井下信息,因声波在金属钻柱中传播的速度比在钻杆钻井液中的要大,可以传播的弹性波信号的主频(0.4～2 kHz)更高,因而具有更强的载波能力。国外的研究人员已经在这个领域展开了大量的研究工作[10,11,12,13,14,15,16,17,18],Halliburton公司于2000年声称其声波遥测系统(ATS)进入商业化,2004年公布了部分现场试验结果。与国外的研究相比,国内在钻柱声波遥测方面的研究刚刚起步。以往的随钻声波遥测理论的研究基本采用有限差分法,笔者采用传递矩阵法对钻柱声传播特性进行数值模拟,并将所得结果与利用有限差分法得到的结果进行对比分析。 1 密度ac2kac2k2cos 钻柱模型由多根长度分别为d1,d2,…,dN的钻杆和接箍连接而成(图1)。采用传递矩阵法进行分析时,首先考虑钻柱的中间任一段存在势函数Φ,其表达式为 Φ=Φiejkx+Φre-jkx. (1) 式中,Φi和Φr分别为入射波和反射波势函数的幅度;k为纵波波数。式(1)中的势函数均略去了公因子ejωt,其中j表示虚数单位,ω表示声波的角频率。 把钻杆截面的法向位移u和法向作用力F用矩阵形式可以表达为 [uF]=A(x)[Φi+ΦrΦi?Φr].(2)[uF]=A(x)[Φi+ΦrΦi-Φr].(2) 其中 A(x)=[?ksin(kx)ρac2k2cos(kx)jkcos(kx)ρac2k2jsin(kx)].A(x)=[-ksin(kx)jkcos(kx)ρac2k2cos(kx)ρac2k2jsin(kx)]. 式中,ρ为密度;c为纵波波速;a为横截面积。 把x=0和x=dN代入式(2),并推导可以得到 [uF]x=d?N=MN[uF]x=0+,(3)[uxFx]x=0?=MN?1[uxFx]x=d+N?1,(4)[uF]x=dΝ-=ΜΝ[uF]x=0+,(3)[uxFx]x=0-=ΜΝ-1[uxFx]x=dΝ-1+,(4) 继续推导可以得到 [uF]x=?(h?d1)?=M1[uF]x=?h+.(5)[uF]x=-(h-d1)-=Μ1[uF]x=-h+.(5) 其中 MN=A(dN)A?1(0),h=∑i=1N?1di.ΜΝ=A(dΝ)A-1(0),h=∑i=1Ν-1di. 考虑到分界面处的法向位移和法向作用力连续的条件,即有 [uF]x=d+=[uF]x=d?.(6)[uF]x=d+=[uF]x=d-.(6) 于是,可以得到 [uF]x=d+N=M[uF]x=?h?.(7)[uF]x=dΝ+=Μ[uF]x=-h-.(7) 其中 M=MNMN-1…M1. 对整个钻柱来说,用ΦI和ΦR表示入射波和反射波势函数的幅度,ΦT来表示透射波势函数的幅度,则定义反射系数R和透射系数T分别为 R=ΦR/ΦI,T=ΦT/ΦI. 假设 M=[m11m21m12m22]Μ=[m11m12m21m22] ,经推导可得到 [jk?ρac2k2jm11k+m12ρac2k2jm21k+m22ρac2k2]?[RT]=[jkρac2k2].(8)[jkjm11k+m12ρac2k2-ρac2k2jm21k+m22ρac2k2]?[RΤ]=[jkρac2k2].(8) 因此,从式(8)中可以求解反射系数和透射系数,从而得到钻柱中声传播的频率响应曲线。 2 计算值的结论 分别采用传递矩阵法和有限差分法对不同长度的钻柱模型在脉冲函数激励下的频谱进行了模拟计算,使用的钻杆和接箍的参数如表1所示。 2.19 通带内的传播 由10个接箍、9根钻杆组成的钻柱的总长度为82.75 m,模拟计算得到的脉冲函数激励下的频谱曲线如图2所示。 可以看出:除了第一个阻带外,其余的阻带都有效地阻止了带内信号的传播;相比而言,通带内的信号却可以几乎不衰减地进行传播。图2(c)所示的是图2(a)和图2(b)中第五个通带的放大视图,可以看到,通带内还存在很多小的谱峰,每个通带内都由间距大约为20Hz的一系列谱峰构成。这里可以给出一个定性的解释:谱峰数量与钻柱的钻杆根数有关,即存在几根钻杆就有几个谱峰。对于这种现象的解释是,每个通带的宽度由为了匹配每根附加的半个波长而引起的