CN103821499B - 用于油井动液面深度检测的声音信号处理方法 - Google Patents
用于油井动液面深度检测的声音信号处理方法 Download PDFInfo
- Publication number
- CN103821499B CN103821499B CN201410059106.9A CN201410059106A CN103821499B CN 103821499 B CN103821499 B CN 103821499B CN 201410059106 A CN201410059106 A CN 201410059106A CN 103821499 B CN103821499 B CN 103821499B
- Authority
- CN
- China
- Prior art keywords
- oil well
- difference
- fluid level
- fft
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于油井动液面深度检测的声音信号处理方法,通过系统参数设置之后采集油井声场信号,将获取的声音信号进行窗函数处理、快速傅里叶变换(FFT)、频谱信号处理等手段得到油井声场的各阶共振频率,从而计算各相邻两阶共振频率差值的平均值,根据管柱声场模型采用共振频率差值法推出深度计算公式,计算出油井动液面深度并显示出来。其显著效果是:本发明采用管柱声场模型推出的共振频率差值法模型,计算油井动液面的深度,硬件搭建简单,体积小,能够有效地避免油井中声场信号检测过程中发生的整体波形移频,不但减小了体积,还提升了油井动液面检测的精度。
Description
技术领域
本发明涉及到油井动液面深度检测技术领域,具体地说,是一种用于油井动液面深度检测的声音信号处理方法。
背景技术
在石油开采的过程中,通过检测油井动液面深度,能够科学地了解油井供应能力,确定抽油泵的沉没深度、油层压力,分析能量衰减的异常原因等,从而合理安排采油工艺,使油井产油率最大化。因此,油井动液面深度的检测在油田开发中显得十分重要。
油井动液面深度检测的一种新方法是管柱声场模型法,该方法是在井口利用声源发生器产生的声波沿油管、套管内的环行空间向下传播,遇到液面产生反射波,从而在管柱内形成声场,利用井口拾音器接收声场波,并对数据处理之后利用声场模型计算液面的深度。
然而,接收到的声场波中通常存在各种干扰,如果不进行专门的声音信号处理,将会对检测精度产生很大的影响,从而影响生产,造成损失。
发明内容
针对现有技术的不足,本发明的目的是提供对接收到的声音信号进行处理,实现油井动液面深度检测,并提高检测精度的声音信号处理方法。
为达到上述目的,本发明表述一种用于油井动液面深度检测的声音信号处理方法,其关键在于按照如下步骤进行:
步骤1:系统参数设置,包括设置声场信号采集的采样频率、窗函数种类、快速傅里叶变换(FFT)点数和AD输入量程;
步骤2:从管柱声场模型中采集油井声场信号的时域序列信号;
步骤3:窗函数处理,即从步骤2获取的时域序列信号中截取一个时间片段,然后将其进行周期延拓处理,得到虚拟的无限长的时域信号;
步骤4:快速傅里叶变换(FFT),即将步骤3中得到的时域信号进行快速傅里叶变换(FFT),将时域序列转换成对应的频域序列;
步骤5:频谱信号处理,即将步骤4中得到的频域波形进行求导处理,找到导数为零的各个点所对应的频率,即为油井声场的各阶共振频率;
步骤6:计算各相邻两阶共振频率之间差值的平均值Δf,即根据步骤5计算出各相邻两阶共振频率之间的差值,然后采用中值滤波和均值滤波算法得到各相邻两阶共振频率之间差值的平均值Δf;
步骤7:油井动液面深度计算,即根据步骤6得到的平均差值Δf,按照计算出油井动液面的深度l,其中,c为油井中的声音传播速度,Δf为各相邻两阶共振频率之间差值的平均值,d表示管柱的直径。
作为进一步描述,所述步骤3中的窗函数为矩形窗、汉宁窗或高斯窗,其中:
矩形窗函数为:
汉宁窗函数为:
高斯窗函数为:
作为进一步描述,所述步骤6具体按照以下步骤进行:
步骤6-1:将步骤5中得到的K阶共振频率点中的相邻两个点分别求差值,得到K-1个共振频率差值数据;
步骤6-2:将K-1个共振频率差值数据进行升序排序,然后采用中值滤波技术取出中间的J个差值数据;
步骤6-3:采用均值滤波技术对J个差值数据求平均值,得到共振频率平均差值Δf。
本发明的显著效果是:采用管柱声场模型推出的共振频率差值模型,计算油井动液面的深度,硬件搭建简单,体积小,能够有效地避免油井中声场信号检测过程中发生的整体波形移频,不但减小了体积,还提升了油井动液面检测的精度。
附图说明
图1是本发明的原理框图;
图2是步骤4中时域序列变换为频域序列波形图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
参见附图1,一种用于油井动液面深度检测的声音信号处理方法,其具体步骤为:
首先进入步骤1:系统参数设置,包括设置声场信号采集的采样频率、窗函数种类、快速傅里叶变换(FFT)点数和AD输入量程;
其中,在实施过程中采样频率和快速傅里叶变换(FFT)点数可以自由的输入,但一般默认采样频率为2000HZ,可根据需要进行修改,快速傅里叶变换点数一般默认为1024点,可根据需要进行修改;窗函数种类提供三种选择,分别为矩形窗、汉宁窗和指数窗,默认为矩形窗;AD输入量程同样提供四种选择,分别为±12V档、±5V档、±2V档和±200mV档,在使用过程中可根据检测到的声音信号增益进行合适的选择;
然后进入步骤2:从管柱声场模型中采集油井声场信号的时域序列信号,即配合MP420E采集卡、前端拾音器、抗混叠滤波器使用,负责将经过抗混叠滤波器滤波之后的油井声场信号采集到固定的内存中,以供软件其它部分使用;
然后进入步骤3:窗函数处理,即从步骤2获取的时域序列信号中截取一个时间片段,然后将其进行周期延拓处理,得到虚拟的无限长的时域信号;所述窗函数为矩形窗、汉宁窗或高斯窗,其中:
矩形窗函数为:
汉宁窗函数为:
高斯窗函数为:
其中各个窗口的长度M为1000个采样点;
然后进入步骤4:快速傅里叶变换(FFT),即将步骤3中得到的时域信号进行快速傅里叶变换(FFT),将时域序列转换成对应的频域序列;
其中离散傅里叶变换为:
其中N为快速傅里叶变换的点数,WN为旋转因子,呈现周期性、可约性和对称型,具体表述为:
周期性
对称性
可约性
利用旋转因子WN的周期性、对称性和可约性即可在离散傅里叶变换的基础上大大简化计算步骤,实现快速傅里叶变换(FFT),具体的快速傅里叶变换(FFT)之后的参考数据如图2所示;
然后进入步骤5:频谱信号处理,即将步骤4中得到的频域波形进行求导处理,找到导数为零的各个点所对应的频率,去除波谷的点,即为油井声场的各阶共振频率;
然后进入步骤6:计算各相邻两阶共振频率之间差值的平均值Δf,即根据步骤5计算出各相邻两阶共振频率之间的差值,然后采用中值滤波和均值滤波算法得到各相邻两阶共振频率之间差值的平均值Δf;
其具体步骤为:
步骤6-1:将步骤5中得到的K阶共振频率点中的相邻两个点分别求差值,得到K-1个共振频率差值数据;
步骤6-2:将K-1个共振频率差值数据进行升序排序,然后采用中值滤波技术取出中间的J个差值数据,通常J的取值为5;
步骤6-3:采用均值滤波技术对J个差值数据求平均值,得到共振频率平均差值Δf。
最后进入步骤7:油井动液面深度计算,即根据步骤6得到的各相邻两阶共振频率差值的平均值Δf,采用根据管柱声场模型推出的共振频率差值公式,计算出油井动液面的深度l,其具体步骤为:
其中管柱声场模型为:
其中,c为油井中的声音传播速度,fn表示第n阶共振频率,d表示管柱的直径,l表示油井的深度。
根据管柱声场模型则可以推出
则Δf为
则由管柱声场模型推出的油井动液面深度计算公式为:
该方法主要适用于管柱声场模型中声音信号的处理,提高计算的准确性。
Claims (1)
1.一种用于油井动液面深度检测的声音信号处理方法,其特征在于:按照以下步骤进行:
步骤1:系统参数设置,包括设置声场信号采集的采样频率、窗函数种类、快速傅里叶变换(FFT)点数和AD输入量程;
步骤2:从管柱声场模型中采集油井声场信号的时域序列信号;
步骤3:窗函数处理,即从步骤2获取的时域序列信号中截取一个时间片段,然后将其进行周期延拓处理,得到虚拟的无限长的时域信号;
步骤4:快速傅里叶变换(FFT),即将步骤3中得到的时域信号进行快速傅里叶变换(FFT),将时域序列转换成对应的频域序列;
步骤5:频谱信号处理,即将步骤4中得到的频域波形进行求导处理,找到导数为零的各个点所对应的频率,即为油井声场的各阶共振频率;
步骤6:计算各相邻两阶共振频率之间差值的平均值Δf,即根据步骤5计算出各相邻两阶共振频率之间的差值,然后采用中值滤波和均值滤波算法得到各相邻两阶共振频率之间差值的平均值Δf;
步骤7:油井动液面深度计算,即根据步骤6得到的平均差值Δf,按照计算出油井动液面的深度l,其中,c为油井中的声音传播速度,Δf为各相邻两阶共振频率之间差值的平均值,d表示管柱的直径;
所述步骤6具体按照以下步骤进行:
步骤6-1:将步骤5中得到的K阶共振频率点中的相邻两个点分别求差值,得到K-1个共振频率差值数据;
步骤6-2:将K-1个共振频率差值数据进行升序排序,然后采用中值滤波技术取出中间的J个差值数据;
步骤6-3:采用均值滤波技术对J个差值数据求平均值,得到共振频率平均差值Δf。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410059106.9A CN103821499B (zh) | 2014-02-20 | 2014-02-20 | 用于油井动液面深度检测的声音信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410059106.9A CN103821499B (zh) | 2014-02-20 | 2014-02-20 | 用于油井动液面深度检测的声音信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103821499A CN103821499A (zh) | 2014-05-28 |
CN103821499B true CN103821499B (zh) | 2017-03-29 |
Family
ID=50756761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410059106.9A Expired - Fee Related CN103821499B (zh) | 2014-02-20 | 2014-02-20 | 用于油井动液面深度检测的声音信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103821499B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105735971B (zh) * | 2014-12-11 | 2018-11-23 | 安徽惠洲地质安全研究院股份有限公司 | 一种基于弹性波的钻孔深度检测系统及其检测方法 |
CN105822289B (zh) * | 2016-03-25 | 2019-03-19 | 重庆科技学院 | 用于油井动液面检测的频率估算方法 |
CN107013206B (zh) * | 2017-05-11 | 2020-05-22 | 重庆科技学院 | 一种油井动液面测量信号处理方法 |
CN108416282B (zh) * | 2018-02-28 | 2021-06-04 | 西安石油大学 | 一种基于油管接箍的井下动液面回波信号声速提取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1250133A (zh) * | 1999-09-30 | 2000-04-12 | 石油大学(华东) | 用超声脉冲反射法测量射孔孔眼深度的方法 |
CN101050702A (zh) * | 2007-04-20 | 2007-10-10 | 中国科学院声学研究所 | 一种过流式共振声谱多相流动态检测的测量装置及方法 |
CN102146791A (zh) * | 2011-05-24 | 2011-08-10 | 中国科学院声学研究所 | 一种油井动液面测量方法及其测量装置 |
CN202374230U (zh) * | 2011-12-28 | 2012-08-08 | 西安科技大学 | 一种油气井液面测试用滤波电路模块 |
CN203008893U (zh) * | 2013-01-17 | 2013-06-19 | 东北石油大学 | 自动控制动液面检测仪 |
CN103321633A (zh) * | 2013-06-25 | 2013-09-25 | 重庆科技学院 | 油井动液面深度检测方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060018190A1 (en) * | 2004-07-21 | 2006-01-26 | Schlumberger Technology Corporation | Tool for Measuring Perforation Tunnel Depth |
-
2014
- 2014-02-20 CN CN201410059106.9A patent/CN103821499B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1250133A (zh) * | 1999-09-30 | 2000-04-12 | 石油大学(华东) | 用超声脉冲反射法测量射孔孔眼深度的方法 |
CN101050702A (zh) * | 2007-04-20 | 2007-10-10 | 中国科学院声学研究所 | 一种过流式共振声谱多相流动态检测的测量装置及方法 |
CN102146791A (zh) * | 2011-05-24 | 2011-08-10 | 中国科学院声学研究所 | 一种油井动液面测量方法及其测量装置 |
CN202374230U (zh) * | 2011-12-28 | 2012-08-08 | 西安科技大学 | 一种油气井液面测试用滤波电路模块 |
CN203008893U (zh) * | 2013-01-17 | 2013-06-19 | 东北石油大学 | 自动控制动液面检测仪 |
CN103321633A (zh) * | 2013-06-25 | 2013-09-25 | 重庆科技学院 | 油井动液面深度检测方法及装置 |
Non-Patent Citations (2)
Title |
---|
窗函数用于频谱分析的研究;周新星等;《微计算机信息》;20091225;第25卷(第36期);184-185,173 * |
频谱分析中窗函数的研究;李杭生等;《微计算机信息》;20080405;第24卷(第10期);272-273 * |
Also Published As
Publication number | Publication date |
---|---|
CN103821499A (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10472957B2 (en) | Entropy based multiphase flow detection | |
CN103821499B (zh) | 用于油井动液面深度检测的声音信号处理方法 | |
CN102116868B (zh) | 一种地震波分解方法 | |
Goryunova et al. | Application of signal analysis for diagnostics | |
CN102692650B (zh) | 一种具有假频压制功能的井筒波分离方法 | |
CN102313904B (zh) | 利用地震信息自动划分地层旋回的方法 | |
CN103321633B (zh) | 油井动液面深度检测方法及装置 | |
CN101551465A (zh) | 一种自适应识别和消除地震勘探单频干扰的方法 | |
CN104597502A (zh) | 一种新的石油地震勘探数据去噪方法 | |
CN108252708B (zh) | 一种油井动液面识别方法 | |
RU2014136576A (ru) | Система и способ вертикального сейсмического профилирования с представлением разведочных данных в виде комбинации параметризованных компрессионного, сдвигового и дисперсивного волновых полей | |
CN107143323A (zh) | 基于welch多段平均功率谱法的油井动液面检测方法 | |
CN108416282B (zh) | 一种基于油管接箍的井下动液面回波信号声速提取方法 | |
CN103674235A (zh) | 基于短时傅里叶变换的单一频率报警声音特征检测方法 | |
CN109164492B (zh) | 一种提取套管井地层声波速度的方法 | |
Wang et al. | Research on comprehensive and effective acoustic signal processing methods for caculating downhole liquid level depth | |
He et al. | Downsampling-based synchrosqueezing transform and its applications on large-scale vibration data | |
US8306762B2 (en) | Systems and methods for analysis of downhole data | |
CN104778342B (zh) | 一种基于小波奇异熵的心音特征提取方法 | |
Xu et al. | A novel method for extracting maximum kurtosis component and its applications in rolling bearing fault diagnosis | |
CN101907726B (zh) | 一种自动识别和消除地震勘探工业电干扰的方法 | |
CN107036706B (zh) | 一种套管振动井口监听检测设备 | |
Dang et al. | Partly duffing oscillator stochastic resonance method and its application on mechanical fault diagnosis | |
CN114859404A (zh) | 超采样地震波形匹配方法及装置 | |
CN116955941B (zh) | 一种随钻测量连续波信号去噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170329 Termination date: 20180220 |