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.

Wax deposition modeling in oil-water stratified pipe flow

Wax deposition in oil-water stratified flow is commonly encountered onshore and offshore oil production pipe systems,and typically reduces transportation capacity of oil.The accurate predicted model of wax deposition has becomes an indispensable approach to design effective remediation strategies.However,a reliable mechanistic model for wax deposition prediction in oil-water two-phase stratified pipe flow is lacking to validate the deposition process.In this work,a three-dimensional(axial,radial,and angular)robust wax deposit model for oil-water stratified circular pipe flow was developed.The model of formation of a gel deposit based on the first principles of rheology was developed,associated with the results obtained from hydrodynamics and heat/mass transfer simulations.The predictions for wax deposition are found to compare satisfactorily with experimental data with two different oils for single phase and four different water cuts for oil-water stratified pipe flow.It can be seen from the wax gelation mechanism that an increase in water cut can help to reduce the wall/oil-deposit interface shear stress,thereby leading to an increase in the degree of gelation as well as the deposit rate.Furthermore,a local deposit analysis in the circumferential direction was conducted,for water cut 75%and total flow rate 5 m3/h,which provided insights to understand that the thickness on pipe wall was roughly uniformly distributed locates near the top of the pipe and the nearer the position gets close to two points,where the oil-water interface contacts the inner wall,the deposition thickness quickly dropped to 0.It was attributed to the fact that a roughly uniformly thickness far away from the oil-water interface contact the inner wall resulted in the slowly changes temperature along the circumferential pipe wall wetted by oil.

Numerical Analysis of Heavy Oil-Water Flow and Leak Detection in Vertical Pipeline

Pipeline is a conventional, efficient and economic way for oil transportations. The use of a good system for detecting and locating leaks in pipeline contribute significantly to operational safety and cost saving in petroleum industry. This paper aims to study the heavy oil-water flow in vertical ducts including leakage. A transient numerical analysis, using the ANSYS-CFX? 11.0 commercial software is performed. The mathematical modeling considers the effect of drag and gravitational forces between the phases and turbulent flow. Mass flow rate of the phases in the leaking orifice, the pressure drop as a function of the time and the velocity distributions are presented and discussed. We can conclude that volumetric fraction of phases and fluid mixture velocity affect pressure drop and mass flow rate at the leak hole.

Ultrasonic method for measuring water holdup of low velocity and high-water-cut oil-water two-phase flow

Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration（of oil droplets） in oil-water two-phase flow,which makes it difficult to measure water holdup in oil wells.In this study,we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in lowvelocity and high water-cut conditions.First,we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling.Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor.Based on the results,we then investigate the effects of oildroplet diameter and distribution on the ultrasonic field.To further understand the measurement characteristics of the ultrasonic sensor,we perform a flow loop test on vertical upward oilwater two-phase flow and measure the responses of the optimized ultrasonic sensor.The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow（D OS/W flow）,but the resolution is favorable for dispersed oil in water flow（D O/W flow） and very fine dispersed oil in water flow（VFD O/W flow）.This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.

Enhancing Heavy Crude Oil Flow in Pipelines through Heating-Induced Viscosity Reduction in the Petroleum Industry

The process of transporting crude oil across pipelines is one of the most critical aspects of the midstream petroleum industry.In the present experimental work,the effect of temperature,pressure drop,and pipe diameter on the flow rate of heavy crude oil have been assessed.Moreover,the total discharge and energy losses have been evaluated in order to demonstrate the improvements potentially achievable by using solar heating method replacing pipe,and adjusting the value of the initial pressure difference.Crude oil of API=20 has been used for the experiments,with the studied pipelines sections connecting the separator unit to the storage tank operating at a temperature of 25℃-100℃,pressure drop of 3,4,5,and 6 kg/cm^(2),and with pipe diameter of 4,6,and 8 in.The results show that on increasing the temperature and/or the pressure drop,the flow rate through the pipeline becomes higher,thus raising the total pumping energy(as the pipe diameter increase),while energy losses increase from the last separator to the storage tank in the field.A pipe diameter increase can also produce a growth of the total pumping energy(i.e.,energy losses increase).The results of the present analysis suggest that employing an optimal temperature(50℃)is needed to ensure good performance.

Partial phase flow rate measurements for stratified oil-water flow in horizontal wells

To accurately measure and evaluate the oil-water production profile of horizontal wells, a dynamic measurement experiment of oil-water two-phase flow in horizontal wells and numerical simulation were combined to establish a method for measuring the partial phase flow rate of oil-water two-phase stratified flow in horizontal wells. An experimental work was performed in horizontal oil-water two-phase flow simulation well using combination production logging tool including mini-capacitance sensor and mini-spinner. The combination tool provides a recording of holdup and velocity profiles at five different heights of the borehole cross-section. The effect of total flow rate and water-cut on the response of spinner and capacitive sensor at five measured positions were investigated. The capacitance water holdup interpolation imaging algorithm was used to determine the local fluid property and oil-water interface height, and the measured local fluid speed was combined with the numerical simulation result to establish an optimal calculation model for obtaining the partial phase flow rate of the oil-water two-phase stratified flow in the horizontal well. The calculated flow rates of five measured points are basically consistent with the experimental data, the total flow rate and water holdup from calculation are in agreement with the set values in the experiment too, suggesting that the method has high accuracy.

Simulation of Oil-Water Flow in Shale Oil Reservoirs Based on Smooth Particle Hydrodynamics

A Smooth Particle Hydrodynamics(SPH)method is employed to simulate the two-phase flow of oil and water in a reservoir.It is shown that,in comparison to the classical finite difference approach,this method is more stable and effective at capturing the complex evolution of this category of two-phase flows.The influence of several smooth functions is explored and it is concluded that the Gaussian function is the best one.After 200 days,the block water cutoff for the Gaussian function is 0.3,whereas the other functions have a block water cutoff of 0.8.The effect of various injection ratios on real reservoir production is explored.When 14 and 8 m^(3)/day is employed,the water breakthrough time is 130 and 170 days,respectively,and the block produces 9246 m^(3) and 6338 m^(3) of oil cumulatively over 400 days.

Multifractal analysis of inclined oil-water countercurrent flow

Characterizing countercurrent flow structures in an inclined oil-water two-phase flow from one-dimensional measurement is of great importance for model building and sensor design.Firstly,we conducted oil-water two-phase flow experiments in an inclined pipe to measure the conductance signals of three typical water-dominated oil-water flow patterns in inclined flow,i.e.,dispersion oil-in-water pseudo-slug flow （PS）,dispersion oil-in-water countercurrent flow （CT）,and transitional flow （TF）.In pseudo-slug flow,countercurrent flow and transitional flow,oil is completely dispersed in water.Then we used magnitude and sign decomposition analysis and multifractal analysis to reveal levels of complexity in different flow patterns.We found that the PS and CT flow patterns both exhibited high complexity and obvious multifractal dynamic behavior,but the magnitude scaling exponent and singularity of the CT flow pattern were less than those of the PS flow pattern; and the TF flow pattern exhibited low complexity and almost monofractal behavior,and its magnitude scaling was close to random behavior.Meanwhile,at short time scales,all sign series of two-phase flow patterns exhibited very similar strong positive correlation; at high time scales,the scaling analysis of sign series showed different anti-correlated behavior.Furthermore,with an increase in oil flow rate,the flow structure became regular,which could be reflected by the decrease in the width of spectrum and the difference in dimensions.The results suggested that different oil-water flow patterns exhibited different nonlinear features,and the varying levels of complexity could well characterize the fluid dynamics underlying different oil-water flow patterns.

Simulation of Oil-Water Flow in a Shale Reservoir Using a Radial Basis Function

Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.

Momentum-Dependent Symmetry Potential from Nuclear Collective Flows in Heavy Ion Collisions at Intermediate Energies

Within the isospin-dependent quantum molecular dynamics model, we investigate the nuclear collective flows produced in semi-central 197 Au＋197 Au collisions at intermediate energies. The neutron proton differential flows and difference of neutron proton collective flows are sensitive to the momentum-dependent symmetry potential. This sensitivity is less affected by both the isoscalar part of nuclear equation of state and in-medium nucleon- nucleon cross sections. Moreover, this sensitivity becomes pronounced with increasing the rapidity cut.

Heavy Metals Removal from Swine Wastewater Using Constructed Wetlands with Horizontal Sub-Surface Flow

The removal efficiency of Cu and Zn from swine wastewater was evaluated as effected by three variables: the hydraulic retention time (HRT) (24, 48, 72 and 96 hours), two different plant species (Typha domingensis Pers. and Eleocharis cellulosa) and two different sizes of filter media (5 and 15 mm) using a horizontal sub-surface flow constructed wetland. From the results, a significant difference was observed in the removal efficiency of Cu and Zn with respect to different hydraulic retention times. The best results were obtained in the HRT of 96 hours for Zn where 96% removal of Zn with Typha domingensis Pers. specie with gravel of 15 mm (experimental unit 6) was achieved. For Cu, at 72 hours of HRT, the efficiency was nearly 100% in five of the six study units (1, 2, 3, 5 and 6). In contrast, in experimental unit 4 with gravel of 15 mm and without plants, only 86% Cu removal was achieved.

Heavy-organic particle deposition from petroleum fluid flow in oil wells and pipelines

Suspended asphaltenic heavy organic particles in petroleum fluids may stick to the inner walls of oil wells and pipelines. This is the major reason for fouling and arterial blockage in the petroleum industry. This report is devoted the study of the mechanism of migration of suspended heavy organic particles towards the walls in oil-producing wells and pipelines. In this report we present a detailed analytical model for the heavy organics suspended particle deposition coefficient corresponding to petroleum fluids flow production conditions in oil wells. We predict the rate of particle deposition during various turbulent flow regimes. The turbulent boundary layer theory and the concepts of mass transfer are utilized to model and calculate the particle deposition rates on the walls of flowing conduits. The developed model accounts for the eddy diffusivity, and Brownian diffusivity as well as for inertial effects. The analysis presented in this paper shows that rates of particle deposition （during petroleum fluid production） on the walls of the flowing channel due solely to diffusion effects are small. It is also shown that deposition rates decrease with increasing particle size. However, when the process is momentum controlled （large particle sizes） higher deposition rates are expected.

Influence of Chiral Mean Field on Kaon In-plane Flow in Heavy Ion Collisions

The influence of the chiral mean field on the in-plane flow in heavy ion collisions at SIS energy is investigated within covariant kaon dynamics. For the kaon mesons inside the nuclear medium a quasi-particle picture including scalar and vector fields is adopted and compared to the standard treatment with a static potential. It is confirmed that a Lorentz force from spatial component of the vector field provides an important contribution to the in-medium kaon dynamics and strongly counterbalances the influence of the vector potential on the in-plane flow. The calculated results show that the new FOPI data can be reasonably described using the Brown & Rho parametrization, which partly takes into account the correction of higher order contributions in the chiral expansion. This indicates that one can abstract the information on the kaon potential in a nuclear medium from the analysis of the in-plane flow.

Thermo-Hydrodynamics of Core-Annular Flow of Water, Heavy Oil and Air Using CFX

The transport of heavy and ultra-viscous oil employing the core-flow technique has been increasing recently, because it provides a greater reduction of the pressure drop during the flow. In this context, the effect of temperature and the presence of gas on the thermo-hydrodynamics of a three-phase water-heavy oil-air flow in a horizontal pipe under the influence of gravity and drag forces, using the commercial software ANSYS CFX?, have been evaluated. The standard κ ? ε turbulence model, the mixture model for heavy oil-water system and the particle model for heavy oil-gas and water-gas systems, were adopted. Results of velocity, volume fraction, pressure and temperature fields of the phases present along the pipe are presented and discussed. It has been found that the presence of the air phase and the variation in the temperature affect the behavior of annular flow and pressure drop.

Prediction of pressure gradient and hold-up in horizontal liquid-liquid pipe flow

This paper aims to propose correlations to predict pressure gradient,friction factor and fluid phase hold-up in liquid-liquid horizontal pipe flow.To develop the correlations,experiments are conducted using high viscous oils(202 and 630 mPa⋅s)in a steel pipe of length 11.25 m and length-to-diameter ratio of 708.In addition,the experimental data from the literature comprising wide range of flow and fluid properties is analyzed.For the analysis,the liquid-liquid pipe flow data is categorized into two as:stratified and dispersed.The existing friction factor correlations are modified to incorporate the effects of viscosity of the oil phase,interfacial curvature(contact/wetting angle-in lieu of material of the pipe)and fluid phase fraction.In the two-fluid model of stratified flow,the wall stress and interfacial stress correlations are substituted with superficial velocities of fluids and superficial Reynolds numbers of fluid phases replacing fluid phase velocities and fluid Reynolds numbers.Similarly,for dispersed flow,an effective Reynolds number is described as the sum of superficial Reynolds number of oil and water phases.Substituting the generally employed mean or mixture Reynolds number with the effective Reynolds number into the existing single-phase turbulent flow friction factor correlation,an effective friction factor for oil-water flow is proposed.Employing the proposed correlations,the pressure gradient across the oil-water flow and hold-up volume fraction are predicted with significant reduction in error compared with that of conventionally employed correlations.The average error and standard deviation values of−7.06%,20.72%and 0.31%,18.79%are found for stratified flow and dispersed flow respectively.

倾斜管不同黏度气液两相段塞流持液率实验及预测方法

准确预测持液率可为井筒中的流型判断和压降预测提供重要依据,段塞流为稠油井筒中最常见流型,井筒高黏流体将加重气液两相滑脱,导致现有持液率模型应用于高黏流体的预测精度变差。为此,基于多相管流实验平台,在内径为60 mm的测试管段内开展段塞流持液率实验,获取了不同黏度、不同倾角段塞流流型及持液率数据,研究了黏度对持液率和流型转变的影响规律,并结合实验现象和理论推导,提出了倾斜管不同黏度气液两相段塞流持液率新模型。研究结果表明:液相黏度的增加,会增加液相与管壁之间的黏滞力,导致持液率增加;黏度对持液率的影响会改变段塞流与其他流型的转变界限。以Kora液塞区持液率关系式为基础,建立了用混相黏度代替液相黏度的倾斜管气液两相段塞流持液率新模型,并利用实验数据和文献数据进行了验证,证实新模型具有更高的精度。该研究可为预测稠油井井筒压降提供技术支持。

超稠油油藏热化学驱油水渗流特征

针对超稠油油藏热化学驱开发过程中不同温度区域内的油水渗流特征不明的问题,利用微观可视化实验和一维物理模拟实验,定量研究了不同温度下热水和驱油剂对驱油效率的影响及相对渗透率的变化规律,分析了热水和驱油剂驱油的致效机理和交互作用。实验结果表明:温度为70℃时,高温驱油剂驱的油相相对渗透率增大,水相相对渗透率变化较小;温度为150℃时,热水和驱油剂的协同增效作用更显著,热水驱转高温驱油剂驱和直接高温驱油剂驱的油相、水相相对渗透率均明显增大;温度超过200℃后,驱油剂在高温限制下驱油作用减弱,热水对驱油效率的提升大幅增加。研究表明:不同温度下,热水驱和高温驱油剂驱均可提高驱油效率;随着温度升高,热水对提高驱油效率的作用不断增大,驱油剂对驱油效率的贡献先增大后减小;热化学驱通过热水、驱油剂在不同温度区域的接替驱油和协同作用,能够实现超稠油油藏效益开发。该研究可为热化学驱提高超稠油油藏采收率提供参考。

基于试验模拟的滑坡泥石流灾害链风险监测预警

[目的]开展以元阳梯田滑坡泥石流地质灾害链为基础的风险监测预警研究,为改善降雨诱发型滑坡泥石流灾害链的预警精度,提高防灾减灾工作针对性提供科学依据。[方法]采用降雨模型模拟试验来探讨元阳梯田滑坡—泥石流灾害链破坏过程,通过对降雨量、体积含水率和坡体变形破坏过程进行实时监测,分析不同降雨条件下灾害链风险预警阈值。试验设置长时间中雨(50 mm/h)、长时间暴雨(120 mm/h)、循环间歇暴雨和干旱突发型暴雨4种降雨工况。[结果]在降雨持续作用下,灾害链发展过程为:地表冲刷、蠕滑、快速滑动、泥石流运移及运移结束。坡体破坏最先从左侧坡脚开始,右侧坡脚稍后。不同降雨工况下坡体破坏形式不同:暴雨工况下坡体中部会产生一条横向裂缝,而形成上下两个滑块;中雨工况下坡体则陆续形成多个滑块;干旱突发暴雨工况下泥石流形成时间最短,运移距离最远,流通和堆积范围相比其他降雨工况更广。[结论]滑坡泥石流灾害链预警监测指标主要为累积降雨量、降雨持续时间、含水率变化和坡体变形发展现场监测。这些监测参数受降雨条件、土体性质、坡角和外部动力条件等诸多因素影响,因此在实际应用中需要综合考虑这些因素。

外源秸秆对污染土壤氧化还原过程水分散性胶体态重金属的影响

土壤胶体尤其是水分散性胶体,作为一种重要的污染物载体,在重金属吸附、迁移以及生物吸收过程起着重要作用。通过土壤培养试验探讨了外源添加水稻秸秆在氧化还原波动条件下如何影响土壤液相中及水分散性胶体中重金属的分布。结果表明,外源秸秆增加厌氧过程液相中溶解性有机碳(DOC)、砷(As)、铁(Fe)、锰(Mn)、钙(Ca)、钾(K)、硅(Si)、铝(Al)、镁(Mg)等浓度,降低了厌氧过程氧化还原电位(Eh)及铜(Cu)、铅(Pb)等浓度,提高了好氧过程Pb浓度。利用非对称流场流分馏-紫外可见-电感耦合等离子体-质谱联用技术(AF4-UV-ICP-MS),测得水分散性胶体颗粒主要分布在0.3~3 kDa、3~40 kDa和130 kDa~450 nm三个粒径范围,各粒径颗粒的组成有所差异,主要含有机质、无机黏土矿物和铁矿物等。外源秸秆促进液相中Fe和As由胶体态向溶解态转化,促进镉(Cd)和Cu由溶解态向胶体态转化。本研究有助于揭示农业活动影响重金属迁移转化及有效性的界面机制。

某工程泄水建筑物泄洪雾化范围计算分析

为了减少泄洪雾化带来的工程影响和损失,笔者结合工程实例,根据理论公式对表孔溢洪洞、泄洪冲沙洞在不同水位下的雾化范围进行估算,并依据估算结果提出针对性工程设计建议。结果表明:表孔溢洪洞在校核洪水位下雾化影响范围最大;该工程的厂区虽位于强暴雨区横向影响范围以外,但处于雾流降雨区范围之内,设计时需考虑雾化对厂区电气设备带来的不利影响;出口边坡局部位于强暴雨区和雾流降雨区范围内,需对边坡进行加固处理并采取表面防护措施,确保边坡运行安全稳定。计算成果可为工程布置、设计提供理论支撑。