CN106840516B - 一种基于多项式拟合的压力计温漂标定方法 - Google Patents
一种基于多项式拟合的压力计温漂标定方法 Download PDFInfo
- Publication number
- CN106840516B CN106840516B CN201710121066.XA CN201710121066A CN106840516B CN 106840516 B CN106840516 B CN 106840516B CN 201710121066 A CN201710121066 A CN 201710121066A CN 106840516 B CN106840516 B CN 106840516B
- Authority
- CN
- China
- Prior art keywords
- pressure
- temperature
- dependent variable
- formula
- calibration point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000001419 dependent effect Effects 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
Abstract
本发明涉及测量领域,特别涉及一种基于多项式拟合的压力计温漂标定方法,包括步骤:S1,建立压力计温度‑温度应变量最小二乘多项式拟合公式,多项式最高项次数由温度应变量个数决定;S2,建立压力标定点压力应变量‑温度最小二乘多项式拟合公式组,公式组包含的公式个数由压力标定点个数决定,各个公式的多项式最高项次数由温度个数决定;S3,输入温度应变量,计算当前温度和所述当前温度下各压力标定点对应的压力应变量;S4,建立当前温度下压力‑压力应变量最小二乘多项式拟合公式;S5,输入当前压力应变量,得到压力标定值。该算法相比传统抛物插值法多项式项数可调,更灵活;计算公式更少,计算效率更为高效。
Description
技术领域
本发明涉及测量领域,特别涉及一种基于多项式拟合的压力计温漂标定方法。
背景技术
在压力测量领域,尤其是油气田井下的压力测量领域,通常采用压阻式传感器或者石英传感器进行前端压力的采集。上述传感器因采用半导体器件而具有明显的温度漂移现象,所以都需要进行温度漂移校正。
目前对压力传感器的输出特性曲线,一般都是采用最小二乘法进行线性拟合,得到压力的多项式表达式,且多项式的系数一般都是固定的,例如4项或者5项,然后标定4~5组温度,得到每一组温度下面的压力输出特性曲线。对于没有标定到的温度点,就需要采用一定的方法进行插值。目前一般采用的都是抛物插值法,对两个标定的温度点见的温度和对应的压力应变量进行插值,然后计算出对应的压力多项式计算公式,用以计算此温度下的压力值。
然而,抛物线插值法插值后产生的计算公式较多,所需要的计算空间和存储空间大,计算过程需要查表获得相对应的温度值,计算效率低下。此外,抛物线插值法插值得到的压力应变量漂移曲线不光滑不连续,插值的温度间隔会影响压力温漂校正精度,不连续点处容易造成测量过程中的出现压力非点,影响测量精度。
发明内容
本发明的目的在于克服现有抛物插值法计算效率低和抛物插值得到的压力应变量漂移曲线不连续不光滑影响压力温漂校正精度,进而影响测量精度的问题,提供一种基于纵横多项式拟合的压力计温漂标定方法。
为了实现上述发明目的,本发明提供了以下技术方案:
一种基于多项式拟合的压力计温漂标定方法,包括:
S1,建立压力计温度-温度应变量最小二乘多项式拟合公式,多项式最高项次数由温度应变量个数决定;
S2,建立压力标定点压力应变量-温度最小二乘多项式拟合公式组,公式组包含的公式个数由压力标定点个数决定,各个公式的多项式最高项次数由温度个数决定;
S3,输入温度应变量,计算当前温度和所述当前温度下各压力标定点对应的压力应变量;
S4,建立当前温度下压力-压力应变量最小二乘多项式拟合公式;
S5,输入当前压力应变量,得到压力标定值。
进一步地,步骤S1中多项式最高项次数为3~5次。
进一步地,所述步骤S2中压力标定点个数为5~7个。
进一步地,所述步骤S2中多项式最高项次数为3~5次。
与现有技术相比,本发明的有益效果:
本发明通过纵横两次多项式拟合的方法实现压力计温漂标定,不限定标定时的温度组数和压力点数量,该方法可以动态根据温度组数和压力标定点个数确定压力计温度-温度应变量,压力应变量-温度多项式拟合关系,具有更灵活的特点;纵横两次多项式拟合方法进行压力温漂校正采用的计算公式相比抛物线插值法更少,计算效率更为高效;纵横两次多项式拟合相比抛物线插值法可有效避免温度插值出现的跳变或者非点、败点,有效提高温漂校正精度,进而提高压力计测量精度。
附图说明:
图1为本发明基于多项式拟合的压力计温漂标定方法流程图;
图2为0MPa时压力应变量-温度抛物插值曲线和多项式拟合曲线对比图;
图3为2MPa时压力应变量-温度抛物插值曲线和多项式拟合曲线对比图;
图4为本发明应用基于多项式拟合的压力计温漂标定方法温漂校正效果图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
一种基于多项式拟合的压力计温漂标定方法,如图1所示,包括:
S1,建立压力计温度-温度应变量最小二乘多项式拟合公式,多项式最高项次数由温度应变量个数决定;具体的,选用5组温度-温度应变量数据求解得到最小二乘条件下的多项式温度-温度应变量多项式拟合公式
T=λ(t)=a0t5+b0t4+c0t3+d0t2+et+g (1)
其中,a、b、c、d、e、g均为常数,表1中选取20℃、50℃、80℃、110℃、140℃的温度及其对应的温度应变量求解(1)中参数。
S2,建立压力标定点压力应变量-温度最小二乘多项式拟合公式组,公式组包含的公式个数由压力标定点个数决定,各个公式的多项式最高项次数由温度个数决定;
具体地,采用(1)式求得5个温度值,通过所述5个温度值及其对应的压力应变量建立压力应变量-温度多项式拟合关系。
进一步地,如表1所述,根据0MPa时的5组压力应变量-温度数据建立0MPa时的压力应变量-温度多项式拟合关系:
σ(0,T)=m0T5+n0T4+o0T3+p0T2+q0T+r0 (2)
同理,分别求解出10MPa、20MPa、30MPa、40MPa时的压力应变量-温度多项式拟合关系:
σ(10,T)=m10T5+n10T4+o10T3+p10T2+q10T+r10 (3)
σ(20,T)=m20T5+n20T4+o20T3+p20T2+q20T+r20 (4)
σ(30,T)=m30T5+n30T4+o30T3+p30T2+q30T+r30 (5)
σ(40,T)=m40T5+n40T4+o40T3+p40T2+q40T+r40 (6)
(2)~(5)即为压力标定点压力应变量-温度最小二乘多项式拟合公式组。
S3,输入温度应变量,计算当前温度和所述当前温度下各压力标定点对应的压力应变量;
具体地,实际温度应变量为t0,代入(1)式得到当前标定温度T0,T0代入(2)~(5)式求得当前温度下压力标定点压力应变量-温度最小二乘多项式拟合公式组,得到当前温度T0时标定点10MPa、20MPa、30MPa、40MPa对应的压力应变量δ(10,T0)、δ(20,T0)、δ(30,T0)、δ(40,T0)。
S4,建立当前温度下压力-压力应变量最小二乘多项式拟合公式;
具体地,根据所述标定点10MPa、20MPa、30MPa、40MPa以及对应的压力应变量δ(10,T0)、δ(20,T0)、δ(30,T0)、δ(40,T0)建立当前温度下压力-压力应变量最小而成多项式拟合公式:
F(δ)=Aδ5+Bδ4+Cδ3+Dδ2+Eδ+G (7)
S5,将当前压力应变量δ0代入(7),得到压力标定值F(δ0)。
表1纵横多项式拟合温漂标定参数表
0MPa | 10MPa | 20MPa | 30MPa | 40MPa | 温度应变量 | |
20℃ | x00 | x01 | x02 | x03 | x04 | t0 |
50℃ | x10 | x11 | x12 | x13 | x14 | t1 |
80℃ | x20 | x21 | x22 | x23 | x24 | t2 |
110℃ | x30 | x31 | x32 | x33 | x34 | t3 |
140℃ | x40 | x41 | x42 | x43 | x44 | t4 |
进一步地,图2为0MPa时压力应变量-温度抛物插值曲线和多项式拟合曲线对比图,图3为2MPa时压力应变量-温度抛物插值曲线和多项式拟合曲线对比图。从图2、图3可以清楚地看到,抛物插值的曲线存在不连续的情况,在对应的温度处,抛物插值给出的测量值会出现跳变或者非点、败点,影响测量效果。纵横多项式拟合方法给出的温漂插值数据是一直连续光滑的,给出的修正数据更平滑和准确。
图4为本发明应用基于多项式拟合的压力计温漂标定方法温漂校正效果图。图4种,坐标轴上显示有压力和温度范围。压力曲线分为2段,左段部分,温度从约25℃变化到约13℃,压力强度恒定为0MPa;右段部分,温度从13℃变化到43℃,压力强度恒定为5MPa。压力强度差异是由于外加在活塞压力计的压力引起的,与温漂算法无关。该曲线说明通过温漂校正后,在大温度范围变化的情况下,压力值稳定性非常好,说明了此方法的温漂校正效果的准确性和高精度。
Claims (4)
1.一种基于多项式拟合的压力计温漂标定方法,其特征在于,包括:
S1,建立压力计温度-温度应变量最小二乘多项式拟合公式,多项式最高项次数由温度应变量个数决定;
S2,建立压力标定点压力应变量-温度最小二乘多项式拟合公式组,公式组包含的公式个数由压力标定点个数决定,各个公式的多项式最高项次数由温度个数决定;
S3,输入温度应变量,计算当前温度和所述当前温度下各压力标定点对应的压力应变量;
S4,根据所述当前温度下各压力标定点对应的压力应变量,建立当前温度下压力-压力应变量最小二乘多项式拟合公式;
S5,输入当前压力应变量,并将输入的当前压力应变量代入当前温度下压力-压力应变量最小二乘多项式拟合公式,计算出对应温度下的压力值,即得到压力标定值。
2.根据权利要求1所述的标定方法,其特征在于,步骤S1中多项式最高项次数为3~5次。
3.根据权利要求2所述的标定方法,其特征在于,所述步骤S2中压力标定点个数为5~7个。
4.根据权利要求3所述的标定方法,其特征在于,所述步骤S2中多项式最高项次数为3~5次。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710121066.XA CN106840516B (zh) | 2017-03-02 | 2017-03-02 | 一种基于多项式拟合的压力计温漂标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710121066.XA CN106840516B (zh) | 2017-03-02 | 2017-03-02 | 一种基于多项式拟合的压力计温漂标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106840516A CN106840516A (zh) | 2017-06-13 |
CN106840516B true CN106840516B (zh) | 2019-03-15 |
Family
ID=59138070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710121066.XA Active CN106840516B (zh) | 2017-03-02 | 2017-03-02 | 一种基于多项式拟合的压力计温漂标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106840516B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108593202B (zh) * | 2018-04-28 | 2021-05-28 | 中国石油天然气集团有限公司 | 一种用于扭矩测量的标定方法系统 |
CN110553786B (zh) * | 2019-10-11 | 2021-09-24 | 北京七星华创流量计有限公司 | 压力传感器的补偿方法和系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2937727B1 (fr) * | 2008-10-24 | 2010-12-31 | Snecma | Procede pour installer et proteger un capteur sur un substrat |
CN102032974B (zh) * | 2010-11-15 | 2012-06-27 | 浙江工商大学 | 一种压力传感器温度补偿方法 |
CN102353481B (zh) * | 2011-06-30 | 2013-05-08 | 华南理工大学 | 基于二维正交函数的压力传感器温度和压力互补的方法及装置 |
CN103162901A (zh) * | 2013-03-28 | 2013-06-19 | 北京国浩传感器技术研究院(普通合伙) | 一种压力传感器的多温度点非线性校准方法 |
CN103968998A (zh) * | 2014-04-12 | 2014-08-06 | 沈阳仪表科学研究院有限公司 | 一种通过上位机对扩散硅传感器进行压力标定方法 |
CN105092145B (zh) * | 2015-05-27 | 2017-12-26 | 重庆川仪自动化股份有限公司 | 应用于智能变送器上的温压补偿线性修正方法 |
CN106248296A (zh) * | 2016-09-29 | 2016-12-21 | 胡海峰 | 压力变送器的多变元、变温标定方法 |
-
2017
- 2017-03-02 CN CN201710121066.XA patent/CN106840516B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106840516A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101839737B (zh) | 质量流量计及控制器以及质量流量计及控制器系统 | |
CN103743427B (zh) | 一种线性校准方法 | |
CN104236652B (zh) | 用于确定流动气体的至少一个气体参数的方法 | |
CN103837169B (zh) | 用于磁电编码器的自校正装置和方法以及磁电编码器 | |
CN101858811B (zh) | 高精度压力传感器信号补偿方法 | |
CN101846572B (zh) | 一种减小压力传感器基本误差的方法 | |
CN106840516B (zh) | 一种基于多项式拟合的压力计温漂标定方法 | |
CN106199483B (zh) | 一种数字万用表的非线性误差补偿方法 | |
CN106248296A (zh) | 压力变送器的多变元、变温标定方法 | |
CN104697659B (zh) | 一种采用ntc热敏电阻的温度测量装置与方法 | |
CN105527056A (zh) | 一种基于温度参考的压力补偿校准方法 | |
CN105572191A (zh) | 一种电化学气体传感器的压力补偿方法 | |
CN103884391A (zh) | 一种双反馈气体流量传感器 | |
CN104614125A (zh) | 一种压力计零位漂移解决方法 | |
CN201788042U (zh) | 数字式压力传感器 | |
CN115931178A (zh) | 一种基于bjt的智能温度传感器的校正方法 | |
RU2587644C1 (ru) | Способ измерения радиальных зазоров между торцами лопаток рабочего колеса и статорной оболочкой турбомашины | |
CN105628297A (zh) | 一种9mm压电激励振动筒压力传感器的拟合方法 | |
CN106644193B (zh) | 一种压强值的测定方法及系统 | |
CN106546267A (zh) | 一种速率陀螺的零位温度标定方法 | |
CN105301154B (zh) | 一种基于温度补偿的通用气路控制系统和方法 | |
CN105571590A (zh) | 一种融合补偿方法 | |
JP2005190461A (ja) | 多変数伝送器及び多変数伝送器の演算処理方法 | |
KR101397566B1 (ko) | 3차원 그린함수와 가중계수를 사용한 열응력 계산방법 | |
CN106371955B (zh) | 一种确定ipmi协议中线性传感器计算公式系数的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240117 Address after: No. 80 Puhe Road, Xindu Industrial East District, Chengdu City, Sichuan Province, 610500 Patentee after: Chengdu Wanji Petroleum Machinery Manufacturing Co.,Ltd. Address before: 610225, No. 24, Section 1, Xuefu Road, Southwest Economic Development Zone, Chengdu, Sichuan Patentee before: CHENGDU University OF INFORMATION TECHNOLOGY |
|
TR01 | Transfer of patent right |