CN105912505B - 一种基于幂律函数的气体状态描述方法 - Google Patents

一种基于幂律函数的气体状态描述方法 Download PDF

Info

Publication number
CN105912505B
CN105912505B CN201610273190.3A CN201610273190A CN105912505B CN 105912505 B CN105912505 B CN 105912505B CN 201610273190 A CN201610273190 A CN 201610273190A CN 105912505 B CN105912505 B CN 105912505B
Authority
CN
China
Prior art keywords
gas
power
law
state
law function
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.)
Expired - Fee Related
Application number
CN201610273190.3A
Other languages
English (en)
Other versions
CN105912505A (zh
Inventor
陈文�
梁英杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610273190.3A priority Critical patent/CN105912505B/zh
Publication of CN105912505A publication Critical patent/CN105912505A/zh
Application granted granted Critical
Publication of CN105912505B publication Critical patent/CN105912505B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明公开了一种基于幂律函数的气体状态描述方法，选定某种气体作为研究对象，确定该气体的温度作为试验条件，并获取该温度下气体的位力系数及其对应的昂内斯气体状态方程，采用幂律函数描述气体状态，估计出幂律函数中的参数值，验证该气体的幂律状态方程，并将满足幂律状态方程的实际气体归类为幂律气体。本发明用于描述气体状态的幂律状态方程仅包含幂律函数的阶数和非整数阶量纲系数两个参数，其中幂律函数的阶数为任意实数，刻画了实际气体偏离理想气体的程度，当温度一定时，幂律状态方程描述了压强和体积之间的幂律关系。该发明有广泛的工程应用前景，可用于试验数据拟合、气体状态定量分析仪器的设计、浮空气囊泄露的检测、石油伴生气日产量的测量等。

Description

一种基于幂律函数的气体状态描述方法
技术领域
本发明涉及热力学领域，具体涉及一种基于幂律函数的气体状态描述方法。
背景技术
气体状态方程是热力学工程应用的基础，广泛应用在如采矿、制造和珠宝鉴定等领域。大量研究发现，当温度较低、压力较高时，气体的行为都不同程度地偏离理想气体状态方程所假设的前提条件，特别是难液化的气体，如氮气和氦气等(见文献：BottinB.Thermodynamic properties of arbitrary perfect gas mixtures at low pressuresand high temperatures[J].Progress in Aerospace Sciences,2000,36(36):547-579.)。造成偏离的主要原因是气体分子间的引力和分子本身的体积，当压力较低时，分子间的引力起主要作用；当压力较高时，分子的体积影响较大，分子间斥力的影响不能被忽略。
低温和高压技术的快速发展要求建立更符合工程的实际气体状态方程。目前，存在有上百种气体状态方程，包括著名的范德华气体状态方程和昂内斯气体状态方程(见文献：Czepirski L,J.Virial-type thermal equation of gas-solid adsorption[J].Chemical Engineering Science,1989,44(4):797-801.)。荷兰著名物理学家昂内斯在研究氢、氦等沸点很低气体的液化时，提出的一种描述实际气体的状态方程，被称为昂内斯气体状态方程，该方程是体积或压强展开的级数。从理论上讲，昂内斯气体状态方程可以有无穷个级数展开。但在实际应用时，通常仅保留无穷级数的前4项来描述实际气体的状态。以体积展开为例，昂内斯气体状态方程的表达式为：
其中R为摩尔气体常数，P压强，V体积，T温度，B，C，D分别为第一、第二和第三位力系数，位力系数通常由试验确定。一般情况下，温度越低，位力系数的数值越大，气体越偏离理想状态。理想气体状态方程和范德华气体状态方程是昂内斯气体状态方程的特例。
这些模型的主要特点是纯经验和半经验型，参数较多，且均具有一定的适用范围。在国内外，已有多项基于理想气体状态方程如范德华气体状态方程和昂内斯气体状态方程的专利技术，如专利CN201410436242.5“一种石油伴生气日产量测量方法”，利用理想气体状态方程导出的公式，计算伴生气的日产量；CN201410436242.5“基于理想气体状态方程的刚性物体体积测量仪器及方法”，采用理想气体状态方程，测量刚性物体的体积；US2003/0045228A1Ventilation system and method，基于理想气体状态方程设计通风系统；CN201310739123.2“一种浮空气囊的泄漏检测方法”，根据理想气体状态方程建立计算泄露气体体积的方程；US1994/5307668A1Gas density meter and method，基于昂内斯气体状态方程设计气体密度测量仪；CN201310163362.8“同步跟踪测定气固吸附过程微热量变化的装置及测定方法”，根据上述三种气体状态方程，计算计量腔中气体的摩尔数。
其中，范德华气体状态方程通过两个矫正因子，分别考虑了气体分子间的引力和斥力，但未考虑矫正因子与温度的关系，仅适用于描述中低压力下气体的行为。而昂内斯气体状态方程是范德华气体状态方程的推广，其中位力系数是温度的函数，与实际气体状态的行为相吻合。由昂内斯气体状态方程可知，在恒定温度下，压强和体积的乘积是体积的幂级数。但昂内斯气体状态方程参数较多，不方便使用。
综上，理想气体状态方程和范德华气体状态方程不能够准确刻画实际气体。而现有的实际气体状态方程方法，如昂内斯气体状态方程，参数较多，实际应用程序较复杂。工程应用需要一种简单有效地确定实际气体状态的工程方法。
发明内容
发明目的：本发明的目的在于针对现有技术的不足，提供一种基于幂律函数的气体状态描述方法，通过两个参数即能准确刻画实际气体偏离理想气体的程度。
技术方案：本发明提供了一种基于幂律函数的气体状态描述方法，包括以下步骤：
(1)选定某种气体作为研究对象，确定该气体的温度作为试验条件，并获取该温度下气体的位力系数及其对应的昂内斯气体状态方程；
(2)采用体积或压强的幂律函数描述热力学量之间的关系；
(3)通过步骤(2)中的幂律函数拟合步骤(1)中的昂内斯气体状态方程，得出幂律函数中的参数值；
(4)将得到的幂律函数描绘成曲线与昂内斯气体状态方程的曲线进行比较，验证该气体的幂律状态方程，并将满足幂律状态方程的实际气体归类为幂律气体。
进一步，步骤(2)当气体的压强和温度一定时，采用体积的幂律函数描述热力学量之间的关系：
PVα＝qRT
其中，α为幂律函数的阶数，直接刻画实际气体的状态特征，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
步骤(3)通过步骤(2)中的幂律函数拟合步骤(1)中的昂内斯气体状态方程，得出幂律函数中参数α和q的值。
进一步，步骤(2)当气体的体积和温度一定时，采用压强的幂律函数描述热力学量之间的关系：
PαV＝qRT
其中，α为幂律函数的阶数，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
步骤(3)通过步骤(2)中的幂律函数拟合步骤(1)中的昂内斯气体状态方程，得出幂律函数中参数α和q的值。
进一步，步骤(3)采用最小二乘法确定幂律函数的阶数α和非整数阶量纲的系数q。
进一步，变换该气体不同的试验条件，重复步骤(2)和步骤(3)，确定幂律状态方程。
有益效果：本发明采用幂律函数确定热力学量之间的关系，然后结合实际气体的试验条件和试验数据，估计该幂律函数中的参数值，最终确定该气体的幂律方程用于描述该气体的状态。得到的幂律状态方程仅包含幂律函数的阶数和非整数阶量纲系数两个参数，其中幂律函数的阶数为任意实数，刻画了实际气体偏离理想气体的程度，阶数越小，实际气体偏离理想气体的程度越高。当温度一定时，幂律状态方程描述了压强和体积之间的幂律关系，温度越低，幂律函数的阶数越小，气体的实际状态越偏离理想气体。本发明幂律状态方程是理想气体状态方程的推广，满足幂律状态方程的实际气体为幂律气体。该发明有广泛的工程应用前景，可用于试验数据拟合、气体状态定量分析仪器的设计、浮空气囊泄露的检测、石油伴生气日产量的测量等。
附图说明
图1为本发明一种基于幂律函数的气体状态描述方法的流程图；
图2为温度为80K时，幂律气体状态方程描述氮气N2压强和体积的关系图；
图3为温度为120K时，幂律气体状态方程描述氮气N2压强和体积的关系图；
图4为温度为323K时，幂律气体状态方程描述四氟甲烷CF4压强和体积的关系图；
图5为温度为373K时，幂律气体状态方程描述四氟甲烷CF4压强和体积的关系图。
具体实施方式
下面对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
实施例：
实施例1：(1)本实施例选择氮气N2作为研究对象。表1给出了温度为80K和120K两种试验条件下，氮气N2的位力系数值：
表1 不同温度下，氮气N2的位力系数
将上述位力系数对应的昂内斯气体状态方程作为试验数据，分别为：
(2)采用体积的幂律函数描述气体状态方程：
PVα＝qRT
其中，α为幂律函数的阶数，直接刻画实际气体的状态特征，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
(3)使用步骤(2)中的幂律函数拟合上述两个昂内斯气体状态方程，分别通过最小二乘法确定步骤(2)中两个温度下幂律函数的阶数α和非整数阶量纲的系数q，表2给出了两个温度下，幂律气体状态方程参数的值。由表2可知，幂律函数的阶数均不为1，表明氮气N2在这两个温度下的状态不是理想气体，此外，温度越低，α的值越小，表明N2的状态越偏离实际气体，该结果与昂内斯气体状态方程位力系数反映的结果一致。
表2 不同温度下，氮气N2幂律气体状态方程参数
(4)两个温度即为不同的试验条件。采用表1和表2的数值，描绘幂律气体状态方程和昂内斯气体状态方程的曲线，分别见图2和图3。由图2和图3可见，两种温度下压强和体积的关系均近似为幂律函数，幂律气体状态方程与昂内斯气体状态方程的结果吻合。由此可知，氮气N2为幂律气体，其幂律气体状态方程可以用较少的两个参数准确描述氮气N2状态方程中压强和体积的幂律关系，幂律阶数刻画了氮气N2偏离理想气体的程度。
实施例2：(1)本实施例选择四氟甲烷CF4作为研究对象。表3给出了温度为323K和373K两种试验条件下，四氟甲烷CF4的位力系数值。
表3 不同温度下，四氟甲烷CF4的位力系数
将上述位力系数对应的昂内斯气体状态方程作为试验数据，分别为：
(2)采用体积的幂律函数描述气体状态方程：
PVα＝qRT
其中，α为幂律函数的阶数，直接刻画实际气体的状态特征，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
(3)使用步骤(2)中的幂律函数拟合上述两个昂内斯气体状态方程，分别通过最小二乘法确定步骤(2)中两个温度下幂律函数的阶数α和非整数阶量纲的系数q，表4给出了两个温度下，幂律气体状态方程参数的值。由表4可知，幂律函数的阶数均不为1，表明四氟甲烷CF4在这两个温度下的状态不是理想气体，此外，温度越低，α的值越小，表明四氟甲烷CF4的状态越偏离实际气体，该结果与昂内斯气体状态方程位力系数反映的结果一致。
表4 不同温度下，四氟甲烷CF4幂律气体状态方程参数的值
(4)两个温度即为不同的试验条件。采用表3和表4的数值，描绘幂律气体状态方程和昂内斯气体状态方程的曲线，分别见图4和图5。由图4和图5可见，两种温度下压强和体积的关系均近似为幂律函数，幂律气体状态方程与昂内斯气体状态方程的结果吻合。由此可知，四氟甲烷CF4为幂律气体，其幂律气体状态方程可以用较少的2参数准确描述四氟甲烷CF4状态方程中压强和体积的幂律关系，幂律阶数刻画了四氟甲烷CF4偏离理想气体的程度。
本发明采用幂律函数确定热力学量之间的关系，然后结合实际气体的试验条件和试验数据，估计该幂律气体状态方程参数的值，最终确定该气体的幂律状态方程。该发明的特点是描述气体状态的方程仅包含两个参数，即幂律函数的阶数和非整数阶量纲系数。其中幂律函数的阶数为任意实数，刻画了实际气体偏离理想气体的程度。当温度一定时，幂律状态方程描述了压强和体积之间的幂律关系，温度越低，幂律函数的阶数越小，气体的实际状态越偏离理想气体。幂律状态方程是理想气体状态方程的推广，满足幂律状态方程的实际气体为幂律气体。该发明有广泛的工程应用前景，可用于试验数据拟合、气体状态定量分析仪器的设计、浮空气囊泄露的检测、石油伴生气日产量的测量等。

Claims (4)

1.一种基于幂律函数的气体状态描述方法，其特征在于：包括以下步骤：
(1)选定某种气体作为研究对象，确定该气体的温度作为试验条件，并获取该温度下气体的位力系数及其对应的昂内斯气体状态方程；
(2)采用体积或压强的幂律函数描述热力学量之间的关系：
当气体的压强和温度一定时，采用体积的幂律函数描述热力学量之间的关系：
PVα＝qRT
其中，α为幂律函数的阶数，直接刻画实际气体的状态特征，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
当气体的体积和温度一定时，采用压强的幂律函数描述热力学量之间的关系：
PαV＝qRT
其中，α为幂律函数的阶数，q为非整数阶量纲的系数，R为摩尔气体常数，P为压强，V为体积，T为温度；
(3)通过步骤(2)中的幂律函数拟合步骤(1)中的昂内斯气体状态方程，得出幂律函数中的参数值；
(4)将得到的幂律函数描绘成曲线与昂内斯气体状态方程的曲线进行比较，验证该气体的幂律状态方程，并将满足幂律状态方程的实际气体归类为幂律气体。
2.根据权利要求1所述的基于幂律函数的气体状态描述方法，其特征在于：步骤(3)通过步骤(2)中的幂律函数拟合步骤(1)中的昂内斯气体状态方程，得出幂律函数中参数α和q的值。
3.根据权利要求1所述的基于幂律函数的气体状态描述方法，其特征在于：步骤(3)采用最小二乘法确定幂律函数的阶数α和非整数阶量纲的系数q。
4.根据权利要求1所述的基于幂律函数的气体状态描述方法，其特征在于：变换该气体不同的试验条件，重复步骤(2)和步骤(3)，确定幂律状态方程。
CN201610273190.3A 2016-04-28 2016-04-28 一种基于幂律函数的气体状态描述方法 Expired - Fee Related CN105912505B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610273190.3A CN105912505B (zh) 2016-04-28 2016-04-28 一种基于幂律函数的气体状态描述方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610273190.3A CN105912505B (zh) 2016-04-28 2016-04-28 一种基于幂律函数的气体状态描述方法

Publications (2)

Publication Number Publication Date
CN105912505A CN105912505A (zh) 2016-08-31
CN105912505B true CN105912505B (zh) 2018-05-11

Family

ID=56753053

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610273190.3A Expired - Fee Related CN105912505B (zh) 2016-04-28 2016-04-28 一种基于幂律函数的气体状态描述方法

Country Status (1)

Country Link
CN (1) CN105912505B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112086135B (zh) * 2020-08-21 2021-12-28 北京信息科技大学 状态分析方法、装置和设备
CN115600052B (zh) * 2022-12-15 2023-09-12 昆明理工大学 一种确定幂律水泥浆液稠度系数的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102693331A (zh) * 2011-03-25 2012-09-26 鸿富锦精密工业（深圳）有限公司 非线性标的的实验设计系统及方法
CN105138727A (zh) * 2015-07-23 2015-12-09 长江勘测规划设计研究有限责任公司 地下电站基于长垂直封闭母线自然热压的建模方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102693331A (zh) * 2011-03-25 2012-09-26 鸿富锦精密工业（深圳）有限公司 非线性标的的实验设计系统及方法
CN105138727A (zh) * 2015-07-23 2015-12-09 长江勘测规划设计研究有限责任公司 地下电站基于长垂直封闭母线自然热压的建模方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Fractional derivative anomalous diffusion equation modeling prime number distribution;CHEN Wen等;《Fractional Calculus ＆ Applied Analysis》;20150523;第18卷(第3期);第789-798页 *
Gas solubilities in ionic liquids using a generic van der W aals equation of state;Yokozeki A等;《Journal of Supercritical Fluids》;20101231;第55卷(第2期);第846-851页 *
T hermodynamic properties of arbitrary perfect gas mixtures at low pressures and high temperatures;Bottin B;《Progress in Aerospace Sciences》;20001031;第36卷(第7期);第547-579页 *
Virial-type thermal equation of gas-solid adsorption;C zepirski L等;《Chemical Engineering Science》;19891231;第44卷(第4期);第797-801页 *
全程压力下气体状态方程新模型;戴浩 等;《湖北第二师范学院学报》;20110820;第28卷(第8期);第10-12页 *
实际气体状态方程的讨论;聂跃光;《川北教育学院院刊》;19870615;第20-26页 *

Also Published As

Publication number Publication date
CN105912505A (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
Munkejord et al. CO2 transport: Data and models–A review
Murata et al. A new determination method of absolute adsorption isotherm of supercritical gases under high pressure with a special relevance to density-functional theory study
Camacho et al. Adsorption equilibrium of carbon dioxide and nitrogen on the MIL-53 (Al) metal organic framework
CA2753410C (en) Methods and systems for optimizing carbon dioxide sequestration operations
Nazeri et al. Measured densities and derived thermodynamic properties of CO2-rich mixtures in gas, liquid and supercritical phases from 273 K to 423 K and pressures up to 126 MPa
CN110489927B (zh) 页岩吸附气吸附相密度模型的构建及绝对吸附量计算方法
Anwar et al. Carbon Dioxide Thermodynamic Properties Handbook: Covering Temperatures from-20° to 250° C and Pressures Up to 1000 Bar
CN105912505B (zh) 一种基于幂律函数的气体状态描述方法
Farzaneh-Gord et al. Numerical procedures for natural gas accurate thermodynamic properties calculation
Sun et al. CO2-ECBM in coal nanostructure: modelling and simulation
Wareing et al. High pressure CO2 CCS pipelines: Comparing dispersion models with multiple experimental datasets
Baltadjiev et al. An investigation of real gas effects in supercritical CO2 centrifugal compressors
Huo et al. Subgrid-scale models for large-eddy simulation of supercritical combustion
Vishnyakov et al. Statistical geometry of cavities in a metastable confined fluid
Takbiri-Borujeni et al. Effect of kerogen type and maturity on performance of carbon dioxide storage in shale
Jarrahian et al. A simple correlation to estimate natural gas thermal conductivity
Rani et al. Comparison of void volume for volumetric adsorption studies on shales from India
Aakenes et al. Frictional pressure drop for two-phase flow of carbon dioxide in a tube: Comparison between models and experimental data
İbrahimoğlu et al. Determination of liquid phase range of matters by graphic-analytical method
Lype et al. Progress in International Research on Thermodynamic and Transport Properties
Yang et al. Molecular modeling of interfacial properties of the hydrogen+ water+ decane mixture in three-phase equilibrium
Deng-Feng et al. Error analyses in volumetric method for measuring methane and CO2 adsorption on coal
Kazemi et al. Flow of multicomponent gases in carbon-based organic nanopores
Yong et al. Study of CO2 fluid density calculation model based on grayscale image
Barsotti Capillary condensation in shale: A narrative review

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180511

CF01 Termination of patent right due to non-payment of annual fee