一种压缩因子的获取方法及系统的制作方法
【专利摘要】本发明公开了一种压缩因子的获取方法及系统，首先获得预设固定体积V的待测工质，在处于不同的预设固定质量时的(Ti，pi)状态点下对应的所述待测工质的实验压缩因子Zi，即通过实验方法获得离散压缩因子，然后根据离散压缩因子构建训练样本集，最后根据训练样本集和LS?SVM方法构建得到压缩因子预测模型，该压缩因子预测模型即所需求解的状态点(T，p)与求解对应的压缩因子Z之间的关系式。本发明首先通过实验方法获得离散压缩因子，然后通过对离散压缩因子进行模型预测得到压缩因子预测模型，根据该压缩因子预测模型即可获得连续压缩因子。因此，本发明可完成任意工质在不同温度、不同压力下压缩因子的获取。
【专利说明】
一种压缩因子的获取方法及系统
技术领域
[0001] 本发明涉及压缩因子技术领域，更具体的说，涉及一种压缩因子的获取方法及系 统。
【背景技术】
[0002] 压缩因子(compressibility factor，简记为Z)是工质重要的热物性参数之一，常 被用来计算其他物性参数，由于压缩因子Z随温度T、压力p的变化而变化，因此，在进行工程 计算时，需要获得不同（T，p)状态下的压缩因子Z。
[0003] 采用对应态定律求解压缩因子是目前常用的方法。该方法具体为:根据待测工质 的温度Τ、压力ρ以及对应的临界温度Τ。和临界压力ρ。，得到对比温度T r和对比压力pr，然后 根据对应态定律，得到压缩因子Z与对比温度Tr、对比压力p r的关系式Z = Zc(HprJr),最后 通过Z对于Tr、p r的关系图一一通用压缩因子图，查询对应的压缩因子Z值，其中，Z。为临界压 缩因子，其取值依据主观经验确定，通常为〇. 27。
[0004] 可以看出，采用对应态定律求解压缩因子Z属于经验方法，对于通用压缩因子图的 高温高压区域，不同物质会不同程度的偏离曲线，且临界压缩因子Z。的取值本身也是主要 经验确定的，因此，获得的压缩因子Z存在较大误差。与此同时，对于需要事先已知的待测工 质的临界温度P。和临界温度T。，在实际操作中，对于新型工质或混合工质而言，临界温度ρ。 和临界温度T。的数值往往是缺失的。因此，如何提供一种压缩因子的获取方法及系统，以实 现对任意工质在不同温度、不同压力下压缩因子的获取是本领域技术人员亟待解决的技术 问题。
【发明内容】
[0005] 有鉴于此，本发明公开了一种压缩因子的获取方法及系统，以对任意工质在不同 温度、不同压力下压缩因子的获取。
[0006] -种压缩因子的获取方法，包括：
[0007] 获得预设固定体积V的待测工质，在处于不同的预设固定质量时的（T1,Pl)状态点 下对应的所述待测工质的实验压缩因子Z 1，其中，T1为第i个温度，Pl为第i个压力，I < i <n， η为正整数；
[0008] 根据各个状态点下对应的各个实验压缩因子，以及公式（1)，构建压缩因子预测模 型的训练样本集S;
[0009] S={((Tl,Pl),Zl),((T2,P2),Z2),---,((Tn,Pn),Zn)} (D；
[0010] 对所述训练样本集S采用最小二乘支持向量机LS-SVM方法构建得到公式(2)所表 示的所述压缩因子预测模型；
[0011；
[0012]九屮，Z 73 /7Γ边ff测丄谀2伯囚子，a i为拉格朗日乘子，b为偏置量，
η，Zj为第j个(Tj，pj)状态点下对应的压缩因子，Tj为第j个状态点对应的温度，pj为第j个状 态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（T，p)为所需求解的状 态点。
[0013] 优选的，所述获得预设固定体积V的待测工质，在处于不同的预设固定质量时的 (T1, P1)状态点下对应的所述待测工质的实验压缩因子Z1的步骤在压缩因子测量装置中执 行；
[0014] 所述压缩因子测量装置包括:压力容器、温度传感器、压力传感器、加热装置、温控 装置、真空栗和计算机;所述加热装置设置在所述压力容器的外部且与所述温控装置连接， 所述真空栗与所述压力容器通过阀门连接，所述压力容器分别通过所述温度传感器、所述 压力传感器与所述计算机连接，其中，所述压力容器设置有加注口；
[0015] 所述步骤包括：
[0016] 将所述待测工质通过所述加注口加注进容积为V的所述压力容器，直至所述压力 容器内的所述待测工质的质量为第一预设固定质量m，所述第一预设固定质量m为所述不同 的预设固定质量中的一种；
[0017] 打开所述阀门，利用所述真空栗将所述压力容器内抽真空；
[0018]当对所述压力容器真空抽取完毕后，关闭所述阀门；
[0019] 通过所述温控装置对所述压力容器进行加热，使所述压力容器内的所述待测工质 达到过热态；
[0020] 在所述待测工质达到过热态后，利用所述温度传感器和所述压力传感器测量所述 压力容器在继续加热过程中nl个温度和压力的对应变化量(Tnp 1),并输出至所述计算机记 录保存；
[0021] 将所述待测工质的相对分子量M、通用气体常数R、所述第一预设固定质量m、所述 预设固定体积V、温度T1和压力 ?1带入公式(3)所示的实际气体状态方程，计算得到nl个(T1, P1)状态点下对应的所述待测工质的实验压缩因子Z1;
[0022]
[0023] 改变加注进所述压力容器内所述待测工质的质量，再次得到n2个(T1, Pl)状态点下 对应的所述待测工质的实验压缩因子21，如此反复，直至得到η个(T 1, Pl)状态点下对应的所 述待测工质的实验压缩因子Zh
[0024] 优选的，所述对所述训练样本集S采用最小二乘支持向量机LS-SVM方法构建得到 公式(2)所表示的所述压缩因子预测模型包括：
[0025] 将所述训练样本集S代入公式(4)，构造带有约束的优化问题；
[0026]
[0027]
[0028]
[0029] 式中，J( ω，e，b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，γ为惩罚系 数，史(.〇是核空间映射函数；
[0030] 所述核空间映射函数满足公式（5):
[0031]
[0032]式中，K(Xl，Xj)为核函数；
[0033]根据拉格朗日乘子α，将所述优化问题转化为对应的无约束优化函数L，所述无约 束优化函数L的表达式参见公式(6);
[0034]
[0035] 根据所述无约束优化函数L求解优化问题，得到最优拉格朗日乘子α = (α?，α?，...， an)7 ;
[0036] 根据所述最优拉格朗日乘子a和所述训练样本集S、所述核空间映射函数P(X)以及 所述核函数，得到权向量ω的表达式参见公式(7)，并通过公式(8)计算得到偏置量b的值；
[0039] 将所述公式(7)和所述公式(8)代入预设置函数，得到压缩因子Z与温度T和压力p 的预测模型，即所述压缩因子的预测模型，所述预设置函数的表达式参见公式(9);
[0037]
[0038]
[0040]
[0041] 优选的，所述核函数K(Xl，Xj)为径向基核函数，所述径向基核函数的表达式参见公 式(10);
[0042]
[0043] 式中，Xi=(Ti，pi)，Xj=(Tj，pj)，〇为核函数宽度。
[0044] -种压缩因子的获取系统，包括：
[0045] 获得单元，用于获得预设固定体积V的待测工质，在处于不同的预设固定质量时的 (T1,Pl)状态点下对应的所述待测工质的实验压缩因子Z 1，其中，T1为第i个温度，Pl为第i个 压力，1<1<11，11为正整数；
[0046] 第一构建单元，用于根据各个状态点下对应的各个实验压缩因子，以及公式（1)， 构建压缩因子预测模型的训练样本集S;
[0047] S={((Tl,Pl),Zl),((T2,P2),Z2),---,((Tn,Pn),Zn)} (1);
[0048] 第二构建单元，用于对所述训练样本集S采用最小二乘支持向量机LS-SVM方法构 建得到公式(2)所表示的所述压缩因子预测模型；
< j < η，Zj为第j个(Tj，pj)状态点下对应的压缩因子，Tj为第j个状态点对应的温度，pj为第j 个状态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（T，p)为所需求解 的状态点。
[0051] 优选的，所述获得单元在压缩因子测量装置中执行；
[0052] 所述压缩因子测量装置包括:压力容器、温度传感器、压力传感器、加热装置、温控 装置、真空栗和计算机;所述加热装置设置在所述压力容器的外部且与所述温控装置连接， 所述真空栗与所述压力容器通过阀门连接，所述压力容器分别通过所述温度传感器、所述 压力传感器与所述计算机连接，其中，所述压力容器设置有加注口；
[0053] 所述获得单元包括：
[0054] 加注子单元，用于将所述待测工质通过所述加注口加注进容积为V的所述压力容 器，直至所述压力容器内的所述待测工质的质量为第一预设固定质量m，所述第一预设固定 质量m为所述不同的预设固定质量中的一种；
[0055] 抽真空子单元，用于打开所述阀门，利用所述真空栗将所述压力容器内抽真空；
[0056] 阀门关闭子单元，用于当对所述压力容器真空抽取完毕后，关闭所述阀门；
[0057] 加热子单元，用于通过所述温控装置对所述压力容器进行加热，使所述压力容器 内的所述待测工质达到过热态；
[0058]温度和压力采集子单元，用于在所述待测工质达到过热态后，利用所述温度传感 器和所述压力传感器测量所述压力容器在继续加热过程中nl个温度和压力的对应变化量 (T1，Pl )，并输出至所述计算机记录保存；
[0059] 实验压缩因子获取子单元，用于将所述待测工质的相对分子量M、通用气体常数R、 所述第一预设固定质量m、所述预设固定体积V、温度T 1和压力?1带入公式(3)所示的实际气 体状态方程，计算得到nl个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1;
[0060]
[0061]重复子单元，用于改变加注进所述压力容器内所述待测工质的质量，再次得到n2 个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1，如此反复，直至得到η个(T1, Pl) 状态点下对应的所述待测工质的实验压缩因子Zh [0062]优选的，所述第二构建单元包括：
[0063] 优化问题构造子单元，用于将所述训练样本集S代入公式(4)，构造带有约束的优 化问题；
[0064]
[0065]
[0066] ei>0,i = l,2,...l
[0067]式中，J( ω ,e,b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，γ为惩罚系 数，供(T)是核空间映射函数；
[0068] 4太炫卓l、m_射函撒? H R:
[0069]
[0070]式中，K(Xl，Xj)为核函数；
[0071 ]无约束优化函数构建子单元，用于根据拉格朗日乘子α，将所述优化问题转化为对 应的无约束优化函数L，所述无约束优化函数L的表达式参见公式(6);
[0072]
[0073] 最优拉格朗日乘子获取子单元，用于根据所述无约束优化函数L求解优化问题，得 到最优拉格朗日乘子α = (α?，α?，. . .;
[0074] 权向量和偏置量获取子单元，用于根据所述最优拉格朗日乘子α和所述训练样本 集S、所述核空间映射函数^(X)以及所述核函数，得到权向量ω的表达式参见公式(7)，并通 过公式(8)计算得到偏置量b的值；
[0075]
[0076] 尸i
[0077] 压缩因子的预测模型获取子单元，用于将所述公式(7)和所述公式(8)代入预设置 函数，得到压缩因子Z与温度T和压力p的预测模型，即所述压缩因子的预测模型，所述预设 置函数的表达式参见公式(9);
[0078]
[0079] 优选的，所述核函数K(Xl，Xj)为径向基核函数，所述径向基核函数的表达式参见公 式(10);
[0080]
[0081] 式中，Xi=(Ti，pi)，Xj=(Tj，pj)，〇为核函数宽度。
[0082] 从上述的技术方案可以看出，本发明公开了一种压缩因子的获取方法及系统，首 先获得预设固定体积V的待测工质，在处于不同的预设固定质量时的η个(T 1, Pl)状态点下对 应的所述待测工质的实验压缩因子Z1，即通过实验方法获得离散压缩因子，然后根据离散 压缩因子构建训练样本集，最后根据训练样本集和LS-SVM方法构建得到压缩因子预测模 型，该压缩因子预测模型即所需求解的状态点（Τ，ρ)与求解对应的压缩因子Z之间的关系 式。可以看出，本发明首先通过实验方法获得离散压缩因子，然后通过对离散压缩因子进行 模型预测得到压缩因子预测模型，根据该压缩因子预测模型即可获得连续压缩因子。因此， 本发明可完成任意工质在不同温度、不同压力下压缩因子的获取，从而解决了现有技术中 的难题。
【附图说明】
[0083] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 公开的附图获得其他的附图。
[0084] 图1为本发明实施例公开的一种压缩因子的获取方法流程图；
[0085] 图2为本发明实施例公开的一种压缩因子测量装置的结构示意图；
[0086] 图3为本发明实施例公开的一种压缩因子的获取系统的结构示意图。
【具体实施方式】
[0087] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0088] 本发明实施例公开了一种压缩因子的获取方法及系统，以对任意工质在不同温 度、不同压力下压缩因子的获取。
[0089] 参见图1，本发明实施例公开的一种压缩因子的获取方法流程图，该方法包括：
[0090] 步骤S11、获得预设固定体积V的待测工质，在处于不同的预设固定质量时的（T1, P1)状态点下对应的所述待测工质的实验压缩因子Z1;
[0091] 其中，Ti为第i个温度，Pi为第i个压力，η为正整数，I < i <n。
[0092] 需要说明的是，本步骤中的实验压缩因子Z1是通过实验方法测得的离散压缩因 子，每个状态点对应一个实验压缩因子。
[0093] 由实际气体状态方程
可以得到压缩因子Z的表达式
从该表达式中可以看出，测量获得待测工质的温度T、压力P、密度P(质量m、 体积V)，即可计算得到该待测工质相应的压缩因子Z。
[0094] 在高温、高压，尤其是超临界状态下，由于密封性等问题，待测工质变化的质量和 体积往往会带来测量误差。因此，本发明为避免这种误差，采用待测工质固定质量和固定体 积的方式，通过向容积V的容器内加入固定质量的待测工质，通过加热的方式使待测工质达 到过热态和超临界状态，在待测工质达到过热态后，获取容器内温度T和压力p的变化。
[0095] 步骤S12、根据各个状态点下对应的各个实验压缩因子，以及公式（1)，构建压缩因 子预测模型的训练样本集S;
[0096] S={((Tl,Pl),Zl),((T2,P2),Z2),---,((Tn,Pn),Zn)} (1);
[0097] 步骤S13、对所述训练样本集S采用最小二乘支持向量机方法构建得到公式(2)所 表示的所述压缩因子预测模型；
< j < η，Zj为第j个(Tj，pj)状态点下对应的压缩因子，Tj为第j个状态点对应的温度，pj为第j 个状态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（T，p)为所需求解 的状态点。
[0?00] 需要说明的是，支持向量机（Support Vector Machine，SVM)是与相关的学习算法 有关的监督学习模型，其建立在统计学习理论的VC维理论和结构风险最小原理基础上，根 据有限的样本信息在模型的复杂性(即对特定训练样本的学习精度)和学习能力（即无错误 地识别任意样本的能力)之间寻求最佳折中，以求获得最好的推广能力，用于分类和回归分 析。
[0101 ]最小二乘支持向量机（Least Squares Support Vector Machine，LS_SVM)是改进 的SVM，在本发明中，LS-SVM是以实验数据为基础，拟合出压缩因子Z对于压强p、温度T的预 测模型。
[0102] 可以看出，本发明通过实验方法测得离散压缩因子，相对于其它压缩因子的求解 方法，本发明采用对压缩因子进行实际测量的方法，从而避免了经验公式引起的误差，尤其 对于高温高压区域。
[0103] 另外，本发明采用预测模型方法获得连续压缩因子，从而解决了实验只能测得特 定压强P、温度T下的压缩因子Z的问题，即本发明可以得到压缩因子Z在任意压强p、温度T下 的连续值。与此同时，将LS-SVM作为预测模型方法，相对于其他预测模型方法，LS-SVM能够 更好的处理小样本问题（由于实验测得的数据有限，相对于全部压缩因子值来说，实验数据 是稀少的，在这种情况下，就需要预测方法具有良好的小样本问题处理能力)和非线性问题 (压缩因子Z随温度T和压力p而变化，压缩因子Z与温度T和压力p的函数关系是非线性的。所 以挖掘这种非线性关系是预测模型必须具备的能力）。
[0104] 综上可以看出，本发明公开了一种压缩因子的获取方法，首先获得预设固定体积V 的待测工质，在处于不同的预设固定质量时的(T1, Pl)状态点下对应的所述待测工质的实验 压缩因子21，即通过实验方法获得离散压缩因子，然后根据离散压缩因子构建训练样本集， 最后根据训练样本集和LS-SVM方法构建得到压缩因子预测模型，该压缩因子预测模型即所 需求解的状态点（T，p)与求解对应的压缩因子Z之间的关系式。可以看出，本发明首先通过 实验方法获得离散压缩因子，然后通过对离散压缩因子进行模型预测得到压缩因子预测模 型，根据该压缩因子预测模型即可获得连续压缩因子。因此，本发明可完成任意工质在不同 温度、不同压力下压缩因子的获取，从而解决了现有技术中的难题。
[0105] 通过进一步分析可知，步骤Sll中所用容积V的容器内除了待测工质外，还存在空 气，所以获得的压力是待测工质与容器内空气的总压力，并不是待测工质的实际压力。对于 理想气体，根据道尔顿分压定律可以求得"组分a"的分压力。 ρ： η
[0106] 由公式^ = =?可以得到?3 =以旭，其中，Pa为组分a的压力，P总为所有组分 P& /2总 的总压力，na为组分a的物质的量，η总为总物质的量。
[0107] 道尔顿分压定律的基础是各组分压缩因子恒定且等于1。对于实际气体，压缩因子 Z是随着温度T和压力ρ的变化而变化的，若要计算实际气体的压力，需要先得到各组分在各 个状态下的压缩因子Ζ，这显然与本发明的目的相悖。因此，本发明为排除其他组分(空气） 对待测工质压力测量的影响，在加热进行之前对容器抽真空，从而使测量得到的压力为待 测工质的实际压力。
[0108] 参见图2,本发明实施例公开的一种压缩因子测量装置的结构示意图，包括:压力 容器11、温度传感器12、压力传感器13、加热装置14、温控装置15、真空栗16和计算机17;
[0109] 其中，加热装置14设置在压力容器11的外部且与温控装置15连接，真空栗16与压 力容器11通过阀门18连接，压力容器11分别通过温度传感器12、压力传感器13与计算机17 连接，其中，压力容器11设置有加注口 111;
[0110] 从加注口 111向压力容器11内加注预设固定质量的待测工质，通过温控装置15调 节加热装置14的加热温度，使加热装置14对整个压力容器11进行加热，并利用温度传感器 12和压力传感器13测量压力容器11内的温度和压力，并将测量得到的温度和压力输出至计 算机17记录保存。
[0111] 需要说明的是，图1公开的实施例中的通过实验方法测得的离散压缩因子的过程 在图2中的压缩因子测量装置中执行，因此，为进一步优化上述实施例，步骤Sll包括：
[0112] 步骤Slll、将所述待测工质通过加注口 111加注进容积为V的压力容器11，直至压 力容器11内的所述待测工质的质量为第一预设固定质量m，所述第一预设固定质量m为不同 的预设固定质量中的一种；
[0113] 步骤S112、打开阀门18,利用真空栗16将压力容器11内抽真空；
[0114] 步骤S113、当对压力容器11真空抽取完毕后，关闭阀门111;
[0115] 步骤S114、通过温控装置15对压力容器11进行加热，使压力容器11内的所述待测 工质达到过热态；
[0116] 步骤S115、在所述待测工质达到过热态后，利用温度传感器12和压力传感器13测 量压力容器11在继续加热过程中nl个温度和压力的对应变化量(Ti， Pi)，并输出至计算机17 记录保存；
[0117] 步骤S116、将所述待测工质的相对分子量M、通用气体常数R(通常R=8.314)、所述 第一预设固定质量m、所述预设固定体积V、温度T 1和压力？1带入公式(3)所示的实际气体状 态方程，计算得到nl个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1;
[0118]
[0119] 步骤SllYj又父那汪进压力容器11内所述待测工质的质量，重复步骤Slll~步骤 S116,再次得到n2个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1，如此反复改 变待测工质的质量，直至得到η个(T 1, Pl)状态点下对应的所述待测工质的实验压缩因子Zh
[0120] 可以看出，本发明公开的测量方法及压缩因子测量装置的设计保证了测量过程中 密封性与其他参量测量的准确性，进而保证最终得到的压缩因子Z的精度。
[0121] 由于实验每次测量是在特定m、V(密度一定）下进行，得到是等密度线上温度Τ、压 力P的对应变化，如果想要得到具体的某个(T，p)状态点上对应的压缩因子Z是无法实现的， 即无法人为设定具体的某个（T，p)状态点下对应的压缩因子Z，因此，无法应用于实际。显 然，仅通过实验手段测量压缩因子Z是不能满足工程计算需要的，故提出了基于LS-SVM的Z 的预测方法，基于实验压缩因子Z1的实验数据，建立温Τ-ρ平面上的计算模型，通过该模型 便可以得到任意(T，p)状态点下对应的压缩因子Z。压缩因子Z的预测就是要建立压缩因子Z 与温度T和压力ρ的函数关系:Z = f (T，ρ)。
[0122] 因此，为进一步优化上述实施例，步骤Sl 3具体包括：
[0123] 步骤S131、将所述训练样本集S代入公式(4)，构造带有约束的优化问题；
[0124]
[0125]
[0126]
[0127]式中，J( ω，e，b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，γ为惩罚系 数，P(X)是核空间映射函数；
[0128] 所述核空间映射函数满足公式(5):
[0129]
[0130]式中，K(Xl，Xj)为核函数
[0131] 步骤S132、根据拉格朗日乘子α，将所述优化问题转化为对应的无约束优化函数L， 所述无约束优化函数L的表达式参见公式(6);
[0132]
[0133] 步骤S133、根据所述无约束优化函数L求解优化问题，得到最优拉格朗日乘子α = (αχ,αι, . . . ,?)7 ；
[0134] 步骤S134、根据所述最优拉格朗日乘子α和所述训练样本集S、所述核空间映射函 数切以)以及所述核函数，得到权向量ω的表达式参见公式(7)，并通过公式(8)计算得到偏 置量b的值；
[0135]
[0136]
[0139] 需要说明的是，上述实施例中的核函数有多种形式，一般情况下采用径向基核函数。[0140] 因此，本发明中的核函数K(Xl，Xj)优选为径向基核函数，所述径向基核函数的表达
[0137] 步骤S135、将所述公式(7)和所述公式(8)代入预设置函数，得到压缩因子Z与温度 T和压力ρ的预测模型，即所述压缩因子的预测模型，所述预设置函数的表达式参见公式 (9)；
[0138 式参见公式(10);
[0141]
[0142] 式中，Xi=(Ti，pi)，Xj=(Tj，pj)，〇为核函数宽度。
[0143] 与上述方法实施例相对应，本发明还公开了一种压缩因子的获取系统。
[0144] 参见图3,本发明实施例公开的一种压缩因子的获取系统的结构示意图，包括：
[0145] 获得单元21，用于获得预设固定体积V的待测工质，在处于不同的预设固定质量时 的(T1,Pl)状态点下对应的所述待测工质的实验压缩因子Z 1，其中，T1为第i个温度，Pl为第i 个压力，ISi <n，n为正整数；
[0146] 需要说明的是，实验压缩因子21是通过实验方法测得的离散压缩因子。
[0147] 由实际气体状态方程
可以得到压缩因子Z的表达式
从该表达式中可以看出，测量获得待测工质的温度T、压力p、密度P(质量m、 体积V)，即可计算得到该待测工质相应的压缩因子Z。
[0148] 在高温、高压，尤其是超临界状态下，由于密封性等问题，待测工质变化的质量和 体积往往会带来测量误差。因此，本发明为避免这种误差，采用待测工质固定质量和固定体 积的方式，通过向容积V的容器内加入固定质量的待测工质，通过加热的方式使待测工质达 到过热态和超临界状态，在待测工质达到过热态后，获取容器内温度T和压力p的变化。
[0149] 第一构建单元22,用于根据各个状态点下对应的各个实验压缩因子，以及公式 (1 )，构建压缩因子预测模型的训练样本集S;
[0150] S={((Tl,Pl),Zl),((T2,P2),Z2),---,((Tn,Pn),Zn)} (1)；
[0151] 第二构建单元23,用于对所述训练样本集S采用最小二乘支持向量机LS-SVM方法 构建得到公式(2)所表示的所述压缩因子预测模型；
[0152]
；
[0153]式中，Z为所述待测工质的压缩因子，αι为拉格朗日乘子，b为偏置量，
K(xi，xj)为核函数，Xi= (Ti，pi)，xj = (Tj，pj)，1 < j < η，Zj为第j个(Tj，pj)状态点下对应的压缩因子，Tj为第j个状态点对应的温度，pj为第j 个状态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（T，p)为所需求解 的状态点。
[0154] 可以看出，本发明通过实验方法测得离散压缩因子，相对于其它压缩因子的求解 过程，本发明采用对压缩因子进行实际测量，从而避免了经验公式引起的误差，尤其对于高 温高压区域。
[0155] 另外，本发明采用预测模型方法获得连续压缩因子，从而解决了实验只能测得特 定压强Ρ、温度T下的压缩因子Z的问题，即本发明可以得到压缩因子Z在任意压强ρ、温度T下 的连续值。与此同时，将LS-SVM作为预测模型方法，相对于其他预测模型方法，LS-SVM能够 更好的处理小样本问题（由于实验测得的数据有限，相对于全部压缩因子值来说，实验数据 是稀少的，在这种情况下，就需要预测方法具有良好的小样本问题处理能力)和非线性问题 (压缩因子Z随温度T和压力P而变化，压缩因子Z与温度T和压力P的函数关系是非线性的。所 以挖掘这种非线性关系是预测模型必须具备的能力）。
[0156] 综上可以看出，本发明公开了一种压缩因子的获取系统，首先获得预设固定体积V 的待测工质，在处于不同的预设固定质量时的(T1, Pl)状态点下对应的所述待测工质的实验 压缩因子21，即通过实验方法获得离散压缩因子，然后根据离散压缩因子构建训练样本集， 最后根据训练样本集和LS-SVM方法构建得到压缩因子预测模型，该压缩因子预测模型即所 需求解的状态点（T，p)与求解对应的压缩因子Z之间的关系式。可以看出，本发明首先通过 实验方法获得离散压缩因子，然后通过对离散压缩因子进行模型预测得到压缩因子预测模 型，根据该压缩因子预测模型即可获得连续压缩因子。因此，本发明可完成任意工质在不同 温度、不同压力下压缩因子的获取，从而解决了现有技术中的难题。
[0157] 通过进一步分析可知，获得单元21中所用容积V的容器内除了待测工质外，还存在 空气，所以获得的压力是待测工质与容器内空气的总压力，并不是待测工质的实际压力。对 于理想气体，根据道尔顿分压定律可以求得"组分a"的分压力。
[0158] 由公
可以得到?3 = ^ P总，其中，Pa为组分a的压力，P总为所有组分 的总压力，na为组分a的物质的量，η总为总物质的量。
[0159] 道尔顿分压定律的基础是各组分压缩因子恒定且等于1。对于实际气体，压缩因子 Z是随着温度T和压力ρ的变化而变化的，若要计算实际气体的压力，需要先得到各组分在各 个状态下的压缩因子Ζ，这显然与本发明的目的相悖。因此，本发明为排除其他组分(空气） 对待测工质压力测量的影响，在加热进行之前对容器抽真空，从而使测量得到的压力为待 测工质的实际压力。
[0160] 需要说明的是，获得单元21在图2示出的压缩因子测量装置中执行，压缩因子测量 装置的具体组成及各组成部分的工作原理请参见方法实施例对应部分。
[0161] 获得单元21包括：
[0162] 加注子单元211，用于将所述待测工质通过所述加注口加注进容积为V的所述压力 容器，直至所述压力容器内的所述待测工质的质量为第一预设固定质量m，所述第一预设固 定质量m为所述不同的预设固定质量中的一种；
[0163] 抽真空子单元212,用于打开所述阀门，利用所述真空栗将所述压力容器内抽真 空；
[0164]阀门关闭子单元213，用于当对所述压力容器真空抽取完毕后，关闭所述阀门；
[0165] 加热子单元214,用于通过所述温控装置对所述压力容器进行加热，使所述压力容 器内的所述待测工质达到过热态；
[0166] 温度和压力采集子单元215,用于在所述待测工质达到过热态后，利用所述温度传 感器和所述压力传感器测量所述压力容器在继续加热过程中nl个温度和压力的对应变化 量(Ti，pi)，并输出至所述计算机记录保存；
[0167] 实验压缩因子获取子单元216,用于将所述待测工质的相对分子量M、通用气体常 数R、所述第一预设固定质量m、所述预设固定体积V、温度T 1和压力？1带入公式(3)所示的实 际气体状态方程，计算得到nl个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1;
[0168]
[0169] 重复子单元217,用于改变加注进所述压力容器内所述待测工质的质量，再次得到 n2个(T1, Pl)状态点下对应的所述待测工质的实验压缩因子Z1，如此反复，直至得到η个(T1, P1)状态点下对应的所述待测工质的实验压缩因子
[0170] 可以看出，本发明公开的测量方法及压缩因子测量装置的设计保证了测量过程中 密封性与其他参量测量的准确性，进而保证最终得到的压缩因子Z的精度。
[0171] 由于实验每次测量是在特定m、V(密度一定）下进行，得到是等密度线上温度Τ、压 力P的对应变化，如果想要得到具体的某个(T，p)状态点上对应的压缩因子Z是无法实现的， 即无法人为设定具体的某个（Τ，ρ)状态点下对应的压缩因子Ζ，因此，无法应用于实际。显 然，仅通过实验手段测量压缩因子Z是不能满足工程计算需要的，故提出了基于LS-SVM的Z 的预测方法，基于实验压缩因子Z1的实验数据，建立温Τ-ρ平面上的计算模型，通过该模型 便可以得到任意(Τ，ρ)状态点下对应的压缩因子Ζ。压缩因子Z的预测就是要建立压缩因子Z 与温度T和压力ρ的函数关系:Z = f (T，ρ)。
[0172] 因此，为进一步优化上述实施例，第二构建单23元包括：
[0173] 优化问题构造子单元231，用于将所述训练样本集S代入公式(4)，构造带有约束的 优化问题·
[0174]
[0175]
[0176]
[0177] 式中，J( ω，e，b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，γ为惩罚系 数，糾>·)是核空间映射函数；
[0178] 所述核空间映射函数满足公式(5):
[0179]
[0180] 式中，K(Xi，Xj)为核函数；
[0181 ]无约束优化函数构建子单元232，用于根据拉格朗日乘子α，将所述优化问题转化 为对应的无约束优化函数L，所述无约束优化函数L的表达式参见公式(6);
[0182]
[0183] 最优拉格朗日乘子获取子单元233,用于根据所述无约束优化函数L求解优化问 题，得到最优拉格朗日乘子0 = ... ,?/ ;
[0184]权向量和偏置量获取子单元234,用于根据所述最优拉格朗日乘子α和所述训练样 本集S、所述核空间映射函数^(Λ〇以及所述核函数，得到权向量ω的表达式参见公式(7 )，并 通过公式(8)计算得到偏置量b的值；
[0185]
[0186]
[0187] 压缩因子的预测模型获取子单元235,用于将所述公式(7)和所述公式(8)代入预 设置函数，得到压缩因子Z与温度T和压力p的预测模型，即所述压缩因子的预测模型，所述 预设置函数的表达式参见公式(9):
[0188]
[0189] 需要说明的是，上述实施例中的核函数有多种形式，一般情况下采用径向基核函 数。
[0190] 因此，本发明中的核函数K(Xl，Xj)优选为径向基核函数，所述径向基核函数的表达 式参见公式(10);
[0191]
[0192] 後函数宽度。
[0193] 需要说明的是，系统实施例各组成部分的具体工作原理，请参见方法实施例对应 部分，本发明在此不做赘述。
[0194] 最后，还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将 一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作 之间存在任何这种实际的关系或者顺序。而且，术语"包括"、"包含"或者其任何其他变体意 在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那 些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者 设备所固有的要素。在没有更多限制的情况下，由语句"包括一个……"限定的要素，并不排 除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0195] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他 实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
[0196] 对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。 对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的 一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明 将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一 致的最宽的范围。
【主权项】
1. 一种压缩因子的获取方法，其特征在于，包括： 获得预设固定体积V的待测工质，在处于不同的预设固定质量时的（Τι，ρι)状态点下对 应的所述待测工质的实验压缩因子Zi，其中，Τι为第i个溫度，Pi为第i个压力，1含i含η，η为 正整数； 根据各个状态点下对应的各个实验压缩因子，W及公式（1)，构建压缩因子预测模型的 训练样本集S; S= {( (Τ?，ρ?)，Ζ?)，（（Τ2，Ρ2)，Ζ2)，...，（（Τη, Ρη)，Ζη)} (1); 对所述训练样本集S采用最小二乘支持向量机LS-SVM方法构建得到公式(2)所表示的 所述压缩因子预测模型；式中，Ζ为所述待测工质的压缩因子，α 1为拉格朗日乘子，b为偏置量，，K(xi，xj)为核函数，xi=(Ti，pi)，xj=(Tj，pj)，l < j < η，Zj为第j个(Tj，Pj)状态点下对应的压缩因子，Tj为第j个状态点对应的溫度，Pj为第j 个状态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（Τ，ρ)为所需求解 的状态点。2. 根据权利要求1所述的获取方法，其特征在于，所述获得预设固定体积V的待测工质， 在处于不同的预设固定质量时的（Τι，Ρι)状态点下对应的所述待测工质的实验压缩因子Zi 的步骤在压缩因子测量装置中执行； 所述压缩因子测量装置包括：压力容器、溫度传感器、压力传感器、加热装置、溫控装 置、真空累和计算机;所述加热装置设置在所述压力容器的外部且与所述溫控装置连接，所 述真空累与所述压力容器通过阀口连接，所述压力容器分别通过所述溫度传感器、所述压 力传感器与所述计算机连接，其中，所述压力容器设置有加注口； 所述步骤包括： 将所述待测工质通过所述加注口加注进容积为V的所述压力容器，直至所述压力容器 内的所述待测工质的质量为第一预设固定质量m，所述第一预设固定质量m为所述不同的预 设固定质量中的一种； 打开所述阀口，利用所述真空累将所述压力容器内抽真空； 当对所述压力容器真空抽取完毕后，关闭所述阀口； 通过所述溫控装置对所述压力容器进行加热，使所述压力容器内的所述待测工质达到 过热态； 在所述待测工质达到过热态后，利用所述溫度传感器和所述压力传感器测量所述压力 容器在继续加热过程中nl个溫度和压力的对应变化量(Τι，ρι)，并输出至所述计算机记录保 存； 将所述待测工质的相对分子量Μ、通用气体常数R、所述第一预设固定质量m、所述预设 固定体积V、溫度Τι和压力Pi带入公式(3)所示的实际气体状态方程，计算得到nl个(Τι，ρι) 状态点下对应的所述待测工质的实验压缩因子Zl;I 改变加注进所述压力容器内所述待测工质的质量，再次得到n2个(Τι，ρι)状态点下对应 的所述待测工质的实验压缩因子Zi，如此反复，直至得至Ijn个(Τι，ρι)状态点下对应的所述待 测工质的实验压缩因子Zl。3. 根据权利要求1所述的获取方法，其特征在于，所述对所述训练样本集S采用最小二 乘支持向量机LS-SVM方法构建得到公式(2)所表示的所述压缩因子预测模型包括： 将所述训练样本集S代入公式(4)，构造带有约束的优化问题；式中，J( ω，e，b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，丫为惩罚系数， 巧約是核空间映射函数； 所述核空间映射函数满足公式巧）：式中，K(xi，xj)为核函数； 根据拉格朗日乘子〇,将所述优化问题转化为对应的无约束优化函数L，所述无约束优 化函数L的表达式参见公式(6);根据所述无约束优化函数L求解优化问题，得到最优拉格朗日乘子α=(αι，αι，...， 口 η)'; 根据所述最优拉格朗日乘子α和所述训练样本集S、所述核空间映射函数餐(句W及所述 核函数，得到权向量ω的表达式参见公式(7)，并通过公式(8)计算得到偏置量b的值；将所述公式(7)和所述公式(8)代入预设置函数，得到压缩因子Z与溫度T和压力P的预 测模型，即所述压缩因子的预测模型，所述预设置函数的表达式参见公式(9);己4. 根据权利要求1所述的获取方法，其特征在于，所述核函数K(xi，xj)为径向基核函数， 所述径向基核函数的表达式参见公式(10);式中，Xi=(Ti，Pi)，Xj=(Tj，Pj)，0为核函数宽度。5. -种压缩因子的获取系统，其特征在于，包括： 获得单元，用于获得预设固定体积V的待测工质，在处于不同的预设固定质量时的（Ti， Pi)状态点下对应的所述待测工质的实验压缩因子Zl，其中，Τι为第i个溫度，Pi为第i个压 力，l<i<n，n为正整数； 第一构建单元，用于根据各个状态点下对应的各个实验压缩因子，W及公式（1)，构建 压缩因子预测模型的训练样本集S; S= {( (Τ?，ρ?)，Zl)，（（Τ2，Ρ2)，Ζ2)，...，（（Τη, Pn)，Zn)} (1); 第二构建单元，用于对所述训练样本集s采用最小二乘支持向量机LS-SVM方法构建得 到公式(2)所表示的所述压缩因子预测模型；式中，Z为所述待测工质的压缩因子，α 1为拉格朗日乘子，b为偏置量，，K(Xi，Xj)为核函数，Xi=(Ti，pi)，Xj=(Tj，pj)，l < j < η，Zj为第j个(Tj，Pj)状态点下对应的压缩因子，Tj为第j个状态点对应的溫度，Pj为第j 个状态点对应的压力，e为松弛变量，ej为第j个状态点对应的松弛变量，（Τ，ρ)为所需求解 的状态点。6. 根据权利要求5所述的获取系统，其特征在于，所述获得单元在压缩因子测量装置中 执行； 所述压缩因子测量装置包括：压力容器、溫度传感器、压力传感器、加热装置、溫控装 置、真空累和计算机;所述加热装置设置在所述压力容器的外部且与所述溫控装置连接，所 述真空累与所述压力容器通过阀口连接，所述压力容器分别通过所述溫度传感器、所述压 力传感器与所述计算机连接，其中，所述压力容器设置有加注口； 所述获得单元包括： 加注子单元，用于将所述待测工质通过所述加注口加注进容积为V的所述压力容器，直 至所述压力容器内的所述待测工质的质量为第一预设固定质量m，所述第一预设固定质量m 为所述不同的预设固定质量中的一种； 抽真空子单元，用于打开所述阀口，利用所述真空累将所述压力容器内抽真空； 阀口关闭子单元，用于当对所述压力容器真空抽取完毕后，关闭所述阀口； 加热子单元，用于通过所述溫控装置对所述压力容器进行加热，使所述压力容器内的 所述待测工质达到过热态； 溫度和压力采集子单元，用于在所述待测工质达到过热态后，利用所述溫度传感器和 所述压力传感器测量所述压力容器在继续加热过程中nl个溫度和压力的对应变化量(Ti， Pi )，并输出至所述计算机记录保存； 实验压缩因子获取子单元，用于将所述待测工质的相对分子量M、通用气体常数R、所述 第一预设固定质量m、所述预设固定体积V、溫度Τι和压力Pi带入公式(3)所示的实际气体状 态方程，计算得到nl个(Τι，ρι)状态点下对应的所述待测工质的实验压缩因子Zi;重复子单元，用于改变加注进所述压力容器内所述待测工质的质量，再次得到n2个(Ti， Pi)状态点下对应的所述待测工质的实验压缩因子Zi，如此反复，直至得到η个(Τι，ρι)状态 点下对应的所述待测工质的实验压缩因子Zi。7.根据权利要求5所述的获取系统，其特征在于，所述第二构建单元包括： 优化问题构造子单元，用于将所述训练样本集S代入公式(4)，构造带有约束的优化问 题；式中，J(w，e，b)为损失函数，ω为权向量，e为松弛变量，b为偏置量，丫为惩罚系数， 抑.Y)是核空间映射函数； 所述核空间映射函数满足公式巧）：式中，K(xi，xj)为核函数； 无约束优化函数构建子单元，用于根据拉格朗日乘子〇,将所述优化问题转化为对应的 无约束优化函数L，所述无约束优化函数L的表达式参见公式(6);最优拉格朗日乘子获取子单元，用于根据所述无约束优化函数L求解优化问题，得到最 优拉格朗日乘子口 = (口 1，口 1，...，口η)'; 权向量和偏置量获取子单元，用于根据所述最优拉格朗日乘子α和所述训练样本集S、 所述核空间映射函数辦对W及所述核函数，得到权向量ω的表达式参见公式(7)，并通过公 式(8)计算得到偏置量b的值；压缩因子的预测模型获取子单元，用于将所述公式（7)和所述公式(8)代入预设置函 数，得到压缩因子Z与溫度T和压力P的预测模型，即所述压缩因子的预测模型，所述预设置 函数的表达式参见公式(9);8.根据权利要求5所述的获取系统，其特征在于，所述核函数K(xi，xj)为径向基核函数， 所述径向基核函数的表达式参见公式(10);式中，Xi二（Ti，pi)，Xj二（Tj，pj)，曰为核函数宽度。
【文档编号】G01N33/00GK105842273SQ201610149229
【公开日】2016年8月10日
【申请日】2016年3月16日
【发明人】朱霄珣
【申请人】华北电力大学（保定）