CN111060168A - 一种流量信号采样方法和装置 - Google Patents
一种流量信号采样方法和装置 Download PDFInfo
- Publication number
- CN111060168A CN111060168A CN201911363057.7A CN201911363057A CN111060168A CN 111060168 A CN111060168 A CN 111060168A CN 201911363057 A CN201911363057 A CN 201911363057A CN 111060168 A CN111060168 A CN 111060168A
- Authority
- CN
- China
- Prior art keywords
- signal
- sampling
- low
- pass filtering
- filtering processing
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/56—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using electric or magnetic effects
- G01F1/58—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using electric or magnetic effects by electromagnetic flowmeters
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种流量信号采样方法和装置,用于油田注水井电磁流量计的流量信号采样。该方法包括:步骤S1:获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;步骤S2:利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;步骤S3:利用递推数字滤波算法对经步骤S2进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;步骤S4:将步骤S3得到的稳定的采用信号作为流量信号输出。该方法由单片机实施,用算法实现,能够代替一部分硬件电路,节省占用空间;可以在井下复杂电磁环境下,减少噪声影响,及时准确地测量出流量信号,并满足测量需求的稳定性和鲁棒性。
Description
技术领域
本发明涉及油田注水井计技术领域,具体涉及一种用于油田注水井电磁流量计的流量信号采样方法和装置。
背景技术
油田注水是油田开发过程中向地层补充能量、提高油田采收率的重要手段之一,其作用是弥补原油采出后所造成的地下亏空,保持或提高油层压力,实现油田高产稳产。在实际生产注水过程中,每口注水井都应该连续均衡地完成配注计划,注水不能过多也不能过少,更不能无故停注。因此,需要通过流量计及时准确地了解到当前时刻的注水流量。目前市场上的电磁流量计绝大多数是把检测到的感生电动势直接转化为对应的流量值。
但是,针对油田注水井独特的井下环境,上述方式存在一些不足:
1、井下注水设备有电磁干扰复杂、充满工频噪声等一系列不利于测量流量的因素,会较大影响流量测量值的准确性;
2、电磁流量计电路板本身也具有一定的机械加工糙点,在平常生产环境可能不会有多大影响,但考虑到井下设备无规律抖动、地质震动、注水流体冲击加压等外部应力,可能会造成电路板上测量电极以及其它电子元件抖动,也会造成测量噪声;
3、井下测量环境狭小,各种传感器和测量设备众多,传统电磁流量计电子元件多,占用空间大。
发明内容
本发明实施例的目的在于提供一种用于油田注水井电磁流量计的流量信号采样方法,用于减少噪声影响,提高准确度,并有助于减少硬件电路。本发明实施例的目的还在于提供相应的装置。
第一方面,提供一种流量信号采样方法,用于油田注水井电磁流量计的流量信号采样,该方法包括:
步骤S1:获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
步骤S2:利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
步骤S3:利用递推数字滤波算法对经步骤S2进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
步骤S4:将步骤S3得到的稳定的采用信号作为流量信号输出。
第二方面,提供一种流量信号采样装置,用于油田注水井电磁流量计的流量信号采样,其特征在于,该装置包括:
输入模块,用于获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
低通滤波处理模块,用于利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
递推数字滤波处理模块,用于利用递推数字滤波算法对经低通滤波处理模块进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
输出模块,用于将递推数字滤波处理模块得到的稳定的采用信号作为流量信号输出。
第三方面,提供一种单片机,包括处理器和存储器,其特征在于,所述存储器中存储有程序,当所述单片机运行时,所述处理器通过执行所述存储器中存储的程序,使所述单片机执行如第一方面所述的流量信号采样方法。
第四方面,提供一种存储介质,存储有程序,当所述程序被单片机执行时,使所述单片机执行如第二方面所述的流量信号采样方法。
从以上技术方案可以看出,本发明实施例具有以下优点:
本发明的用于油田注水井电磁流量计的流量信号采样方法,由单片机实施,用算法实现,能够代替一部分硬件电路,节省硬件装置的占用空间;可以在电磁干扰和各种噪声干扰的井下复杂电磁环境下,减少噪声影响,及时准确地测量出流量信号,并满足测量需求的稳定性和鲁棒性。
附图说明
为了更清楚地说明本发明实施例技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。
图1是本发明一个实施例提供的一种流量信号采样方法的流程示意图;
图2是本发明一个实施例中低通滤波处理的流程示意图;
图3是本发明一个实施例中递推数字滤波处理的流程示意图;
图4是本发明一个实施例提供的一种流量信号采样装置的结构示意图;
图5是本发明一个实施例提供的一种单片机的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”等是用于区别不同的对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面通过具体实施例,分别进行详细的说明。
请参考图1,本发明的一个实施例,提供一种流量信号采样方法,用于油田注水井电磁流量计的流量信号采样。
该方法可包括:
步骤S1,信号输入:获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
步骤S2,低通滤波处理:利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
步骤S3,递推数字滤波处理:利用递推数字滤波算法对经步骤S2进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
步骤S4,信号输出:将步骤S3得到的稳定的采用信号作为流量信号输出。
一些实施例中,该方法可由单片机实施。其中,步骤S1的信号输入,可以是通过单片机从电磁流量计捕获低频的脉冲信号;步骤S2至S3采用的低通滤波算法和递推滤波算法,可以用单片机C语言编程实现;步骤S4的信号输出可以是通过单片机的TX端口和RX端口来输出。
如图2所示,是一些实施例中,低通滤波处理的流程示意图。
步骤S2的低通滤波处理可包括:按照设定的采样频率对脉冲信号采样;采用FIR(Finite Impulse Response,有限长单位冲激响应)滤波函数对采样得到的信号进行移位、加权和卷积运算,实现低通滤波处理,该低通滤波处理操作用于滤除采样信号中的高频噪声,包括电路噪声和工频噪声等,从而在软件层次上即可达到FIR滤波器硬件电路的效果,使得信号采样结果具有良好的鲁棒性。FIR(Finite Impulse Response,有限长单位冲激响应)滤波函数可以是Matlab软件的内建函数。一种可能的实现方式中,具体流程如下。
(1)设定采样频率为300Hz,可定义低频目标脉冲信号x(k)、20Hz电路噪声n1(k)和50Hz工频噪声n2(k)的各个表达公式:x(k)=1,n1(k)=sin(2*π*k*20/300),n2(k)=sin(2*π*k*50/300),其中k=1,2,3,…,1000为采样时刻。
(2)声明一个能容纳1000个char类型数据的采样数组x[k]。
(3)对脉冲信号采样,可将采样信号用数组x[k]表示。
实际的采样信号中包含了电路噪声和工频噪声,按照步骤(1)里的三个公式,对脉冲信号采样相当于分别对1Hz的低频目标脉冲信号、20Hz电路噪声、50Hz工频噪声进行采样和求和,采样得到的信号存进步骤(2)里声明的采样数组x[k],x[k]的表达式可记为为x[k]=1+sin(2*π*k*20/300)+sin(2*π*k*50/300),其中k=1,2,3,…,1000。
(4)利用Matlab软件的内建滤波函数firl()进行低通滤波。
首先,可利用Matlab软件计算低通滤波的卷积系数h(i),i=1,2,…,N。N即为低通滤波阶数,可表示卷积系数h(i)为一个N点的冲激序列。设定低通滤波算法的通带截止频率fp=50Hz,阻带起始频率fst=50Hz,抽样频fs=300Hz,阻带衰减不少于-50dB。将这些初始参数输入Matlab软件自带的内建函数firl(),即可求出低通滤波的卷积系数h(i)。
然后,可利用下述公式求出经过低通滤波处理后的输出脉冲信号y(k):
该公式中,是先将采样信号x(k)进行移位,即x(k-i);然后与卷积系数h(i)求积加权,即h(i)*x(k-i);最后对积求和,进行卷积运算,即算出脉冲信号y(k)。计算得到的脉冲信号y(k)已经滤除了20Hz的电路噪声和50Hz的工频噪声,仅含有低频目标脉冲信号x(k),这就完成了对采样信号的低通滤波处理。
如图3所示,是一些实施例中,递推数字滤波处理的流程示意图。
步骤S3的递推数字滤波处理可包括:对于经步骤S2进行低通滤波处理处理后的信号,利用单片机中断技术实时更新信号数值,通过依次对M个最新的信号数值求平均值,得到稳定的采样信号。其中,M为正整数,例如M可以等于10。该步骤对信号进行递推求平均值,可使输出结果更加线性稳定,具有良好的稳定性。一种可能的实现方式中,具体流程如下:
(1)声明一个能存放十个char类型数据的递推求平均值数组buf[];
(2)声明一个用于求和的char类型变量sum,sum表示数组buf[]里十个数的和;
(3)单片机每发生一次中断,就把当前发来的信号的数值存进数组buf[]里;
(4)发生十次中断后,对数组buf[]进行求平均值avg,其公式为avg=sum/10,将求得的数值输出,如此反复进行,输出稳定的采样信号。
该稳定的采用信号作为流量信号输出,可用于转化为流量值。
以上,本发明实施例公开了一种流量信号采样方法,该方法是一种专门用于油田注水井电磁流量计在井下复杂电磁环境下及时准确测量出流量信号的采样算法。该流量信号采样方法取得了以下技术效果:
(1)该方法的低通滤波处理步骤通过对采样信号进行低通滤波处理,可有效滤除电磁干扰和高频噪声,提高采样结果准确率,保证较好的数据线性相位;
(2)该方法的递推数字滤波处理步骤利用单片机中断技术,可以实时更新采样数据,代替传统的暴力算法,克服单片机存储空间的限制,解放芯片算力;
(3)该方法的递推数字滤波步骤可以去除因电路噪声、水压突变等原因引起的计数漂移,使流量信号的测量值更加线性稳定贴合实际;
(4)该方法仅在软件层次上即可达到低通滤波硬件电路的效果,代替掉相当多的电子元件,大大减小电路板板上面积,使电路整体设计更加符合井下狭小的测量环境,使得流量信号测量过程更加灵活高效,节省生产成本。
请参考图4,本发明的一个实施例,还提供一种流量信号采样装置,用于油田注水井电磁流量计的流量信号采样,该装置包括:
输入模块41,用于获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
低通滤波处理模块42,用于利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
递推数字滤波处理模块43,用于利用递推数字滤波算法对经低通滤波处理模块进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
输出模块44,用于将递推数字滤波处理模块得到的稳定的采用信号作为流量信号输出。
可选的,所述低通滤波处理模块42,具体用于:按照设定的采样频率对脉冲信号采样;采用FIR滤波函数对采样得到的信号进行移位、加权和卷积运算,实现低通滤波处理,该低通滤波处理操作用于滤除采样信号中的电路噪声和工频噪声。
可选的,所述递推数字滤波处理模块43,具体用于:对于经所述低通滤波处理模块进行低通滤波处理处理后的信号,利用单片机中断技术实时更新信号数值,通过依次对M个最新的信号数值求平均值,得到稳定的采样信号。
请参考图5,本发明的一个实施例,还提供一种单片机50,包括处理器51和存储器52,所述存储器52中存储有程序,当所述单片机50运行时,所述处理器51通过执行所述存储器52中存储的程序,使所述单片机50执行如图1实施例所述的流量信号采样方法。
本发明的一个实施例,还提供一种存储介质,存储有程序,当所述程序被单片机执行时,使所述单片机执行如图1实施例所述的流量信号采样方法。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
上述实施例仅用以说明本发明的技术方案,而非对其限制;本领域的普通技术人员应当理解:其依然可以对上述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种流量信号采样方法,用于油田注水井电磁流量计的流量信号采样,其特征在于,该方法包括:
步骤S1:获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
步骤S2:利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
步骤S3:利用递推数字滤波算法对经步骤S2进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
步骤S4:将步骤S3得到的稳定的采用信号作为流量信号输出。
2.根据权利要求1所述的方法,其特征在于,步骤S2包括:
按照设定的采样频率对脉冲信号采样;
采用FIR滤波函数对采样得到的信号进行移位、加权和卷积运算,实现低通滤波处理,该低通滤波处理操作用于滤除采样信号中的电路噪声和工频噪声。
3.根据权利要求1所述的方法,其特征在于,步骤S3包括:
对于经步骤S2进行低通滤波处理处理后的信号,利用单片机中断技术实时更新信号数值,通过依次对M个最新的信号数值求平均值,得到稳定的采样信号。
4.一种流量信号采样装置,用于油田注水井电磁流量计的流量信号采样,其特征在于,该装置包括:
输入模块,用于获取输入的脉冲信号,该脉冲信号是电磁流量计在测量注水流量过程中生成的;
低通滤波处理模块,用于利用低通滤波算法对获取的脉冲信号进行采样,并对采样得到的信号进行低通滤波处理;
递推数字滤波处理模块,用于利用递推数字滤波算法对经低通滤波处理模块进行低通滤波处理处理后的信号进行递推求平均值,得到稳定的采样信号;
输出模块,用于将递推数字滤波处理模块得到的稳定的采用信号作为流量信号输出。
5.根据权利要求4所述的装置,其特征在于,
所述低通滤波处理模块,具体用于:按照设定的采样频率对脉冲信号采样;采用FIR滤波函数对采样得到的信号进行移位、加权和卷积运算,实现低通滤波处理,该低通滤波处理操作用于滤除采样信号中的电路噪声和工频噪声。
6.根据权利要求4所述的装置,其特征在于,
所述递推数字滤波处理模块,具体用于:对于经所述低通滤波处理模块进行低通滤波处理处理后的信号,利用单片机中断技术实时更新信号数值,通过依次对M个最新的信号数值求平均值,得到稳定的采样信号。
7.一种单片机,包括处理器和存储器,其特征在于,所述存储器中存储有程序,当所述单片机运行时,所述处理器通过执行所述存储器中存储的程序,使所述单片机执行如权利要求1至3中任一所述的流量信号采样方法。
8.一种存储介质,存储有程序,当所述程序被单片机执行时,使所述单片机执行如权利要求1至3中任一所述的流量信号采样方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911363057.7A CN111060168A (zh) | 2019-12-26 | 2019-12-26 | 一种流量信号采样方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911363057.7A CN111060168A (zh) | 2019-12-26 | 2019-12-26 | 一种流量信号采样方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111060168A true CN111060168A (zh) | 2020-04-24 |
Family
ID=70303688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911363057.7A Pending CN111060168A (zh) | 2019-12-26 | 2019-12-26 | 一种流量信号采样方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111060168A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1714277A (zh) * | 2001-04-27 | 2005-12-28 | 米克罗利斯公司 | 对质量流控制器和质量流量计的输出滤波的系统和方法 |
US20070095123A1 (en) * | 2005-11-02 | 2007-05-03 | Ishida Co., Ltd. | Noise eliminating apparatus, weight measurement apparatus, method of eliminating noise, and method of designing digital filter |
CN101203732A (zh) * | 2005-05-20 | 2008-06-18 | 微动公司 | 电子计量器和快速从科里奥利流量计信号中确定多相流体的质量分数的方法 |
CN101268341A (zh) * | 2005-08-18 | 2008-09-17 | 微动公司 | 用于处理用于流量计中的多相流动材料的传感器信号的流量计电子器件和方法 |
CN101410698A (zh) * | 2006-03-14 | 2009-04-15 | 罗斯芒特公司 | 用于测量过程流体流量的磁流量计及方法 |
CN106595783A (zh) * | 2016-12-08 | 2017-04-26 | 中船重工环境工程有限公司 | 一种电磁流量计的非满管检测电路 |
US20190294649A1 (en) * | 2016-11-11 | 2019-09-26 | Oxford University Innovation Limited | Method and system for tracking sinusoidal wave parameters from a received signal that includes noise |
-
2019
- 2019-12-26 CN CN201911363057.7A patent/CN111060168A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1714277A (zh) * | 2001-04-27 | 2005-12-28 | 米克罗利斯公司 | 对质量流控制器和质量流量计的输出滤波的系统和方法 |
CN101203732A (zh) * | 2005-05-20 | 2008-06-18 | 微动公司 | 电子计量器和快速从科里奥利流量计信号中确定多相流体的质量分数的方法 |
CN101268341A (zh) * | 2005-08-18 | 2008-09-17 | 微动公司 | 用于处理用于流量计中的多相流动材料的传感器信号的流量计电子器件和方法 |
US20070095123A1 (en) * | 2005-11-02 | 2007-05-03 | Ishida Co., Ltd. | Noise eliminating apparatus, weight measurement apparatus, method of eliminating noise, and method of designing digital filter |
CN101410698A (zh) * | 2006-03-14 | 2009-04-15 | 罗斯芒特公司 | 用于测量过程流体流量的磁流量计及方法 |
US20190294649A1 (en) * | 2016-11-11 | 2019-09-26 | Oxford University Innovation Limited | Method and system for tracking sinusoidal wave parameters from a received signal that includes noise |
CN106595783A (zh) * | 2016-12-08 | 2017-04-26 | 中船重工环境工程有限公司 | 一种电磁流量计的非满管检测电路 |
Non-Patent Citations (1)
Title |
---|
张理恒 等: "基于DSP的电磁流量计信号处理", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102620729B (zh) | 一种机抖激光陀螺惯性测量单元数字滤波器设计方法 | |
CN105553443A (zh) | 一种强噪声条件下微弱信号提取及数字处理系统 | |
CN102332894A (zh) | 数字滤波器 | |
CN204989893U (zh) | 一种强噪声条件下微弱信号提取及数字处理系统 | |
JP5490828B2 (ja) | 線形系係数推定方法及びそれを用いた集積回路、タッチパネルシステム、及び電子機器 | |
CN105262503B (zh) | 一种基于群时延校准的多径时延产生装置及方法 | |
CN102435849A (zh) | 使用数字滤波器提高频率测量精度的方法 | |
CN116013240B (zh) | 基于变分模态分解与数字滤波的钢管混凝土信号降噪方法 | |
CN112487604A (zh) | 海洋重力仪输出数据长时间非线性漂移补偿方法 | |
CN104280047A (zh) | 一种多传感器融合的陀螺漂移滤波系统及方法 | |
CN111060168A (zh) | 一种流量信号采样方法和装置 | |
Bogdan et al. | LabVIEW modeling and simulation, of the low-pass and high-pass analog filters | |
CN105227158B (zh) | 一种高精度信号实时滤波器 | |
CN113054949A (zh) | 一种水击压力波信号的滤波方法、装置和设备 | |
CN106153029A (zh) | 二频机抖激光陀螺抖动信号抵消装置 | |
CN108832908A (zh) | 基于fpga的多路高速滤波器实现方法 | |
CN107636583A (zh) | 确定触摸位置的方法和触摸控制芯片 | |
CN110082642B (zh) | 基于全相位微分滤波的电网工况故障时刻检测方法及装置 | |
CN102664724B (zh) | 数字通信系统中的数字信号的处理方法和装置 | |
CN102455707A (zh) | 一种基于无限冲击响应数字滤波技术的舰船姿态传感器 | |
CN207518557U (zh) | 一种激光陀螺小抖动稳频无延迟带通滤波器 | |
CN104198205A (zh) | 一种水下机器人多普勒故障随机共振检测装置及检测方法 | |
CN117033910B (zh) | 一种海面高精度信号的处理方法及系统 | |
CN109214092A (zh) | 基于希尔伯特黄变换的四分量钻孔应变数据异常提取方法 | |
CN106551686A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200424 |
|
RJ01 | Rejection of invention patent application after publication |