CN106485087B - 一种针对压力降法定量计算泄漏率的修正方法 - Google Patents
一种针对压力降法定量计算泄漏率的修正方法 Download PDFInfo
- Publication number
- CN106485087B CN106485087B CN201610913288.0A CN201610913288A CN106485087B CN 106485087 B CN106485087 B CN 106485087B CN 201610913288 A CN201610913288 A CN 201610913288A CN 106485087 B CN106485087 B CN 106485087B
- Authority
- CN
- China
- Prior art keywords
- leakage
- drawdown
- curve
- density
- pressure
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Abstract
本发明涉及一种针对压力降法定量计算泄漏率的修正方法,其中压力降法使用压力替代密度进行相对年泄漏率的计算,该计算基于理想气态方程中压力与密度成线性关系,所述方法包括以下步骤:先根据绝缘气体压力与密度变化关系实际曲线作出其在一段区域内的拟合曲线;判断拟合曲线与绝缘气体实际曲线的拟合程度;之后分别计算拟合曲线和绝缘气体理想曲线的相对密度降;最后计算拟合曲线和绝缘气体理想曲线的相对密度降之比,即为压力降法计算相对泄漏率的修正系数。与现有技术相比,本发明具有使压力降法定量计算泄漏率准确率提升等优点。
Description
技术领域
本发明涉及一种泄漏率计算方法的修正方法,尤其是涉及一种针对压力降法定量计算泄漏率的修正方法。
背景技术
泄漏率是气体绝缘设备的重要监督和管理指标,目前相关各类技术标准都规定每个气室或隔室的相对年泄漏率必须小于0.5%。常用定量泄漏检漏技术对气室进行相对年泄漏率的计算,定量计算泄漏率的试验方法包括:扣罩法、挂瓶法、局部包扎法、压力降法等。其中,压力降法无需仪器设备,方法简便直观,适用于电气设备在运行期间进行泄漏率的计算。然而,压力降法的定义及其存在的近似计算误差还不为电力行业的运行和检修人员所熟悉,如何减小压力降法计算带来的误差,使计算所得泄漏率更接近实际值是当前技术人员亟需解决的问题。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种针对压力降法定量计算泄漏率的修正方法。该修正方法提高了压力降法在泄漏率检测中的计算精度和判断准确性。
本发明的目的可以通过以下技术方案来实现:
一种针对压力降法定量计算泄漏率的修正方法,其中压力降法使用压力替代密度进行相对年泄漏率的计算,该计算基于理想气态方程中压力与密度成线性关系,
所述方法包括以下步骤:
S1、根据绝缘气体压力与密度变化关系实际曲线作出其在一段区域内的拟合曲线;
S2、判断拟合曲线与绝缘气体实际曲线的拟合程度,若拟合程度符合要求则进行步骤S3,否则重新拟合;
S3、分别计算拟合曲线和绝缘气体理想曲线的相对密度降;
S4、计算拟合曲线和绝缘气体理想曲线的相对密度降之比,即为压力降法计算相对泄漏率的修正系数。
所述的绝缘气体为六氟化硫气体。
所述泄漏率计算公式为
式中m1为泄漏前的质量,单位为kg;
m2为泄漏后的质量,单位为kg;
Δt为漏气至m2的时间间隔,单位为月。
对于同一个气室或隔室体积为一定时,泄漏率计算公式变为
式中ρ1为泄漏前的密度,单位为kg/m3;
ρ2为泄漏后的密度,单位为kg/m3;
Δt为漏气至ρ2的时间间隔,单位为月。
所述压力降法定量计算泄漏率公式为
式中,p1为泄漏前的相对压力,单位为MPa;
p2泄漏后的相对压力,单位为MPa;
Δt为漏气至p2的时间间隔,单位为月。
所述的步骤S1中采用取点拟合法作出拟合曲线。
所述的步骤S2中判断拟合程度大于98%即可进行步骤S3,否则需重新拟合。
步骤S3中相对密度降计算公式为x1为泄漏前的密度值,x2为泄漏后的密度值。
该方法得到的相对密度降的比值为p1为泄漏前的相对压力,B为拟合曲线在压力轴上的截距。
与现有技术相比,本发明提出了针对压力降法定量计算泄漏率的修正方法,通过拟合曲线,得到了修正系数的计算公式;使得压力降法计算泄漏率既具有无需仪器设备,方法简便直观,适用于电气设备在运行期间进行泄漏率的计算的优势的同时,准确率又有所提升。
附图说明
图1为理想气体与六氟化硫气体的压力与密度变化关系图(t=20℃);
图2为0.4-0.8MPa之间理想气体与六氟化硫气体的压力与密度变化关系图(t=20℃);
图3为拟合曲线与六氟化硫气体的压力与密度变化关系图(t=20℃)。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
实施例
泄漏率的定义
相对年泄漏率为一定时间间隔内测定的漏气量与总充气量之比,以年漏气百分率表示,值得注意的是该定义中所指漏气量和总充气量都应为质量,根据上述定义,相对年泄漏率的计算公式可表示为式(1)
式中m1为泄漏前的质量,单位为kg
m2为泄漏后的质量,单位为kg
Δt为漏气至m2的时间间隔,单位为月
对于同一个气室或隔室体积为一定时,式(1)可变形为式(2)
式中ρ1为泄漏前的密度,单位为kg/m3;
ρ2为泄漏后的密度,单位为kg/m3;
Δt为漏气至ρ2的时间间隔,单位为月;
压力降法
GB11023-1989高压开关设备六氟化硫气体密封试验方法中给出了使用压力降法计算相对年泄漏率的公式,如式(3)
式中:p1为泄漏前的相对压力(表压),单位为MPa;
p2泄漏后的相对压力(表压),单位为MPa;
Δt为漏气至p2的时间间隔,单位为月;
误差分析
比较式(2)和式(3)可以发现,压力降法使用压力替代密度进行相对年泄漏率的计算。这一近似计算的依据是由于理想气态方程中压力与密度成线性关系。然而,当气压增大时,六氟化硫气体并非理想气体,工程上常采用下列经验公式(4)计算六氟化硫气体的状态参数。
式中:p——六氟化硫气体的绝对压力,单位为MPa;
ρ——六氟化硫气体的密度,单位为kg/m3;
T——六氟化硫气体的温度,单位为K;
根据六氟化硫的理想气体方程和式(4),使用matlab将理想气体和六氟化硫气体的状态参数之间的关系绘于图1,并将压力在0.4-0.8MPa(GIS等开关设备常用压力范围)之间的区域放大至图2,结合图1和图2可知,在压力降一定的情况下,使用压力降替代质量降或密度降计算相对年泄漏率,将使计算结果比实际值偏小。为了得到更加精确的泄漏率,为设备提供可靠的评价数据,应对压力降法的计算结果进行修正。
修正方法
从上图2中可以发现,选取的GIS等开关设备的压力范围(0.4-0.8MPa)内,六氟化硫气体的状态方程近似为直线,选取点(25,0.4042),(30,0.4802),(35,0.5545),(40,0.6273),(45,0.6985),(50,0.7681),(55,0.8362),使用matlab将该直线方程拟合,得到直线方程为y=0.01440*x+0.04820,相关系数达到0.9996,拟合程度较好,该直线代表的压力与密度的关系如图3所示,结合拟合程度与图2可说明,拟合曲线与六氟化硫气体的状态参数曲线非常吻合。
对于相同的压力降,即已知泄漏发生前的压力p1和泄漏发生后的压力p2,使用拟合曲线方程P=0.01440*ρ+0.04820得到的相对密度降和六氟化硫气体理想方程P=0.0167*ρ得到的相对密度降之比,即为压力降法计算相对年泄漏率的修正系数。通过该修正方法得到的相对密度降的比值为
从得到的修正系数可以看出,修正系数仅与泄漏发生前的压力有关,针对GIS等开关设备的一般压力区间0.4-0.8MPa,将常用的系数列于表1。从表1中的数据可以看出,当电气设备泄漏发生前的压力越小,使用压力降法计算相对年泄漏率的误差越大。尤其值得注意的是,对于设备表压处于0.4-0.5MPa之间时,压力降法的相关系数为1.1左右,即方法误差达到10%左右。因此,在使用压力降法计算泄漏率时应稍加注意方法本身的误差。当然,对于处于其他压力下的设备,可以参考本文4中提出的修正方法,得到修正公式或系数。
表1 GIS设备常用压力的修正系数表
泄漏前的绝对压力,MPa | 0.40 | 0.45 | 0.50 | 0.55 | 0.60 | 0.65 | 0.70 | 0.75 | 0.80 |
压力降法的修正系数 | 1.14 | 1.12 | 1.10 | 1.09 | 1.09 | 1.08 | 1.07 | 1.07 | 1.06 |
修正系数仅与泄漏前的压力有关,设备泄漏发生前的压力越小,所得修正系数越大,压力降法计算的相对年泄漏率误差越大。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (4)
1.一种针对压力降法定量计算泄漏率的修正方法,其中压力降法使用压力替代密度进行相对年泄漏率的计算,该计算基于理想气态方程中压力与密度成线性关系,其特征在于,
所述方法包括以下步骤:
S1、根据绝缘气体压力与密度变化关系实际曲线采用取点拟合法作出其在一段区域内的拟合曲线,所述的绝缘气体为六氟化硫气体;
S2、判断拟合曲线与绝缘气体实际曲线的拟合程度,若拟合程度大于98%则进行步骤S3,否则重新拟合;
S3、分别计算拟合曲线和绝缘气体理想曲线的相对密度降,相对密度降计算公式为x1为泄漏前的密度值,x2为泄漏后的密度值;
S4、计算拟合曲线和绝缘气体理想曲线的相对密度降之比,即为压力降法计算相对泄漏率的修正系数,相对密度降的比值为p1为泄漏前的相对压力,B为拟合曲线在压力轴上的截距。
2.根据权利要求1所述的一种针对压力降法定量计算泄漏率的修正方法,其特征在于,
所述泄漏率计算公式为
式中m1为泄漏前的质量,单位为kg;
m2为泄漏后的质量,单位为kg;
Δt为漏气至m2的时间间隔,单位为月。
3.根据权利要求2所述的一种针对压力降法定量计算泄漏率的修正方法,其特征在于,对于同一个气室或隔室体积为一定时,泄漏率计算公式变为
式中ρ1为泄漏前的密度,单位为kg/m3;
ρ2为泄漏后的密度,单位为kg/m3;
Δt为漏气至ρ2的时间间隔,单位为月。
4.根据权利要求1所述的一种针对压力降法定量计算泄漏率的修正方法,其特征在于,
所述压力降法定量计算泄漏率公式为
式中,p1为泄漏前的相对压力,单位为MPa;
p2泄漏后的相对压力,单位为MPa;
Δt为漏气至p2的时间间隔,单位为月。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610913288.0A CN106485087B (zh) | 2016-10-20 | 2016-10-20 | 一种针对压力降法定量计算泄漏率的修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610913288.0A CN106485087B (zh) | 2016-10-20 | 2016-10-20 | 一种针对压力降法定量计算泄漏率的修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106485087A CN106485087A (zh) | 2017-03-08 |
CN106485087B true CN106485087B (zh) | 2019-03-22 |
Family
ID=58270250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610913288.0A Active CN106485087B (zh) | 2016-10-20 | 2016-10-20 | 一种针对压力降法定量计算泄漏率的修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106485087B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107835236A (zh) * | 2017-10-30 | 2018-03-23 | 中国电力科学研究院有限公司 | 一种充气型电力设备气压远程监控方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102042945A (zh) * | 2010-11-03 | 2011-05-04 | 北京航空航天大学 | 一种测量密闭式齿轮箱油雾浓度的方法 |
CN102420021A (zh) * | 2011-11-15 | 2012-04-18 | 中广核工程有限公司 | 一种核电站安全壳泄漏率监测方法及系统 |
CN102735408A (zh) * | 2012-02-28 | 2012-10-17 | 宝山钢铁股份有限公司 | 一种基于标准压力变化曲线的密闭容器密封性能检测方法 |
CN104077443A (zh) * | 2014-06-25 | 2014-10-01 | 重庆科技学院 | 一种气液联动球阀最大压降速率设置方法 |
CN104131811A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气股份有限公司 | 一种气井标况下体积泄漏速率获取方法及装置 |
CN105608331A (zh) * | 2016-01-26 | 2016-05-25 | 上海燃气工程设计研究有限公司 | 一种随高程变化的附加压头的计算方法 |
-
2016
- 2016-10-20 CN CN201610913288.0A patent/CN106485087B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102042945A (zh) * | 2010-11-03 | 2011-05-04 | 北京航空航天大学 | 一种测量密闭式齿轮箱油雾浓度的方法 |
CN102420021A (zh) * | 2011-11-15 | 2012-04-18 | 中广核工程有限公司 | 一种核电站安全壳泄漏率监测方法及系统 |
CN102735408A (zh) * | 2012-02-28 | 2012-10-17 | 宝山钢铁股份有限公司 | 一种基于标准压力变化曲线的密闭容器密封性能检测方法 |
CN104077443A (zh) * | 2014-06-25 | 2014-10-01 | 重庆科技学院 | 一种气液联动球阀最大压降速率设置方法 |
CN104131811A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气股份有限公司 | 一种气井标况下体积泄漏速率获取方法及装置 |
CN105608331A (zh) * | 2016-01-26 | 2016-05-25 | 上海燃气工程设计研究有限公司 | 一种随高程变化的附加压头的计算方法 |
Non-Patent Citations (4)
Title |
---|
安全壳泄漏率在线监测系统原理及数据分析;章春伟等;《核安全》;20140630;第13卷(第2期);全文 |
浮环密封流场数值模拟及泄漏率公式的修正;陈志等;《四川大学学报(工程科学版)》;20160131;第48卷(第1期);第2.1-3.3节 |
破前漏技术中气体的泄漏率;张征明等;《清华大学学报(自然科学版)》;20081231;第48卷(第9期);全文 |
高压燃气管道破裂的定量风险分析;赖建波等;《天津大学学报》;20070531;第40卷(第5期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106485087A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2856365B1 (en) | System and method to generate molecular formula distributions beyond a predetermined threshold for a petroleum stream | |
CN105628511B (zh) | 一种高温合金蠕变寿命预测方法 | |
CN102184292B (zh) | 服从指数分布的电子产品可靠性预计模型修正方法 | |
CN102564696A (zh) | 便携式真空规校准系统及方法 | |
CN102589803A (zh) | 一种便携式多功能真空校准系统及方法 | |
CN102735408A (zh) | 一种基于标准压力变化曲线的密闭容器密封性能检测方法 | |
CN106485087B (zh) | 一种针对压力降法定量计算泄漏率的修正方法 | |
CN103399203A (zh) | 一种基于复合迭代算法的谐波参数高精度估计方法 | |
CN104133404A (zh) | 一种信号处理方法及装置 | |
CN110261698A (zh) | 基于mmc换流阀应用工况的金属化膜电容器可靠性测评方法 | |
NO20180810A1 (en) | System and method for identifying and recovering from a temporary sensor failure | |
Tian et al. | Fault diagnosis of gas pressure regulators based on CEEMDAN and feature clustering | |
CN117195608B (zh) | 任意应力分布下裂纹最深点处应力强度因子的计算方法 | |
CN113570475A (zh) | 智慧管网能量计量数据跨网共享方法及系统 | |
CN202442842U (zh) | 一种便携式真空规校准装置 | |
CN117236069B (zh) | 任意应力分布下裂纹自由表面处应力强度因子的计算方法 | |
CN109358276A (zh) | 空气净距计算方法及装置 | |
CN110361648B (zh) | 一种六氟化硫气体密度继电器校验方法及系统 | |
CN109307158B (zh) | 一种确定管道线路泄漏的方法和装置 | |
CN108197342B (zh) | 一种航空发动机空气系统非设计性漏气计算方法 | |
CN105548882B (zh) | 基于线性凝聚度的发电机空载特性线性段辨识方法 | |
CN208350586U (zh) | 一种水汽渗透性能测试系统 | |
CN104408540B (zh) | 一种汽轮机改造节能量的评估方法 | |
CN113421674B (zh) | 高温堆冷试期间一回路泄漏率计算方法及装置 | |
CN110531608A (zh) | 基于冗余设计的高可靠电子设备定量fmeca分析方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |