《光电信息工程论文.docx》由会员分享,可在线阅读,更多相关《光电信息工程论文.docx(35页珍藏版)》请在三一办公上搜索。
1、光电信息工程论文目: 用声波幅度来确定产液的油水含油率院:电子工程学院部:光电子技术系业:光电信息工程目录摘要3ABSTRACT4前言5第一章绪论61.1声波测井技术的发展现状61.2早期测量油井含水率主要方法61.3声波测井技术在油井测量中遇到的问题71.4未来测井的发展方向7第二章力学网络的分析与计算82.1换能器的位移方程82.2换能器所受力对位移矢量的影响92.2.1薄圆环换能器所受力对质点运动的影响92.2.2求解压电换能器发出的声场对压电换能器的反作用10第三章 换能器电学网络的分析与计算113.1薄圆环压电换能器的物理效应113.2换能器机电等效图分析11第四章脉冲信号的声源冲击
2、响应和电-声传输函数144.1激励电压为脉冲信号144.2接收器声-电冲击响应和声-电传输函数16第五章液体和固体中平面波的反射系数205.1液-固界面的声波传播205.2滑行波作为首波接受的条件205.3求解液-固界面的反射系数21第六章利用反射系数求油井的含油率256.1求解当接收器首先接收到反射波的临界角256.2求解液-固界面的反射系数26结论30致谢31参考文献32摘要石油产业发展的初期的油井一般都是垂直井,在油井内近似认为是油和水的混合 液,我们可以用电容法含水率来测井,它是一种测量井内流体持水率的重要方法,利用 油气同水的介电特性差异测定水的含量。而随着油井的不断开采,出现越来越
3、多的水平 井,电容法测量会出现很大的误差。声波测井的应用,使得误差得到减小,薄圆环压电 陶瓷的换能器是近年发展起来的新型换能器件,它具有体积小、重量轻、精度和分辨率 高、频响高等优点,在垂直井中得到了很好的应用。分析了声波测井油井中薄圆环压电换能器的暂态响应,研究并建立了其等效电路,并且求出了电驱动信号与换能器的瞬时运动之间的传输函数。利用留数 定理,可以对任何电驱动信号,解析地求出换能器的运动在时域与频域的暂态 函数。推导出液-固界面的反射系数,当换能器处于不同的介质时,反射系数是 不同的,因此建立起反射系数与油井含油率的对应关系。这样就很容易的利用薄 圆环换能器测得水平井换能器所在位置的含
4、油率。关键词:声波测井,压电换能器,传输函数,反射系数ABSTRACTOil industry development, the original Wells early commonly in oil Wells is vertical Wells, that is within approximation of oil and water mixture, we can use capacitive moisture content to logging, it is a kind of measurement in well hold water rate of fluid import
5、ant method, using the dielectric properties with water oil determination of water content of differences. With the continuous mining Wells, appear more and more horizontal well, capacitive measurement may appear in large error. Acoustic logging application, make errors, thin rings have reduced the p
6、iezoelectric transducers is developed in recent years of new transducer parts, it has small volume, light weight, precision and high resolution, the advantages of high frequency response in vertical Wells, get a good application.This paper analyses the transient response of thin shell transducer, de
7、velops its equivalent circuit and determines the relation between an electric-driving signal and a radiated acoustic signal. By using the residue principle, the transient function s of the temporal and frequency domains for the transducer are determined analytically for any electrical-driving signal
8、. Deduces liquid-solid interface of reflection coefficient, when in different medium transducer, reflection coefficient is different, so set up reflection coefficient and the corresponding relations opu oil Wells. This is easy to utilize thin rings transducer horizontal location measurement transduc
9、er the opu.Key word : sonic-logging, piezoelectric-transducer, transient function, reflection coefficient前言声波测井技术起始于20世纪40年代,经过70多年的发展,并且在计算机产业和 数字处理技术迅速发展的带动下,加快了对声波测井技术的不断优化,可以或得更多的 油井声学性质,提高了声波测井在油井中测量的准确度。声波测井是地球物理测井技术 中一种发展最快应用最广泛的现代测井方法之一。该技术采用了声波的速度及幅度在岩 石、岩层孔隙、含油气水等介质中传播时的全波记录得出的不同物理地质特征来研究和
10、 解决地质问题,进而发现油气、煤、金属与非金属、放射性、地热、地下水等矿产资源。 近年来已扩展到工程地质、灾害地质、生态环境、考古研究等应用领域。声波测井已经 成为地球物理测井科学的重要领域,是发展最快和应用最广泛的现代测井方法之一。根据目前国内外的声波测井技术的研究现状,基于声波测井井壁反射波的简单工程 模式,本文在利用声波测井传输网络模型计算声源的电-声冲击响应和接收器的声-电转 换传输函数。利用数字信号处理理论计算了油井中声场传播模型,对电压驱动信号激励 薄圆环压电换能器产生声信号作为点声源,分析计算得出了纵波、横波和反射波的信号 模型,计算并推导出薄圆环换能器和接收器传输函数,建立声波
11、测井传输网络之后,着 重研究分析了液-固界面的反射系数在不同的介质中的变化情况,根据反射系数与面积 的对应关系,确定油井中的油水含油率。在本论文中,我通过查阅资料和向法老师请教,自己也做了一些理论研究,着重分 析了建立的薄圆环压电换能器的机电网络等效图,得出了声源的电-声冲击响应和传输 函数和收器的声-电转换传输函数。着重研究分析了液-固界面的反射系数在不同的介质 中的变化情况,得出反射系数与面积的对应关系,并且用matlab仿真软件确定了反射系 数与油井含油率的对应曲线图。第一章绪论1.1声波测井技术的发展现状测井技术于1927年起源于法国斯仑贝谢公司,1939年中国开始将测井技术应用于 石
12、油工业。80多年来,测井技术经历了从模拟测井到数字测井、数控测井、成像测井 的发展历程,目前已发展成为十大石油学科之一。它被广泛应用于油气田的勘探与开发 的全过程,为石油地质和工程技术人员寻找和评价油气层提供了重要的资料和数据,解 决了一系列的地质难题。同时,测井技术还是进行煤炭、金属等矿产资源勘探的重要手 段,并被扩展到工程及其它领域。近几年来,我国斥巨资直接请斯仑贝谢测井公司服务,并引进了几十套包括MAC在 内的ECLIPS2000系统。借鉴国外在声波测井仪器研制方面的成功技术,开发新一代偶极 声波测井仪器是提高我国测井装备水平的重要途径。石油勘探技术发展迅速。主要表现 为:地面记录系统向
13、高性能大型复合型方向发展;声、电、核、磁等各系列的井下仪器 全面向成像化方向发展,尤其是核磁成像测井技术,发展特别迅速;测井资料处理解释 技术向解决实际问题的个性化方向发展;测井软件技术则向大型综合性方向发展。为了 更加准确的测量油水的含油率,声波幅度测井作为一种较为常用的测量数据的方法,也 得到了很广的运用。1.2早期测量油井含水率主要方法目前生产井主要通过测量井内持水率及流体密度来识别流体。测量混合流体持水率 的仪器主要有电容式含水仪、微波含水仪、开关电导式含水仪以及压差密度计等。这类 测量所得的数据可以识别井内流体的类型,从而确定油井的含油率。电容法含水率测井是利用油气同水的的介电特性差
14、异测定水的含量。由于油气等碳 氢化合物与水具有不同的介电常数。水的介电常数为80,油气的介电常数1.04.0,因 而具有较高的分辨率,仪器是对所处位置的流体进行采样,然后在仪器的取样室内油水 靠重力分离,该室内置一电容电极,电极与取样室外壳构成圆柱状电容,通过测量圆柱 状电容器电容量的变化就可以得到持水率,确定含油率。微波含水率测井是利用高频电磁波的谐振状态来测量原油中的水分。高频电磁波在 含水原油中传播时,其波长随含水率的不同而不同,并引起谐振回路频率的变化,致使 改变谐振回路和晶振回路之间的谐振状态,根据这种变化可以测定原油中的含水率。开关电导法测井是以电导为基础,由于油的电导趋于零,水的
15、电导较大,当油、水 以各种不同的百分含量流过探头的电极时,探头都会有相应的输出电压,这样得出各种 百分比的对应值,从而确定油水的含油率。压差密度计又称密度梯压计,利用两个0.6m的压敏元件,测量井筒内流体柱两点 间的压力差值,而测出的压力梯度正比于流体密度,这样就可以测得流体的密度,确定 含油率。1.3声波测井技术在油井测量中遇到的问题在油田勘探开发中,测井技术是确定和评价油层的重要手段,也是解决一系列地质 问题的有效途径,但是油井计量主要存在以下问题:油井产液量普遍波动较大,且没有规律可循,低产油井还存在产液间歇现象,间歇 时间长短不一,短时计量很难得出真实的产量,采用两相计量分离器配玻璃管
16、量油的方 法已难以适应。油井产量差别大、波动大导致不同油井伴生气的产量差别和波动也较大, 孔板配双波纹管差压计等测气方法已不能适用。由于采用高液量生产和油层压裂等措 施,增加了油井产液中的砂和杂质,含水率高又使仪表易结垢,齿轮流量计、腰轮流量 计已不适用。由于油井产液含水率高,有大量游离水存在,因此,人工取样随机性很强, 很难取准,无法得到真实、稳定、准确的含水率值。在声波测井发展到目前这种状态,无论是研究新的测井方法,还是结合其他学科, 对现有的技术进行应用拓展和改进,都是极为艰难的。现在声波测井面临以下难题:如何摆脱信息量小的制约是一个关键的问题,再一个就是如何在复杂的介质条件下 加大对测
17、井问题的解决,使得地质情况对测量的影响降低到最小,以及如何使得采集的 测井信息更加真实的反映油井的原始状态。提高测井的效率也是需要解决的一个问题。1.4未来测井的发展方向未来测井的方向是要提高测井技术,发展井间测井新技术,能够提高测井的纵向分 辨率,并且发展综合技术,提高测井技术综合应用能力,以及要提高测井综合信息采集 能力,并且提高油气采收率测井技术,同时提高测井的深度做好剩余油的检测。这都是 未来声波测井的发展要求及方向。光电信息工程论文第二章力学网络的分析与计算2.1换能器的位移方程薄圆环压电换能器的剖面图如图2-1所示,为了研究我们引入柱面坐标系对 换能器的工作原理进行分析与计算。当平
18、均半径为,,壁厚为lt沿着半径方向进行极化,设它的密度为 p。由于 薄圆环压电换能器的壁较薄,因而我们可以近似地认为 g r( r是换能器的半径)。b由于它的质点位移具有轴对称性,切向应力为零,所以在薄圆环内不能形成应力波,若此换能器的外壁自由,则可得简化方程:(2-1)(2-2)其中T、七分别是薄圆环切向和轴向的正应力,ur和uz分别是相应径向和轴向的质点位移分量。r Z(2-3)(2-4)设0,乙r分别由数字1,2,3代表,则与径向极化有关的压电方程为1:S = Set + d E111 131 3D =d T+8tE331 133 3式中S1是薄圆环在0方向上的伸缩应变量,S是换能器材料
19、的柔顺系数,e r和d 分别是换能器的介电常数和压电常数,E3是电场的径向分量。2.2换能器所受力对位移矢量的影响2.2.1薄圆环换能器所受力对质点运动的影响当薄圆环处于耦合液体中,并使它的两端与周围耦合液体隔离,只让其外壁 与液体相接触,当换能器在液体中径向振动时,他会交替引起液体的膨胀与收缩 运动,在声波测井仪中,声源通常是放在耦合液体中,因而会向外辐射声波。同 时,换能器也处于自己的辐射声波中。所以,它受声场的反作用力,该反作用力与薄圆环换能器的运动状态有关,用 F表示,除此之外,薄圆环的振动还会引起 摩擦力阻Ff,它与质点的位移速度和液体的摩擦力阻 七成正比,方向与振动方 向相反,如下
20、:m(Rk2r: + jkrb1 + k 2 r 21 + k 2 r 2bb+ jX、r dtdU(2-5)式中R = 2Kr HU p ,p和。分别是耦合液体的密度和声速,R和X分别b m mm mr r是辐射阻和辐射抗。F j =- R ( 2-6)其中,Rm与液体粘滞系数成正比。薄圆环换能器声辐射面积 A = 2兀rH,则它所受外力:(2-7)F = F + F = _(R + R + jX )dUrUrrb由于薄圆环换能器的轴对称性,质点位移和伸缩变量的关系为:(2-8)由于薄圆环换能器高度H 0, D = 0三种情况:1)D 0 ,等式s3+ as 2 + bs + C = 0有三
21、个不相等的实根*、%、 声源为过阻尼模式,声源的冲击响应为:h (t)= h (t) + h (t)+、(t) = K est + K2g + K3es3t对参数D 0,D 0,有一个实根s1 = a和两个复根s23=a8八声源为振荡模式,声源的 冲击响应为h (t)= K e-at + 2 K e-Pt cos 0t-。)(4-6)其中,a = a (x + y)3、3( x y)24( x + y a)27(x + y)2 + 3o 24d .J心 2 (2x + 2 y + + 1) + X y (3 + a)V 326 = arctanJ 2+3(x+ y)2+2a(x+ y) 4据+
22、 y )+ 3加载中心频率为w3(x - y)J =2兀4兀4.2接收器声-电冲击响应和声-电传输函数在发射声波中,我们可以利用机电等效电路来分析,实现力学和电学的结合, 如图4-1所示,左边是力学网络,右边是电学网络。虬3 mr3 Rm3 叫 J N3AjI IIrvwx图4-1接收器的机电等效电路图薄圆环接收器的冲击响应的拉普拉斯变换表达式:H (s)= -sd(4-7)-s 3 - as 2 + bs + c对(4-7)式进行傅里叶变换可得接收器的传输函数表达式:H (jw)=jwdjw3 C02 + jwb + c(4-8)对于三次方程s3 + as2 + bs+c = 有 实用数学手
23、册4:2a 3 ab273(4-9)1) D V 0时,有三个不相等的实根S 1,S 2 ,为换能器的过阻尼模式x+y 5S3 =2 - a4 ad 1227 C + y2 a + 3 (x +18 a (x + y3(x + y ) 10 a ,侦、s ) = H 湛 ( ) + H12 123(x + y )d36 a (x + y )d(,+ H 13( S)(4-10):T d 24 a 2 - 18 a (x + y )(4-11)h (t) = h (t) + h (t) + h (t)3111213=K est + K es t + K est112232)D = 0时且;2=(
24、P33)A 0时,三个实根中有两个相等为换能器的临界阻尼模式S=-(x+y)2 a + 3 (% + y ), K11 =9(% +y)K2112 (x + y )- 4 a27 (x + y)2-222 a + 3(x + y9(x + y )G)二sd(s - s )(s - s )2 H11 ()+ H 21 ts) + h (t) + h(t)+ h(t)112122(4-12 )H 22 (s )=K11 e叩 +( K 21t + K 22)心2t(4-13 )3) D 0时,有一个实根 和两个复根% ,为换能器的振荡模式s = -as P + j s3 P 3 j气A (x +
25、y )2B - G - y) 2P a + A3 3w = :3 B-adC 9A2 +w 2d睥2D + 3A 2 w(9中arct an w 2 + 3 P A w(P + 3 A )(4-14)H 3( s)=(s s )G ds )G s )H11 (s )+ H12 + H 21 3H (s )与H (5 )互为共轴J. Z匕 _Lh = Ce 一喝 + De -阶 cos(o 1 + 中)(4-15 )3接收器的中心频率为L =0)/271 =必(x- y)/47i光电信息工程论文第五章液体和固体中平面波的反射系数5.1液-固界面的声波传播反射系数:有一列平面简谐波从一种介质进入另
26、一种介质,则会在界面处发生反 射现象。假设入射波的振幅A0,反射波的振幅A1,则定义反射系数R = A1。A 0斯奈尔定律:声波的反射和折射遵循一定的规律,即斯奈尔定律:6(5-1)VVV打0= 70n = 一0Ti Zn VpsV为介质中油气层纵波速度、钢管纵波速度、钢管横波速度、界面),气为入射角、钢管纵波折射角、钢管横波折射角;斯奈尔定律是式中V,V上的相速度;0描述波在弹性分界面上发生反射、折射后波的传播方向的定律,可解释为入射波、反射 波和折射波沿分界面视速度相等。5.2滑行波作为首波接受的条件在井中离反射源一定距离的接收点可接收到的波有直达波、反射波和滑行波等。在 井中放置一居中发
27、射器T和接收器R,发射器与接收探头距离(源距)为L,井眼半径 为a,油气层和钢管速度分别为*油气和*钢,见图5-1。图5-1井内声波传播示意图依照几何声学理论,首先,比较直接波和反射传播时间,由于直达波和反射波在油 气中的传播速度都一样,而直达波的传播路径比反射波传播路径都短,因此对于直达波和反射波而言,直达波先到达接收器。但是我们可以在发射器和接收器之间加上和一种 材料,使得直达波不能通过。其次,比较直达波与滑行纵波达到接收器的时间,只要选 择合适的源距就能使反射波比滑行纵波先到达接收器。6直达波达到接收器所需要的时间为:( L、2(5-2)W2 + T_ k 2 JDV油气设第一临界角为0
28、,滑行纵波经过井眼,沿井壁滑行后折回井内,最后到达接收器所需 要的时间为:2 aL 一 2 a tan 0p V 油 cos 0V要使反射波先于滑行纵波到达接收器,应使TD T,即: 1p!/ r、2li a 2 + k 2 )2 aL 一 2 a t an 0ZT +V油气V油气cos 0V钢由于sin 0=V油气,因此发射器和接收器之间最大源距为:V钢L = 2 a 钢油气卜钢-V油气(5-3)(5-4)(5-5)5.3求解液-固界面的反射系数一般情况下,固体中有纵波和横波两种形式,当声波投射到液体-固体分界面上时 引起反射和折射,在固体中同时产生折射的纵波和横波,而流体中只可能有反射纵波
29、。设液体中单频入射声波的位移电势写作:(5-6) 其中,A一入射声波的振幅;0 一入射角;k =s.c,c 一液体中波数和波速。见图5-2。G, Z)= Aejk(xsin0-zcos0)图5-2液-固界面的声波的反射和折射在液-固界面处S 3为钢管折射横波,p为油层的入射波,p为油层的反射波,则:1P2为钢管的折射纵波,e j( t - kx sin 0 + kz cos 0 )e j ( t - kx sin 0 - kz cos 0 )(sin 0 ) I- cos 0 I p /e j ot-kpx sin 0 p + kpz cos 0 )e j (o t - k x sin 0,.
30、 + k z cos 0 )R一反射波的反射系数,R2、R3 一折射波的折射系数。 由界面应力平衡条件和界面垂直位移连续条件得:其中,sin0 + R sin0 = R sin0 - R cos012 p 3 s(5-7)一 cos0 + R cos0 = -R cos0 一 R sin012 p 3 s(5-8)对于油层一侧的液体而言,根据介质的本构方程推导应力表达式,由本构方程:T = c: S,S=V U(5-9)s得:T=c:VsU-T - CCC000 一1111212TCCC0002121112TCCC0003=121211T000C00444T0000C0544T00000C1-
31、6J44 -1ad x000aa zaa yaa y0aa z0aa x在XZ平面上,uy= ay= 在液体中 c 120aa zaa yaa x0UxUyUz=0,气4 = 0代入液体的cJ】:Tc0111T0c211T003=T004T005T006J000000c001100000000000a 云0a0a za0aa za x0Ta u(0 )1 xa x00u(0)z_对于入射波,可得:T(0) = C11a u(。)ITc11(jk cos 0 )(- cos 0 )= - jkc cos 2 0(5-10)对于反射波,可得:a u(1) c z11 a zc11R (-1jkco
32、s0 )cos 0 = - jR. kccos 2 0(5-11)11可得:000T(1)=3c、12c、11c、12c、12c、12c、11000c、0044T 3=120000c、44000000aaxaazaazaaxc、44aazaaxu (1)x0u(1)za u(1)xa xT(2) = -jR2k、11a u(1)11 a z一 2c、sin 20 )44p(5-12)(5-13)T(3) = - 1 jR k sin 29323 s(5-14)由声波在液-固界面应力平衡得:T(0) + T(i)= T(2) + T(3)kc cos2 9 + kc R cos2 9= R k C、一 2c、sim9 )1 R k sin291111 12 p 1144 p 2 3 s s、+ c、)1112(5-15)(5-16)联立式(5-7)、(5-8)和(5-16)解得反射系数为:r - z + z cos 2 20 + z sin 2 29 z + z cos 2 20 + z、sin 2 20,(5-17)第六章利用反射系数求油井的含油率通过查声学手册得到以下数据:石油的密度P油800 kg/ m3水的密度P水1000 kg/ m3钢的密度7900 m/s声波在钢中的速度D钢5000 m/s声波在水中的速度u水1500 m/s