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.

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.

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技术能够可视化心腔内血流模式,对左室流出道峰值流速的测量具有高度可重复性,且与超声心动图的测量结果具有显著的一致性。

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.

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

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

Liquid-liquid two-phase flow in a wire-embedded concentric microchannel: Flow pattern and mass transfer performance

In this work, flow pattern and mass transfer of liquid-liquid two-phase flow in a wire-embedded concentric microchannel are studied using toluene-water system. Droplet flow, slug flow, oval flow and annular flow are observed in the wire-embedded concentric microchannel. The effects of embedded wires and physical properties on flow patterns are investigated. The embedded wire insert is conducive to the formation of annular flow. The flow pattern distribution regions are distinguished by the Caaq(capillary number)±We_(org)(Weber number) flow pattern map. When Weorg<0.001, slug flow is the main flow pattern, and when Weorg>0.1, annular flow is the main flow pattern. Oval flow and droplet flow are between We_(org)= 0.001-0.1, and oval flow is transformed into droplet flow with the increase of Caaq. The effect of flow rate, phase ratio, initial acetic acid concentration, insert shape and flow patterns on mass transfers are studied. Mass transfer process is enhanced under annular flow conditions, the volumetric mass transfer coefficient is up to 0.36 s^(-1) because of the high interfacial area and interface renewal rate of annular flow.

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

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

V Flow技术测量颈动脉壁面剪应力的一致性研究

目的评价血流向量成像(V Flow)技术在测量健康成年人颈动脉壁面剪应力(WSS)中的一致性。方法于2021年2月至2021年3月招募健康成年志愿者20人,由2名不同年资的超声医师使用配备V Flow功能的Mindray Resona 7超声仪和3~9 MHz线阵探头进行双侧颈动脉扫查,分别采集双侧颈总动脉远段、颈总动脉分叉至颈内动脉起始部两段动脉的动态V Flow图像,测量两侧颈总动脉远段的近心端、远心端、分叉处及颈内动脉起始部的前、后壁的WSS,使用组内相关系数(ICC)和Bland-Altman图评估组内一致性及组间一致性。结果双侧颈动脉前、后壁的4个不同节段之间WSS值差异均有统计学意义(P<0.05)。高年资超声医师2次测量结果的一致性结果显示,左侧颈总动脉远段的远心端一致性极好(ICC 0.779),右侧颈总动脉远段的近心端(ICC 0.605)、远心端(ICC 0.585)、颈内动脉起始部(ICC 0.457)、左侧颈总动脉分叉处(ICC 0.606)及颈内动脉起始部(ICC 0.702)一致性均较好;不同年资超声医师的测量结果显示,仅右侧颈总动脉分叉处(ICC 0.486)及左侧颈总动脉远段的远心端(ICC 0.576)一致性较好。结论V Flow技术可显示不同位点间颈动脉WSS之间的差异,其组内一致性较好,但在不同年资超声医师间存在一定的差异。

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.

基于Moldflow软件对薄壁件某塑料喷嘴的注塑方案优化研究

喷嘴注塑产品的质量和性能很大程度上由注塑方案和工艺参数所决定,不当的注塑方案和工艺参数会导致充填不足、熔接线过多、翘曲变形量大等问题。为了可以减少生产周期,节约生产成本,利用三维软件UG对喷嘴进行建模,基于模流分析软件Moldflow软件对喷嘴注塑产品进行分析和优化,通过比较点浇口和侧浇口两种方案的分析结果,最终确定该产品点浇口的方案更好。并通过设计正交实验,得到优化后的工艺参数。

Experimental Investigation of Regular or Wavy Two-Phase Flow in a Manifold

An experimental study was conducted to investigate the properties of stratified regular or wavy two-phase flow in two parallel separators located after a manifold.A total of 103 experiments with various gas and liquid velocity combinations in three inlet pipes were conducted,including 77 groups of outlet pipe resistance symmetry and 26 groups of outlet pipe resistance asymmetry trials.The experimental results have revealed that when the gas-liquid flow rate is low,the degree of uneven splitting is high,and“extreme”conditions are attained.When the superficial gas velocity is greater than that established in the extreme case,the direction of the liquid-phase displacement is reversed,while that of the gas remains unchanged.Thus,the degree of gas phase bias tends to be mitigated with an increase in the gas velocity,while the uneven splitting degree of liquid approaches 10%.Finally,varying the gas-phase outlet pipe resistance is shown to effectively change the gas-liquid two-phase flow distribution.

Simulation of gas–liquid two-phase flow in a flow-focusing microchannel with the lattice Boltzmann method

A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.

Comprehensive modeling of frictional pressure drop during carbon dioxide two-phase flow inside channels using intelligent and conventional methods

Environmentally friendly nature of CO_(2),associated with its safety and high efficiency,has made it a widely used working fluid in heat exchangers.Since CO_(2)has strange thermophysical features,specific models are required to estimate its two-phase characteristics,particularly frictional pressure drop(FPD).Herein,a widespread dataset,comprising 1195 experimental samples for two-phase FPD of CO_(2)was adopted from 10 sources to fulfill this requirement.The literature correlations failed to provide satisfactory precisions and exhibited the average absolute relative errors(AAREs)between 29.29% and 67.69% from the analyzed data.By inspiring the theoretical method of Lockhart and Martinelli,three intelligent FPD models were presented,among which the Gaussian process regression approach surpassed the others with AARE and R^(2)values of 5.48% and 98.80%,respectively in the test stage.A novel simple correlation was also derived based on the least square fitting method,which yielded opportune predictions with AARE of 19.76% for all data.The truthfulness of the newly proposed models was assessed through a variety of statistical and visual analyses,and the results affirmed their high reliability over a broad range of conditions,channel sizes and flow patterns.Furthermore,the novel models performed favorably in describing the physical attitudes corresponding to two-phase FPD of CO_(2).Eventually,the importance of operating factors in controlling the FPD was discussed through a sensitivity analysis.

Experiments on two-phase flow in hydraulic jump on pebbled rough bed:Part 1–Turbulence properties and particle chord time and length

This study reported and discussed turbulence characteristics,such as turbulence intensity,correlation time scales,and advective length scales.The characteristic air–water time scale,including the particle chord time and length and their probability density functions(PDFs),was investigated.The results demonstrated that turbulence intensity was relatively greater on a rough bed in the roller length,whereas further downstream,the decay rate was higher.In addition,the relationship between turbulence intensity and dimensionless bubble count rate reflected an increase in turbulence intensity associated with the number of entrained particles.Triple decomposition analysis(TDA)was performed to determine the contributions of slow and fast turbulent components.The TDA results indicated that,regardless of bed type and inflow conditions,the sum of the band-pass(T'_(u))and high-pass(T″_(u))filtered turbulence intensities was equal to the turbulence intensity of the raw signal data(T_(u)).T″_(u) highlighted a higher turbulence intensity and larger vorticities on the rough bed for an identical inflow Froude number.Additional TDA results were presented in terms of the interfacial velocity,auto-and cross-correlation time scales,and longitudinal advection length scale,with the effects of low-and high-frequency signal components on each highlighted parameter.The analysis of the air chord time indicated an increase in the proportion of small bubbles moving downstream.The second part of this research focused on the basic properties of particle grouping and clustering.

Experiments on two-phase flow in hydraulic jump on pebbled rough bed:Part 2–Bubble clustering

A survey on bubble clustering in air–water flow processes may provide significant insights into turbulent two-phaseflow.These processes have been studied in plunging jets,dropshafts,and hydraulic jumps on a smooth bed.As a first attempt,this study examined the bubble clustering process in hydraulic jumps on a pebbled rough bed using experimental data for 1.70<Fr_(1)<2.84(with Fr_(1) denoting the inflow Froude number).The basic properties of particle grouping and clustering,including the number of clusters,the dimensionless number of clusters per second,the percentage of clustered bubbles,and the number of bubbles per cluster,were analyzed based on two criteria.For both criteria,the maximum cluster count rate was greater on the rough bed than on the smooth bed,suggesting greater interactions between turbulence and bubbly flow on the rough bed.The results were consistent with the longitudinal distribution of the interfacial velocity using one of the criteria.In addition,the clustering process was analyzed using a different approach:the interparticle arrival time of bubbles.The comparison showed that the bubbly flow structure had a greater density of bubbles per unitflux on the rough bed than on the smooth bed.Bed roughness was the dominant parameter close to the jump toe.Further downstream,Fr_(1) predominated.Thus,the rate of bubble density decreased more rapidly for the hydraulic jump with the lowest Fr_(1).

Dynamic behavior of outburst two-phase flow in a coal mine T-shaped roadway:The formation of impact airflow and its disaster-causing effect

The study of the dynamic disaster mechanism of coal and gas outburst two-phase flow is crucial for improving disaster reduction and rescue ability of coal mine outburst accidents.An outburst test in a T-shaped roadway was conducted using a self-developed large-scale outburst dynamic disaster test system.We investigated the release characteristics of main energy sources in coal seam,and obtained the dynamic characteristics of outburst two-phase flow in a roadway.Additionally,we established a formation model for outburst impact flow and a model for its flow in a bifurcated structure.The results indicate that the outburst process exhibits pulse characteristics,and the rapid destruction process of coal seam and the blocking state of gas flow are the main causes of the pulse phenomenon.The outburst energy is released in stages,and the elastic potential energy is released in the vertical direction before the horizontal direction.In a straight roadway,the impact force oscillates along the roadway.With an increase in the solid–gas ratio,the two-phase flow impact force gradually increases,and the disaster range extends from the middle of the roadway to the coal seam.In the area near the coal seam,the disaster caused by the two-phase flow impact is characterized by intermittent recovery.In a bifurcated roadway,the effect of impact airflow on impact dynamic disaster is much higher than that of two-phase flow,and the impact force tends to weaken with increasing solid-gas ratio.The impact force is asymmetrically distributed;it is higher on the left of the bifurcated roadway.With an increase in the solid-gas ratio,the static pressure rapidly decreases,and the bifurcated structure accelerates the attenuation of static pressure.Moreover,secondary acceleration is observed when the shock wave moves along the T-shaped roadway,indicating that the bifurcated structure increases the shock wave velocity.

Simulation of Gas-Water Two-Phase Flow in Tight Gas Reservoirs Considering the Gas Slip Effect

A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.

基于MoldFlow太阳能储能系统下盖成型模拟及应用

以太阳能储能系统下盖模流分析为例,反复修改流道和浇口截面尺寸与长度以及主流道位置,确定流道尺寸。运用MoldFlow对设计好的流道进行充填分析,以使成型达到流动平衡,并进行冷却、翘曲分析,预测成型缺陷,提出解决措施,为模具设计提供浇注系统和冷却系统方案和注射成型工艺参数,并应用于模具设计和实际生产。

Transient pressure analysis of polymer flooding fractured wells with oil-water two-phase flow

The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.