发明名称 一种基于示功图的有杆泵抽油井的多故障诊断方法
摘要 本发明公开了一种基于示功图的有杆泵抽油井的多故障诊断方法,涉及石油生产技术领域。将采集的地面示功图转换为井下泵功图,根据油田生产中已有的示功图数据建立各故障类型的标准特征库,对于每一个不变曲线矩特征向量的特征值,由区间数据的形式进行表示,然后计算所采集样本与标准特征库中各故障类型的关联度,判断所采集样本可能具有的故障类型,提高了故障诊断的可信性。最后将所采集样本可能具有的故障类型进行各种多个故障的组合,通过计算各个组合的F指标值来确定所具有的多故障类型。该方法原理简单,计算复杂性小,容易实现,诊断的准确度高。
申请公布号 CN106121622A 申请公布日期 2016.11.16
申请号 CN201610597512.X 申请日期 2016.07.27
申请人 渤海大学 发明人 李琨;韩莹;佘东生;魏泽飞;杨一柳;于震
分类号 E21B44/00(2006.01)I;E21B47/008(2012.01)I;G06F19/00(2011.01)I 主分类号 E21B44/00(2006.01)I
代理机构 西安铭泽知识产权代理事务所(普通合伙) 61223 代理人 俞晓明
主权项 一种基于示功图的有杆泵抽油井的多故障诊断方法,其特征在于,包括:采集地面示功图,并将所述地面示功图转化为井下泵功图;对所述井下泵功图进行归一化处理,归一化区间为[0,1],提取归一化后的井下泵功图的不变曲线矩特征向量,并对提取的不变曲线矩特征向量进行修正;由修正后的不变曲线矩特征向量建立所述井下泵功图的物元模型,如下:<img file="FDA0001061013010000011.GIF" wi="1302" he="342" />其中,δ<sub>1</sub>,δ<sub>2</sub>,...,δ<sub>7</sub>分别为7个所述修正后的不变曲线矩特征向量,v<sub>1</sub>,v<sub>2</sub>,...,v<sub>7</sub>分别为7个所述修正后的不变曲线矩特征向量的取值;由油田生产中所记录的已知不同故障类型的地面示功图样本建立标准特征库,并建立每一种故障类型的物元模型,如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>M</mi><mi>j</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>W</mi><mi>j</mi></msub></mtd><mtd><msub><mi>&delta;</mi><mn>1</mn></msub></mtd><mtd><mrow><mo>&lsqb;</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>1</mn><mi>a</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>1</mn><mi>b</mi></mrow></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><msub><mi>&delta;</mi><mn>2</mn></msub></mtd><mtd><mrow><mo>&lsqb;</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>2</mn><mi>a</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>2</mn><mi>b</mi></mrow></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><msub><mi>&delta;</mi><mn>7</mn></msub></mtd><mtd><mrow><mo>&lsqb;</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>7</mn><mi>a</mi></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mn>7</mn><mi>b</mi></mrow></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001061013010000012.GIF" wi="1334" he="358" /></maths>其中,j=1,2,...,T,T表示已知故障类型的数量,W<sub>j</sub>表示第j种故障类型;[v<sub>j1a</sub>,v<sub>j1b</sub>]表示所述标准特征库中第j种故障类型下的第1个修正后的不变曲线矩特征向量的取值区间,[v<sub>j2a</sub>,v<sub>j2b</sub>]表示所述标准特征库中第j种故障类型下的第2个修正后的不变曲线矩特征向量的取值区间;依次类推,[v<sub>j7a</sub>,v<sub>j7b</sub>]表示所述标准特征库中第j种故障类型下的第7个修正后的不变曲线矩特征向量的取值区间;计算所述井下泵功图与所述标准特征库中各故障类型的关联度,计算公式如下:<maths num="0002"><math><![CDATA[<mrow><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msub><mi>v</mi><mi>k</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi><mi>a</mi></mrow></msub><mo>+</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi><mi>b</mi></mrow></msub></mrow><mn>2</mn></mfrac><mo>|</mo><mo>-</mo><mfrac><mrow><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi><mi>b</mi></mrow></msub><mo>-</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi><mi>a</mi></mrow></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001061013010000021.GIF" wi="1348" he="168" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>&xi;</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><mi>X</mi><mo>)</mo></mrow><mo>-</mo><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>v</mi><mi>k</mi></msub><mo>&NotElement;</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mfrac><mrow><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>|</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>|</mo></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>v</mi><mi>k</mi></msub><mo>&Element;</mo><msub><mi>v</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001061013010000022.GIF" wi="1384" he="367" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>&lambda;</mi><mrow><mo>(</mo><msub><mi>W</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>7</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><msub><mi>&xi;</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001061013010000023.GIF" wi="1172" he="157" /></maths>其中,k=1,2,...,7,v<sub>k</sub>表示所述井下泵功图第k个修正后的不变曲线矩特征向量的取值,v<sub>jk</sub>表示所述标准特征库中第j种故障类型下第k个修正后的不变曲线矩特征向量的取值,v<sub>jka</sub>和v<sub>jkb</sub>分别表示所述标准特征库中第j种故障类型下第k个修正后的不变曲线矩特征向量取值的最小值和最大值,|v<sub>jk</sub>|表示所述标准特征库中第j种故障类型下第k个修正后的不变曲线矩特征向量取值范围的大小;节域X=[0,20];λ(W<sub>j</sub>)表示所述井下泵功图与所述标准特征库中第j种故障类型的关联度,ε(v<sub>k</sub>,v<sub>jk</sub>)表示所述井下泵功图第k个修正后的不变曲线矩特征向量与所述标准特征库中第j种故障类型下第k个修正后的不变曲线矩特征向量的距离,ξ<sub>jk</sub>(v<sub>k</sub>)表示第k个修正后的不变曲线矩特征向量下所述井下泵功图与所述标准特征库中第j种故障类型的关联函数,ε(v<sub>k</sub>,X)表示所述井下泵功图第k个修正后的不变曲线矩特征向量与节域的距离;对λ(W<sub>j</sub>)≥0的各故障类型进行多故障组合的分类,将一组多个故障的组合看作一类,每一个所述多个故障的组合中均至少包括λ(W<sub>j</sub>)最大的故障类型以及所述井下泵功图;对λ(W<sub>j</sub>)≥0的故障类型的各修正后的不变曲线矩特征向量进行基于“无偏转换”的区间数据的均值估计,并计算所述均值估计的无偏估计值;根据所述均值估计的无偏估计值以及以下公式计算每一组多个故障的组合的F指标值,并确定具有最小F指标值的多个故障组合为所述井下泵功图所具有的多故障类型:<img file="FDA0001061013010000031.GIF" wi="1685" he="318" />其中,k=1,2,...,num,num表示每一组多个故障组合中故障类型的个数;expt<sub>k</sub>表示第k个故障类型下进行实验的个数,即为第k个故障类型下修正后的不变曲线矩特征向量的个数,expt<sub>k</sub>=7,expt为不同故障类型下所有实验的个数,有:<img file="FDA0001061013010000032.GIF" wi="438" he="159" /><img file="FDA0001061013010000033.GIF" wi="105" he="94" />表示第k个故障类型下第<img file="FDA00010610130100000314.GIF" wi="31" he="47" />个修正后的不变曲线矩特征向量的区间数据的均值估计的无偏估计值,<img file="FDA00010610130100000315.GIF" wi="421" he="67" /><img file="FDA0001061013010000034.GIF" wi="101" he="103" />表示第k个故障类型下所有修正后的不变曲线矩特征向量的区间数据的均值估计的无偏估计值的平均值,有:<img file="FDA0001061013010000035.GIF" wi="565" he="183" /><img file="FDA0001061013010000036.GIF" wi="103" he="111" />表示所有故障类型下的所有修正后的不变曲线矩特征向量的区间数据的均值估计的无偏估计值的总平均值,有:<img file="FDA0001061013010000037.GIF" wi="598" he="174" />sum<sub>1</sub>表示<img file="FDA0001061013010000038.GIF" wi="104" he="103" />与<img file="FDA0001061013010000039.GIF" wi="108" he="109" />的差异之和,有:<img file="FDA00010610130100000310.GIF" wi="694" he="143" />Sum<sub>2</sub>表示<img file="FDA00010610130100000311.GIF" wi="106" he="95" />与<img file="FDA00010610130100000312.GIF" wi="100" he="103" />的差异之和,有:<img file="FDA00010610130100000313.GIF" wi="662" he="159" />
地址 121013 辽宁省锦州市松山新区科技路19号渤海大学