A unified fractional flow framework for predicting the liquid holdup in two-phase pipe flows

Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper we apply the fractional flow theory to multiphase flow in pipes and present a unified modeling framework for predicting the fluid phase volume fractions over a broad range of pipe flow conditions.Compared to existing methods and correlations,this new framework provides a simple,approximate,and efficient way to estimate the phase volume fraction in two-phase pipe flow without invoking flow patterns.Notably,existing correlations for estimating phase volume fraction can be transformed and expressed under this modeling framework.Different fractional flow models are applicable to different flow conditions,and they demonstrate good agreement against experimental data within 5%errors when compared with an experimental database comprising of 2754 data groups from 14literature sources,covering various pipe geometries,flow patterns,fluid properties and flow inclinations.The gas void fraction predicted by the framework developed in this work can be used as inputs to reliably model the hydraulic and thermal behaviors of two-phase pipe flows.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Identifying the enhancement mechanism of Al/MoO_(3) reactive multilayered films on the ignition ability of semiconductor bridge using a one-dimensional gas-solid two-phase flow model

Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.

4D-Flow MRI在肥厚型心肌病左室流出道血流评估中的价值探索

目的 探索四维血流(four-dimensional flow,4D-Flow)磁共振成像(magnetic resonance imaging,MRI)技术在左心室腔内应用的可行性。材料与方法 本研究为前瞻性、横断面研究,纳入2022年8月至2023年1月于我院接受心脏MRI检查的21例肥厚型心肌病患者,采用3.0 T MRI扫描仪进行二维血流(tow-dimensional flow,2D-Flow)及4D-Flow成像,收集患者一周内进行的超声心动图检查结果。采用组内相关系数(inter-class correlation coefficient,ICC)、变异系数(coefficients of variation,COV)及Bland-Altman分析比较2D-Flow、4D-Flow评估左室流出道峰值流速的可重复性及一致性,并通过Pearson相关性分析探究二者与超声心动图测量结果的关系。结果 2D-Flow和4D-Flow观察者内/观察者间的ICC分别为0.999/0.999和0.995/0.992,COV分别为0.5%/0.5%和2.4%/2.6%。4D-Flow与超声心动图的测量结果呈中度相关,相关系数r值为0.574(P=0.006),但一致性较差,ICC为0.375(P=0.013)。2D-Flow与4D-Flow和超声心动图间无显著的一致性及相关性。结论 4D-Flow技术能够可视化心腔内血流模式,对左室流出道峰值流速的测量具有高度可重复性,且与超声心动图的测量结果具有显著的一致性。

针对无线控制环境的OpenFlow流表申请优化机制

在无线控制环境下,受限于控制通道的带宽和时延,OpenFlow交换机通过Packet-In消息申请流表的过程耗时较长,这个时段内交换机可能收到同一条流的大量后续报文,并分别产生与首包报文完全相同的Packet-In消息发往控制器,明显浪费网络资源。同时,针对每条流独立发送一个Packet-In报文的方式容易导致控制通道负载以小包为主、包头开销占比过大、网络资源利用率不高的问题。因此,针对无线控制条件下报文粒度Packet-In消息触发和离散化Packet-In消息发送机制存在的问题,提出了一种基于流粒度实现Packet-In消息触发和Packet-In消息聚合发送的方案,并通过仿真试验,验证了方案的效果。

Numerical modeling and parametric sensitivity analysis of heat transfer and two-phase oil and water flow characteristics in horizontal and inclined flowlines using OpenFOAM

Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.

基于 Moldflow 的汽车中控台框架翘曲变形分析及优化

以某汽车中控台框架为研究对象,测量试模样品发现其翘曲变形量超过了装配要求。通过Moldflow软件模拟了该塑件实际的注塑过程,翘曲变形量的模拟值与实测平均值的最大误差为5.98%,发现该塑件翘曲变形的主要因素为冷却不均和收缩不均。本文在原物料中添加质量分数为25%的玻璃纤维以及优化工艺参数后,翘曲变形量的模拟值与初始方案相比降低了86.22%。试模验证表明,优化后的翘曲变形量模拟值与实测平均值的最大误差为4.35%,证明了Moldflow软件模拟分析的准确性。试模后各检测点的最大翘曲变形量降到了1.6 mm以下,较优化之前降低了80%以上,为类似大型复杂注塑件的翘曲变形分析及优化提供了思路。

基于Flow Simulation的某发动机涡轮压气机流场与效率分析

涡轮机内部流场对涡轮增压器的性能和效率有着重要影响,采用SolidWorks Flow Simulation模块对某发动机涡轮压气机侧流场和压气效率分析。在六种不同空气体积流量工况下,体积流量为0.29时,压气机效率最高,达到80%左右。模型的建模和流体分析均在SolidWorks环境下,分析效率高,为涡轮增压器设计和优化提供了支撑。

Numerical Simulation of a Two-Phase Flow with Low Permeability anda Start-Up Pressure Gradient

A new numerical model for low-permeability reservoirs is developed.The model incorporates the nonlinear characteristics of oil-water two-phase flows while taking into account the initiation pressure gradient.Related numerical solutions are obtained using a finite difference method.The correctness of the method is demonstrated using a two-dimensional inhomogeneous low permeability example.Then,the differences in the cumulative oil and water production are investigated for different starting water saturations.It is shown that when the initial water saturation grows,the water content of the block continues to rise and the cumulative oil production gradually decreases.

iFlow彩色血流编码成像技术在下肢动脉硬化闭塞症诊断中的应用价值

目的研究iFlow彩色血流编码成像技术在下肢动脉硬化闭塞症(LEASO)诊断中的应用价值。方法选择2022年3月至2023年10月期间确诊的106例LEASO患者作为本研究的LEASO组,以一般资料与LEASO组匹配且无动脉病变的80例志愿者作为对照组。两组受试者均进行数字减影血管造影(DSA)并采用iFlow彩色血流编码成像技术检测股骨头区域和踝关节区域达峰时间(TTP)、计算踝关节区域与股骨头区域TTP的差值,测量踝肱指数(ABI)。结果两组研究对象年龄、性别、体质量指数、吸烟史、高血压病史、糖尿病史、冠心病史、股骨头区域TTP的比较,差异无统计学意义(P>0.05);LEASO组踝关节区域TTP及TTP差值均高于对照组,差异有统计学意义(P<0.05);LEASO组中不同Rutherford分类患者股骨头区域TTP的比较以及左侧病变患者与右侧病变患者股骨头区域TTP、踝关节区域TTP、TTP差值的比较,差异无统计学意义(P>0.05),Rutherford分类越高,踝关节区域TTP及TTP差值越低(P<0.05);经Pearson检验,LEASO患者的踝关节区域TTP、TTP差值与ABI呈负相关(P<0.05);经受试者工作特征(ROC)曲线分析,踝关节区域TTP、TTP差值对LEASO具有诊断效能;经Delong检验,TTP差值诊断的ROC曲线下面积高于踝关节区域TTP(P<0.05)。结论iFlow彩色血流编码成像技术测定踝关节区TTP及TTP差值是诊断LEASO的量化指标。

袋型阻尼密封动力学特性双控制体Bulk Flow模型

为快速准确预测袋型阻尼密封泄漏特性和动力学特性,针对传统单控制体Bulk Flow模型预测精度低、无法预测交叉动力系数的问题,提出了袋型阻尼密封双控制体Bulk Flow模型和动力学特性数值预测方法,并开发了计算程序。首先,依据边界层理论,将袋型密封腔室划分为两个控制体,推导了控制体的连续性、周向动量和能量方程,引入Swamee-Jain和Takahashi方程,计算流体-壁面间和流体-流体间的周向黏性摩擦力;其次,采用牛顿-拉夫森算法和摄动分析法分别求解0阶和1阶控制方程,获得各刚度、阻尼动力特性系数;然后,通过与袋型阻尼密封泄漏量和动力特性系数的实验值、单控制体Bulk Flow模型和非定常计算流体动力学(CFD)数值结果进行比较,验证了模型和方法的准确性和可靠性;最后,研究了转子转速(10 000、15 000、20 000 r/min)和预旋比(0.067、0.724、0.997)对袋型阻尼密封动力学特性的影响。结果表明:所发展的模型和方法具有计算速度快、预测精度高(泄漏量预测误差小于6%,动力特性系数预测误差小于38%)的优点;转子转速和进口预旋的增大均会导致袋型阻尼密封有效阻尼显著减小,穿越频率显著增大,易诱发轴系失稳。

基于Flow-3D的桩林坝不同布置形式对坝基的冲刷侵蚀

桩林坝是一种以桩为基础的竖向格子坝,既有良好的受力条件,又相对节约材料,同时拥有“拦粗排细”的功效,现有桩林坝的研究,多在冲击荷载、拦截固源和调控泥石流性能等方面构建不同布置形式的桩林坝模型。以杨家沟为例,通过野外现场踏勘和数值仿真计算,基于Flow-3D软件开展数值模拟研究,建立不同布置形式(单排桩、双排桩、桩间距、排间距)的桩林坝模型,探讨不同影响因素下桩林坝坝基冲刷侵蚀规律。结果表明:桩林坝的整体设计中桩林不宜过密,桩间距和排间距以2.5 m为宜;单排桩的稳定性较弱,在泥石流冲击下极易发生剪切破坏,且对上游的“固源”作用远不如多排桩,但单排桩的冲刷深度及侵蚀规模均小于多排桩。研究成果对桩林坝的布置形式及冲刷机制提供数据支撑,并为震后泥石流防治工程设计提供参考。

肝硬化患者4D flow MRI血流动力学参数与中医证型的相关性研究

目的探讨肝硬化患者4D flow MRI血流动力学参数与中医证型的相关性。方法将118例肝硬化患者依据中医辨证分为肝气郁结证、湿热蕴结证、肝肾阴虚证、脾肾阳虚证、瘀血阻络证5个证型,所有患者均行门静脉4D flow MRI检查,统计不同证型肝硬化患者分布情况,观察门静脉系统(门静脉主干、脾静脉、肠系膜上静脉)的血流动力学参数,包括血流量、流速、壁剪切力等,比较不同证型患者门静脉血流动力学参数差异。结果肝硬化代偿期以肝气郁结证、湿热蕴结证为主,肝硬化失代偿期以脾肾阳虚、瘀血络阻证为主;A级以肝气郁结证、湿热蕴结证为主,B、C级以瘀血络阻证为主。瘀血络阻证肝硬化患者门静脉主干及脾静脉血流量明显高于肝气郁结证、湿热蕴结证患者(P<0.05);脾肾阳虚证门静脉主干血流量明显高于肝气郁结证肝硬化患者(P<0.05);瘀血络阻证肝硬化患者门静脉主干流速及剪切力较肝气郁结证和湿热蕴结证低。结论肝硬化患者中医辨证分型与门静脉血流动力学参数具有一定相关性,4D flow MRI可为肝硬化的中医辨证提供血流动力学参考。

Numerical Simulation of Oil and Gas Two-Phase Flow in Deep Condensate Gas Reservoirs in Bohai Buried Hills

The BZ19-6 gas field is characterized by high temperature and high pressure (HTHP), high condensate content, little difference between the formation pressure and dew point pressure, and large amount of reverse condensate liquid. During the early stage of depletion development, the production gas-oil ratio (GOR) and production capacity remain relatively stable, which is inconsistent with the conventional reverse condensate seepage law. In view of the static and dynamic conflict in development and production, indoor high-temperature and high-pressure PVT experiment was carried out to reveal the mist-like condensation phenomenon of fluids in the BZ19-6 formation. And the seepage characteristics of condensate gas reservoirs with various degrees of depletion under the condition of HTHP were analyzed based on production performance. The change rule of fluid phase state was analyzed in response to the characterization difficulties of the seepage mechanism. The fluid state was described using the miscible mechanism. And the interphase permeability interpolation coefficient was introduced based on interfacial tension. By doing so, the accurate characterization of the “single-phase flow of condensate gas-near-miscible mist-like quasi single-phase flow-oil-gas two-phase flow” during the development process was achieved. Then the accurate fitting of key indicators for oilfield development was completed, and the distribution law of formation pressure and the law of condensate oil precipitation under different reservoir conditions are obtained. Based on research results, the regulation strategy of variable flow rate production was developed. Currently, the work system has been optimized for 11 wells, achieving a “zero increase” in the GOS of the gas field and an annual oil increase of 22,000 cubic meters.

Beautiful Flow Plus与FUjiⅨ玻璃离子水门汀在牙颈部楔状缺损治疗中的应用效果比较

目的:比较Beautiful Flow Plus与FUjiⅨ玻璃离子水门汀在牙颈部楔状缺损治疗中的应用效果。方法:选取牙颈部楔状缺损患者198例(567颗患牙),根据修复时充填材料不同将患者分为观察组101例289颗患牙(Beautiful Flow Plus修复治疗)和对照组97例278颗患牙(FUjiⅨ玻璃离子水门汀修复治疗)。比较两组患者治疗后1、2年填充体状态及修复成功率。比较两组患者治疗后2年咀嚼功能、美观度评分。记录治疗期间和治疗后2年内患者不良事件发生情况。结果:两组治疗后1年患牙填充体状态良好率及修复成功率比较差异无统计学意义(均P>0.05)。治疗后2年,观察组色泽协调性良好及无边缘着色的牙齿比例高于对照组(均P<0.05)。两组其余充填体状态良好率及修复成功率比较差异无统计学意义(均P>0.05)。治疗后2年,两组干燥咀嚼物X50数值较治疗前降低,且观察组低于对照组(均P<0.05)。治疗后2年,两组患者和医生美观度评分较治疗前升高,且观察组高于对照组(均P<0.05)。在治疗期间和治疗后2年内,两组不良事件总发生率比较差异无统计学意义(P<0.05)。结论:与FUjiⅨ玻璃离子水门汀比较,Beautiful Flow Plus在牙颈部楔状缺损治疗中具有较好的远期美学效果,且有助于改善患者咀嚼功能。