CN114689460B - 一种在线密度自动测量的原油密度值校正方法 - Google Patents

一种在线密度自动测量的原油密度值校正方法 Download PDF

Info

Publication number
CN114689460B
CN114689460B CN202210611918.4A CN202210611918A CN114689460B CN 114689460 B CN114689460 B CN 114689460B CN 202210611918 A CN202210611918 A CN 202210611918A CN 114689460 B CN114689460 B CN 114689460B
Authority
CN
China
Prior art keywords
density
online
crude oil
measurement
fitting
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
Application number
CN202210611918.4A
Other languages
English (en)
Other versions
CN114689460A (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.)
Second Oil Transportation Division Of Changqing Oilfield Branch Of Cnpc
Northwest University
Original Assignee
Second Oil Transportation Division Of Changqing Oilfield Branch Of Cnpc
Northwest University
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 Second Oil Transportation Division Of Changqing Oilfield Branch Of Cnpc, Northwest University filed Critical Second Oil Transportation Division Of Changqing Oilfield Branch Of Cnpc
Priority to CN202210611918.4A priority Critical patent/CN114689460B/zh
Publication of CN114689460A publication Critical patent/CN114689460A/zh
Application granted granted Critical
Publication of CN114689460B publication Critical patent/CN114689460B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N35/00—Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Mathematical Optimization (AREA)
  • Pathology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种在线密度自动测量的原油密度值校正方法，属于原油密度值校正技术领域。该方法包括：获取原油在线密度测量的第一测量时间、第一测量温度以及在线视密度；同时对原油密度进行人工测量，分别获取原油的第二测量时间、第二测量温度以及人工视密度；选取第一测量时间等于第二测量时间时的相同样品数据；分别计算得到人工标密和在线标密；并计算得到校正标密。本发明建立了二步法的在线原油密度校正方法，第一步计算是在线视密求取在线标密，进行温度校正，第二步是以在线标密为自变量进行多项式计算校正标密，进行系统误差的校正。本发明解决了在线密度仪所测视密进行标密度换算时不准确的问题。

Description

一种在线密度自动测量的原油密度值校正方法
技术领域
本发明属于原油密度值校正技术领域，尤其涉及种在线密度自动测量的原油密度值校正方法。
背景技术
原油密度是水和各种有机组份混合液在一定温度和压力下单位体积的质量，因受温度和压力影响，通常把20°C，一个大气压下的密度值称为标准密度，其他温度下测的密度值称为视密度。原油密度是原油的基本属性参数，在石油科学和产业中都具有重要作用。地层原油密度分布与地层、构造的关系，可以分析原油的赋存规律以及构造控制规律，对于认识石油的“生储盖圈运保”的成油理论具有重要作用。原油密度不仅仅可以估算石油储量，并可建立随时间动态变化规律，对于开采具有重要的指导价值。原油密度同时也是原油分类的重要依据，决定了后续原油冶炼工艺选取。
原油密度获取主要有两种方式，一种是通过其他参数间接求取原油密度，但间接求取方式受参数值准确度影响大，另一种是直接测量原油密度。现今测量原油密度主要依据油水两相介质的物理学特点进行测量，依据油水两相介质不相容的密度法，依据油水介电常数不同的电容法和射频法，依据力学理论的方法，如基于阿基米德原理的方法、液体柱压力理论的方法，依据声学传播理论的超声波测量法，依据物质对放射性衰减理论的𝛾射线方法，依据振动频率理论的方法等，方法虽然多，但测量方式无非两种，一是人工测量，二是在线自动测量，人工测量工作效率低难以满足现今石油管网密布的测量要求，在线自动测方式虽然可以时时测量，但受系统影响测量的密度值与人工测量值之间存在明显的差异性，既充分利用在线自动测量高效的特点，又能保证测量值准确就显得极为重要，因此在线密度测量是当前密度测量发展的趋势，而在线密度值的校准工作是亟待解决的首要问题。
发明内容
针对现有技术中的上述不足，本发明提供的种在线密度自动测量的原油密度值校正方法，解决了在线密度测量与人工测量之间标密度计算存在差异的问题，把在线密度测量影响因素归于温度和系统误差影响，通过二步法计算进行温度和系统误差校正，使在线密度测量计算的标密达到人工测量的精度。
为了达到上述发明目的，本发明采用的技术方案为：一种在线密度自动测量的原油密度值校正方法，包括以下步骤：
S1、在原油管道上安装在线密度测量仪，并利用在线密度测量仪对原油进行在线测量，分别获取原油在线密度测量的第一测量时间、第一测量温度以及在线视密度；
S2、同时对原油密度进行人工测量，分别获取原油的第二测量时间、第二测量温度以及人工视密度；
S3、选取第一测量时间等于第二测量时间时的相同样品；
S4、根据第二测量温度和人工视密度，获取相应视密度和温度的温度补偿系数，并对相同样品人工测量数据进行计算，得到人工标密；
S5、根据第一测量温度和在线视密度，获取相应视密度和温度的温度补偿系数，并对相同样品在线测量数据进行计算，得到在线标密；
S6、预设拟合多项式的系数，选择拟合多项式的拟合幂次，以在线标密作为自变量，并依据选取的拟合幂次计算得到校正标密；
S7、利用在线标密进行温度校正，以及利用校正标密进行系统误差校正，完成对原油密度值的校正。
本发明的有益效果是：本发明通过安装不同类型的在线密度性，确定了适合于原油密度测理仪器，通过长时间的人工与在线密度测量的对比与分析，建立了二步法的在线原油密度校正方法，第一步计算是在线视密求取在线标密，进行温度校正，第二步是以在线标密为自变量进行多项式计算校正标密，进行系统误差的校正。本发明解决了在线密度仪所测视密进行标密度换算时不准确的问题，把在线密度仪影响归于温度和系统误差影响，通过二步法的计算可以使在线密度仪所测视密度计算的校正标密达到人工标密数值，且误差极小，可实现在线密度仪时时测量，解决人工测量时人力和时间不足的问题，为在线密度测量的推广奠定理论和计算基础。
进一步地，所述人工标密的表达式如下：
Figure 691818DEST_PATH_IMAGE001
其中，
Figure 771770DEST_PATH_IMAGE002
表示20摄氏度的人工标密，
Figure 390970DEST_PATH_IMAGE003
表示某一温度T时的人工视密 度，r表示温度补偿系数。
上述进一步方案的有益效果是：通过计算把人工视密转化为标密，校正了人工测量的温度影响，计算的人工标密作为后续在线测量系统误差校正的多项式函数系数求取的真实值。
再进一步地，所述在线标密的表达式如下：
Figure 923583DEST_PATH_IMAGE004
其中，
Figure 388062DEST_PATH_IMAGE005
表示20摄氏度的在线标密，
Figure 322520DEST_PATH_IMAGE006
表示某一温度T时的在线视密 度，r表示温度补偿系数。
上述进一步方案的有益效果是：通过计算把在线测量的视密转化为在线标密，校正了在线测量的温度影响，并以在线标密作为后续在线测量系统误差校正的多项式函数的自变量。
再进一步地，所述校正标密的表达式如下：
Figure 112621DEST_PATH_IMAGE007
其中，
Figure 132530DEST_PATH_IMAGE008
表示校正标密，i表示拟合多项式的拟合幂次，n表示拟合多项式的 最高拟合幂次，ai表示拟合多项式的系数，
Figure 400700DEST_PATH_IMAGE009
表示第i次拟合幂次的校正标密。
上述进一步方案的有益效果是：通过选定的多项式幂次函数，把在线测量的在线标密作为函数自变量，通过计算获得校正标密，此步计算校正在线测量的系统误差，校正标密作为多项式系数求取的计算值。
再进一步地，所述拟合多项式系数的确定方法如下：
A1、最小化所有相同样品数据的人工标密和校正标密的误差；
A2、基于最小化误差，以最小二乘法进行求取，得到拟合多项式的系
上述进一步方案的有益效果是：通过最小二乘法求取的多项式系数，可以保障所有校正标密与人工标密尽可能一致。
再进一步地，所述最小化所有相同样品数据的人工标密和校正标密的误差的条件为：
Figure 691130DEST_PATH_IMAGE010
其中，E表示人工标密和校正标密的误差，
Figure 652132DEST_PATH_IMAGE011
表示校正标密，
Figure 893758DEST_PATH_IMAGE012
表示 20摄氏度的人工标密。
上述进一步方案的有益效果是：确保人工标密与经过温度和系统误差校正的在线测量的校正标密的误差最小，以便求取的多项式函数的计算的校正标密与人工标密相一致，且误差小。
再进一步地，以最小二乘法进行求取的表达式如下：
Figure 965619DEST_PATH_IMAGE013
Figure 140248DEST_PATH_IMAGE014
其中，A表示待求拟合多项式系数列向量，P表示人工标密列向量，an表示拟合多项 式的系数，T表示矩阵转置，n表示拟合多项式系数的最高幂次，m表示标密测量样品总数量， G-1表示广义逆，-1表示矩阵逆，D表示在线测量密度仪在线标密所形成的系数矩阵，
Figure 740994DEST_PATH_IMAGE015
表示 第k个在线样品的在线标密的i拟合幂次，
Figure 735495DEST_PATH_IMAGE016
，
Figure 345468DEST_PATH_IMAGE017
，k表示在线样品的次 序，i表示拟合多项式的拟合幂次，
Figure 109024DEST_PATH_IMAGE018
表示第m个人工标密。
上述进一步方案的有益效果是：在最小二乘原则下求取多项式函数的系数，确保由多项式函数求取的校正标密，作为最终标密输出能达到人工标密的精度，为在线测量推广和直接求取标密奠定基础。
附图说明
图1为本发明的方法流程图。
图2为本实施例中标密-密度差与温度的关系示意图。
图3为本实施例中在线密度标密与人工标密的关系图。
图4为本实施例中一次多项式拟合计算结果与误差分析示意图。
图5为本实施例中一次多项式拟合的函数图。
图6为本实施例中二次多项式拟合计算结果与误差分析示意图。
图7为本实施例二次多项式拟合的函数图。
图8为本实施例中三次多项式拟合计算结果与误差分析示意图。
图9为本实施例中三次多项式拟合的函数图。
图10为本实施例中四次多项式拟合计算结果与误差分析示意图。
图11为本实施例中四次多项式拟合的函数图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
本发明通过对比人工测量与在线测量的差异，对于在线测量的误差原因进行了探讨，并对提出了一种“二步法”的在线密度校正方法，主要认识如下：
1）在线密度测量主要受温度和系统影响，且在线密度测量与人工测量密度的差异主要是由于在线密度测量系统误差没有校正造成的。
2）通过二步法进行在线密度的校正，第一步校正温度影响，第二步校正系统误差，校正后的标密值与人工标密相近，其误差在±0.004之内，具有很高的精度。
3）不同系数和幂次的多项式拟合函数图形具有高度一致性表明在线密度系统误差具有很好的稳定性，不因拟合多项式的不同而有明显差异。
基于上述问题，如图1所示，本发明提供了一种在线密度自动测量的原油密度值校正方法，其实现方法如下：
S1、在原油管道上安装在线密度测量仪，并利用在线密度测量仪对原油进行在线测量，分别获取原油在线密度测量的第一测量时间、第一测量温度以及在线视密度；
S2、同时对原油密度进行人工测量，分别获取原油的第二测量时间、第二测量温度以及人工视密度；
S3、选取第一测量时间等于第二测量时间时的相同样品；
S4、根据第二测量温度和人工视密度，获取相应视密度和温度的温度补偿系数，并对相同样品人工测量数据进行计算，得到人工标密：
Figure 880671DEST_PATH_IMAGE019
（1）
其中，
Figure 395092DEST_PATH_IMAGE020
表示20摄氏度的人工标密，
Figure 808755DEST_PATH_IMAGE021
表示某一温度T时的人工视密 度，r表示温度补偿系数；
S5、根据第一测量温度和在线视密度，获取相应视密度和温度的温度补偿系数，并对相同样品在线测量数据进行计算，得到在线标密；
Figure 426819DEST_PATH_IMAGE004
（2）
其中，
Figure 900525DEST_PATH_IMAGE022
表示20摄氏度的在线标密，
Figure 604039DEST_PATH_IMAGE023
表示某一温度T时的在线视密 度，r表示温度补偿系数；
S6、预设拟合多项式的系数，选择拟合多项式的拟合幂次，以在线标密作为自变量，并依据选取的拟合幂次计算得到校正标密：
Figure 555815DEST_PATH_IMAGE024
（3）
其中，
Figure 293963DEST_PATH_IMAGE025
表示校正标密，i表示拟合多项式的拟合幂次，n表示拟合多项式的 最高拟合幂次，ai表示拟合多项式的系数，
Figure 407413DEST_PATH_IMAGE026
表示第i次拟合幂次的校正标密；
所述拟合多项式系数的确定方法如下：
A1、最小化所有相同样品数据的人工标密和校正标密的误差；
A2、基于最小化误差，以最小二乘法进行求取，得到拟合多项式的系数；
保证所有相同样品的人工标密与校正标密误差最小，即最小化所有相同样品数据的人工标密和校正标密的误差的条件为：
Figure 863802DEST_PATH_IMAGE027
（4）
其中，E表示人工标密和校正标密的误差，
Figure 353689DEST_PATH_IMAGE028
表示校正标密，
Figure 946345DEST_PATH_IMAGE029
表示20 摄氏度的人工标密；
以最小二乘法进行求取的表达式如下：
Figure 263319DEST_PATH_IMAGE030
Figure 941425DEST_PATH_IMAGE031
其中，A表示待求拟合多项式系数列向量，P表示人工标密列向量，an表示拟合多项 式的系数，T表示矩阵转置，n表示拟合多项式系数的最高幂次，m表示标密测量样品总数量， G-1表示广义逆，-1表示矩阵逆，D表示在线测量密度仪在线标密所形成的系数矩阵，
Figure 969424DEST_PATH_IMAGE032
表示 第k个在线样品的在线标密的i拟合幂次，
Figure 416585DEST_PATH_IMAGE033
，
Figure 402996DEST_PATH_IMAGE034
，k表示在线样品的次 序，i表示拟合多项式的拟合幂次，
Figure 568398DEST_PATH_IMAGE035
表示第m个人工标密，i=1时上标不标注；
本实施例中，基于选择是相同时间的样品进行的测量，所以人工测量与在线测量是对应的，有人工测量样品与在线测量样品值相对应，P与D中的标密是相同时刻对相同性质原油进行的测量，P是通过人工取样的，D是管道内的，可认为是不同的样品。
S7、利用在线标密进行温度校正，以及利用校正标密进行系统误差校正，完成对原油密度值的校正。
本发明根据二步计算法，实现了对在线测量视密进行标密求取的温度和系统误差校正，即按照公式（2）和公式（3）的在线密度测量的原油密度值校正方法，分别进行温度和误差校正。
本发明通过分析人工测量与在线测量的差异及影响因素，认为在线密度测量的影响主要是温度和测量系统两个因素，在线测量与人工测量的密度差异原因是在线测量未进行系统误差的校正。本发明建立了两步法在线测量密度校正方法，第一步通过温度补偿系数，计算得到在线标密校正温度影响；第二步以温度校正后的标密作为初值（即以在线标密作为初值），进行多项式拟合得到系统校正后的密度值，即校正标密，以系统校正密度作为计算值，以人工测量密度作为真实值，保证计算值与真实值的差最小以最小二乘的方式求解，获取多项式拟合函数系数，得到第二步系统误差校正表达式。通过实际在线密度和人工测量实验，进行了一次、二次、三次、四次的拟合校正，校正后的标密与人工测量值相近，其误差范围在±0.004以内，虽然拟合幂次不同，函数表达式不同，但其拟合段的函数图形具有高度一致性，进一步说明在线密度仪的系统误差具有稳定性。本发明建立的二步法密度校准方法，对于原油在线密度仪的校准有直接的指导价值，对于进一步的基于原油密度的研究、生产和交易具有重要意义。
下面对本发明作进一步的说明。
本实施例中，针对人工原油标密码的计算，人工原油密度受人工因素、原油性质、仪器因素、测量方法、环境因素的影响，像人工因素、仪器因素和测量方法的影响只要注意规范，其影响是可以避免的，而原油性质和环境因素的温度是直接影响原油密度的，其公式为：
Figure 665667DEST_PATH_IMAGE036
（1）
其中，
Figure 967335DEST_PATH_IMAGE037
表示20摄氏度的人工标密，
Figure 859068DEST_PATH_IMAGE038
表示某一温度T时的人工视密度， r表示温度补偿系数，不同的
Figure 511766DEST_PATH_IMAGE039
具有不同温度补偿系数，公式（1）表明标密码直接受原油 性质和温度的影响。
本实施例中，针对在线密仪标密的计算，现今在线密度仪主要有基于振动频率的音叉密度仪和振动密度仪，基于液柱压力的压差密度仪，基于放射性衰减的放射衰减密度仪以及基于科里奥利力的振荡的质量流量仪等，在线密度仪测量的是某一温度下的视密度，现今未有针对在线密度仪建立的视密度与标密度计算的公式，因为温度对原油的热胀冷缩的物理影响不因测试方法的改变而改变，因此在线密度仪的标密计算依然采用公式（1）方法，即有：
Figure 147147DEST_PATH_IMAGE040
（2）
其中，
Figure 568901DEST_PATH_IMAGE041
表示20摄氏度的在线标密，
Figure 631535DEST_PATH_IMAGE042
表示某一温度T时的在线视密 度，r表示温度补偿系数。
但由公式（1）和公式（2）对同一样品计算的标密不同。由人工计算的标密是严格执行规范测视密，然后再计算标密的，应为原油标密标准值，也可称为真值，而在线密度仪计算标密与人工不同，说明直接套用公式计算标密是存在问题，其原因主要在于影响因素没有进行合适的剔除或校正。人工测量时受人工因素、仪器因素、测量方法、环境因素和原油性质的影响，在测量时按测量方式的相应规范把人工因素、仪器因素和测量方法的影响尽可能进行了消除，而在线密度测量仪其影响因素同样包括仪器因素、测量方法、环境因素和原油性质的影响，在公式（2）只体现了环境因素和原油性质，并未对仪器因素和测量方法带来的影响进行消除，因此才造成的两个标密的不同。
本实施例中，针对在线密测量仪校正的方式，仪器因素和测量方法影响造成了人工标密与在线密度仪标密的差异，而仪器因素和测量方法的影响难以用一个量值去衡量。在线密度仪是一个独立的系统，其仪器因素和测量方法影响可以归并为系统因素，可认为系统因素对密度有一个系统固有影响值，因此在线密度标密的校正，把在线密度仪的标密做为测量值，人工标密做为真实值，测量值与真实值之间的差异影响是系统因素，只要建立起测量值与真实值，即在线密度仪的标密与人工标密的函数关系，使两者之间最小，即校正了系统因素误差。可以用多项式函数来表述校正标密，有公式（3）：
Figure 284713DEST_PATH_IMAGE043
（3）
其中，
Figure 723784DEST_PATH_IMAGE044
表示校正标密，i表示拟合多项式的拟合幂次，n表示拟合多项式的最 高拟合幂次，ai表示拟合多项式的系数，
Figure 265624DEST_PATH_IMAGE045
表示第i次拟合幂次的校正标密，
Figure 499159DEST_PATH_IMAGE046
为公 式（2）计算结果，欲想校正系统性误差，需要使（3）式计算结果与公式（1）计算结果差值最 小，即有：
Figure 860871DEST_PATH_IMAGE047
（4）
其中，E为校正标密与人工标密误差2范数，此时方程（4）变为拟合多项式系数ai为 函数，只需要使
Figure 838054DEST_PATH_IMAGE048
最小，即可求取
Figure 234400DEST_PATH_IMAGE049
值，采用最小二乘方法进行求解，即有：
Figure 638837DEST_PATH_IMAGE050
（5）
其中，A表示待求拟合多项式系数列向量，
Figure 753423DEST_PATH_IMAGE051
，
Figure 268718DEST_PATH_IMAGE052
表示人工标密列向 量，
Figure 519571DEST_PATH_IMAGE053
，m表示人工标密测样数量，n表示拟合多项式系数的数量，T 表示矩阵转置，-1表示矩阵逆，G-1表示广义逆， D表示在线测量密度仪在线标密所形成的系 数矩阵：
Figure 861953DEST_PATH_IMAGE054
（6）
公式（6）中的所有密度均为公式（2）计算的在线标密值，D表示在线测量密度仪在 线标密所形成的系数矩阵，
Figure 198256DEST_PATH_IMAGE055
表示第k个在线样品的在线标密的i拟合幂次，
Figure 517242DEST_PATH_IMAGE056
，
Figure 888181DEST_PATH_IMAGE057
，k表示在线样品的次序，i表示拟合多项式的拟合幂次，
Figure 634420DEST_PATH_IMAGE058
表示第m个人工 标密，i=1时上标不标注。
本实施例中，安装压差式密度在线测量仪，测量管道内原油的温度和视密，同时采用人工方法按照相关国标规定测量同时间管道内原油的温度和视密，运用公式（1）计算了人工测量的标密值，运用公式（2）计算在线密度仪的标密值，如图2所示，图2左图为压差式在线密度仪的标密与视密差与温度关系，图2右图为人工测量的相应样品的标密与视密差与温度关系，标密-视密差与温度均成线性关系，其斜率是温度补偿系数。
本实施例中，在线标密与人工标密的关系如图3，若两个标密一致应在图中分布于过原点的斜率为45度的直线上，从图中可以看出在线标密与人工标密明显不一致，此不一致的原因主要受在线密度仪的系统性误差影响。
本实施例中，以在线标密做为理论值的输入，采用公式（3）计算校正标密，以人工标密做为实际值，使校正标密和人工标密满足公式（4），采用公式（5）进行求解，分别做了公式（3）的一次、二次多项式、三次、四次多项式的拟合。图4为一次多项式拟合结果，原先明显离散度高的压差标密通过拟合校正其系统误差后变得集中且与人工标密相接近，并计算了拟合标密与人工标密的误差，误差值为±0.004之内。
本实施例中，一次多项式拟合的函数如图5所示，其公式为：
Figure 458019DEST_PATH_IMAGE059
（7）
其中，
Figure 580696DEST_PATH_IMAGE060
。二次多项式拟合结果图如图6所示，原先明 显离散度高的压差标密通过拟合校正其系统误差后变得集中且与人工标密相接近，并计算 了拟合标密与人工标密的误差，误差值为±0.004之内。
本实施例中，二次多项式拟合的函数如图7所示，其公式为：
Figure 540562DEST_PATH_IMAGE061
（8）
其中，
Figure 723282DEST_PATH_IMAGE062
。三次多项式拟合结果图如图8所示， 原先明显离散度高的压差标密通过拟合校正其系统误差后变得集中且与人工标密相接近， 并计算了拟合标密与人工标密的误差，误差值也在±0.004之内。
本实施例中，三次多项式拟合的函数如图9所示，其公式为：
Figure 34177DEST_PATH_IMAGE063
（9）
其中：a0=19.54188156，a1=-68.5853107，a2=-83.71289905，a3=-34.02360877四次多项式拟合结果图如图10所示，原先明显离散度高的压差标密通过拟合校正其系统误差后变得集中且与人工标密相接近，并计算了拟合标密与人工标密的误差，误差值也在±0.004之内。
本实施例中，四次多项式拟合的函数如图11所示，其公式为：
Figure 694966DEST_PATH_IMAGE064
（10）
其中： a0=66.43060154，a1=-295.76429788，a2=496.45905804，a3=-367.29587613，a4=100.90909091，从公式（7）、公式（8）、公式（9）和公式（10）可以看出虽然拟合函数的次数不同，拟合系数不同，但图5、图7、图9和图11的拟合函数的图形具有很高的一致性，说明虽然拟合函数不同，但系统误差是固有的且一致的，不因拟合函数不同而不同。

Claims (4)

1.一种在线密度自动测量的原油密度值校正方法，其特征在于，包括以下步骤：
S1、在原油管道上安装在线密度测量仪，并利用在线密度测量仪对原油进行在线测量，分别获取原油在线密度测量的第一测量时间、第一测量温度以及在线视密度；
S2、同时对原油密度进行人工测量，分别获取原油的第二测量时间、第二测量温度以及人工视密度；
S3、选取第一测量时间等于第二测量时间时的相同样品；
S4、根据第二测量温度和人工视密度，获取相应视密度和温度的温度补偿系数，并对相同样品人工测量数据进行计算，得到人工标密；
S5、根据第一测量温度和在线视密度，获取相应视密度和温度的温度补偿系数，并对相同样品在线测量数据进行计算，得到在线标密；
S6、预设拟合多项式的系数，选择拟合多项式的拟合幂次，以在线标密作为自变量，并依据选取的拟合幂次计算得到校正标密；
所述拟合多项式系数的确定方法如下：
A1、最小化所有相同样品数据的人工标密和校正标密的误差；
A2、基于最小化误差，以最小二乘法进行求取，得到拟合多项式的系数；
所述最小化所有相同样品数据的人工标密和校正标密的误差的条件为：
Figure 62370DEST_PATH_IMAGE001
其中，E表示人工标密和校正标密的误差，
Figure 879017DEST_PATH_IMAGE002
表示校正标密，
Figure 943925DEST_PATH_IMAGE003
表示20摄氏度的人工标密；
以最小二乘法进行求取的表达式如下：
Figure 28162DEST_PATH_IMAGE004
Figure 520323DEST_PATH_IMAGE005
其中，A表示待求拟合多项式系数列向量，P表示人工标密列向量，an表示拟合多项式的系数，T表示矩阵转置，n表示拟合多项式系数的最高幂次，m表示标密测量样品总数量，G-1表示广义逆，-1表示矩阵逆，D表示在线测量密度仪在线标密所形成的系数矩阵，
Figure 824266DEST_PATH_IMAGE006
表示第k个在线样品的在线标密的i拟合幂次，
Figure 161706DEST_PATH_IMAGE007
，
Figure 572221DEST_PATH_IMAGE008
，k表示在线样品的次序，i表示拟合多项式的拟合幂次，
Figure 969704DEST_PATH_IMAGE009
表示第m个人工标密；
S7、利用在线标密进行温度校正，以及利用校正标密进行系统误差校正，完成对原油密度值的校正。
2.根据权利要求1所述的在线密度自动测量的原油密度值校正方法，其特征在于，所述人工标密的表达式如下：
Figure 760943DEST_PATH_IMAGE010
其中，
Figure 902074DEST_PATH_IMAGE011
表示20摄氏度的人工标密，
Figure 960904DEST_PATH_IMAGE012
表示某一温度T时的人工视密度，r表示温度补偿系数。
3.根据权利要求1所述的在线密度自动测量的原油密度值校正方法，其特征在于，所述在线标密的表达式如下：
Figure 794867DEST_PATH_IMAGE013
其中，
Figure 807823DEST_PATH_IMAGE014
表示20摄氏度的在线标密，
Figure 752645DEST_PATH_IMAGE015
表示某一温度T时的在线视密度，r表示温度补偿系数。
4.根据权利要求1所述的在线密度自动测量的原油密度值校正方法，其特征在于，所述校正标密的表达式如下：
Figure 668911DEST_PATH_IMAGE016
其中，
Figure 673776DEST_PATH_IMAGE017
表示校正标密，i表示拟合多项式的拟合幂次，n表示拟合多项式的最高拟合幂次，ai表示拟合多项式的系数，
Figure 439606DEST_PATH_IMAGE018
表示第i次拟合幂次的校正标密。
CN202210611918.4A 2022-06-01 2022-06-01 一种在线密度自动测量的原油密度值校正方法 Active CN114689460B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210611918.4A CN114689460B (zh) 2022-06-01 2022-06-01 一种在线密度自动测量的原油密度值校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210611918.4A CN114689460B (zh) 2022-06-01 2022-06-01 一种在线密度自动测量的原油密度值校正方法

Publications (2)

Publication Number Publication Date
CN114689460A CN114689460A (zh) 2022-07-01
CN114689460B true CN114689460B (zh) 2022-08-23

Family

ID=82131158

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210611918.4A Active CN114689460B (zh) 2022-06-01 2022-06-01 一种在线密度自动测量的原油密度值校正方法

Country Status (1)

Country Link
CN (1) CN114689460B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2170909A (en) * 1985-02-08 1986-08-13 Spectra Tek Uk Limited Apparatus and method for monitoring crude oil
CN1030980A (zh) * 1987-05-06 1989-02-08 富士胶片公司 密度校正方法及装置
US5687100A (en) * 1996-07-16 1997-11-11 Micro Motion, Inc. Vibrating tube densimeter
WO2009003963A2 (de) * 2007-06-30 2009-01-08 Endress+Hauser Flowtec Ag Messsystem für ein in einer prozessleitung strömendes medium
WO2010134911A1 (en) * 2009-05-20 2010-11-25 Halliburton Energy Services, Inc. Determining fluid density
CN103674772A (zh) * 2013-11-27 2014-03-26 安徽秋田电子有限公司 柴油密度的温度补偿方法
CN109387456A (zh) * 2017-08-09 2019-02-26 中国石油化工股份有限公司 一种原油密度测量方法
CN109642864A (zh) * 2016-07-27 2019-04-16 环球油品有限责任公司 使用多个传感器进行密度测量的方法
WO2019096576A1 (de) * 2017-11-20 2019-05-23 Endress+Hauser Flowtec Ag Verfahren zum signalisieren einer standardfrequenz eines dichtemessers, welcher mindestens ein schwingfähiges messrohr zum führen eines mediums aufweist

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5355278B2 (ja) * 2009-07-29 2013-11-27 京都電子工業株式会社 振動式密度計の校正パラメータ決定方法及び密度算出方法
CN103063599B (zh) * 2011-10-18 2015-08-26 中国石油化工股份有限公司 由透射红外光谱预测原油密度的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2170909A (en) * 1985-02-08 1986-08-13 Spectra Tek Uk Limited Apparatus and method for monitoring crude oil
CN1030980A (zh) * 1987-05-06 1989-02-08 富士胶片公司 密度校正方法及装置
US5687100A (en) * 1996-07-16 1997-11-11 Micro Motion, Inc. Vibrating tube densimeter
WO2009003963A2 (de) * 2007-06-30 2009-01-08 Endress+Hauser Flowtec Ag Messsystem für ein in einer prozessleitung strömendes medium
WO2010134911A1 (en) * 2009-05-20 2010-11-25 Halliburton Energy Services, Inc. Determining fluid density
CN103674772A (zh) * 2013-11-27 2014-03-26 安徽秋田电子有限公司 柴油密度的温度补偿方法
CN109642864A (zh) * 2016-07-27 2019-04-16 环球油品有限责任公司 使用多个传感器进行密度测量的方法
CN109387456A (zh) * 2017-08-09 2019-02-26 中国石油化工股份有限公司 一种原油密度测量方法
WO2019096576A1 (de) * 2017-11-20 2019-05-23 Endress+Hauser Flowtec Ag Verfahren zum signalisieren einer standardfrequenz eines dichtemessers, welcher mindestens ein schwingfähiges messrohr zum führen eines mediums aufweist
CN111566467A (zh) * 2017-11-20 2020-08-21 恩德斯+豪斯流量技术股份有限公司 用于发信号通知具有至少一个用于传导介质的振动测量管的密度计的标准频率的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
U形振动管法测定原油密度的影响因素;郝晨新等;《理化检验(物理分册)》;20161231(第04期);第253-256页 *
原油动态计量的误差分析与减少误差的方法;冯海东等;《油气储运》;19961231;第15卷(第11期);第35-37页 *
石油标准密度简便换算方法研究与实践;张宁 等;《广州化工》;20200229;第48卷(第03期);第122-124页 *
自动密度仪测定原油密度及注意事项;李权斌;《理化检验(物理分册)》;20061231;第42卷(第03期);第134-136页 *

Also Published As

Publication number Publication date
CN114689460A (zh) 2022-07-01

Similar Documents

Publication Publication Date Title
CN110514391B (zh) 一种五孔压力探针测量三维流场不确定度评定方法
CN110514390B (zh) 一种三孔压力探针测量二维流场不确定度评定方法
CN1128989C (zh) 用来校准压差流体流量测量系统的方法
CN106096223A (zh) 一种五孔探针数据处理方法
CN114459573A (zh) 超声波燃气表的新型校准模型、校准方法及超声波燃气表
CN115389567A (zh) 水质电导率传感器的温度补偿算法
CN114689460B (zh) 一种在线密度自动测量的原油密度值校正方法
CN110553631A (zh) 一种关于水位流量关系的水位测量系列误差分析方法
CN106525180B (zh) 天然气压缩因子计算方法
Spitzer et al. pH and electrolytic conductivity as parameters to characterize bioethanol
Folgerø et al. Uncertainty analysis of multiphase flow meters used for allocation measurements: Field experiences and future challenges
CN110609248A (zh) 一种基于多参考样品的同轴谐振腔校准方法及系统
Montecinos et al. An ADER-type scheme for a class of equations arising from the water-wave theory
CA2432403A1 (en) Method for testing liquid hydrocarbons using refractive index
Eremin Error Estimation of Indirect Measurement Results for Certain Quantities
US7349819B2 (en) Method for providing robustness to the solution of the cubic equation of state
Sysoev Analysis of the duration of metrologically reliable operation of measuring instruments by stochastic extrapolation techniques
CN111380587A (zh) 一种对采用蒙特卡罗法测量立式罐容量的结果误差补偿方法
Vetter Quantifying measurement uncertainty in analytical chemistry–A simplified practical approach
CN117553938B (zh) 一种高低温环境下温度二次仪表的误差修正方法
Ponomarev et al. An estimate of the errors in measuring the thermal properties of solids
Scott Distribution of soil conductivity and its relation to underground corrosion
RU2682540C9 (ru) Способ настройки измерительного канала расхода среды с сужающим устройством
Teplukh et al. Design of linear capillary measuring transducers for low gas flow rates
CN117890454A (zh) 工业在线质谱分析方法

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