CN103924963B - 一种示功仪采样率自动切换方法 - Google Patents
一种示功仪采样率自动切换方法 Download PDFInfo
- Publication number
- CN103924963B CN103924963B CN201410169545.5A CN201410169545A CN103924963B CN 103924963 B CN103924963 B CN 103924963B CN 201410169545 A CN201410169545 A CN 201410169545A CN 103924963 B CN103924963 B CN 103924963B
- Authority
- CN
- China
- Prior art keywords
- indicator
- frequency
- actual stroke
- group
- stroke
- 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.)
- Active
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种示功仪采样率自动切换方法,根据所述示功图估计冲程频率计算出示功图实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出的该组采样数据实际冲程周期采集点数与示功仪预设的单周期采集点数范围值进行对比,适时的调整示功仪的采样率,以达到精确的采样数据形成精确功图的同时,避免采样点过多导致的示功仪的功耗过大,影响示功仪使用寿命的问题。
Description
技术领域
本发明涉及油井测量领域,尤其涉及一种示功仪采样率自动切换方法。
背景技术
抽油机井抽油是石油生产中最重要的工艺过程之一。及时准确获知的油井运行状态可考察油井单位时间的产量、跟踪掌握地下油藏区块的动态情况、预测和评价油井区块的开发潜力、总体分析和把握全局产能情况、合理制定油田生产方案。
油井示功图是油井运行状态的重要表征,可以全面反映井下抽油泵运行状况以及原油的开采情况。示功图用来反映深井泵工作情况好坏,可由专门的仪器测出,并在坐标纸上进行制图,图上被封闭的线所围成的面积表示驴头在依次往复运动中抽油泵所做的功。通过分析示功图,可以了解油井动态及抽油装置的各项参数选择是否合理、反映出深井泵在井下工作中的异常现象,以及可以通过分析获取抽油机井的工作制度是否合理等等信息。
传统对示功图的测量是一个月一次的人工测量,但是这种人工的测量方式,工作量大而繁重,且难以及时、准确、全面直观、连续地反应生产动态。鉴于传统人工测量方式的上述问题,现有技术发展了可以进行远程监控的无线示功仪。现有的无线示功仪对一个冲程周期的功图信号进行采样,然后通过无线网络将采集到的数据上传至上位机,通过上位机绘制功图。根据实际需要,示功仪每次所采集的数据应该是正好包括抽油机一个冲程,因此示功仪采样周期的精确性决定采集的功图信号的完整性,是整个设备可靠性的关键。一般而言,对于某个抽油井一个冲程周期的采样点数不应少于150点,确保所绘制的示功图的精确性。
对于周期信号的周期估计常用的经典方法有傅里叶变换法,自相关法,同态滤波法和AMDF(平均幅度差函数法)。
傅里叶变换法利用频谱分析计算信号的频率,该方法只适合分析正弦信号,且参数估计存在偏差。对于复杂的周期信号,单周期信号频谱是连续的,而多周期信号频谱是离散,对于单周期信号频谱的离散采样,即存在大量的谐波分量,因此傅里叶变换法得到的最强频谱点可能只是信号的谐波频率。
自相关法利用信号的自相关函数求解信号周期,而有限采样信号的自相关函数是阶梯减小的,给周期估计带来不便,同时常会导致估计偏差。
AMDF方法由Ross等人提出,该方法利用最小值间隔确定周期,对于有限长度信号,该方法的输出最小值是阶梯增加的,给求解带来不便。一般而言,AMDF和估计精度较自相关法略好,但实际仿真表明,估计精度在一定的程度上和周期信号的样式有关。
解卷积同态系统是求取信号周期的另一种方法,利用卷积同态滤波得到的低时和高时部分,然后恢复单个周期信号和冲击函数序列,冲击函数周期即为信号周期。由于需要通过复倒频谱的时域滤波,滤波截断窗选取会对信号周期估计产生影响,时域截断滤波的选取不容易控制,受噪音影响较大,因此该方法估计性能差。
而上述四种方法对于功图的冲程周期的估计并不适用,这是因为我们要分析的信号频率非常低。目前,在役的几乎全部的有杆泵抽油机的冲次均在1.5次/分~12次/分之间,也就是示功图采集到的信号的周期一般在5~40秒(频率0.025~0.2赫兹)之间。另外,油田规定,对于功图信号周期的测量,一般要求精确到小数点后一位,也就是需要精确到0.1秒。而在0.025~0.2赫兹之间,在频率上计算的任何微小的偏差与当前频率转换为周期后,则误差会被大幅度放大,无法满足油田提出的0.1秒的要求。而要达到要求的时间精度,信号采集的时间就会长到我们难以接受的时间,显然是不现实的。
而对于这些频率极低的信号来说,现有技术中通常利用两种方法进行采样周期估计,一种是通过过零比较,把原来的类似于正弦波的信号变成方波,然后通过测量方波脉宽的方法来取得信号周期,由于示功图测量的加速度信号上边有些脉冲甚至高于周期正弦波信号,通过这种过零比较的方法很难按正确的信号周期求解出来。还有一种方法是确定过零点后,拟合出一条穿过过零点的直线,通过最小二乘法求得直线的斜率,然后根据直线斜率求取信号周期,由于示功图上行下行受力不一致,所以信号的图并不是完全正弦波,上行下行波峰是有些偏心的,因此采用这种方法也很难求得正确的周期。
但是,在现有的示功仪中,通常仅采用频域分析的方法计算抽油井的冲程周期,而由于数据长度以及频域分析阶数对分辨率存在一定影响,导致该冲程周期作为采样周期进行采样并不精确,而当实际的冲程周期中采样点过少时,不能得到精确的采样数据形成精确的功图,当采样点过多又会使得示功仪的功耗过大,影响示功仪的使用寿命。
发明内容
有鉴于此,为解决上述技术问题,本发明提供了一种示功仪采样率自动切换方法。
本发明采用的技术手段如下:一种示功仪采样率自动切换方法,其特征在于,包括:
以初始采样率进行数据采样,得到一组采样数据;
根据该组采样数据进行快速傅氏变换运算,计算出示功图估计冲程频率;
根据所述示功图估计冲程频率计算出示功图实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出该组采样数据实际冲程周期的采集点数;
根据所述实际冲程周期的采集点数判断所述初始采样率是否合适,当采样率不合适时,选择合适的采样率。
进一步,所述以初始采样率进行数据采样,得到一组采样数据包括:
获取示功仪的实际冲程周期值范围,以初始采样率N进行数据采样,得到一组采样数据{Ai},i∈[0,n*N-1],其中n为正整数,且n为示功仪实际冲程周期值范围中最大实际冲程周期值的至少两倍。
进一步,根据该组采样数据进行快速傅氏变换运算,计算出示功图估计冲程频率包括:
对采样数据{Ai}进行快速傅氏变换,得到一组谱线;
设该组谱线中第M个谱线为幅值最大的谱线,Δf为该组谱线中任意两相邻谱线之间的间隔,则示功图估计冲程频率为Fk=M*Δf。
进一步,当该组谱线存在多个幅值相同的最大幅值时,以频点最小的谱线为第M个谱线。
进一步,根据所述示功图估计冲程频率计算出示功图实际冲程的最大可能频率和最小可能频率包括:
根据估计示功图实际冲程频率Fk,则示功图实际冲程的最小可能频率为Fk-Δf/2,示功图实际冲程的最大可能频率为Fk+Δf/2。
进一步,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出该组采样数据实际冲程周期的采集点数包括:
设L为示功图实际冲程周期的采集点数,则
计算两个相邻示功图实际冲程周期的幅值差的平方的单点平均值为
计算Δ(L)最小时L的值L1,作为示功图实际冲程周期的采集点数。
进一步,根据所述实际冲程周期的采集点数判断所述初始采样率是否合适,当采样率不合适时,选择合适的采样率包括:
比较示功图实际冲程周期的采集点数L1与示功仪预设的单周期采集点数范围值,当L1不属于所述示功仪预设的单周期采集点数范围值时,调整示功仪的采样率,以使示功图实际冲程周期的采集点数L1属于所述示功仪预设的单周期采集点数范围值。
进一步,所述示功仪预设的单周期采集点数范围值为150至300。
本发明所提供的示功仪采样率自动切换方法,根据所述示功图估计冲程频率计算出示功图实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出的该组采样数据实际冲程周期采集点数与示功仪预设的单周期采集点数范围值进行对比,适时的调整示功仪的采样率,以达到精确的采样数据形成精确功图的同时,避免采样点过多导致的示功仪的功耗过大,影响示功仪使用寿命的问题。
附图说明
图1为本发明一种示功仪采样率自动切换方法流程示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
本发明是基于以下考虑实现的:
对于固定的抽油机井而言,其实际冲程周期在一段时间内是不会发生变化的,换而言之,即相邻两实际冲程周期可默认为相同,并且示功图的采样信号为频率极低的周期信号。对于已知初始采样率的抽油机而言,可基于上述客观条件对抽油机实际冲程频率范围进行估计后,通过精确估计示功图实际冲程周期的采集点数调整示功仪的采样率。
本发明提供了一种示功仪采样率自动切换方法,包括:
以初始采样率进行数据采样,得到一组采样数据;具体的,获取示功仪的实际冲程周期值范围,以初始采样率N进行数据采样,得到一组采样数据{Ai},i∈[0,n*N-1],其中n为正整数,且n为示功仪实际冲程周期值范围中最大实际冲程周期值的至少两倍;
根据该组采样数据进行快速傅氏变换运算,计算出估计冲程频率;具体的,对采样数据{Ai}进行快速傅氏变换,得到一组谱线;设该组谱线中第M个谱线为幅值最大的谱线,Δf为该组谱线中任意两相邻谱线之间的间隔,则估计冲程频率为Fk=M*Δf,其中,当该组谱线存在多个幅值相同的最大幅值时,以频点最小的谱线为第M个谱线;
根据所述估计冲程频率计算出实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出该组采样数据实际冲程周期的采集点数;具体的,首先,根据估计示功图实际冲程频率Fk,则示功图实际冲程的最小可能频率为Fk-Δf/2,示功图实际冲程的最大可能频率为Fk+Δf/2;然后,设L为示功图实际冲程周期的采集点数,则计算两个相邻示功图实际冲程周期的幅值差的平方的单点平均值为计算Δ(L)最小时L的值L1,作为示功图实际冲程周期的采集点数;
根据所述实际冲程周期的采集点数判断所述初始采样率是否合适,当采样率不合适时,选择合适的采样率;具体的,可通过比较示功图实际冲程周期的采集点数L1与示功仪预设的单周期采集点数范围值,当L1不属于所述示功仪预设的单周期采集点数范围值时,调整示功仪的采样率,以使示功图实际冲程周期的采集点数L1属于所述示功仪预设的单周期采集点数范围值。
作为本发明的典型实施例,假定全部的功图周期在5~40秒之间,因此,此时需要采集80秒的数据(因为需要进行两个周期比较,所以至少两个周期),假设采样率为N,这样我们得到一组数据序列:
{Ai},i∈[0,80*N-1]
首先需要对此序列进行FFT变换,计算完毕后,求取幅值最大的第M个谱线,用M乘以谱线的间隔Δf,则可以得到功图信号的频率为Fk=M*Δf。由于每条谱线的间隔为Δf,则会对测量频率造成Δf/2的频率偏差,也就是实际信号频率可能在[Fk-Δf/2,Fk+Δf/2]之间。取倒数后可以得到信号的周期范围在由于信号的频率是小于1Hz的,所以转换成周期时间后,周期时间偏差会被放大,所以还需要进行二次比较提高精度。换句话说,通过这种算法求得了一个功图周期的范围此数据乘以采样率N,则可以得到组成一个周期的采样点个数可能的范围在之间。
如果L为周期的采集点数,那么,则有且L为正整数,两个相邻周期的幅值差的平方的单点平均值就为求取当Δ(L)最小时L的值L1,此值即为该信号周期的采集点数,为该功图的周期点数。此时,那么信号周期的时间长度就为点数除以采样率,即周期为T1=L1/N。
由于油田规定一个合理的功图的采样点数应该在150~300点之间,过少的点数会导致功图描述不正确,而太多的点数则会使存储区白白浪费,基于此,设定求得的功图周期在5~10秒之间时,采样率为30SPS;周期在10~20秒之间时,采样率为15SPS;周期在20~40秒之间时,采样率为7.5SPS,这样,就可以确保每个周期的点数控制在150~300之间。当确定完示功图实际冲程周期的采集点数L1后,可计算示功图的实际冲程周期T1=L1/N,如果周期在5~10秒之间,则可以直接使用数据;如果周期在10~20秒之间,则对原数据进行滤波后再进行抽值运算,使输出采样率变为15SPS;如果周期在20~40秒之间,则对原数据进行两次滤波和抽值运算,使输出采样率变为7.5SPS。
采用本发明提供的示功仪采样率自动切换方法,根据所述示功图估计冲程频率计算出示功图实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出的该组采样数据实际冲程周期采集点数与示功仪预设的单周期采集点数范围值进行对比,适时的调整示功仪的采样率,以达到精确的采样数据形成精确功图的同时,避免采样点过多导致的示功仪的功耗过大,影响示功仪使用寿命的问题。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (7)
1.一种示功仪采样率自动切换方法,其特征在于,包括:
以初始采样率进行数据采样,得到一组采样数据;
根据该组采样数据进行快速傅氏变换运算,计算出估计冲程频率;
根据所述估计冲程频率计算出实际冲程的最大可能频率和最小可能频率,使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出该组采样数据实际冲程周期的采集点数;
根据所述实际冲程周期的采集点数判断所述初始采样率是否合适,当采样率不合适时,选择合适的采样率;其中,
使用最小二乘法拟合周期函数同时使用最大和最小可能频率计算出该组采样数据实际冲程周期的采集点数包括:
设L为示功图实际冲程周期的采集点数,则
计算两个相邻示功图实际冲程周期的幅值差的平方的单点平均值为
计算Δ(L)最小时L的值L1,作为示功图实际冲程周期的采集点数;
其中,N为初始采样率,Fk为估计冲程频率,Ai为采样数据{Ai},i∈[0,n*N-1]中的第i个,Δf为对采样数据{Ai}进行快速傅氏变换所得到的一组谱线中任意两相邻谱线之间的间隔;
其中,Fk-Δf/2为示功图实际冲程的最小可能频率,Fk+Δf/2为示功图实际冲程的最大可能频率。
2.根据权利要求1所述的方法,其特征在于,所述以初始采样率进行数据采样,得到一组采样数据包括:
获取示功仪的实际冲程周期值范围,以初始采样率N进行数据采样,得到一组采样数据{Ai},i∈[0,n*N-1],其中n为正整数,且n为示功仪实际冲程周期值范围中最大实际冲程周期值的至少两倍。
3.根据权利要求2所述的方法,其特征在于,根据该组采样数据进行快速傅氏变换运算,计算出估计冲程频率包括:
对采样数据{Ai}进行快速傅氏变换,得到一组谱线;
设该组谱线中第M个谱线为幅值最大的谱线,Δf为该组谱线中任意两相邻谱线之间的间隔,则估计冲程频率为Fk=M*Δf。
4.根据权利要求3所述的方法,其特征在于,当该组谱线存在多个幅值相同的最大幅值时,以该多个幅值相同的最大幅值中的频点最小的谱线为第M个谱线。
5.根据权利要求4所述的方法,其特征在于,根据所述估计冲程频率计算出实际冲程的最大可能频率和最小可能频率包括:
根据估计冲程频率Fk,则示功图实际冲程的最小可能频率为Fk-Δf/2,示功图实际冲程的最大可能频率为Fk+Δf/2。
6.根据权利要求1所述的方法,其特征在于,根据所述实际冲程周期的采集点数判断所述初始采样率是否合适,当采样率不合适时,选择合适的采样率包括:
比较示功图实际冲程周期的采集点数L1与示功仪预设的单周期采集点数范围值,当L1不属于所述示功仪预设的单周期采集点数范围值时,调整示功仪的采样率,以使示功图实际冲程周期的采集点数L1属于所述示功仪预设的单周期采集点数范围值。
7.根据权利要求6所述的方法,其特征在于,所述示功仪预设的单周期采集点数范围值为150至300。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410169545.5A CN103924963B (zh) | 2014-04-25 | 2014-04-25 | 一种示功仪采样率自动切换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410169545.5A CN103924963B (zh) | 2014-04-25 | 2014-04-25 | 一种示功仪采样率自动切换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103924963A CN103924963A (zh) | 2014-07-16 |
CN103924963B true CN103924963B (zh) | 2017-02-15 |
Family
ID=51143329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410169545.5A Active CN103924963B (zh) | 2014-04-25 | 2014-04-25 | 一种示功仪采样率自动切换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103924963B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104564021A (zh) * | 2014-04-21 | 2015-04-29 | 山东天工石油装备有限公司 | 油井示功图数据采集的一种方法 |
CN107544334B (zh) * | 2017-09-18 | 2019-07-19 | 吉林大学 | 自动调整采样率的数据采集卡 |
CN107576988B (zh) * | 2017-09-29 | 2021-06-18 | 中国科学院地质与地球物理研究所 | 一种星载磁强计采样率自动切换方法及装置 |
CN110185437A (zh) * | 2019-07-22 | 2019-08-30 | 南京南大数字科技有限责任公司 | 一种适用于低冲次变频油井的示功图采集方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4509901A (en) * | 1983-04-18 | 1985-04-09 | Fmc Corporation | Method and apparatus for detecting problems in sucker-rod well pumps |
CN2558767Y (zh) * | 2001-12-28 | 2003-07-02 | 石油大学(北京) | 抽油机间抽与功率切换控制装置 |
CN101305187B (zh) * | 2005-10-13 | 2010-12-08 | 井泵技术有限公司 | 井下流体产量优化系统及方法 |
CN102072788A (zh) * | 2009-11-19 | 2011-05-25 | 上海神开石油化工装备股份有限公司 | 无线示功仪及示功图测量方法 |
CN101950169A (zh) * | 2010-09-06 | 2011-01-19 | 北京时代科仪新能源科技有限公司 | 抽油机控制器、控制方法及系统 |
CN102288816A (zh) * | 2011-10-23 | 2011-12-21 | 黑龙江大学 | 抽油机多用途能效分析仪及测量方法 |
CN102966346A (zh) * | 2012-11-09 | 2013-03-13 | 西安冠林智能科技有限公司 | 基于Zigbee网络的无线低功耗油井示功仪 |
CN103293990B (zh) * | 2013-05-16 | 2016-03-23 | 贵州航天林泉电机有限公司 | 一种永磁同步电机直驱式抽油机控制器及控制方法 |
CN203241817U (zh) * | 2013-05-29 | 2013-10-16 | 中国石油天然气集团公司 | 抽油机运行参数在线智能分析及控制系统 |
CN103573252A (zh) * | 2013-11-05 | 2014-02-12 | 上海理工大学 | 无线示功仪、无线示功仪系统及其示功图测量方法 |
-
2014
- 2014-04-25 CN CN201410169545.5A patent/CN103924963B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103924963A (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100573381C (zh) | 用于过程控制系统的诊断设备 | |
US8898017B2 (en) | Automated hydrocarbon reservoir pressure estimation | |
CN103924963B (zh) | 一种示功仪采样率自动切换方法 | |
EP2494145B1 (en) | Automated hydrocarbon reservoir pressure estimation | |
CN102087787B (zh) | 短时交通状态预测装置及预测方法 | |
CN103995178A (zh) | 一种基于时频聚集特性准则s变换的电压暂降检测方法 | |
CN104406680B (zh) | 一种电力变压器表面测点振动加速度信号特征的提取方法 | |
CN104931040B (zh) | 基于机器学习的北斗ⅱ代导航系统电力铁塔变形监测设备安装和调试方法 | |
US20190316942A1 (en) | Methodologies and apparatus for the recognition of production tests stability | |
CN102425409B (zh) | 一种井斜和工具面测量装置及垂直钻井设备 | |
CN103149276A (zh) | 确定模型水轮机转轮叶片空蚀起始点的方法 | |
CN106250904A (zh) | 基于改进s变换的电能扰动分析仪及分类方法 | |
CN104198811B (zh) | 低频信号频率测量方法和装置 | |
CN103995147B (zh) | 一种适用于声学多普勒流速仪的数据后处理系统与应用 | |
CN111651728A (zh) | 一种辨识实测风速非平稳性与非平稳度的方法 | |
CN105569639A (zh) | 一种基于地面示功图计算油井动液面的方法 | |
US10832152B2 (en) | Wet well prediction using real time data | |
US8306762B2 (en) | Systems and methods for analysis of downhole data | |
CN201974940U (zh) | 短时交通状态预测装置 | |
Zhang et al. | Mathematical model of time difference for Coriolis flow sensor output signals under gas-liquid two-phase flow | |
KR102019075B1 (ko) | 스펙트럼의 피로손상도 계산 최적화 시스템 및 그 방법 | |
CN110486009A (zh) | 一种无限大地层的参数自动反求方法及系统 | |
CN104777356A (zh) | 一种基于神经网络的实时高精度谐波检测方法 | |
CN205721844U (zh) | 基于改进s变换的电能扰动分析仪 | |
CN104808643A (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 | ||
CB02 | Change of applicant information |
Address after: 100085 Beijing city Haidian District on the seven street, Huizhong building 710 Applicant after: BEIJING WILL CREATE TECHNOLOGY CO., LTD. Address before: 100085 Beijing city Haidian District on the seven street, Huizhong building room 710 Applicant before: Beijing Beetech Technology Inc. |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |