气体同位素质谱与色谱结合分析天然气组成.docVIP

 气体引起圆柱谐振腔共振 频率偏移的数值模拟 徐德龙1 ,2 , 王秀明2 , 宋延杰1 , 丛健生1 , 陈德华1 ( 1 . 大庆石油学院 地球科学学院 ,黑龙江 大庆 163318 ; 2 . 中国科学院 声学研究所 ,北京 100080 ) 摘 要 :基于共振声谱法 ,研究油中存在气体时引起圆柱谐振腔共振频率偏移的问题 ; 求取单相油存在时圆柱谐振 腔的声势与共振频率的表达式 ; 采用等效速度和等效密度的方法 ,研究气体存在时谐振腔共振频率的偏移 ,并得出理论 模拟曲线. 结果表明 : 随着圆柱谐振腔内气体积分数的增大 ,谐振腔的共振频率明显降低 ; 当气体体积分数增大到 30 % 左右时 ,谐振腔内气的影响是主要的. 此时谐振腔的共振频率与谐振腔内完全充满气体时的共振频率差别不大. 该结果 可解决油井内多相流问题 ,特别是对气体流量的测量有重要意义. 关 键 词 :气体 ; 共振频率 ; 多相流 ; 谐振腔 ; 数值模拟 中图分类号 : P631 文献标识码 : A 文章编号 :1000 1891 ( 2007) 04 0016 03 0 引言 在油井中 ,井内流体绝大多数是多相流. 为确定油井内的多相流流动特征 ,许多理论和方法被用于油 井多相流的研究[ 1 - 3 ] . 共振声谱法被用于高精度测量液体的流速[ 4 ] . 1998 年 , Ha r ri s J M 首次提出利用 共振声谱法测量岩石声学特性[ 5 ,6 ] . 徐德龙等在理论方面对存在扰动体时球形和圆柱形谐振腔共振频率 偏移的模拟计算做过研究[ 7 , 8 ] ;在实验方面 ,讨论了利用共振声谱求取小扰动体声速的 2 种方法[ 9 ] . 王秀 明等曾利用数值模拟方法研究了套管井中声波的共振谱问题 ,并讨论了套管井中水泥胶结状态对共振声 谱的影响[ 10 ] . 李铭等对声频谱测井方法进行了研究 ,并讨论了该方法在确定储层参数 、识别油气等方面 的应用[ 11 ] . 利用共振声谱法进行油气两相流测量的研究还未见报道 . 笔者尝试利用共振声谱法原理研究 气液两相存在的多相流测量问题 . 首先研究单相油存在时圆柱谐振腔的共振频率测量原理 ,然后采用等 效速度和等效密度的方法研究气体存在引起谐振腔共振频率的偏移问题 ,并给出气体体积分数不同时的 理论模拟曲线. 该研究避免了烦琐的理论推导 ,更有利于模拟计算和现场应用的实现. 1 共振频率 圆柱谐振腔模型见图 1 . 设圆柱坐标系 ( r ,θ, z ) 的坐标原点是圆柱谐振腔 的底面圆心 , z 轴与圆柱谐振腔的中心轴重合. 圆柱的半径为 a , 上 、下底面为 z = 0和 z = L . 谐振腔内介质的密度为ρ, 声速为 c. 圆柱谐振腔的边侧面为刚性 , 上 、下底面为应力自由 . 设圆柱谐振腔的声 势为 . 根据数学物理方程并结合边界条件 , 可以得出此时谐振腔声势解析式 : ω z z ωr r 6 φm ( mθ) Zn ( r ,θ, z) = R m , m , n = 0 , 1 , 2 , , c c 图 1 圆柱体谐振腔 模型示意 ω ωr r πl c z z 其中 ,φm ( mθ) , Zn 分 别 为θ, z , r 方 向的 本 征模 式. ωz = , R , m c c L 收稿日期 :2006 03 21 ; 审稿人 :蔺景龙 ; 编辑 : 陆雅玲 基金项目 :国家自然科学基金 ( ;黑龙江省自然科学基金 ( A2004 - 3) ;大庆油田有限责任公司研究生基金 ( 2006) 作者简介 :徐德龙 ( 1976 - ) , 男 ,博士生 ,讲师 ,主要从事岩石物理声学测量 、地球物理方法原理及声空化方面的研究. πcαmn ωr = , l = 0 , 1 , 2 , . 由此可得圆柱谐振腔共振频率 f mn 的表达式 : a 2 2 α c l mn ( 1) f mn = + , 2 L a 式中 :αmn 是 J m (πr) 的第 n 个零点值 ,J m 是 m 阶柱贝塞尔函数 . 扰动体在所建立的柱坐标系中是轴对称的 , 所以声势函数和共振频率应与角度θ无关 , 即为 0 . 此时 沿 z 轴方向传播的声波共振频率的最低模式为 l = 1 和 n = 0 . 共振频率 f 0 与本征模式的表达式分别为 c ( 2) f 0 = , 2 L n = Cn si n ω J 0 z z ωr r ( 3) , c c 式中 :J 0 为零阶柱贝塞尔函数 ;