CN107389787B - 一种基于频域的气侵检测信号处理方法 - Google Patents
一种基于频域的气侵检测信号处理方法 Download PDFInfo
- Publication number
- CN107389787B CN107389787B CN201710630541.6A CN201710630541A CN107389787B CN 107389787 B CN107389787 B CN 107389787B CN 201710630541 A CN201710630541 A CN 201710630541A CN 107389787 B CN107389787 B CN 107389787B
- Authority
- CN
- China
- Prior art keywords
- signal
- flow velocity
- check signal
- incursion check
- gas incursion
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 15
- 238000005553 drilling Methods 0.000 claims abstract description 49
- 239000012530 fluid Substances 0.000 claims abstract description 40
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 238000002474 experimental method Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 abstract description 15
- 238000000605 extraction Methods 0.000 abstract description 6
- 239000007789 gas Substances 0.000 description 47
- 238000000034 method Methods 0.000 description 21
- 238000001514 detection method Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 238000002604 ultrasonography Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 239000002245 particle Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 230000011514 reflex Effects 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 239000002002 slurry Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
- G01N29/028—Analysing fluids by measuring mechanical or acoustic impedance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/018—Impedance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02433—Gases in liquids, e.g. bubbles, foams
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于频域的气侵检测信号处理方法，先采用高采样率对气侵检测信号进行A/D转换采集，然后进行累加平均抽取、低通滤波，其中，低通滤波用滤波器依据发射信号的频率范围估算来进行参数设置，这样提高了信号的信杂比，对噪声也有很好地滤除，为后面的流速估算、含气量估算提供了良好的条件；最后，通过功率谱得到频偏来计算钻井液的流速，再通过总功率值与钻井液流速，计算出钻井液的含气量。与现有技术的气侵检测信号处理方法相比，处理方法简单，并能够准确地钻井液的流速和含气量。
Description
技术领域
本发明属于气侵检测技术领域，更为具体地讲，涉及一种基于频域的气侵检测信号处理方法。
背景技术
1、气侵检测方法
气侵检测方法按照传感器安装位置可分为井口检测、海水段检测和随钻检测。井口检测主要包含泥浆池液面增量检测法和泥浆进出口流量差法；海水段检测有压差检测法和声波检测法；随钻检测是将传感器安装于钻柱，通过检测钻井液性能变化判断是否产生气侵。
超声波多普勒检测方法是基于多普勒效应获取信号进行处理的声波检测方法，其基本原理为：当管道内随钻井液运动的反射物与超声波传感器之间存在相对运动时，接收到的反射信号会产生多普勒频偏，当传感器与钻井液中的反射物相互靠近时，接收到的反射信号会产生正频偏；当传感器与钻井液中的反射物相互远离时，接收到的反射信号会产生负频偏。那么反射信号的频偏值就反映了钻井液的流速信息。
此外，接收到的反射信号的强度取决于反射面的大小和反射体的个数。当超声波在钻井液中传播时，遇到了钻井液和气泡交界面时，超声波就会反射回来，被接收传感器接收。那么反射信号的幅值反映了钻井液中的含气量信息。
2、气侵检测信号的获取及其特点
气侵检测信号是由超声波接收传感器接收到的回波信号经过硬件电路处理后输出的待处理信号，气侵检测信号产生框图如图1所示。
首先激励信号产生电路产生激励信号驱动超声波发射传感器向管道中发射超声波信号，超声波信号在钻井液中传播时遇到反射体将信号反射回来，由超声波接收传感器接收信号，但是接收到的回波信号不仅包含了钻井液中气泡的反射信号，还有超声波沿管壁传播接收到的信号和各种干扰信号。
由于管道中的气泡的大小和分布都是随机的，所以接收到的反射信号幅度和相位也具有随机性，反射信号的频偏受到管道流体流速的影响。
若发射信号为u0＝U0sinω0t，那么超声波接收传感器接收到的回波信号表示如下：
由单个气泡反射回的信号其幅值Ui和相位的大小具有随机性，产生的频偏Δωi受到气泡位置和运动速度的影响，超声波接收传感器接收到的气泡的反射信号是每个气泡反射的信号的叠加，叠加得到的频率和幅值分别包含了流速和气泡含量信息。
在实际应用中，超声波发射、接收传感器安装位置固定，沿管壁传播接收到的信号也基本不变，噪声信号N(t)包含了各种干扰信号。
硬件电路对回波信号ur进行检波放大后，得到信号uout表示如下
uout信号是由多个不同频率Δωi、不同幅值Ui和相位的信号的叠加，以及包含了环境中的噪声信号N(t)，其波形如图2所示。uout信号即需要进行处理的气侵检测信号即待处理信号
发明内容
本发明的目的在于克服现有技术的不足，提出一种基于频域的气侵检测信号处理方法，以准确地获得钻井液的流速和含气量。
为实现上述发明目的，本发明基于频域的气侵检测信号处理方法，其特征在于，包括以下步骤：
(1)、气侵检测信号的预处理
使用高采样率的数据采集卡对气侵检测信号做A/D转换采集，得到数字信号进行累加平均抽取，以降低数据冗余，减少后面数据处理的运算量，并且提高信号的信噪比；
然后进行低通滤波，以滤除气侵检测信号中的各种噪声，其中，低通滤波时，根据反射信号的频率范围估算，来设置低通滤波器的参数；
(2)、流速估算
对预处理后的气侵检测信号做功率谱估计，然后根据功率谱计算出频偏值，依据流速与频偏值的对应关系，得到钻井液的流速；
其中，所述流速与频偏值的对应关系，采用实测流速值与气侵检测信号的频偏值进行分析，得到它们之间的对应关系。
(3)、含气量估算
对步骤(2)得到的功率谱的谱值进行积分，得到气侵检测信号的总功率值，然后根据步骤(2)计算出的钻井液流速，依据总功率值与流速、含气量之间的关系，得到钻井液的含气量；
其中，所述总功率值与流速、含气量之间的关系，采用实测不同流速和含气量下获取的气侵检测信号的总功率值进行分析得到。
本发明基于频域的气侵检测信号处理方法，先采用高采样率对气侵检测信号进行A/D转换采集，然后进行累加平均抽取、低通滤波，其中，低通滤波用滤波器依据发射信号的频率范围估算来进行参数设置，这样提高了信号的信杂比，对噪声也有很好地滤除，为后面的流速估算、含气量估算提供了良好的条件；最后，通过功率谱得到频偏来计算钻井液的流速，再通过总功率值与钻井液流速，计算出钻井液的含气量。与现有技术的气侵检测信号处理方法相比，处理方法简单，并能够准确地钻井液的流速和含气量。
附图说明
图1是气侵检测信号的产生原理框图；
图2是气侵检测信号一具体实例波形图；
图3是本发明基于频域的气侵检测信号处理方法一具体实施方式流程图；
图4是气泡上升示意图；
图5是频偏值随时间变化图；
图6是采集到信号抽取前后比较图；
图7是气侵检测信号及其功率谱一具体实例图；
图8是实验装置示意图；
图9是流速与频偏均方根关系图；
图10是不同含气量实验的功率值结果图；
图11是流速和含气量与功率的拟合结果图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。
一、信号处理过程
在本实施例中，如图3所示，本发明对气侵检测信号的处理过程分为三个部分，包括气侵检测信号的预处理、流速估算和含气量估算。
步骤S1：气侵检测信号的预处理
使用高采样率的数据采集卡对气侵检测信号做A/D转换采集，得到数字信号进行累加平均抽取，高采样率采集信号可以减少A/D转换带来的误差，以降低数据冗余，减少后面数据处理的运算量，并且提高信号的信噪比。
然后进行低通滤波，以滤除气侵检测信号中的各种噪声，其中，低通滤波时，根据反射信号的频率范围估算，来设置低通滤波器的参数。
下面对气侵检测信号的预处理作详细说明。
1.1)、信号的累加平均抽取
由于气侵检测信号包含噪声，信号本身存在很多毛刺，若是将信号直接用低采样率做A/D转换采集，那么采集得到的信号与原信号相比就会存在很大的误差，造成流速和含气量的计算不准确，所以本发明使用高采样率的数据采集卡对气侵检测信号做A/D转换采集，得到数字信号即高采样率采集信号，然后对高采样率采集做累加平均抽取，在减小误差的同时减少数据冗余，降低数据处理运算量，并能够提高信噪比。
对反射信号的频率范围可做如下估计：
若超声波的发射频率为f0，在钻井液中的波速为c。设单个小气泡在环空管中以速度v匀速上升，如图4所示，小气泡上升到A点时，超声波接收传感器能够接收到气泡从此处反射回的信号，然后逐渐上升，小气泡经过B点后，由于与超声波接收传感器距离较远，反射信号将在泥浆中衰减不会被超声波接收传感器接收。由于超声波传感器的尺寸与环空管的直径相比要小得多，所以可将紧挨着的两个传感器即超声波发射、接收传感器的安装位置看作一个点C。
如图4所示，A点和B点应该有个对称性，即AC＝BC，那么应将气泡上升过程分两段来看：
在AO段，气泡和超声波波束之间的夹角θ在逐渐增大，θ随时间的变化可表示为
气泡接收并反射的超声波频率为f1，接收传感器接收到的超声波频率为f2，那么AO段的频偏为：
在OB段，θ随时间的变化可表示为：
OB段的频偏为：
公式(4)和公式(6)的Δfd表达式一样，但是θt的取值范围不一样，可得出单个气泡上升过程中，回波信号的频偏先正后负，并且Δfd的取值由c、v、θt、f0确定。
例如假设初始夹角θ0＝π/6，气泡到管壁的距离为r＝0.1m，上升速度为2m/s。超声波发射频率为f0＝64kHz。超声波在钻井液中的传播速度为1850m/s。通过MATLAB仿真得到频偏曲线如图5所示。
从图5可看出，气泡在通过超声波可检测范围时，多普勒频偏绝对值先增大再减小。且变化速度越来越快。图5中的频偏值为正时，是气泡在AO段的上升过程时得到的频偏；频偏值为负时，是气泡在OB段上升过程得到的频偏。当θ0＝π/6时，经计算频偏绝对值最大为120Hz；当θ0趋于0时，计算得到频偏绝对值最大约为138Hz。实验中获得的信号是通过管道中很多反射粒子所反射信号的叠加，这些信号所产生的频偏范围与单个粒子引起的频偏范围一致。
累加平均算法利用噪声信号为随机过程其均值为零这一统计特性，对同一个样本点进行多次采样后计算其均值，使噪声大大降低。实验中对信号做A/D转换采集的采样率远大于有效信号的频偏范围，所以对信号做抽取不会造成有用信号在频谱上的混叠，降低了数据量并且信号的幅度几乎和抽取前没发生变化。若信号频率为f，采样频率为fs，当采样频率远大于信号频率时，可以将相邻m个采样点的采样值看做同一采样点的采样值，对每m个点进行累加平均，则得到累加平均后的采样信号为：
累加平均后，等效于信号采样频率降低为fs/m，并且经过m次累加信噪比为：
累加平均提高了信噪比，降低了噪声。取一组采集到的实验数据进行抽取前后的对比如图6所示。其中采样率fs＝80kHz，抽取倍数m＝40。
从图6的两组波形可观察出信噪比得到了明显的改善，采样率降低了m倍，去掉了冗余信号，减少了数据处理时间。
1.2)、信号的低通滤波
对信号做累加平均抽取后提高了信噪比，但仍然不能将噪声完全去除。所以使用一个低通滤波器对信号进行滤波是十分必要的。
选择FIR低通滤波器，那么需要设置的参数就有滤波器截止频率、滤波器窗函数以及滤波器阶数。滤波器的截止频率设置可根据有用信号的频率范围估算结果来确定。考虑一定余量，将截止频率设置为稍微大于对信号估计的最大频偏值。例如若是计算出最大频偏值为138Hz，由于管道中有的反射粒子运动速度较快，产生的频偏值也会大一些，那么可将截止频率设置为200Hz。由于得到的信号是随机信号，并且频谱比较复杂，选择旁瓣峰值较小并且衰减速度较快的汉宁窗。滤波器阶数过高会增大计算量，并且产生较大延时，对于N阶的FIR滤波器的延时是N/(2fs)，可将阶数设置为50阶。
预处理气侵检测信号用x(n)表示，n表示采样点数。
步骤S2：流速估算
在本发明中，对预处理后的气侵检测信号做功率谱估计，然后根据功率谱计算出频偏值，依据流速与频偏值的对应关系，得到钻井液的流速；其中，所述流速与频偏值的对应关系，采用实测流速值与气侵检测信号的频偏值进行分析，得到它们之间的对应关系。
下面对钻井液的流速估计作详细说明。
2.1)、钻井液流速估算原理
根据多普勒效应可知，信号的频偏值能够反映管道中的钻井液流速值，对于单个反射体而言，由公式(4)和公式(6)可得管道内反射体流速和频偏的关系如下，其中超声波发射频率f0远大于频偏值Δfd：
管道内的实际情况是在任意时刻都有多个角度不同、大小不同的气泡反射信号，所以有用信号的频谱不是表现出一个单一的频点，而是表现了多个气泡的反射信号带来的不同频偏值。所以流速不能根据公式(9)直接计算，需要根据信号计算出一个能够反映整体流速的平均频偏进行估计。
其中，为反射颗粒或气泡的多普勒频偏的平均值；Ni(Δfdi)表示产生多普勒频偏值Δfdi的粒子数；∑iNi(Δfdi)表示所有能够反射超声波的粒子个数；Δfdi表示任一颗粒产生的多普勒频偏值。
2.2)、信号的功率谱估计
信号的功率谱密度定义为自相关函数的傅里叶变换。经典谱估计法中的周期图法是最早提出的算法，可用FFT实现。气侵检测信号x(n)的功率谱计算公式为：
为了能够直接度量信号功率大小，在本实施例中，将周期图法做了一些修改，计算过程为将x(n)的傅里叶变换的模乘以2/N得到单边幅度谱，然后再取平方，得到的功率谱即表示每赫兹的瓦特数。
其中，N表示将预处理气侵检测信号x(n)的长度，公式(12)的气侵检测信号x(n)的功率谱用频率f表示为S(f)，以便于频偏的计算。
2.3)、计算信号频偏值
频偏估计法对平均频率的定义如下
在检测环境中，当管道内的气泡产生的频偏值越大时，除了少部分靠近管壁的反射体外，通常这些反射体距离超声波传感器的位置都会较远，那么经过在泥浆中的传播衰减后传感器接收到的信号能量就较小，所以采用下式求取信号的频偏均方根值，相当于对分子上的S(f)乘上了一个系数w＝f，再开方得到frms。
根据公式(12)得到功率谱后计算信号的频偏均方根值，积分下限是0Hz，积分上限根据有用信号的频率估计范围来确定，将有效频带的功率谱积分。得到的功率谱是离散的，需要将公式(14)离散化，积分变为求和。假设计算出信号的最大频偏在138Hz，可将积分上下限设为0～150Hz，那么频偏均方根的计算可表示为：
上式中Δf为功率谱的频率分辨力
fs为信号抽取后的实际采样率，N为数据长度。
2.4)、估算流速
采集不同流速条件下的实验数据，对这些数据进行分析得到它们对应的频偏均方根值，寻找出频偏均方根值与实验中流量计测得的实际流速值之间的关系。通过对大量的实验数据分析，得到实验中流速v与频偏均方根frms为三次拟合关系，表示如下：
frms＝av3+bv2+cv+d (17)
具体的拟合系数a、b、c、d由实验确定，不同密度的钻井液(泥浆)实验对应不同的拟合系数。那么在实际测量中，在某一钻井液(泥浆)密度下的流速可根据检测信号的频偏均方根计算得到。
步骤S3：含气量估算
对步骤S2得到的功率谱的谱值进行积分，得到气侵检测信号的总功率值，然后根据步骤S2计算出的钻井液流速，依据总功率值与流速、含气量之间的关系，得到钻井液的含气量；
其中，所述总功率值与流速、含气量之间的关系，采用实测不同流速和含气量下获取的气侵检测信号的总功率值进行分析得到。
3.1)、钻井液含气量估算原理
超声波在碰到声阻抗有差别的界面就会发生反射作用，两种介质的声阻抗差别越大，反射信号能量就越强。钻井液的声阻抗远大于气泡的声阻抗，当钻井液中含气时，超声波在钻井液中传输遇到气泡时基本上会发生全反射。
对于钻井液中含气量的分析，由于管道内含气量可以通过气泡的数量和大小反映，即游离气含量可以表示成气泡数量与气泡平均大小的乘积。而在检测过程中，超声波通过气泡表面反射信号。钻井液中气体含量越多，那么全反射面就越大，意味着反射的能量越多，超声波传感器接收到的回波信号也就越强。
此外，考虑超声波在泥浆中的衰减问题，在含气量越多时，管道内的气泡增多，并且这些气泡随机地分布在管道内各个区域，那么在管壁附近区域内的气泡也会相应增多，超声波在钻井液中传播时遇到管壁附近的气泡就会被反射回来，缩短了它的传播距离，衰减也会相应变小。
从上面的分析可知，含气量越多，接收到的回波信号幅度也会越大，那么检波过后的有用信号幅度也会越大，计算出的信号功率值也会越大。同时在实验中发现，流速加快后，也使得有用信号幅度增大，应该是流速加快影响了气泡在管道中的分布，导致气泡更靠近管壁或者增大了反射面积，所以信号的功率与含气量和流速都有关系。
3.2)、计算信号总功率
在钻井液流速估算部分求取了信号的功率谱S(f)，那么对功率谱沿频率轴积分，便得到了信号的总功率，表示如下：
与频偏均方根的计算相同，将上式离散化，积分变为累加求和，求和的频带范围与计算频偏均方根使用的频带范围一致。若取0～150Hz，那么信号总功率值的计算如下式
公式(19)中的Δf与公式(15)中的Δf一样。
假设抽取后的采样频率fs＝2000Hz，每组数据采集点数N＝4000，那么Δf＝fs/N＝2000/4000＝0.5Hz，0～150Hz的范围在功率谱上对应301个离散点，对如图7中的功率谱在0～150Hz范围进行累加求和即得到总功率值。
3.2)、估算含气量
采集不同流速和不同含气量条件下的实验数据，对这些数据进行分析得到每组数据对应的总功率值，寻找出总功率值与计算出的流速和实际的含气量之间的关系。通过对大量的实验数据分析，对流速v、含气量g和信号功率p做曲面拟合，可得到它们之间有如下关系
p＝k00+k10g+k01v+k20g2+k11gv+k02v2 (20)
具体的拟合系数k00、k10、k01、k20、k11、k02由实验确定，不同密度的泥浆实验对应不同的系数。那么在实际的含气量估算中，在某一泥浆密度下的含气量可根据计算出的流速和信号总功率得到。
二、实验结果
1.实验装置
实验装置由空压机、泥浆储存和搅拌桶、活塞装置、高压气瓶、主体循环装置、循环泵及其控制设备6部分组成。首先空压机提供压力将搅拌桶中的泥浆注入主体循环装置。泥浆的循环流动在封闭的主体循环装置中进行，配备循环泵保持泥浆在管中稳定的流动状态，并通过改变泵的转速控制管道中的泥浆速度。使用高压氮气瓶连接活塞装置可向管道内加压，天然气罐连接活塞装置可向管道内注入天然气模拟气侵。装置示意图如图8所示。
2.流速分析结果
实验中通过逐级改变变频器频率来改变流体流速，采集不同流速条件下的实验数据。实验中使用的油基钻井液密度为1.39g/mL。表1是密度为1.39g/mL的泥浆实验中变频器频率与流速的变化关系表。
表1
计算每组实验数据的频偏均方根值如表2所示，表2是密度为1.39g/mL的泥浆实验频偏均方根结果
| 变频器频率 | 15Hz | 18Hz | 21Hz | 24Hz | 27Hz | 30Hz | 33Hz |
| 频偏均方根(Hz) | 3.327 | 3.974 | 5.980 | 11.208 | 17.128 | 21.037 | 27.200 |
| 变频器频率 | 36Hz | 39Hz | 42Hz | 45Hz | 48Hz | 50Hz | |
| 频偏均方根(Hz) | 34.615 | 43.811 | 50.419 | 53.945 | 57.518 | 58.753 |
表2
将表2中的频偏均方根值与表1中对应的流速结合分析，拟合流速和频偏均方根之间的关系，得到的结果如图9所示，拟合曲线表示如下
f＝-105.26v3+219.27v2-70.166v+9.5965 (21)
其中f为频偏均方根，单位Hz、v为泥浆实测流速值，单位m/s。
不同实验条件下得到的拟合系数会有差别，但泥浆流速和频偏之间是三次拟合关系。在泥浆流速计算中，通过计算采集的随机信号的频偏均方根，将得到的频偏均方根值代入与流速的函数关系反解即得到流速值。
3.含气量分析结果
实验使用的泥浆密度为1.39g/mL，含气量分别为0.2％、1.25％、3％，变频器频率从16Hz开始，逐次增加4Hz进行实验，计算采集数据的功率值结果如图10所示。
从图11中可以看出，在不同含气量下的功率值之间的差别较大，在相同流速条件下，即变频器频率一定时，含气量越多，获取到的信号的功率值越大；同时，含气量一定时，功率值也在流速增加时呈逐渐增大趋势。将图10中的功率值和变频器对应的流速值、含气量做曲面拟合，拟合结果如图11所示。
图11中的拟合结果表达式为
p＝0.131-0.0666g-0.2549v+0.02267g2+0.2622gv+0.2726v2 (22)
式中p表示信号的功率，单位W；g为泥浆含气量(％)；v为泥浆流速，单位m/s。
在不同的实验条件下得到的系数会有差别。在气侵检测中，得出钻井液流速后，根据当前流速值下信号的功率值的大小即可判断出含气量的多少。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
Claims (5)
1.一种基于频域的气侵检测信号处理方法，其特征在于，包括以下步骤：
(1)、气侵检测信号的预处理
使用高采样率的数据采集卡对气侵检测信号做A/D转换采集，得到数字信号进行累加平均抽取，以降低数据冗余，减少后面数据处理的运算量，并且提高信号的信噪比；
然后进行低通滤波，以滤除气侵检测信号中的各种噪声，其中，低通滤波时，根据反射信号的频率范围估算，来设置低通滤波器的参数；
(2)、流速估算
对预处理后的气侵检测信号做功率谱估计，然后根据功率谱计算出频偏值，依据流速与频偏值的对应关系，得到钻井液的流速；
其中，所述流速与频偏值的对应关系，采用实测流速值与气侵检测信号的频偏值进行分析，得到它们之间的对应关系；
(3)、含气量估算
对步骤(2)得到的功率谱的谱值进行积分，得到气侵检测信号的总功率值，然后根据步骤(2)计算出的钻井液流速，依据总功率值与流速、含气量之间的关系，得到钻井液的含气量；
其中，所述总功率值与流速、含气量之间的关系，采用实测不同流速和含气量下获取的气侵检测信号的总功率值进行分析得到。
2.根据权利要求1所述的气侵检测信号处理方法，其特征在于，所述的功率谱估计表示每赫兹的瓦特数：
其中，x(n)表示预处理气侵检测信号，n表示采样点数，N表示将预处理气侵检测信号x(n)的长度。
3.根据权利要求2所述的气侵检测信号处理方法，其特征在于，所述的频偏值为频偏均方根值frms：
其中，S(f)为用频率f表示的功率谱。
4.根据权利要求3所述的气侵检测信号处理方法，其特征在于，所述流速与频偏值的对应关系为流速v与频偏均方根frms的三次拟合关系：
frms＝av3+bv2+cv+d
具体的拟合系数a、b、c、d由实验确定，不同密度的钻井液实验对应不同的拟合系数。
5.根据权利要求4所述的气侵检测信号处理方法，其特征在于，所述总功率值p与流速v、含气量g之间的关系为：
p＝k00+k10g+k01v+k20g2+k11gv+k02v2
具体的拟合系数k00、k10、k01、k20、k11、k02由实验确定，不同密度的钻井液实验对应不同的系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710630541.6A CN107389787B (zh) | 2017-07-28 | 2017-07-28 | 一种基于频域的气侵检测信号处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710630541.6A CN107389787B (zh) | 2017-07-28 | 2017-07-28 | 一种基于频域的气侵检测信号处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107389787A CN107389787A (zh) | 2017-11-24 |
| CN107389787B true CN107389787B (zh) | 2019-09-24 |
Family
ID=60342132
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710630541.6A Expired - Fee Related CN107389787B (zh) | 2017-07-28 | 2017-07-28 | 一种基于频域的气侵检测信号处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107389787B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110185433B (zh) * | 2019-05-16 | 2022-09-16 | 中国海洋石油集团有限公司 | 一种基于谱特征分析方法的隔水管气侵监测装置及方法 |
| CN112067874A (zh) * | 2019-06-11 | 2020-12-11 | 广州汽车集团股份有限公司 | 一种适用于车辆传感器信号的滤波方法、装置和存储介质 |
| CN113153263B (zh) * | 2021-04-26 | 2024-05-10 | 中国石油天然气集团有限公司 | 一种高噪声背景下井下溢流多普勒气侵监测装置和方法 |
| CN115126963B (zh) * | 2022-06-21 | 2022-12-30 | 安徽省特种设备检测院 | 一种内检测器的检测信号处理方法及系统 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1907230A (zh) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | 频谱多普勒血流速度的自动检测方法 |
| CN103032064A (zh) * | 2013-01-11 | 2013-04-10 | 西南石油大学 | 一种钻井过程中气侵位置检测的方法及装置 |
| CN103291284A (zh) * | 2013-05-14 | 2013-09-11 | 中国海洋石油总公司 | 基于随钻环空压力测量的井筒气侵早期监测方法及装置 |
| CN103811016A (zh) * | 2014-01-16 | 2014-05-21 | 浙江工业大学 | 一种基于周期图法的冲床噪声功率谱估计改进方法 |
| CN104677437A (zh) * | 2015-02-12 | 2015-06-03 | 延安大学 | 一种超声波液相流量精确化测量方法 |
| CN105041302A (zh) * | 2015-08-26 | 2015-11-11 | 中国海洋石油总公司 | 一种深海钻井中隔水管内泥浆循环模拟实验方法和装置 |
| CN105545285A (zh) * | 2015-10-22 | 2016-05-04 | 中国石油大学(华东) | 基于隔水管气液两相流识别的深水钻井气侵监测方法 |
| CN105890685A (zh) * | 2016-06-28 | 2016-08-24 | 电子科技大学 | 一种基于累积相位差的超声波流量测量装置 |
| CN105933257A (zh) * | 2016-04-28 | 2016-09-07 | 西安电子科技大学 | 一种频谱共享方式下时频重叠信号的信噪比估计方法 |
| CN205644584U (zh) * | 2016-05-09 | 2016-10-12 | 中国海洋石油总公司 | 一种深海钻井气侵检测信号处理电路 |
| CN106194162A (zh) * | 2016-08-30 | 2016-12-07 | 中国石油大学(华东) | 基于环空压差随钻测量的气侵监测装置及监测方法 |
| CN106640051A (zh) * | 2016-12-19 | 2017-05-10 | 中国石油大学(北京) | 隔水管内气侵早期声波监测方法及系统 |
-
2017
- 2017-07-28 CN CN201710630541.6A patent/CN107389787B/zh not_active Expired - Fee Related
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1907230A (zh) * | 2005-08-02 | 2007-02-07 | 深圳迈瑞生物医疗电子股份有限公司 | 频谱多普勒血流速度的自动检测方法 |
| CN103032064A (zh) * | 2013-01-11 | 2013-04-10 | 西南石油大学 | 一种钻井过程中气侵位置检测的方法及装置 |
| CN103291284A (zh) * | 2013-05-14 | 2013-09-11 | 中国海洋石油总公司 | 基于随钻环空压力测量的井筒气侵早期监测方法及装置 |
| CN103811016A (zh) * | 2014-01-16 | 2014-05-21 | 浙江工业大学 | 一种基于周期图法的冲床噪声功率谱估计改进方法 |
| CN104677437A (zh) * | 2015-02-12 | 2015-06-03 | 延安大学 | 一种超声波液相流量精确化测量方法 |
| CN105041302A (zh) * | 2015-08-26 | 2015-11-11 | 中国海洋石油总公司 | 一种深海钻井中隔水管内泥浆循环模拟实验方法和装置 |
| CN105545285A (zh) * | 2015-10-22 | 2016-05-04 | 中国石油大学(华东) | 基于隔水管气液两相流识别的深水钻井气侵监测方法 |
| CN105933257A (zh) * | 2016-04-28 | 2016-09-07 | 西安电子科技大学 | 一种频谱共享方式下时频重叠信号的信噪比估计方法 |
| CN205644584U (zh) * | 2016-05-09 | 2016-10-12 | 中国海洋石油总公司 | 一种深海钻井气侵检测信号处理电路 |
| CN105890685A (zh) * | 2016-06-28 | 2016-08-24 | 电子科技大学 | 一种基于累积相位差的超声波流量测量装置 |
| CN106194162A (zh) * | 2016-08-30 | 2016-12-07 | 中国石油大学(华东) | 基于环空压差随钻测量的气侵监测装置及监测方法 |
| CN106640051A (zh) * | 2016-12-19 | 2017-05-10 | 中国石油大学(北京) | 隔水管内气侵早期声波监测方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| "基于功率谱密度的静电悬浮加速度计分辨率估计方法";刘爽 等;《仪表技术与传感器》;20151231(第10期);第30-34页 * |
| "基于多普勒测量技术的深水隔水管气侵早期监测研究";隋秀香 等;《石油钻探技术》;20140930;第42卷(第5期);第90-94页 * |
| "超声波流速测量的多普勒频移提取方法研究";曹新亮 等;《应用声学》;20170331;第36卷(第2期);第149-154页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107389787A (zh) | 2017-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107389787B (zh) | 一种基于频域的气侵检测信号处理方法 | |
| Blanckaert et al. | Means of noise reduction in acoustic turbulence measurements | |
| CN107389786B (zh) | 一种基于时域的气侵检测信号处理方法 | |
| CN106871980B (zh) | 基于相邻峰值最大差值的气体超声流量计间歇式激励和信号处理方法和系统 | |
| CN103630706B (zh) | 一种在声学多普勒流速剖面仪中获取径向水流速度的方法 | |
| JP2003344131A (ja) | 超音波流速分布及び流量計 | |
| CN108490069A (zh) | 两相流分布参数双模式超声可视化测量装置 | |
| CN101592288B (zh) | 一种管道泄漏辨识方法 | |
| JPH0778437B2 (ja) | 媒体の音響エネルギ伝送特性測定装置 | |
| CN108332063A (zh) | 一种基于互相关的管道泄漏定位方法 | |
| CN104820218B (zh) | 一种基于频域自相关的浅海海底单参数反演方法 | |
| CN103644460A (zh) | 一种输气管线泄漏音波信号的滤波优选方法 | |
| CN105116445A (zh) | 一种水陆检波器地震数据合并处理的方法及装置 | |
| CN106289121A (zh) | 一种变径管等效管长的计算方法 | |
| CN108416282B (zh) | 一种基于油管接箍的井下动液面回波信号声速提取方法 | |
| CN107436451A (zh) | 一种自动计算地震数据光缆耦合噪声强弱程度的振幅谱方法 | |
| CN108252708A (zh) | 一种油井动液面识别方法 | |
| CN105067048B (zh) | 一种抗振动干扰的数字化涡街流量计的测量方法 | |
| CN104181586A (zh) | 一种水陆检波器数据海底反射系数反演的方法 | |
| Wang et al. | Research on comprehensive and effective acoustic signal processing methods for caculating downhole liquid level depth | |
| Morriss et al. | Measurement of velocity profiles in upwards oil/water flow using ultrasonic Doppler velocimetry | |
| CN110108797A (zh) | 利用声阻抗变化信息的介质分界面超声检测方法 | |
| CN102798889B (zh) | 相控震源一致性判定方法 | |
| CN106195648A (zh) | 一种变径管等效管长的实验测试方法 | |
| CN107741263A (zh) | 基于频域信息融合的液位估计方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for 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 |
Granted publication date: 20190924 |