Effects of seepage pressure on the mechanical behaviors and microstructure of sandstone

Surrounding rocks of underground engineering are subjected to long-term seepage pressure,which can deteriorate the mechanical properties and cause serious disasters.In order to understand the impact of seepage pressure on the mechanical property of sandstone,uniaxial compression tests,P-wave velocity measurements,and nuclear magnetic resonance(NMR)tests were conducted on saturated sandstone samples with varied seepage pressures(i.e.0 MPa,3 MPa,4 MPa,5 MPa,6 MPa,7 MPa).The results demonstrate that the mechanical parameters(uniaxial compressive strength,peak strain,elastic modulus,and brittleness index),total energy,elastic strain energy,as well as elastic strain energy ratio,decrease with increasing seepage pressure,while the dissipation energy and dissipation energy ratio increase.Moreover,as seepage pressure increases,the micro-pores gradually transform into meso-pores and macro-pores.This increases the cumulative porosity of sandstone and decreases P-wave velocity.The numerical results indicate that as seepage pressure rises,the number of tensile cracks increases progressively,the angle range of microcracks is basically from 50-120to 80-100,and as a result,the failure mode transforms to the tensile-shear mixed failure mode.Finally,the effects of seepage pressure on mechanical properties were discussed.The results show that decrease in the effective stress and cohesion under the action of seepage pressure could lead to deterioration of strength behaviors of sandstone.

Impact of Osmotic Pressure on Seepage in Shale Oil Reservoirs

Following large-scale volume fracturing in shale oil reservoirs,well shut-in measures are generally employed.Laboratory tests and field trials have underscored the efficacy of fracturing fluid imbibition during the shut-in phase in augmenting shale oil productivity.Unlike conventional reservoirs,shale oil reservoirs exhibit characteristics such as low porosity,low permeability,and rich content of organic matter and clay minerals.Notably,the osmotic pressure effects occurring between high-salinity formation water and low-salinity fracturing fluids are significant.The current understanding of the mobilization patterns of crude oil in micro-pores during the imbibition process remains nebulous,and the mechanisms underpinning osmotic pressure effects are not fully understood.This study introduces a theoretical approach,by which a salt ion migration control equation is derived and a mathematical model for spontaneous imbibition in shale is introduced,which is able to account for both capillary and osmotic pressures.Results indicate that during the spontaneous imbibition of low-salinity fluids,osmotic effects facilitate the migration of external fluids into shale pores,thereby complementing capillary forces in displacing shale oil.When considering both capillary and osmotic pressures,the calculated imbibition depth increases by 12%compared to the case where only capillary forces are present.The salinity difference between the reservoir and the fracturing fluids significantly influences the imbibition depth.Calculations for the shutin phase reveal that the pressure between the matrix and fractures reaches a dynamic equilibrium after 28 days of shut-in.During the production phase,the maximum seepage distance in the target block is approximately 6.02 m.

Effect mechanism of seepage force on the hydraulic fracture propagation

The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.

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.

Experimental study of the effects of a multistage pore-throat structure on the seepage characteristics of sandstones in the Beibuwan Basin:Insights into the flooding mode

To investigate the relationship between grain sizes, seepage capacity, and oil-displacement efficiency in the Liushagang Formation of the Beibuwan Basin, this study identifies the multistage pore-throat structure as a crucial factor through a comparison of oil displacement in microscopic pore-throat experiments. The two-phase flow evaluation method based on the Li-Horne model is utilized to effectively characterize and quantify the seepage characteristics of different reservoirs, closely relating them to the distribution of microscopic pores and throats. It is observed that conglomerate sandstones at different stages exhibit significant heterogeneity and noticeable differences in seepage capacity, highlighting the crucial role played by certain large pore throats in determining seepage capacity and oil displacement efficiency. Furthermore, it was found that the displacement effects of conglomeratic sandstones with strong heterogeneity were inferior to those of conventional homogeneous sandstone, as evidenced by multiple displacement experiments conducted on core samples with varying granularities and flooding systems. Subsequently, core-based experiments on associated gas flooding after water flooding were conducted to address the challenge of achieving satisfactory results in a single displacement mode for reservoirs with significant heterogeneity. The results indicate that the oil recovery rates for associated gas flooding after water flooding increased by 7.3%-16.4% compared with water flooding alone at a gas-oil ratio of approximately 7000 m^(3)/m^(3). Therefore, considering the advantages of gas flooding in terms of seepage capacity, oil exchange ratio, and the potential for two-phase production, gas flooding is recommended as an energy supplement mode for homogeneous reservoirs in the presence of sufficient gas source and appropriate tectonic angle. On the other hand, associated gas flooding after water flooding is suggested to achieve a more favorable development effect compared to a single mode of energy supplementation for strongly heterogeneous sandstone reservoirs.

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.

Experimental investigations on characteristics of water-air two-phase flows during water-fillings in undulation pipelines

Visual and pressurized pipeline systems with single- and multi-undulation layouts were used to study experimentally and analyze theoretically the transient characteristics of water-air two-phase flow during water fillings in undulation pipelines based on the combination action analyses of both the communicating pipe and the gravity of the water-air two-phase flows in the descending pipe. For the single undulation pipeline, the complex two-phase flow-pattern evolutions including full pipe flow and stratified flow for low, medium, high water-filling velocity cases, respectively, lead to a great difference in transient pressure, flow pattern and the water-filling duration. Especially for low and medium water-filling velocity cases, the hydraulic theories related to hydraulic drop and hydraulic jump were employed to investigate the entrapped air pocket evolutions in the descending pipe, and the mechanism of negative pressure at the top of the undulation pipes was analyzed. For the same multi-undulation pipeline, due to the different elevations of the three undulation points along flow direction, namely three different types of pipeline layout, high-medium-low case (high elevation undulation point, medium one, and low one), low-medium-high and high-low-medium ones, their water-filling durations are significantly different, i.e., approximately 80.02 s, 227.34 s and 617.78 s. Meanwhile, there are significant differences in flow patterns in water filling, namely larger entrapped air pockets in three descending pipes for the high-medium-low case, entrapped air pockets in the first two descending pipes and open channel stratified flow in the last one for low-medium-high case, some bubbles in three descending pipes for the high-low-medium case.

Simulation of Two-Phase Flowback Phenomena in Shale Gas Wells

The gas-water two-phaseflow occurring as a result of fracturingfluidflowback phenomena is known to impact significantly the productivity of shale gas well.In this work,this two-phaseflow has been simulated in the framework of a hybrid approach partially relying on the embedded discrete fracture model(EDFM).This model assumes the region outside the stimulated reservoir volume(SRV)as a single-medium while the SRV region itself is described using a double-medium strategy which can account for thefluid exchange between the matrix and the micro-fractures.The shale gas adsorption,desorption,diffusion,gas slippage effect,fracture stress sensitivity,and capillary imbibition have been considered.The shale gas production,pore pressure distribution and water saturation distribution in the reservoir have been simulated.The influences of hydraulic fracture geometry and nonorthogonal hydraulic fractures on gas production have been determined and discussed accordingly.The simulation results show that the daily gas production has an upward and downward trend due to the presence of a large amount of fracturingfluid in the reservoir around the hydraulic fracture.The smaller the angle between the hydraulic fracture and the wellbore,the faster the daily production of shale gas wells decreases,and the lower the cumulative production.Nonplanar fractures can increase the control volume of hydraulic fractures and improve the production of shale gas wells.

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.

A Note on an Order Level Inventory Model with Varying Two-Phased Demand and Time-Proportional Deterioration

The main purpose of this paper is to generalize the effect of two-phased demand and variable deterioration within the EOQ (Economic Order Quantity) framework. The rate of deterioration is a linear function of time. The two-phased demand function states the constant function for a certain period and the quadratic function of time for the rest part of the cycle time. No shortages as well as partial backlogging are allowed to occur. The mathematical expressions are derived for determining the optimal cycle time, order quantity and total cost function. An easy-to-use working procedure is provided to calculate the above quantities. A couple of numerical examples are cited to explain the theoretical results and sensitivity analysis of some selected examples is carried out.

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.

Impact of seepage flow on sediment resuspension by internal solitary waves:parameterization and mechanism

Sediment incipient motion is the first step in sediment resuspension.Previous studies ignored the effect of seepage flow on the mobility of sediment particles and simplified the seabed surface as a rigid boundary.A flume experiment was designed to innovatively divide the seabed into two parts to control the dynamic response of the seabed and control the seepage conditions.In the experiment,the seabed sediments and the amplitude of internal solitary waves(ISWs)were changed to compare and analyze the impact of seepage flow on the sediment resuspension by shoaling ISWs.Moreover,parametric research and verification were carried out.Results indicate that seepage flow can greatly influence fine sand,promote sediment resuspension,and increase the amount of suspension by two times on average.However,seepage flow had a little effect on the suspension of clayey silt and sandy silt.Besides,seepage force was added to the traditional gravity,drag force,and uplift force,and the parameterization of threshold starting shear stress of coarse-grained sediments was developed.The results of this parameterization were verified,and seepage force was critical to parameterization.The threshold starting shear stress was reduced by 54.6%after increasing the seepage force.The physical mechanism of this process corresponded to the vertical reciprocating transient seepage in and out the seabed interface caused by the wave-induced transient excess pore water pressure.This quantitative study on seepage flow for shear stress of coarse-grained sediments induced by ISWs is critical to geohazard assessment.

Sensitivity Analysis of Multi-phase Seepage Parameters Affecting the Clayey Silt Hydrate Reservoir Productivity in the Shenhu Area,South China Sea

Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately predict the productivity of hydrate-bearing sediments(HBS).The multi-phase seepage parameters of HBS include permeability,porosity,which is closely related to permeability,and hydrate saturation,which has a direct impact on hydrate content.Existing research has shown that these multi-phase seepage parameters have a great impact on HBS productivity.Permeability directly affects the transmission of pressure-drop and discharge of methane gas,porosity and initial hydrate saturation affect the amount of hydrate decomposition and transmission process of pressure-drop,and also indirectly affect temperature variation of the reservoir.Considering the spatial heterogeneity of multi-phase seepage parameters,a depressurization production model with layered heterogeneity is established based on the clayey silt hydrate reservoir at W11 station in the Shenhu Sea area of the South China Sea.Tough+Hydrate software was used to calculate the production model;the process of gas production and seepage parameter evolution under different multi-phase seepage conditions were obtained.A sensitivity analysis of the parameters affecting the reservoir productivity was conducted so that:(a)a HBS model with layered heterogeneity can better describe the transmission process of pressure and thermal compensation mechanism of hydrate reservoir;(b)considering the multi-phase seepage parameter heterogeneity,the influence degrees of the parameters on HBS productivity were permeability,porosity and initial hydrate saturation,in order from large to small,and the influence of permeability was significantly greater than that of other parameters;(c)the production potential of the clayey silt reservoir should not only be determined by hydrate content or seepage capacity,but also by the comprehensive effect of the two;and(d)time scales need to be considered when studying the effects of changes in multi-phase seepage parameters on HBS productivity.

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.

A two-phase type-curve method with multiscale fluid transport mechanisms in hydraulically fractured shale reservoirs

The quantitative understanding of hydraulic fracture(HF)properties guides accurate production forecasts and reserve estimation.Type curve is a powerful technique to characterize HF and reservoir properties from flowback and long-term production data.However,two-phase flow of water and hydrocarbon after an HF stimulation together with the complex transport mechanisms in shale nanopores exacerbate the nonlinearity of the transport equation,causing errors in type-curve analysis.Accordingly,we propose a new two-phase type-curve method to estimate HF properties,such as HF volume and permeability of fracture,through the analysis of flowback data of multi-fractured shale wells.The proposed type curve is based on a semianalytical solution that couples the two-phase flow from the matrix with the flow in HF by incorporating matrix influx,slippage effect,stress dependence,and the spatial variation of fluid properties in inorganic and organic pores.For the first time,multiple fluid transport mechanisms are considered into two-phase type-curve analysis for shale reservoirs.We analyze the flowback data from a multi-fractured horizontal well in a shale gas reservoir to verify the field application of the proposed method.The results show that the fracture properties calculated by the type-curve method are in good agreement with the long-time production data.

Groundwater flow through fractured rocks and seepage control in geotechnical engineering: Theories and practices

Groundwater flow through fractured rocks has been recognized as an important issue in many geotechnical engineering practices.Several key aspects of fundamental mechanisms,numerical modeling and engineering applications of flow in fractured rocks are discussed.First,the microscopic mechanisms of fluid flow in fractured rocks,especially under the complex conditions of non-Darcian flow,multiphase flow,rock dissolution,and particle transport,have been revealed through a com-bined effort of visualized experiments and theoretical analysis.Then,laboratory and field methods of characterizing hydraulic properties(e.g.intrinsic permeability,inertial permeability,and unsaturated flow parameters)of fractured rocks in different flow regimes have been proposed.Subsequently,high-performance numerical simulation approaches for large-scale modeling of groundwater flow in frac-tured rocks and aquifers have been developed.Numerical procedures for optimization design of seepage control systems in various settings have also been proposed.Mechanisms of coupled hydro-mechanical processes and control of flow-induced deformation have been discussed.Finally,three case studies are presented to illustrate the applications of the improved theoretical understanding,characterization methods,modeling approaches,and seepage and deformation control strategies to geotechnical engi-neering projects.

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.