发明名称 高转差率电动机驱动的游梁式抽油机平衡判定调整方法
摘要 一种高转差率电动机驱动游梁式抽油机的平衡判定方法,利用电动机的磁链观测方法计算抽油机电动机的转速,结合抽油机动力学分析,进而得出抽油机曲柄的运动规律,使曲柄运动方程和抽油杆柱的波动方程之间解除耦合关系,分别可以对抽油机悬点进行运动学和动力学的分析,分别计算出抽油机悬点的速度、加速度和位移,以及悬点的载荷,计算出抽油机曲柄扭矩,将上、下冲程的最大扭矩进行比较,其比值作为游梁平衡的判定依据,该方法准确度高、抗干扰性强。
申请公布号 CN106522901A 申请公布日期 2017.03.22
申请号 CN201611148754.7 申请日期 2016.12.13
申请人 西安石油大学 发明人 徐向前;周好斌
分类号 E21B43/12(2006.01)I;E21B47/009(2012.01)I 主分类号 E21B43/12(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 弋才富
主权项 高转差率电动机驱动游梁式抽油机平衡判定方法,其特征在于,包括以下步骤:步骤一、计算电动机转速采用三相交流异步电动机定子电压和电流的磁链观测模型,进行了电动机转速的计算,电动机定子电压、电流的磁链观测模型为:ψ<sub>s</sub>=∫(U<sub>s</sub>‑I<sub>s</sub>R<sub>s</sub>)dt式中,ψ<sub>s</sub>为定子磁链,Wb;U<sub>s</sub>为定子电压,V;I<sub>s</sub>为定子电流,A;R<sub>s</sub>为定子电阻,Ω;在ABC坐标系下,定子电压u<sub>a</sub>,u<sub>b</sub>,u<sub>c</sub>和电流i<sub>a</sub>,i<sub>b</sub>,i<sub>c</sub>可以通过坐标变换得到α‑β坐标系下的电压u<sub>α</sub>,u<sub>β</sub>和电流i<sub>α</sub>,i<sub>β</sub>:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><msqrt><mfrac><mn>2</mn><mn>3</mn></mfrac></msqrt><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac></mtd><mtd><mrow><mo>-</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mi>a</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>c</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001179373240000011.GIF" wi="659" he="271" /></maths><maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>i</mi><mi>&alpha;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>&beta;</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><msqrt><mfrac><mn>2</mn><mn>3</mn></mfrac></msqrt><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></mtd><mtd><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac></mtd><mtd><mrow><mo>-</mo><mfrac><msqrt><mn>3</mn></msqrt><mn>2</mn></mfrac></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>i</mi><mi>a</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>b</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>c</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001179373240000012.GIF" wi="629" he="271" /></maths>根据α‑β坐标系下的电压,求出定子电压的幅值u<sub>s</sub>和相位<img file="FDA0001179373240000013.GIF" wi="75" he="48" />如下:<img file="FDA0001179373240000014.GIF" wi="638" he="119" />同样可以得出定子电流的幅值i<sub>s</sub>和相位<img file="FDA0001179373240000015.GIF" wi="76" he="47" />如下:<img file="FDA0001179373240000016.GIF" wi="580" he="135" />α‑β坐标系下,磁链分量ψ<sub>α</sub>和ψ<sub>β</sub>,如下:<maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>&psi;</mi><mi>&alpha;</mi></msub><mo>=</mo><mo>&Integral;</mo><mrow><mo>(</mo><msub><mi>u</mi><mi>&alpha;</mi></msub><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>&alpha;</mi></msub><mo>)</mo></mrow><mi>d</mi><mi>t</mi></mrow></mtd><mtd><mrow><msub><mi>&psi;</mi><mi>&beta;</mi></msub><mo>=</mo><mo>&Integral;</mo><mrow><mo>(</mo><msub><mi>u</mi><mi>&beta;</mi></msub><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>&beta;</mi></msub><mo>)</mo></mrow><mi>d</mi><mi>t</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001179373240000021.GIF" wi="834" he="87" /></maths>定子磁链的幅值ψ<sub>sm</sub>和相位ψ<sub>ψ</sub>,如下:<maths num="0004"><math><![CDATA[<mrow><msub><mi>&psi;</mi><mrow><mi>s</mi><mi>m</mi></mrow></msub><mo>=</mo><msqrt><mrow><msubsup><mi>&psi;</mi><mi>&alpha;</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&psi;</mi><mi>&beta;</mi><mn>2</mn></msubsup></mrow></msqrt><mo>,</mo><msub><mi>&psi;</mi><mi>&psi;</mi></msub><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>s</mi><mi>i</mi><mi>n</mi><mfrac><msub><mi>&psi;</mi><mi>&beta;</mi></msub><msub><mi>&psi;</mi><mrow><mi>s</mi><mi>m</mi></mrow></msub></mfrac></mrow>]]></math><img file="FDA0001179373240000022.GIF" wi="718" he="127" /></maths>由计算得到的定子磁链和定子电流,可以得到异步电动机电磁转矩T<sub>e</sub>,如下:T<sub>e</sub>=p<sub>n</sub>(ψ<sub>α</sub>i<sub>β</sub>‑ψ<sub>β</sub>i<sub>α</sub>)式中,p<sub>n</sub>为极对数;由以上计算得到的电磁转矩、定子磁链,可以计算出电动机转子转差角速度ω<sub>f</sub>,如下:<maths num="0005"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>f</mi></msub><mo>=</mo><msub><mi>R</mi><mi>r</mi></msub><msub><mi>T</mi><mi>e</mi></msub><mo>/</mo><mrow><mo>(</mo><msub><mi>P</mi><mi>n</mi></msub><msubsup><mi>&psi;</mi><mrow><mi>s</mi><mi>m</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001179373240000023.GIF" wi="374" he="70" /></maths>根据ABC坐标系下的电动机定子电压方程,可以计算出电动机定子同步角速度ω<sub>s</sub>,如下:<maths num="0006"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>s</mi></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>u</mi><mi>&beta;</mi></msub><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>&beta;</mi></msub><mo>)</mo><msub><mi>&psi;</mi><mi>&alpha;</mi></msub><mo>-</mo><mo>(</mo><msub><mi>u</mi><mi>&alpha;</mi></msub><mo>-</mo><msub><mi>R</mi><mi>s</mi></msub><msub><mi>i</mi><mi>a</mi></msub><mo>)</mo><msub><mi>&psi;</mi><mi>&beta;</mi></msub></mrow><msubsup><mi>&psi;</mi><mrow><mi>s</mi><mi>m</mi></mrow><mn>2</mn></msubsup></mfrac></mrow>]]></math><img file="FDA0001179373240000024.GIF" wi="694" he="136" /></maths>则,电动机转子机械转速ω,如下:ω=(ω<sub>s</sub>‑ω<sub>f</sub>)/p<sub>n</sub>根据电动机在ABC坐标系下的三相电压电流转换到α‑β坐标系下,利用定子磁链观测的方法,得出了电动机转速;步骤二、曲柄轴净扭矩计算将曲柄作为等效构件,建立等效数学模型,求解减速器曲柄轴净扭矩,曲柄运动规律微分方程如下:<maths num="0007"><math><![CDATA[<mrow><msub><mi>J</mi><mi>e</mi></msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msup><mfrac><mrow><msub><mi>dJ</mi><mi>e</mi></msub></mrow><msub><mi>d</mi><mi>&theta;</mi></msub></mfrac><mo>=</mo><msub><mi>M</mi><mrow><mi>e</mi><mi>d</mi></mrow></msub><mo>-</mo><msub><mi>M</mi><mrow><mi>e</mi><mi>f</mi></mrow></msub></mrow>]]></math><img file="FDA0001179373240000031.GIF" wi="555" he="142" /></maths>式中,M<sub>ed</sub>为等效驱动力矩,N·m;M<sub>ef</sub>为等效阻力矩,N·m;J<sub>e</sub>系统转动惯量,kg·m<sup>2</sup>;θ为t时刻时曲柄所在位置的转角,rad;电动机产生的等效驱动力矩M<sub>ed</sub>,根据电动机的机械特性曲线得出,即<maths num="0008"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>e</mi><mi>d</mi></mrow></msub><mo>=</mo><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><msub><mi>&eta;</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><msub><mi>M</mi><mi>M</mi></msub><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001179373240000032.GIF" wi="469" he="86" /></maths>根据电动机额定参数由下式计算:<maths num="0009"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>e</mi><mi>d</mi></mrow></msub><mo>=</mo><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><msub><mi>&eta;</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mfrac><mrow><mn>2</mn><msub><mi>&lambda;</mi><mi>k</mi></msub><msub><mi>M</mi><mi>H</mi></msub><msub><mi>S</mi><mi>m</mi></msub><msub><mi>&omega;</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mo>-</mo><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow></mrow><mrow><msubsup><mi>S</mi><mi>m</mi><mn>2</mn></msubsup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>&omega;</mi><mn>0</mn></msub><mo>-</mo><msub><mi>i</mi><mrow><mi>M</mi><mi>B</mi></mrow></msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0001179373240000033.GIF" wi="772" he="173" /></maths><maths num="0010"><math><![CDATA[<mrow><msub><mi>S</mi><mi>m</mi></msub><mo>=</mo><mrow><mo>(</mo><mo>(</mo><mrow><msub><mi>n</mi><mn>0</mn></msub><mo>-</mo><msub><mi>n</mi><mi>H</mi></msub></mrow><mo>)</mo><mo>/</mo><msub><mi>n</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>k</mi></msub><mo>+</mo><msqrt><mrow><msubsup><mi>&lambda;</mi><mi>k</mi><mn>2</mn></msubsup><mo>-</mo><mn>1</mn></mrow></msqrt><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001179373240000034.GIF" wi="726" he="84" /></maths>ω<sub>0</sub>=2πn<sub>0</sub>/60式中i<sub>MB</sub>为电动机输出轴到减速箱的传动比;η<sub>MB</sub>为电动机输出轴到减速箱的传动效率;M<sub>M</sub>为电动机在转子角速度为<img file="FDA0001179373240000035.GIF" wi="107" he="87" />的输出扭矩,N·m;λ<sub>k</sub>为电动机最大扭矩与额定扭矩的比;M<sub>H</sub>为额定扭矩,N·m;P<sub>H</sub>为额定功率,kW;n<sub>0</sub>为同步转速,min<sup>‑1</sup>;n<sub>H</sub>为额定转速,min<sup>‑1</sup>,电动机输出扭矩由皮带传递给减速箱得到的扭矩为抽油机减速箱输出净扭矩M<sub>ed</sub>,则减速箱输出轴净扭矩M<sub>ed</sub>的计算公式:M<sub>ed</sub>=T<sub>e</sub>i<sub>p</sub>η<sub>p</sub>i<sub>j</sub>η<sub>j</sub>式中,T<sub>e</sub>为电动机输出扭矩,N·m;η<sub>p</sub>为皮带的传动效率;η<sub>j</sub>为减速箱的传动效率;i<sub>p</sub>为皮带的传动比;i<sub>j</sub>减速箱的传动比;等效阻力矩由下式计算:<maths num="0011"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>e</mi><mi>f</mi></mrow></msub><mo>=</mo><mover><mrow><mi>T</mi><mi>F</mi></mrow><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>P</mi><mi>R</mi><mi>L</mi><mo>-</mo><msub><mi>B</mi><mi>W</mi></msub><mo>)</mo></mrow><msubsup><mi>&eta;</mi><mrow><mi>C</mi><mi>L</mi></mrow><msub><mi>k</mi><mn>1</mn></msub></msubsup><mo>-</mo><msub><mi>M</mi><mi>C</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>0</mn></msub><mo>+</mo><mi>&theta;</mi><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001179373240000036.GIF" wi="910" he="78" /></maths>式中<img file="FDA0001179373240000041.GIF" wi="62" he="63" />为游梁式抽油机扭矩因数,m;PRL为抽油机悬点载荷,N;B<sub>W</sub>为抽油机平衡配重,N;M<sub>C</sub>为曲柄平衡力矩,N·m;τ为平衡重的滞后角,rad;η<sub>CL</sub>为曲柄到悬点的传动效率;k<sub>1</sub>为系数,与悬点速度相关;等效转动惯量J<sub>e</sub>的计算:将由抽油机悬点到电动机转子的所有运动部件的质量和转动惯量依据动能相等的原则转化到减速器的输出轴上,则<maths num="0012"><math><![CDATA[<mrow><msub><mi>J</mi><mi>e</mi></msub><mo>=</mo><msub><mi>&Sigma;J</mi><mi>i</mi></msub><msup><mrow><mo>&lsqb;</mo><mfrac><msub><mi>&omega;</mi><mi>i</mi></msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover></mfrac><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&Sigma;m</mi><mi>i</mi></msub><msup><mrow><mo>&lsqb;</mo><mfrac><msub><mi>v</mi><mi>i</mi></msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover></mfrac><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0001179373240000042.GIF" wi="597" he="190" /></maths>式中,J<sub>i</sub>为第i个部件的转动惯量,kg·m<sup>2</sup>;ω<sub>i</sub>为第i个部件转动速度,1/s;m<sub>i</sub>为第i个部件绕的质量,kg;v<sub>i</sub>为第i个部件的质心运动速度,m/s;根据抽油杆柱的波动方程,用F<sub>i,j</sub>表示节点i在j时刻的抽油杆柱的载荷,用PRL表示悬点载荷,进而根据抽油机电动机的运动规律,利用运动学可以得出悬点的速度、加速度和位移,即可以求解悬点载荷;将悬点载荷、扭矩因数、抽油机平衡配重、曲柄平衡力矩、平衡重的滞后角以及曲柄到悬点的传动效率等参数代入曲柄净扭矩计算公式即可;步骤三、游梁平衡的判定由计算得到的上、下冲程中曲柄轴净扭矩的最大值进行比作为平衡度,即上、下冲程中的扭矩峰值比,如下:<maths num="0013"><math><![CDATA[<mrow><mi>&epsiv;</mi><mo>=</mo><mfrac><mrow><mi>min</mi><mrow><mo>&lsqb;</mo><mrow><msub><mi>T</mi><mrow><mi>max</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>max</mi><mn>2</mn></mrow></msub></mrow><mo>&rsqb;</mo></mrow></mrow><mrow><mi>max</mi><mrow><mo>&lsqb;</mo><mrow><msub><mi>T</mi><mrow><mi>max</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>max</mi><mn>2</mn></mrow></msub></mrow><mo>&rsqb;</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001179373240000051.GIF" wi="469" he="159" /></maths>式中,T<sub>max1</sub>、T<sub>max2</sub>为上、下冲程中的最大的扭矩值。如果,ε大于0.7,则游梁式抽油机处于平衡状态;如果ε小于0.7,则需要调整平衡块位置,使其达到调整后平衡。
地址 710065 陕西省西安市电子二路东段18号