CN109991396A - 一种考虑盲管段影响的成品油管道混油长度的计算方法 - Google Patents

一种考虑盲管段影响的成品油管道混油长度的计算方法 Download PDF

Info

Publication number
CN109991396A
CN109991396A CN201910370608.6A CN201910370608A CN109991396A CN 109991396 A CN109991396 A CN 109991396A CN 201910370608 A CN201910370608 A CN 201910370608A CN 109991396 A CN109991396 A CN 109991396A
Authority
CN
China
Prior art keywords
oil
mixture spread
product
reynolds number
pipeline
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.)
Pending
Application number
CN201910370608.6A
Other languages
English (en)
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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum 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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201910370608.6A priority Critical patent/CN109991396A/zh
Publication of CN109991396A publication Critical patent/CN109991396A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/06Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness specially adapted for measuring length or width of objects while moving
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/26Oils; Viscous liquids; Paints; Inks
    • G01N33/28Oils, i.e. hydrocarbon liquids

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • General Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Pipeline Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种考虑盲管段影响的成品油管道混油长度的计算方法,通过临界雷诺数将A油和B油这两种流体的流动状态区分开来,并判段两种油品的混油特点,从而选择公式来计算混油长度,公式中含拖尾油长度。本发明考虑了输油管道沿程阀室和泵站产生的拖尾油对混油的影响,进一步提高了混油量计算的精确度,更好地指导混油切割作业,节约成本、提高效率。

Description

一种考虑盲管段影响的成品油管道混油长度的计算方法
技术领域
本发明属于成品油顺序输送混油长度计算领域,尤其因阀室或泵站产生的盲管段对混油长度的影响,即成品油长距离输送产生的混油长度计算方法。
背景技术
在成品油顺序输送的过程中会产生混油,在切除混油时往往会忽略拖尾油对油品质量的影响。由于拖尾油的存在会造成相当多的混油。拖尾油形成的原因主要有两种,一种是层流边界层的影响,另一种是死油段或管道起伏段滞留的前行油品。层流底层的传质过程制约着由粘度差异而导致的速度分布变化,进而影响前后两种油品间混油的形成。在顺序输送过程中,前行油品所形成的边界层会慢慢掺混进后行油品,形成拖尾油。
成品油长输管道有许多阀室和泵站,油品在经过这些阀室、泵站时会形成死油段,当后行油品流经充满前行油品的死油段时,首先由于重力作用,后行油品进入死油区,前行油品从中流出,两种油品在接口处形成封闭交换流,这时引起混油的主要因素是对流传递。当死油段前行油品被替换了一半左右后,对流传递的影响逐渐减小,扩散传递开始发挥作用。在湍流扩散作用下,后行油品逐渐混入死油区,前行油品慢慢被替换。在这两个过程中湍流扩散进行的时间要明显大于对流传递,从死油区流出的前行油品会在主管中形成较长的拖尾油。
现有的混油长度计算公式未考虑盲管段产生的拖尾油对混油量的影响。盲管段产生的拖尾油会污染成品油,降低成品油的质量,在输油站场也无法做到精确切割混油段。
发明内容
本发明目的在于克服上述技术缺陷,提供一种更加精确的成品油顺序输送混油长度计算方法,该计算方法能够综合考虑因盲管段所形成的拖尾油对混油长度的影响,使得该式计算结果更加符合工程实际,实现输油管道混油的有效控制、节约成本、提高效率的目的。
本发明的目的是通过以下技术方案来实现的:
一种考虑盲管段影响的成品油管道混油长度的计算方法,通过临界雷诺数将A油和B油这两种流体的流动状态区分开来,并判段两种油品的混油特点,从而选择公式来计算混油长度,公式中含拖尾油长度。
作为优选方式,通过模拟和数据分析后,运用软件对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到混油长度计算公式。
作为优选方式,软件为matlab。
作为优选方式,包括如下步骤:
第一步:收集管道数据,根据管径壁厚计算出管道内径;由运行温度下的前行A油(油品A)与后行B油(油品B)的运动粘度计算出混油粘度;
第二步:根据管道内径、混油粘度、流速计算出混油临界雷诺数及混油雷诺数;
第三步:判断混油临界雷诺数与雷诺数之间的大小关系选择相应的混油长度计算式。
作为优选方式,运用matlab对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到如下式子:
其中T——盲管中油品替换所需时间,s;
V——主管段流速,m/s;
Lm——盲管长度,m。
作为优选方式,对于多个盲管的情况,前一盲管中的前行A油进入主管后势必会对管道中油品密度造成影响,进而影响后续盲管的油品携带速度和拖尾油长度;
利用模拟得到的数据,进行多元非线性拟合得到拖尾油长度C′:
C′=[0.21ln(N)+0.96]TV
其中N——盲管数量,个;
T——盲管中油品替换所需时间,s;
V——主管段流速,m/s。
作为优选方式,只考虑混油段中油品质量不合格的部分,计算公式所计算的混油长度考虑前行油品体积分数大于1%的部分;
在6s-9s之间汽油体积分数汽油体积分数已低于1%,把前行油品浓度为99%-1%的两界面间的距离定义为混油段长度,对拖尾油长度C′进行进一步推导得:
C′=9[0.21ln(N)+0.96]V
其中N——盲管数量,个;
V——主管段流速,m/s。
作为优选方式,当雷诺数大于混油临界雷诺数时混油长度计算公式为:
C=11.75d0.5L0.5Re-0.1+9[0.21ln(N)+0.96]V
当雷诺数小于混油临界雷诺数时混油长度计算公式为:
其中C——混油长度,m;
N——盲管段数量,个;
d——管道内径,m;
L——管道长度,m;
Re——混油雷诺数;
V——流速,m/s;
其中Rej——临界雷诺数;
υA——A油在输送温度下的运动粘度,m2/s;
υB——B油在输送温度下的运动粘度,m2/s;
υ——混油的计算运动粘度,m2/s。
本发明的有益效果是:本发明考虑了输油管道沿程阀室和泵站产生的拖尾油对混油的影响,进一步提高了混油量计算的精确度,更好地指导混油切割作业,节约成本、提高效率。
附图说明
图1是本发明其他两种方法计算所得混油量相比较示意图。
图2是阀室和泵站内由于旁通而形成的盲管段示意图。
图3是本发明有盲管段影响的成品油管道混油长度计算方法的主要计算步骤。
图4是CFD模拟盲管混油的简化模型。
图5是模拟得到的拟合盲管中油品替换所需时间的数据。
图6是模拟得到的盲管数量对拖尾油长度影响的数据。
图7是不同密度下主管道混油量随时间变化曲线图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
图2是阀室和泵站内由于旁通而形成的盲管段示意图,图4为CFD模拟盲管混油的简化模型。本发明提供的一种考虑盲管段影响的成品油管道混油长度的计算方法,通过临界雷诺数将A油和B油这两种流体的流动状态区分开来,并判段两种油品的混油特点,从而选择公式来计算混油长度,公式中含拖尾油长度。如图3所示为一个实施例的流程示意图。
在一个优选实施例中,通过模拟和数据分析后,运用软件对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到混油长度计算公式。
在一个优选实施例中,软件为matlab。
在一个优选实施例中,包括如下步骤:
第一步:收集管道数据,根据管径壁厚计算出管道内径;由运行温度下的前行A油(油品A)与后行B油(油品B)的运动粘度计算出混油粘度;
第二步:根据管道内径、混油粘度、流速计算出混油临界雷诺数及混油雷诺数;
第三步:判断混油临界雷诺数与雷诺数之间的大小关系选择相应的混油长度计算式。
在一个优选实施例中,运用matlab对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到如下式子:
其中T——盲管中油品替换所需时间,s;
V——主管段流速,m/s;e——自然常数,2.718;
Lm——盲管长度,m。
在一个优选实施例中,对于多个盲管的情况,前一盲管中的前行A油进入主管后势必会对管道中油品密度造成影响,进而影响后续盲管的油品携带速度和拖尾油长度;
利用模拟得到的数据,进行多元非线性拟合得到拖尾油长度C′:
C′=[0.21ln(N)+0.96]TV
其中N——盲管数量,个;
T——盲管中油品替换所需时间,s;
V——主管段流速,m/s。
在一个优选实施例中,
在实际计算混油长度时,出于混油切割的需要,通常只考虑混油段中油品质量不合格的部分,出于对混油的定义和现场混油切割的需要,现有的计算公式所计算的混油长度也是考虑前行油品体积分数大于1%的部分,因此在公式推导过程中也只考虑这部分,后面大部分拖尾油对油品质量的影响可以忽略。
由图7数据可知,在6s-9s之间汽油体积分数汽油体积分数已低于1%,把前行油品浓度为99%-1%的两界面间的距离定义为混油段长度,对拖尾油长度C′进行进一步推导得:
C′=9[0.21ln(N)+0.96]V
其中N——盲管数量,个;
V——主管段流速,m/s。
在一个优选实施例中,当雷诺数大于混油临界雷诺数时混油长度计算公式为:
C=11.75d0.5L0.5Re-0.1+9[0.21ln(N)+0.96]V
当雷诺数小于混油临界雷诺数时混油长度计算公式为:
其中C——混油长度,m;
N——盲管段数量,个;
d——管道内径,m;
e——自然常数,2.718;
L——管道长度,m;
Re——混油雷诺数;
V——流速,m/s;
其中Rej——临界雷诺数;
υA——A油在输送温度下的运动粘度,m2/s;
υB——B油在输送温度下的运动粘度,m2/s;
υ——混油的计算运动粘度,m2/s。
本发明一种考虑盲管段影响的成品油管道混油长度的计算方法,对兰成渝管道输送情况进行分析,全线输量为21220.021t,终点混油量为334.109t。如图1所示,扩散理论计算混油量为396.785t,相对误差为18.8%;传统经验计算式计算得到混油量为314.415t,相对误差为5.9%。本发明计算结果为331.067t,相对误差为0.9%。本发明精确度更高。
本发明针对盲管段易形成拖尾油的情形,分析拖尾油的形成机理,将这一工程中实际存在的混油加以考虑,利用CFD软件对不同的流速及盲管长度的情况进行模拟得到混油相关实验数据,运用matlab软件对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到盲管段所产生的拖尾油长度,结合混油长度计算的经验计算式并考虑盲管段影响得到混油长度计算的一种新计算方法,进而使混油长度计算更精确,更好的指导现场的混油切割工作。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,应当指出的是,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:通过临界雷诺数将A油和B油这两种流体的流动状态区分开来,并判段两种油品的混油特点,从而选择公式来计算混油长度,公式中含拖尾油长度。
2.根据权利要求1所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:通过模拟和数据分析后,运用软件对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到混油长度计算公式。
3.根据权利要求2所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:软件为matlab。
4.根据权利要求1所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于,包括如下步骤:
第一步:收集管道数据,根据管径壁厚计算出管道内径;由运行温度下的前行A油与后行B油的运动粘度计算出混油粘度;
第二步:根据管道内径、混油粘度、流速计算出混油临界雷诺数及混油雷诺数;
第三步:判断混油临界雷诺数与雷诺数之间的大小关系选择相应的混油长度计算式。
5.根据权利要求1所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:运用matlab对盲管段油品替换时间、盲管长度和流速进行多元非线性拟合,得到如下式子:
其中T——盲管中油品替换所需时间,s;
V——主管段流速,m/s;
e——自然常数,2.718;
Lm——盲管长度,m。
6.根据权利要求5所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:对于多个盲管的情况,前一盲管中的前行A油进入主管后势必会对管道中油品密度造成影响,进而影响后续盲管的油品携带速度和拖尾油长度;
利用模拟得到的数据,进行多元非线性拟合得到拖尾油长度C′:
C′=[0.21ln(N)+0.96]TV
其中N——盲管数量,个;
T——盲管中油品替换所需时间,s;
V——主管段流速,m/s。
7.根据权利要求5或6所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:
只考虑混油段中油品质量不合格的部分,计算公式所计算的混油长度考虑前行油品体积分数大于1%的部分;
在6s-9s之间汽油体积分数汽油体积分数已低于1%,把前行油品浓度为99%-1%的两界面间的距离定义为混油段长度,对拖尾油长度C′进行进一步推导得:
C′=9[0.21ln(N)+0.96]V
其中N——盲管数量,个;
V——主管段流速,m/s。
8.根据权利要求1所述的一种考虑盲管段影响的成品油管道混油长度的计算方法,其特征在于:当雷诺数大于混油临界雷诺数时混油长度计算公式为:
C=11.75d0.5L0.5Re-0.1+9[0.21ln(N)+0.96]V
当雷诺数小于混油临界雷诺数时混油长度计算公式为:
其中C——混油长度,m;
N——盲管段数量,个;
d——管道内径,m;
e——自然常数,2.718;
L——管道长度,m;
Re——混油雷诺数;
V——流速,m/s;
其中Rej——临界雷诺数;
υA——A油在输送温度下的运动粘度,m2/s;
υB——B油在输送温度下的运动粘度,m2/s;
υ——混油的计算运动粘度,m2/s。
CN201910370608.6A 2019-05-06 2019-05-06 一种考虑盲管段影响的成品油管道混油长度的计算方法 Pending CN109991396A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910370608.6A CN109991396A (zh) 2019-05-06 2019-05-06 一种考虑盲管段影响的成品油管道混油长度的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910370608.6A CN109991396A (zh) 2019-05-06 2019-05-06 一种考虑盲管段影响的成品油管道混油长度的计算方法

Publications (1)

Publication Number Publication Date
CN109991396A true CN109991396A (zh) 2019-07-09

Family

ID=67136061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910370608.6A Pending CN109991396A (zh) 2019-05-06 2019-05-06 一种考虑盲管段影响的成品油管道混油长度的计算方法

Country Status (1)

Country Link
CN (1) CN109991396A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110608368A (zh) * 2019-08-01 2019-12-24 石化盈科信息技术有限责任公司 物料长输线的进厂油头切换时间的预测方法及切换系统
CN114963016A (zh) * 2022-06-30 2022-08-30 国家石油天然气管网集团有限公司 用于确定成品油顺序输送管道混油量的方法及装置
CN116307303A (zh) * 2023-05-24 2023-06-23 中国石油大学(华东) 机理-数据双驱动的成品油管道混油长度预测方法及系统
CN117150931A (zh) * 2023-10-30 2023-12-01 中国石油大学(华东) 基于混合单隐层神经网络的混油长度在线估计方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080277318A1 (en) * 2004-07-21 2008-11-13 Rj Oil Sands Inc. Separation and recovery of bitumen oil from tar sands
CN102606887A (zh) * 2011-01-25 2012-07-25 中国石油天然气集团公司 一种成品油管道顺序输送混油界面位置的判定方法
CN204043699U (zh) * 2014-08-15 2014-12-24 中国石油天然气股份有限公司 一种混油量影响模拟试验装置
CN106408162A (zh) * 2016-08-29 2017-02-15 中国石化销售有限公司华东分公司 一种成品油管道顺序输送批量控制方法
CN107524921A (zh) * 2017-08-09 2017-12-29 中国石油大学(北京) 成品油顺序输送管道混油界面跟踪方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080277318A1 (en) * 2004-07-21 2008-11-13 Rj Oil Sands Inc. Separation and recovery of bitumen oil from tar sands
CN102606887A (zh) * 2011-01-25 2012-07-25 中国石油天然气集团公司 一种成品油管道顺序输送混油界面位置的判定方法
CN204043699U (zh) * 2014-08-15 2014-12-24 中国石油天然气股份有限公司 一种混油量影响模拟试验装置
CN106408162A (zh) * 2016-08-29 2017-02-15 中国石化销售有限公司华东分公司 一种成品油管道顺序输送批量控制方法
CN107524921A (zh) * 2017-08-09 2017-12-29 中国石油大学(北京) 成品油顺序输送管道混油界面跟踪方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宫敬等: "成品油管道的混油计算方法研究", 《管道技术与设备》 *
郑劲: "柴推汽顺序输送混油尾段特性研究", 《中国优秀硕士学位论文全文数据库 工程科技1辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110608368A (zh) * 2019-08-01 2019-12-24 石化盈科信息技术有限责任公司 物料长输线的进厂油头切换时间的预测方法及切换系统
CN110608368B (zh) * 2019-08-01 2021-06-22 石化盈科信息技术有限责任公司 物料长输线的进厂油头切换时间的预测方法及切换系统
CN114963016A (zh) * 2022-06-30 2022-08-30 国家石油天然气管网集团有限公司 用于确定成品油顺序输送管道混油量的方法及装置
CN114963016B (zh) * 2022-06-30 2024-06-04 国家石油天然气管网集团有限公司 用于确定成品油顺序输送管道混油量的方法及装置
CN116307303A (zh) * 2023-05-24 2023-06-23 中国石油大学(华东) 机理-数据双驱动的成品油管道混油长度预测方法及系统
CN116307303B (zh) * 2023-05-24 2023-09-26 中国石油大学(华东) 机理-数据双驱动的成品油管道混油长度预测方法及系统
CN117150931A (zh) * 2023-10-30 2023-12-01 中国石油大学(华东) 基于混合单隐层神经网络的混油长度在线估计方法及系统
CN117150931B (zh) * 2023-10-30 2024-01-30 中国石油大学(华东) 基于混合单隐层神经网络的混油长度在线估计方法及系统

Similar Documents

Publication Publication Date Title
CN109991396A (zh) 一种考虑盲管段影响的成品油管道混油长度的计算方法
Sukhodolov et al. Case study: Effect of submerged aquatic plants on turbulence structure in a lowland river
CN102087130B (zh) 基于cfd技术的多声路超声流量计弯管安装的声路优化方法
CN108304620B (zh) 一种核反应堆堆芯棒束通道内硼扩散过程的计算方法
Yang et al. Performance prediction of erosion in elbows for slurry flow under high internal pressure
CN113006768A (zh) 气井结垢模拟装置及模拟方法
Wang et al. Predicting erosion in a non-Newtonian shear-thinning jet flow with validated CFD models from PIV and PTV measurements
Nizovtseva et al. Simulation of two‐phase air–liquid flows in a closed bioreactor loop: Numerical modeling, experiments, and verification
CN105547885A (zh) 气井冲蚀实验系统及其冲蚀方法和流束调节装置
Su et al. Analysis on residence time distribution curve of continuous casting tundish by combined model
Chen et al. Computational Fluid Dynamics Modeling of Argon–Steel (–Slag) Multiphase Flow in an Ruhrstahl–Heraeus Degasser: A Review of Past Numerical Studies
CN116362153A (zh) 一种基于cfd的变径管成品油顺序输送混油长度的计算方法
CN110728060A (zh) 一种感潮河段船闸开通闸安全运行条件的确定方法
CN205636627U (zh) 用于潮汐河口泥沙物理模型的自动加沙装置
CN205638300U (zh) 一种钻井工况监测系统
RU2013106523A (ru) Эмульгатор и способ определения параметров для эмульгатора
CN208888220U (zh) 一种评价压裂返排液消泡剂动态消泡性能的装置
CN109184796B (zh) 一种正洞掌子面风管管口瓦斯涌出量检测瓦斯浓度的方法
Hashmi Oil film dynamics in aero engine bearing chambers: fundamental investigations and numerical modelling
CN105973319A (zh) 一种控制棒驱动机构排污系统水力特性计算方法
CN202054674U (zh) 一种冷却水加药装置
Cho et al. PREDICTION OF A HEAT TRANSFER TO CO₂ FLOWING IN AN UPWARD PATH AT A SUPERCRITICAL PRESSURE
KR20090092472A (ko) 핵연료 피복관 크러드 모사시험장치 및 방법
CN109817283A (zh) 一个适用于含离子液体体系的联合热力学模型—comso-unifac
Gabl et al. Adapted evaluation criterion for the velocity distribution in front of turbines

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190709