CN112417697A - 混合原油配置方法及装置 - Google Patents
混合原油配置方法及装置 Download PDFInfo
- Publication number
- CN112417697A CN112417697A CN202011336266.5A CN202011336266A CN112417697A CN 112417697 A CN112417697 A CN 112417697A CN 202011336266 A CN202011336266 A CN 202011336266A CN 112417697 A CN112417697 A CN 112417697A
- Authority
- CN
- China
- Prior art keywords
- oil
- mixing
- stirring
- crude oil
- mixed
- 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
Links
- 239000010779 crude oil Substances 0.000 title claims abstract description 89
- 238000002360 preparation method Methods 0.000 title description 3
- 239000003921 oil Substances 0.000 claims abstract description 221
- 238000002156 mixing Methods 0.000 claims abstract description 177
- 238000003756 stirring Methods 0.000 claims abstract description 94
- 238000000034 method Methods 0.000 claims abstract description 82
- 230000008569 process Effects 0.000 claims abstract description 41
- 238000004088 simulation Methods 0.000 claims description 34
- 239000000295 fuel oil Substances 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 9
- 230000009467 reduction Effects 0.000 description 9
- 230000006872 improvement Effects 0.000 description 8
- 239000012895 dilution Substances 0.000 description 6
- 238000010790 dilution Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 239000000084 colloidal system Substances 0.000 description 4
- 238000007865 diluting Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005484 gravity Effects 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000005494 condensation Effects 0.000 description 2
- 239000002283 diesel fuel Substances 0.000 description 2
- 238000003113 dilution method Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003502 gasoline Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 125000000174 L-prolyl group Chemical group [H]N1C([H])([H])C([H])([H])C([H])([H])[C@@]1([H])C(*)=O 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 239000002199 base oil Substances 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000013043 chemical agent Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003085 diluting agent Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Classifications
-
- 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
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本说明书实施例提供混合原油配置方法及装置。所述方法包括：根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油，从而提高混合原油的配置效率。
Description
技术领域
本说明书实施例涉及稠油掺稀技术领域，特别涉及一种混合原油配置方法及装置。
背景技术
在工程现场常需要将两种或两种以上油品进行混合，这种工艺可用于稠油掺稀降粘输送、油品加药改性等原油处理工艺。其中，稠油掺稀工艺是稠油长距离输送常用的方法，通过向稠油中加入粘度较低的稀油，二者根据相似相容原理互相融合，稀释稠油，使稠油中胶质、沥青质浓度下降，增强流动性，从而使稠油粘度下降。
稠油掺稀工艺中，一般常用稀释剂有稀原油、凝析油、柴油等，而且稀释后的混合油利用普通原油管道即可输送。在油罐中运用搅拌器或者在管道中放置混合器即可实现稠油掺稀工艺。
稠油掺稀法操作运行简单，降粘效果显著，掺稀后混合油的凝点比稀油低，长距离输送时沿程摩阻损失极大地降低，使得稠油可实现低能耗长距离输送，不过所掺稀油的密度和粘度直接影响降粘效果，具体影响规律需要取不同稠稀比的混合原油进行多组实验，研究稠稀比、掺稀油密度和稀油的粘度对稠油的降粘规律，这需要在现场多次取样，耗时耗力，受制条件多。
发明内容
本说明书实施例的目的是提供一种混合原油配置方法及装置，以提高混合原油的配置效率。
为解决上述问题，本说明书实施例提供一种混合原油配置方法，所述方法包括：获取工程现场油品混合设备的实际尺寸和运行参数；根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
为解决上述问题，本说明书实施例还提供一种混合原油配置装置，所述装置包括：获取模块，用于获取工程现场油品混合设备的实际尺寸和运行参数；第一模拟模块，用于根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；第二模拟模块，用于根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；获得模块，用于在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
由以上本说明书实施例提供的技术方案可见，本说明书实施例中，可以获取工程现场油品混合设备的实际尺寸和运行参数；根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。本说明书实施例提供的混合原油配置方法，利用至少两种不同的油品，辅以仿真模拟结果，可以在实验室内配制出与工程现场混合状态相同的异种油品混合原油，避免了因实验需要多次返回现场取样的情况，进而提高混合原油的配置效率。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施例一种混合原油配置方法的流程图；
图2为本说明书实施例一种混合原油配置装置的功能模块示意图。
具体实施方式
下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。
在本说明书实施例中，稠油是指相对密度
的原油。稠油是全球石油烃类能源中的重要组成部分。稠油由于粘度高，流动阻力大，不易开采和运输，其特点是含沥青质、胶质较高。为便于稠油的运输和开采，通常需要对稠油进行降粘处理，即降低稠油的粘度。目前常用的降粘方法有加热法、掺稀油法、稠油改质法和化学药剂法等。
其中，掺稀油法主要是利用相似相容性质和稀释原理，有效降低稠油的粘度和密度。其中，稀油可以为相对于稠油来说粘度较低的油品。稀油掺入稠油后，稠油的胶质沥青相对含量降低，稀油中的轻组分使堆积紧密的胶体结构产生一定程度的破坏，宏观表现为稠油粘度降低。稠油掺稀法操作运行简单，降粘效果显著，掺稀后混合油的凝点比稀油低，长距离输送时沿程摩阻损失极大地降低，使得稠油可实现低能耗长距离输送，不过所掺稀油的密度和粘度直接影响降粘效果，具体影响规律需要取不同稠稀比的混合原油进行多组实验，研究稠稀比、掺稀油密度和稀油的粘度对稠油的降粘规律，这需要在现场多次取样，耗时耗力，受制条件多。
考虑到如果通过仿真软件模拟工程现场的不同油品的混合过程，并通过计算表征油品混合状态的指标来确定油品的混合状态，再通过仿真软件模拟实验室不同油品的混合过程，在两者表征油品混合状态的指标相同的情况下，即可将工程现场与实验室不同油品的混合状态关联起来，则有望解决现有技术中需要在现场多次取样，耗时耗力的问题，提高与工程现场混合状态相同的混合原油的配置效率。基于此，本说明书实施例提供了一种混合原油配置方法。
在本说明书实施例中，执行所述混合原油配置方法的主体可以是具有逻辑运算功能的电子设备，所述电子设备可以是服务器或客户端，所述客户端可以为台式电脑、平板电脑、笔记本电脑、工作站等。当然，客户端并不限于上述具有一定实体的电子设备，其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件，该程序软件可以运行于上述电子设备中。
图1为本说明书实施例一种混合原油配置方法的流程图。如图1所示，所述混合原油配置方法可以包括以下步骤。
S110：获取工程现场油品混合设备的实际尺寸和运行参数。
在一些实施例中，不同工程现场油品通常可以包括不同的油品混合设备。所述油品混合设备用于混合多种不同的油品，例如在稠油掺稀过程中药用到所述油品混合设备。所述油品混合设备通常可以包括至少一个稠油入口、至少一个稀油入口和一个混油出口。所述稠油入口用于输出稠油，所述稀油入口用于输入稀油，所述混油出口用于输出经过混合的稠油和稀油。
在一些实施例中，所述油品混合设备可以包括SK型混合器和异径管道。其中，所述SK型混合器又称单螺旋形静态混合器，它的单元是扭转180°或270°的螺旋板，安装入管壳时相邻的螺旋板分别为左旋和右旋。
在一些实施例中，所述油品混合设备通常基于一定的运行参数运行。所述运行参数可以包括稠油入口输入稠油的速度、稀油入口输入稀油的速度、温度等参数。
在一些实施例中，服务器可以采用任何方式获取工程现场油品混合设备的实际尺寸和运行参数。例如，用户可以直接向服务器发送工程现场油品混合设备的实际尺寸和运行参数，服务器可以进行接收；又如除去所述服务器以外的其它电子设备可以向服务器发送工程现场油品混合设备的实际尺寸和运行参数，服务器可以进行接收，在本说明书实施例中，对服务器采用何种方式获取工程现场油品混合设备的实际尺寸和运行参数不作限定。
S120：根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态。
在一些实施例中，可以通过仿真软件来模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度。具体的，可以根据所述工程现场油品混合设备的实际尺寸建立第一仿真模型；将所述运行参数作为模拟参数，使用所述第一仿真模型模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度。
在一些实施例中，可以使用SOLIDWORKS软件和ANSYS软件建立建立与工程现场油品混合设备尺寸一致的第一仿真模型。具体的，可以在SOLIDWORKS软件中输入相应的参数，如油品混合设备设备的形状、尺寸的参数，从而建立模型，再将所述建立的模型导入ANSYS软件中，从而建立所述第一仿真模型。其中，所述SOLIDWORKS软件是一种基于Windows开发的三维CAD系统；所述ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。所述ANSYS软件能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer，NASTRAN，Alogor，I－DEAS，AutoCAD等。
在一些实施例中，所述至少两种预设类别的油品包括至少一种稠油和至少一种稀油。所述稀油可以为粘度较低的油品，例如可以是稀原油、凝析油、柴油等。
在一些实施例中，可以通过在仿真模型中设置油品的密度和油品的动力粘度模拟所述至少两种预设类别的油品。
在实验室研究不同油品混合时，一般需要研究多个混合比例、多个混合温度下不同油品混合油的流变性和物性，混合油样品种类多、实验量大，致使研究工作进展缓慢。因此，可以使用混合不均匀度来表征混合油品的混合状态，从而提高油品混合的研究效率。具体的，某截面上一种油品体积分数的平均值主要取决于不同油品各自的流入流量，在流动状态稳定不变的条件下，截面上一种油品体积分数的平均值基本可以维持稳定，所以，某截面上异种油品混合越不均匀，截面上某点上的体积分数与平均值相差就越大，导致混合不均匀度值越大。根据模型内某截面一种油品的体积分数可得到该截面的混合不均匀度值，其表达式为：
其中，
表示混合不均匀度，ci表示截面上第i个点的一种油品体积分数，i＝1,2…n，
表示截面上的一种油品体积分数平均值，n表示截面上取点的个数。
S130：根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度。
在一些实施例中，可以通过仿真软件来模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度。具体的，可以根据配置混合原油所采用的搅拌设备的形状和尺寸建立第二仿真模型；使用所述第二仿真模型模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度。
在一些实施例中，所述第二仿真模型模拟的至少两种预设类别的油品与所述第一仿真模型模拟的至少两种预设类别的油品相同，即设置相同的油品的密度和油品的动力粘度模拟所述至少两种预设类别的油品。
在一些实施例中，可以模拟所述搅拌设备不同的搅拌速率下配置混合原油的过程，例如分别模拟所述搅拌设备在1000r/min、2000r/min或其他搅拌速率下配置混合原油的过程，并根据混合不均匀度的表达式计算得到不同搅拌速率下混合原油不同流动时间的第二混合不均匀度。
S140：在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
在一些实施例中，可以在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间。例如第一混合均匀度为x，在所述第二混合均匀度也等于x时，可以获得搅拌速率在1000r/min对应的搅拌时间，2000r/min对应的搅拌时间，从而根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
本说明书实施例提供的混合原油配置方法，可以获取工程现场油品混合设备的实际尺寸和运行参数；根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。本说明书实施例提供的混合原油配置方法，利用至少两种不同的油品，辅以仿真模拟结果，可以在实验室内配制出与工程现场混合状态相同的异种油品混合原油，避免了因实验需要多次返回现场取样的情况，进而提高混合原油的配置效率。
为了更好的理解本说明书实施例提供的混合原油配置方法。本说明书实施例还提供一个场景示例。本场景示例中，所用的软件包括SOLIDWORKS和ANSYS 19.2。
在本场景示例中，可以利用两种不同的油为基础油，运用CFD模拟。具体的，运用SOLIDWORKS等建模软件建立与工程现场实际设备尺寸一致的模型，将建立好的模型导入ICEM CFD等网格划分软件进行网格划分并输出网格文件，用Fluent等求解软件读取网格文件并根据工程现场设备运行条件设置相关求解参数进行求解计算，将求解出的数据文件由CFD Post等后处理软件进行后处理导出关键截面一种油品的体积分数计算其混合不均匀度，作为工程现场异种油品混合状态的表征参数。用相同的方法通过模拟确定实验室异种油品不同搅拌速率不同搅拌时间混合状态的表征参数，当二者的表征参数相同时，即可将工程现场与实验室异种油品的混合状态关联起来。
在本场景示例中，模拟工程现场大流量工况下稠稀油的混合状态如下：通过现场调研某油库稠油掺稀设备的实际尺寸和运行状态，可知现场油品混合设备包括两个SK型混合器和异径管道；运用SOLIDWORKS软件根据调研的尺寸建立模型；将建立的模型导入ICEMCFD中设置稠油入口、稀油入口和混油出口，划分网格并导出msh(网格)文件；将msh(网格)文件导入Fluent软件中，在General(通用)里，Solver Type(求解方式)为Pressure-Based(压力基)，SolverVelocity Formulation(求解速度方程)为Absolute(绝对的)，SolverTime(求解模式)为Transient(瞬态)，勾选Gravity(重力)，根据模型特征设置Gravitational Acceleration(重力加速度)中Z为9.81，其余为0；在Models(模型)里，打开Multiphase(多相)，勾选Mixture(混合模型)，并设置Number ofEulerian Phases(欧拉相数)为2，打开Viscous(粘性)，勾选Model(模型)里的k-epsilon(2eqn)，勾选k-epsilonModel里的Realizable(可实现)，勾选Near-Wall Treatment(近壁面处理)里的StandardWall Functions(标准壁面模式)；在Materials(材料库)里，选择fuel-oil-liquid(燃料油)并更改其密度为930kg/m3，更改其动力粘度为0.2kg/(m·s)，作为稠油相，选择gasoil-liquid(汽油)并更改其密度为850kg/m3，更改其动力粘度为0.03kg/(m·s)，作为稀油相；在Boundary Conditions(边界条件)里，设置稀油入口稀油速度为0.71m/s，稠油为0，设置稠油入口稠油速度为2.83m/s，稀油速度为0，设置出口Phase(相)为mixture(混合相)，Flow Rate Weighting(流量加权)为1；调整控制因子，初始化流场，根据网格和模型大小设置时间步长和计算步数，设置每500步保存一次计算结果并开始计算。计算完成后，将数据文件导入至CFD-Post软件中进行后处理，导出不同流动时间出口截面上稠油体积分数，计算不同流动时间出口截面上的混合不均匀度，统计如表1所示，其中，前15s油流从入口处未流至出口。
表1
| 流动时间(s) | 混合不均匀度 |
| 15 | 0.119518771 |
| 17.5 | 0.205335554 |
| 20 | 0.138616104 |
| 22.5 | 0.085294278 |
| 25 | 0.058862380 |
| 27.5 | 0.047652381 |
| 30 | 0.042864437 |
| 32.5 | 0.041517766 |
| 35 | 0.041204589 |
| 37.5 | 0.041168550 |
| 40 | 0.041647057 |
| 42.5 | 0.041678321 |
| 45 | 0.041255947 |
| 47.5 | 0.041709606 |
| 49.34 | 0.041272941 |
| 50 | 0.041274621 |
由表1可知，该工况下混合不均匀度可取为0.041204589，即根据仿真模拟结果，工程现场大流量工况下稠稀油混合的混合不均匀度可降至0.041204589，流动时间为35s。
在本场景示例中，模拟实验室稠稀油的混合状态的过程如下：在实验室调研稠油掺稀搅拌设备的实际尺寸和形状等，运用SOLIDWORKS软件根据调研的尺寸建立模型，将建立的模型导入ICEM CFD中划分网格并导出msh(网格)文件，将msh(网格)文件导入Fluent软件中，在General(通用)里，Solver Type(求解方式)为Pressure-Based(压力基)，SolverVelocity Formulation(求解速度方程)为Absolute(绝对的)，Solver Time(求解模式)为Transient(瞬态)，勾选Gravity(重力)，根据模型特征设置Gravitational Acceleration(重力加速度)中Y为9.81，其余为0；在Models(模型)里，打开Multiphase(多相)，勾选Eulerian(欧拉模型)，并设置Number ofEulerian Phases(欧拉相数)为3，打开Viscous(粘性的)，勾选Model(模型)里的k-epsilon(2eqn)，勾选k-epsilon Model里的RNG，勾选Near-Wall Treatment(近壁面处理)里的StandardWall Functions(标准壁面模式)；在Materials(材料库)里，选择fuel-oil-liquid(燃料油)并更改其密度为930kg/m3，更改其动力粘度为0.2kg/(m·s)，作为稠油相，选择gasoil-liquid(汽油)并更改其密度为850kg/m3，更改其动力粘度为0.03kg/(m·s)，作为稀油相；在Boundary Conditions(边界条件)里，设置wall1(壁面1，即叶轮壁面)为Moving Wall(动壁面)，在Motion(运动)里勾选Relative toAdjacent Cell Zone(相对于相邻网格区)，设置Speed(rpm)(转速(转/分钟))为1000，勾选Rotational(旋转的)，并设置Y(mm)为30，Y为90；为与模拟工程现场的实例作为对照，设置模型空间内下层为稠油，上层为稀油，两层密度不同的原有占据整个容器空间；选择合适的求解方法和控制因子，初始化流场，根据网格和模型大小设置时间步长和计算步数，设置每1000步保存一次计算结果并开始计算。计算完成后，将数据文件导入至CFD-Post软件中进行后处理，导出不同流动时间模型竖向截面上的稠油体积分数，计算不同流动时间出口截面上的混合不均匀度，统计如表2所示。
表2
其中，表2为在搅拌设备的搅拌速率为1000r/min搅拌工况下，混合不均匀度随搅拌时间变化数据。
在本场景示例中，可以将搅拌设备的搅拌速率改为2000r/min，重复实验室稠稀油混合的模拟，得到混合不均匀度随搅拌时间变化如表3所示。
表3
| 搅拌时间(s) | 混合不均匀度 |
| 10 | 0.484840744 |
| 20 | 0.215343605 |
| 30 | 0.104734663 |
| 40 | 0.053682443 |
| 50 | 0.028346447 |
| 60 | 0.015558499 |
| 70 | 0.008907259 |
| 80 | 0.005549907 |
| 90 | 0.003911161 |
| 100 | 0.002887609 |
| 110 | 0.002314298 |
| 120 | 0.002120863 |
模拟工程现场稠稀油混合状态后可知混合不均匀度为0.041204589，分别对照1000r/min和2000r/min时实验室搅拌的稠稀油混合状态模拟，用内插法可得1000r/min时欲使混合不均匀度降至0.041204589需要搅拌约88s，2000r/min时欲使混合不均匀度降至0.041204589需要搅拌约45s。至此，即可将实验室配制稠稀混合油模拟与工程现场稠稀混油混合状态模拟建立关联，从而可以在实验室内配制出与工程现场混合状态一致的异种油品混合原油，避免了因实验需要多次返回现场取样的情况，进而提高混合原油的配置效率。
参阅图2，在软件层面上，本说明书实施例还提供了一种混合原油配置装置，该装置具体可以包括以下的结构模块。
获取模块210，用于获取工程现场油品混合设备的实际尺寸和运行参数；
第一模拟模块220，用于根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；
第二模拟模块230，用于根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；
获得模块240，用于在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
需要说明的是，本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其它实施例的不同之处。尤其，对于装置实施例和设备实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本领域技术人员在阅读本说明书文件之后，可以无需创造性劳动想到将本说明书列举的部分或全部实施例进行任意组合，这些组合也在本说明书公开和保护的范围内。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(AlteraHardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本说明书各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本说明书，本领域普通技术人员知道，本说明书有许多变形和变化而不脱离本说明书的精神，希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。
Claims (10)
1.一种混合原油配置方法，其特征在于，所述方法包括：
获取工程现场油品混合设备的实际尺寸和运行参数；
根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；
根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；
在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
2.根据权利要求1所述的方法，其特征在于，所述工程现场油品混合设备包括SK型混合器和异径管道。
3.根据权利要求1所述的方法，其特征在于，所述至少两种预设类别的油品包括至少一种稠油和至少一种稀油。
4.根据权利要求1所述的方法，其特征在于，所述运行参数至少包括稠油入口输入稠油的速度和稀油入口输入稀油的速度。
5.根据权利要求1所述的方法，其特征在于，通过设置油品的密度和油品的动力粘度模拟所述至少两种预设类别的油品。
6.根据权利要求1所述的方法，其特征在于，所述根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度包括：
根据所述工程现场油品混合设备的实际尺寸建立第一仿真模型；
将所述运行参数作为模拟参数，使用所述第一仿真模型模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度。
7.根据权利要求1所述的方法，其特征在于，所述根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度包括：
根据配置混合原油所采用的搅拌设备的形状和尺寸建立第二仿真模型；
使用所述第二仿真模型模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度。
8.根据权利要求6或7所述的方法，其特征在于，使用SOLIDWORKS软件和ANSYS软件建立第一仿真模型和第二仿真模型。
9.根据权利要求1所述的方法，其特征在于，所述第一混合不均匀度和所述第二混合不均匀度根据以下公式计算：
其中，
表示混合不均匀度，ci表示截面上第i个点的一种油品体积分数，i＝1,2…n，
表示截面上的一种油品体积分数平均值，n表示截面上取点的个数。
10.一种混合原油配置装置，其特征在于，所述装置包括：
获取模块，用于获取工程现场油品混合设备的实际尺寸和运行参数；
第一模拟模块，用于根据所述工程现场油品混合设备的实际尺寸和运行参数，模拟至少两种预设类别的油品的混合过程，得到混合油品不同流动时间的第一混合不均匀度；所述混合不均匀度表征混合油品的混合状态；
第二模拟模块，用于根据配置混合原油所采用的搅拌设备的形状和尺寸模拟所述至少两种预设类别的油品，通过所述搅拌设备不同搅拌速率下配置混合原油的过程，得到混合原油不同流动时间的第二混合不均匀度；
获得模块，用于在所述第二混合均匀度等于所述第一混合均匀度的情况下，获得不同搅拌速率对应的搅拌时间，以便于根据不同搅拌速率对应的搅拌时间配置与工程现场混合状态相同的混合原油。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011336266.5A CN112417697A (zh) | 2020-11-25 | 2020-11-25 | 混合原油配置方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011336266.5A CN112417697A (zh) | 2020-11-25 | 2020-11-25 | 混合原油配置方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112417697A true CN112417697A (zh) | 2021-02-26 |
Family
ID=74842338
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011336266.5A Pending CN112417697A (zh) | 2020-11-25 | 2020-11-25 | 混合原油配置方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112417697A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113406274A (zh) * | 2021-05-28 | 2021-09-17 | 中国人民解放军空军研究院航空兵研究所 | 基于粘度在线检测的混合均匀度分析方法及混合装置 |
| CN114047030A (zh) * | 2021-11-02 | 2022-02-15 | 江苏三沃电子科技有限公司 | 一种助焊剂取样方法和系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101684727A (zh) * | 2008-09-28 | 2010-03-31 | 中国石油化工股份有限公司 | 超深井稠油掺稀比例确定的优化方法及其掺稀混配器 |
| CN202146741U (zh) * | 2011-07-21 | 2012-02-22 | 王顺华 | 原油动态混合器 |
| CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
| US20180348105A1 (en) * | 2016-03-22 | 2018-12-06 | Southwest Petroleum University | Direct method for manufacturing large model fractured core and maintaining original oil-water saturation |
-
2020
- 2020-11-25 CN CN202011336266.5A patent/CN112417697A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101684727A (zh) * | 2008-09-28 | 2010-03-31 | 中国石油化工股份有限公司 | 超深井稠油掺稀比例确定的优化方法及其掺稀混配器 |
| CN202146741U (zh) * | 2011-07-21 | 2012-02-22 | 王顺华 | 原油动态混合器 |
| CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
| US20180348105A1 (en) * | 2016-03-22 | 2018-12-06 | Southwest Petroleum University | Direct method for manufacturing large model fractured core and maintaining original oil-water saturation |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113406274A (zh) * | 2021-05-28 | 2021-09-17 | 中国人民解放军空军研究院航空兵研究所 | 基于粘度在线检测的混合均匀度分析方法及混合装置 |
| CN114047030A (zh) * | 2021-11-02 | 2022-02-15 | 江苏三沃电子科技有限公司 | 一种助焊剂取样方法和系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112417697A (zh) | 混合原油配置方法及装置 | |
| CN103902751B (zh) | 油藏数值模拟的边界处理方法及装置 | |
| Zhang et al. | Analysis of flow and phase interaction characteristics in a gas-liquid two-phase pump | |
| CN104573241B (zh) | 一种rp‑3航空煤油空化的数值模拟方法 | |
| Kruisbrink et al. | An SPH multi‐fluid model based on quasi buoyancy for interface stabilization up to high density ratios and realistic wave speed ratios | |
| Renard | A 2D reservoir streamtube EOR model with periodical automatic regeneration of streamlines | |
| Zou et al. | A numerical study on heat transfer performances of horizontal ground heat exchangers in ground-source heat pumps | |
| CN105893647B (zh) | 一种泡沫沥青混合料搅拌设备结构设计与优化的分析方法 | |
| Hsu et al. | New physics-based 3D hydraulic fracture model | |
| CN109446741A (zh) | 一种混输泵泵腔瞬时温度特性的建模与预测方法 | |
| Wu et al. | A New Projection‐Based Integrally Embedded Discrete Fracture Model and Its Application in Coupled Flow and Geomechanical Simulation for Fractured Reservoirs | |
| Shabani et al. | Fluid flow characterization of liquid–liquid mixing in mixer-settler | |
| CN116844669A (zh) | 碳/碳复合材料cvd装炉构型流场模拟方法、介质及终端 | |
| Hayashi et al. | Void distribution and bubble motion in bubbly flows in a 4× 4 rod bundle. Part II: numerical simulation | |
| Niu et al. | The application of REV-LBM double mesh local refinement algorithm in porous media flow simulation | |
| Cominelli et al. | Deployment of high-resolution reservoir simulator: Methodology & cases | |
| Akanji et al. | A finite-element algorithm for Stokes flow through oil and gas production tubing of uniform diameter | |
| Schlegel et al. | Openfoam-hybrid: A morphology adaptive multifield two-fluid model | |
| Deng et al. | Optimization of blade structural parameters of concrete truck mixer based on discrete element method | |
| Ye et al. | CFD and FEA-Based, 3D metal printing hybrid stage prototype on electric submersible pump ESP system for high-gas wells | |
| Larsson | Modelling of mass transfer phenomena in chemical and biochemical reactor systems using computational fluid dynamics | |
| Suarez et al. | 3D CFD Simulation of ESP Rotary Gas-Separator Performance Under Two-Phase-Flow Condition | |
| McGowen et al. | The effects of fluid preconditioning and test cell design on the measurement of dynamic fluid loss data | |
| Konyukhov et al. | COMPUTER MODELING OF HEAT AND MASS TRANSFER IN THE UNIFIED COMPLEX" OIL RESERVOIR-SYSTEM OF WELLS" | |
| Baiamonte | Discussion of “Hydraulic Model of Transition of Transient to Steady Flows in the Vadose Zone” by Yaguo Zhang, Tonglu Li, Wei Shen, and Yu Wang |
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 |