油气集输中原油处理质量特性的多变量工序能力评价.docx

1、油气集输中原油处理质量特性的多变量工序能力评价2013年油气集输大队技术管理座谈会油气集输中原油处理质量特性的多变量工序能力评价完成单位：集输大队管理办撰 稿 人：马腾祥完成时间：2013年6月21日油气集输中原油处理质量特性的多变量工序能力评价前言随着对油气集输原油含水指标标准的提高，对油气集输原油处理过程的稳定性提出了更高的要求，原油含水指标是否稳定将直接影响油气集输大队的发展及内在质量观念。过程性能指数Pp 、Ppk是QS-9000标准中提出的，能够反映系统当前的实际状态，而并不要求在过程稳态条件下进行计算。Pp和Ppk是对过程能力和管理能力的反映，因其具有不同于过程能力指数Cp 、Cp

2、k的特点，反映了实时过程的性能，利用Pp和Ppk两数值，可对桩西联合站当前的过程性能有更多的了解。在进行工序能力评价时，需对原油集输的多个质量特性进行综合评价。一、常规过程能力指数过程是一组输入转为输出的相互关联或作用的活动，具备三个条件：(1)任何过程必须是有效过程；(2)任何过程必须是稳定受控过程；(3)任何过程都不是孤立的，应形成过程网络。1、过程能力过程能力指在正常状态下，加工产品质量能够满足技术质量标准的能力。过程能力记为：B=6=6s过程能力实际是考核过程的99.73%的能力。用于衡量过程加工产品质量的内在一致性，是稳定受控状态下的最小波动，要求其数值越小越好。它取决于影响质量的因

3、素，如人、机、料、法、环等，即组织的技术水平和管理水平，与产品公差无关。2、过程能力指数过程能力指数指过程能力满足质量要求的程度，记为Cp：Cp=T/B=（TU-TL）/6=（TU-TL）/6s式中：T质量特性的公差范围；TU上公差界限；TL下公差界限；B过程能力；质量特性值总体分布的标准差；s质量特性值样本分布的标准差。3、过程能力指数的计算根据质量特性的类别不同，过程能力指数的计算有以下四种形式：(1)望目值质量特性当X=M时，Cp=T/6=(TU-TL)/6S*注：望目值质量特性当X=M时，Cp=CPU=(TU-X)/3=(TU-X)/3SCp=CPL=(X-TL)/3=(X-TL)/3

4、S(2)望目值质量特性当XM时，此时有偏移量=|X-M|相对偏移量：k=/(T/2)=2/T修正的过程能力指数为：CPK=(1-k)CP=(T-2)/6=(T-2)/6S*注：望目值质量特性当XM时，有CPK=min(CPU，CPL)(3)望小值质量特性CPU=(TU-X)/3=(TU-X)/3s(4)望大值质量特性CPL=(X-TL)/3=(X-TL)/3s4、过程能力的等级评定表1 常规过程能力等级评定表等级特级一级二级三级四级Cp1.671.331.671.001.330.671.000.67评价过高充分尚可不充分不足注：1.认为Cp=1.33最佳，不符合现代高科技对高质量的要求；2.没

5、有考虑产品质量特性重要度分级的不同情况不同对待。表2 产品过程综合能力等级评定表Cp等级特性过程能力指数范围1.671.331.671.001.330.671.00Cp1.33低风险级1.33Cp1.00中等风险级1.00Cp0.67高风险级CpCp1.33理想状态级1.33Cp1.00低风险级1.00Cp0.67中等风险级CpCp1.33能力富余级1.33Cp1.00理想状态级1.00Cp0.67低风险级Cp2m）是矩阵S的特征值，E1，E2，,Em是这些特征值所对应的特征向量。Ei就是各变量在第i个因子上的载荷。计算每个因子的贡献率：Ri=i/i(i=1，2，m) (6)则因子分析的载荷矩

6、阵L为： 1 0 0L=E1，E2，Em 0 0 0 (7) 0 0 m3、综合过程性能指数的计算要全面评价油气集输多个质量特性的综合过程能力，可通过主要素分析，构造相互独立且服从正态分布的新变量，根据各个新变量的过程性能指数及其权重综合计算过程性能指数。（1）对原始数据矩阵X进行标准化处理，消除量纲影响，得到矩阵Xst。（2）对矩阵Xst进行主要素分析，计算出载荷矩阵L，得分矩阵T，计算每个因子的贡献率，并依据累积方差贡献率90%的原则，确定主因素个数l。（3）将各质量特性规范限按照步骤（1）的方法标准化后，转化为对应各因子的规范限。（4）依据Pp和Ppk的计算公式，计算各元素的Pp 、Pp

7、k值。（5）计算综合过程性能指数： lPpr= WiPpi， i=1lPpkr= WiPpk,i i=1其中，l为主因素个数，lWi=i/ n为权重值。 i=1三、原油处理质量特性的多变量工序能力评价随机抽取桩西联合站30天的生产数据作为样本，见表6。表6 桩西联合站原油含水指标检测结果序号进站含水率%一级高效分水器%原油沉降罐%原油稳定率%交油含水率%1739.90.23820.532819.50.18790.4638210.20.16880.444799.80.27920.5757610.50.23850.5168810.30.21790.497769.70.20840.5287410.1

8、0.18810.4197710.20.20910.55108410.40.18900.5111829.90.24940.5312809.70.23890.48138310.10.19880.4314789.80.20790.52157710.40.18920.44168210.20.26930.43177910.10.19890.4718889.80.21870.5819719.90.25880.5320799.80.22920.49218110.10.17880.4622779.90.21780.52238310.00.16790.50248610.10.18840.4625799.80.2

9、7920.4526759.70.22830.56277810.20.11870.55288010.30.19900.4729819.80.24850.49308410.10.17880.48各项指标规范要求为：原油进站含水率80%；一级高效分水器10%；原油沉降罐0.1%；原油稳定率75%；交油含水率0.6%。对表6数据进行标准化处理。计算标准化后原始数据各项指标之间的相关系数，结果如表7所示。表7 标准化后原始数据各项指标之间的相关系数进站含水率一级高效分水器原油沉降罐原油稳定率一级高效分水器0.5181*原油沉降罐0.3114*0.1280原油稳定率0.19100.2859-0.2981交

10、油含水率0.1180-0.2287-0.0163-0.1728注：*表示在=0.05水平下显著相关；*表示在=0.10水平下显著相关。对标准化后数据进行主成分分析，得出载荷矩阵L：0.6034 -0.2805 0.3096 0.0378 -0.67810.6498 0.0581 -0.0643 -0.5721 0.4929L = 0.2241 -0.6612 -0.3600 0.5187 0.33750.3566 0.5901 0.2661 0.6332 0.2299-0.1904 -0.3640 0.8364 -0.0362 0.3610载荷矩阵的特征值为：1.7491、1.3779、1.0

11、133、0.5188、0.3409，各因子方差贡献率分别为：34.98%、27.56%、20.27%、10.38%、6.82%，由于前三个因子累积贡献率为93.18%90%，确定为主成分。将原油进站含水率、一级高效分水器、原油沉降罐、原油稳定率、交油含水率规范限及分布中心标准化后投影到新空间，得到主成分规范限、规格中心、分布中心，其规格限为：PC1（-8.727，7.9143）、PC2（-3.0800，3.6643）、PC3（-2.9932，3.8002）、PC4（-1.6950，2.2853）。规格中心为：-0.4064、0.2923、0.4034、0.2952。由于原始数据经过了标准化处理

12、，分布中心均为0。经检验，四个主成分均符合正态分布，由此得PC1（0，1.3226）、PC2（0，1.1739）、PC3（0，1.0066）、PC4（0，0.7202）。计算PP、PPK值。根据上述公式，求得主成分的PP、PPK值：PP，PC1=2.10，PP，PC2=0.96，PP，PC3=1.12，PP，PC4=0.92；PPK，PC1=2.04，PPK，PC2=0.91，PPK，PC3=1.05，PPK，PC4=0.85。计算综合过程性能指数PP、PPK。首先确定各主成分权重：lWi=i/ n i=1由此得到，W1=37.54%，W2=29.57%，W3=21.75%，W4=11.13%

13、。从而求得：Ppr=1.42；Ppkr=1.36。由于1.33Ppr1.67，说明过程能力充分，且Ppkr与差异不大，说明“管理能力”较好。四、结果与讨论计算原油脱水样本进站含水率、一级高效分水器、原油沉降罐、原油稳定率、交油含水率各质量指标单变量过程性能指数PP、PPK，结果如表8所示。表8 原油脱水样本各质量指标单变量过程性能指数进站含水率一级高效分水器原油沉降罐原油稳定率交油含水率PP1.581.861.691.311.22PPK1.441.761.641.181.09从表8中看出，各单变量过程性能指数计算结果差距较大，不能反映工序的综合过程能力。与之相比，综合过程性能指数充分利用各单变量的质量信息，能够较全面地反映工序的过程能力。本文将过程性能指数与主成分分析相结合，对工序能力进行综合评价，从而有效解决了由于质量特性相关性和多质量特性带来的无法有效计算过程性能指数的问题。