Evaluation of gas wettability and its effects on fluid distribution and fluid flow in porous media

The special gas wettability phenomenon of reservoir rocks has been recognized by more and more researchers.It has a significant effect on efficient development of unconventional reservoirs.First,based on the preferentially gas-covered ability and surface free energy changes,definition and evaluation methods have been established.Second,a method for altering rock wettability and its mechanisms have been studied,surface oriented phenomena of functional groups with low surface energy are the fundamental reason for gas wettability alteration of rock.Third,the effect of gas wettability on the surface energy,electrical properties and dilatability are investigated.Last,the effects of gas wettability on capillary pressure,oil/gas/water distribution and flow are investigated with capillary tubes and etchedglass network models.The gas wettability theory of reservoir rocks has been initially established,which provides theoretical support for the efficient production of unconventional reservoirs and has great significance.

LBM simulation of fluid flow in fractured porous media with permeable matrix

To analyze and depict complicated fluid behaviors in fractured porous media with variably permeable matrix,an integrated discrete computational algorithm is proposed based on lattice Boltzmann method（LBM）.This paper combines with the external force model and statistical material physics to effectively describe the feature changes while the fluid passes through the fractures within the permeable matrix.As an application example,a two dimensional rock sample is reconstructed using the digital image and characterized with different feature values at each LBM grid to distinguish pores,impermeable and permeable matrix by stating its local physical property.Compared with the conventional LBM,the results demonstrate the advantages of proposed algorithm in modeling fluid flow phenomenon in fractured porous media with variably permeable matrix.

Forward modeling of ocean-bottom cable data and wave-mode separation in fluid–solid elastic media with irregular seabed

In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.

Supercritical Fluid Media in Challenges of Substance and Material Dispersion

The paper deals with the problem of material dispersion using supercritical fluid media. At the same time, emphasis is made on modifications(SAS, GAS, SEDS and ASES) of the supercritical fluid anti-solvent method of dispersion. The results of SAS method implementation for dispersion of pure polycarbonate and polycarbonate doped with "CdSe/CdS-core/shell" quantum dots(carried out in the pressure range of 8.0-25.0 MPa at temperatures of 313.15 K and 358.15 K) are submitted. The range of the operating parameters has been established through the example of pure polycarbonate dispersion, which provides the production of nanoparticles with the size range of 10-100 nm. Encapsulation of Cd Se/Cd S quantum dots into polycarbonate using the SAS method has no effect on optical properties of the encapsulated quantum dots. The results of paracetamol dispersion using the SEDS method are presented. The effect of operating conditions of the paracetamol dispersion process on morphology of the obtained product is described. Co-dispersion of ethylene–vinyl acetate copolymers and low-density polyethylene mixtures by SEDS method has been carried out under pressures of 8.0-25.0 MPa at temperatures of 313 K, 323 K, and 333 K. The comparison of melting and crystallization between the resulting copolymer mixtures and mixtures with the same composition obtained by mixing in the liquid melt, has shown that implementation of SEDS results in an increase of crystallinity degree of the polymer mixtures.

Numerical simulation for separation of multi-phase immiscible fluids in porous media

Based on a lattice Boltzmann method and general principles of porous flow, a numerical technique is presented for analysing the separation of multi-phase immiscible fluids in porous media. The total body force acting on fluid particles is modified by axiding relative permeability in Nithiarasu＇s expression with an axiditional surface tension term. As a test of this model, we simulate the phase separation for the case of two immiscible fluids. The numerical results show that the two coupling relative permeability coefficients K12 and K21 have the same magnitude, so the linear flux-forcing relationships satisfy Onsager reciprocity. Phase separation phenomenon is shown with the time evolution of density distribution and bears a strong similarity to the results obtained from other numerical models and the flows in sands. At the same time, the dynamical rules in this model are local, therefore it can be run on massively parallel computers with well computational efficiency.

Hydrothermal fluid circulation in anisotropic permeable media associated with discrete fractures

A proper form of the Rayleigh number, containing the geometric mean of the vertical and horizontal permeabilities was obtained. The critical value for the onset of stable convection was found. The results proved analytically and numerically that anisotropy in permeability resists the initiation of hydrothermal convection. The equivalence between homogeneously anisotropic media and multiply fractured media was also investigated. It was confirmed that multiply fractured models are comparable to anisotropic models as long as they have the same averaged horizontal or vertical permeabilities and other physical parameters.

PARAMETERS INVERSION OF FLUID-SATURATED POROUS MEDIA BASED ON NEURAL NETWORKS

The multi- layers feedforward neural network is used for inversion ofmaterial constants of fluid-saturated porous media. The direct analysis of fluid-saturated porousmedia is carried out with the boundary element method. The dynamic displacement responses obtainedfrom direct analysis for prescribed material parameters constitute the sample sets training neuralnetwork. By virtue of the effective L-M training algorithm and the Tikhonov regularization method aswell as the GCV method for an appropriate selection of regu-larization parameter, the inversemapping from dynamic displacement responses to material constants is performed. Numerical examplesdemonstrate the validity of the neural network method.

PENALTY FINITE ELEMENT METHOD FOR NONLINEAR DYNAMIC RESPONSE OF VISCOUS FLUID-SATURATED BIPHASE POROUS MEDIA

The governing equations as well as boundary land initial conditions for nonlinear dynamic response problems of viscous fluid-saturated biphase porous medium model, based on mixture theory, are presented. With Galerkin weighted residual method the corresponding nonlinear dynamic penalty finite element equation, in which the dependencies of volume fraction and permeation coefficients an deformation are included, is obtained. The iteration solution method of the nonlinear system equation is also discussed. As a numerical example, the dynamic response of a porous medium column under impulsive loading action is analyzed with the developed finite element program. The numerical results demonstrate the efficiency and correctness of the method.

Enduring effect of permeability texture for enhancing accuracy and reducing uncertainty of reservoir fluid flow through porous media

Modeling reservoir permeability is one of the crucial tasks in reservoir simulation studies.Traditionally,it is done by kriging-based methods.More rigorous modeling of the permeability results in more reliable outputs of the reservoir models.Recently,a new category of geostatistical methods has been used for this purpose,namely multiple point statistics(MPS).By this new category of permeability modeling methods,one is able to predict the heterogeneity of the reservoir permeability as a continuous variable.These methods consider the direction of property variation in addition to the distances of known locations of the property.In this study,the reservoir performance of a modified version of the SPE 10 solution project as a pioneer case is used for investigating the efficiency of these methods and paralleling them with the kriging-based one.In this way,the permeability texture concept is introduced by applying some MPS methods.This study is accomplished in the conditions of real reservoir dimensions and velocities for the whole reservoir life.A continuous training image is used as the input of calculation for the permeability modeling.The results show that the detailed permeability of the reservoir as a continuous variable makes the reservoir simulation show the same fluid front movement and flooding behavior of the reservoir similar to the reference case with the same permeability heterogeneity.Some MPS methods enable the reservoir simulation to reproduce the fluid flow complexities such as bypassing and oil trapping during water flooding similar to the reference case.Accordingly,total oil production is predicted with higher accuracy and lower uncertainty.All studied cases are identical except for the permeability texture.Even histograms and variograms of permeabilities for the studied reservoir are quite similar,but the performance of the reservoir shows that kriging-based method results have slightly less accuracy than some MPS methods.Meanwhile,it results in lower uncertainty in outputs for this water flooding case performance.

Numerical modeling of fluid-particle interaction in granular media

Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method （CFD-DEM） approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier-Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

MATHEMATICAL MODEL OF TWO-PHASE FLUID NONLINEAR FLOW IN LOW-PERMEABILITY POROUS MEDIA WITH APPLICATIONS

A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.

Exact Solution of Second Grade Fluid in a Rotating Frame through Porous Media Using Hodograph Transformation Method

In this paper exact solution for a homogenous incompressible, second grade fluid in a rotating frame through porous media has been provided using hodograph-Legendre transformation method. Results are summarised in the form of theorems. Two examples have been taken and streamline patterns are shown for the solutions.

Coscinium fenestratum: Callus and Suspension Cell Culture of the Endangered Medicinal Plant Using Vermicompost Extract and Coelomic Fluid as Plant Tissue Culture Media

In vitro tissue culture of hard woody, endangered, medicinal plant Coscinium fenestratum is most challenging to plant tissue culturists. In the present study, petiole and leaf explants of Coscinium fenestratum were induced to form callus when cultured on vermicompost extract media along with coelomic fluid. Suspension medium was developed using vermicompost extract and coelomic fluid in 3:1 ratio. Phytochemical analysis of the alkaloid berberine was confirmed from callus, suspension cell culture and suspension medium by Thin Layer Chromatography and High Performance Liquid Chromatography. Vermicompost and its extracts with coelomic fluid have shown maximum (100 per cent) response of callus induction. Callus mass enlarged with increasing concentration of coelomic fluid and callus growth was assessed from the biomass. Incubation of culture tubes in dark supported callus development significantly. The Rf value of 0.36 confirmed the presence of berberine by Thin Layer Chromatography. Qualitative analysis confirmed the presence of alkaloid berberine with the retention time of 2.8 minutes similar to that of standard reference sample from Sigma chemicals, USA. The suspension medium turned deep yellow because of the release of the alkaloid. Vermicompost and its extracts along with coelomic fluid have shown the economical approach for micropropagation of economically and medicinally important plants.

THE TRANSIENT ELLIPTIC FLOW OF POWER-LAW FLUID IN FRACTAL POROUS MEDIA

The steady oil production and pressure distribution formulae of vertically fractured well for power-law non-Newtonian fluid were derived on the basis of the elliptic flow model in fractal reservoirs. The corresponding transient flow in fractal reservoirs was studied by numerical differentiation method: the influence of fractal index to transient pressure of vertically fractured well was analyzed. Finally the approximate analytical solution of transient flow was given by average mass conservation law. The study shows that using elliptic flow method to analyze the flow of vertically fractured well is a simple method.

考虑多孔介质渗透率变化的宾汉流体渗透注浆扩散模型

渗透注浆是一种广泛应用于岩土多孔介质加固的方法,但已有基于浆液流变方程和平衡方程的研究,侧重分析浆液本身性质如黏度时变性、滤过效应、自身重力等对渗透注浆扩散范围的影响,尚未考虑浆液在孔隙内的流动阻力对多孔介质渗透率的影响,导致现有理论解析结果与试验测试结果相差较大。因此,基于黏性流体的运动方程、连续方程和平衡方程,引入浆液在孔隙内的流动阻力对多孔介质渗透率的影响,推导获得多孔介质渗透率随介质孔隙特性、注浆施工参数和注浆扩散半径的函数解析表达式,据此建立考虑多孔介质渗透率变化的浆液扩散控制方程。与已有理论结果对比表明,考虑浆液流动阻力影响的理论结果更接近试验结果,可以验证该理论方法的有效性。参数分析结果表明:考虑浆液在孔隙内的流动阻力时,岩土多孔介质渗透率随注浆时间的增加呈现先迅速降低后缓慢上升趋势,随注浆压力的增大呈现逐渐下降趋势;对于工程中常见的符合宾汉流体性质的水泥浆液材料,当注浆时间小于5 min时,多孔介质渗透率变化是影响浆液扩散的主控因素,超过7 min后,浆液黏度时变性成为主控因素。研究成果能为岩土多孔介质渗透注浆工程实际中注浆参数选取和质量控制起到良好的理论指导意义。

粘弹频变固液解耦流体因子地震反演方法研究及应用

地震波在含油气储层中传播时,因为波致流的影响,会发生振幅衰减与弹性参数频散现象,使得地震流体预测描述成为可能。但是,这些现象产生理论不明确和数学表达精度不高,使得含油气识别困难。为此,本研究首先根据地震振幅衰减与弹性参数频散现象,利用Chapman孔裂隙衰减理论考虑波致流影响,引入挤喷流效应,构建粘弹介质固液解耦流体因子。然后,推导构建关于流体因子的反射系数特征方程,并与Zoeppritz方程和Aki近似方程对比分析,验证新方程的精度。最后,基于叠前地震数据的分角度、分频带变化信息,提出粘弹频变固液解耦流体因子反演方法,并利用该流体因子在胜利油田胜北地区进行储层含油气性预测。研究结果表明,新构建的粘弹频变固液解耦流体因子对储层油气具有较强的敏感性,提出的地震反射特征方程和反演方法准确可靠,实际储层流体预测结果与测井解释结果吻合度较高,为复杂储层流体识别提供新的思路和方法。

射孔完井近井地带等效剪切控制体物理模型

为研究地层对含聚液体的剪切降解作用,应用ANSYS软件求解了射孔完井近井地带的压力场和速度场分布,并用MATLAB对压力场和速度场数据进行后处理,划分出10段射孔孔眼控制体,计算各控制体内等压面上的平均流速和各等压面沿流线轨迹到射孔孔壁距离,得出各段孔眼控制体内渗流速度随流线轨迹上渗流距离的变化关系。通过变直径、等直径岩心组合物模实验,测试了各段孔眼控制体内沿流线轨迹不同渗流距离处聚合物溶液的流变性,用流量加权法计算出了整个射孔孔眼控制区域聚合物溶液黏度保留率,按照等效剪切原理确定了射孔孔眼等效剪切控制体的位置及参数,研制了射孔孔眼等效剪切控制体物理模型,该模型在定排量条件下即可准确模拟射孔完井近井地带的变速渗流场对含聚液体的剪切降解作用。

基于孔裂隙储层弹性参数频散衰减机理的有限元数值模拟

饱含流体孔隙裂隙储层介质的弹性参数频变特性常应用于地震勘探和储层评价中。针对传统有效介质理论不准确及计算精度不高的问题,该研究推导了适用于孔隙弹性储层的Biot固结方程;为了研究模型介质渗透率、流体类型和黏度等参数对频变弹性参数的影响机理,开发了相应的有限元数值模拟流程,用于计算由波致流引起的地震波速度频散及衰减。模拟结果表明:1)与现有岩石物理理论模型对比,该数值模拟方法在宽频段模拟一维/二维孔隙裂隙介质的适用性更好;2)渗透率、流体黏度以及流体饱和模式作为岩心重要的物性表征,显著影响地震波弹性响应,渗透率增大,特征频率右移;黏度增大,特征频率左移;水包气模型相对于油包气模型衰减峰值较大,特征频率略高。综上所述,该研究为深入理解地震波在孔裂隙储层中的传播特性提供了准确的预测方法,有助于揭示储层的非均质性和孔隙结构,提高勘探精度和储层评价的可靠性,对储层预测与流体识别的岩石物理解释具有重要的指导意义。

剪切增稠流体在Kevlar复合材料中的应用研究进展

高强度、高韧性、低密度的凯夫拉(Kevlar)材料在防护领域广泛应用,但直接使用Kevlar织物制备防护材料会使其过于笨重。因此,人们研究了轻质灵活的Kevlar复合材料,使用剪切增稠流体来增强其性能。本综述总结了使用不同分散相和分散介质的剪切增稠流体浸渍Kevlar复合材料的研究进展,旨在为防护材料的后续研究提供思路。