Pore-scale Simulation of Water Alternate Gas Injection

We use a three-dimensional mixed-wet random network model representing Berea sandstone to compute displacement paths and relative permeabilities for water alternating gas (WAG) flooding. First we reproduce cycles of water and gas injection observed in previously published experimental studies. We predict the measured oil, water and gas relative permeabilities accurately. We discuss the hysteresis trends in the water and gas relative permeabilities and compare the behavior of water-wet and oil-wet media. We interpret the results in terms of pore-scale displacements. In water-wet media the water relative permeability is lower during water injection in the presence of gas due to an increase in oil/water capillary pressure that causes a decrease in wetting layer conductance. The gas relative permeability is higher for displacement cycles after first gas injection at high gas saturation due to cooperative pore filling, but lower at low saturation due to trapping. In oil-wet media, the water relative permeability remains low until water-filled elements span the system at which point the relative permeability increases rapidly. The gas relative permeability is lower in the presence of water than oil because it is no longer the most non-wetting phase.     

Pore-Scale Flow in Surfactant Flooding

Recovery of oil from the blocks of an initially oil-wet, naturally fractured, reservoir as a result of counter-current flow following introduction of aqueous wettability-altering surfactant into the fracture system is considered, as an example of a practical process in which phenomena acting at the single pore-scale are vital to the economic displacement of oil at the macroscopic scale. A Darcy model for the process is set up, and solutions computed illustrating the recovery rate controlling role of the bulk diffusion of surfactant. A central ingredient of this model is the capillary pressure relation, linking the local values of the pressure difference between the oleic and aqueous phases, the aqueous saturation and the surfactant concentration. Using ideas from single capillary models of oil displacement from oil-wet tubes by wettability-altering surfactant, we speculate that the use of a capillary pressure function, with dependences as assumed, may not adequately represent the Darcy scale consequences of processes acting at the single pore-scale. Multi-scale simulation, resolving both sub-pore and multi-pore flow processes may be necessary to resolve this point. Some general comments are made concerning the issues faced when modelling complex displacement processes in porous media starting from the pore-scale and working upwards.     

Accurate Modelling of Pore-Scale Films and Layers for Three-Phase Flow Processes in Clastic and Carbonate Rocks with Arbitrary Wettability

Three-phase flow is a key process occurring in subsurface reservoirs, for example, during $\text{ CO }_2$ sequestration and enhanced oil recovery techniques such as water alternating gas (WAG) injection. Predicting three-phase flow processes, for example, the increase in oil recovery during WAG, requires a sound understanding of the fundamental flow physics in water- to oil-wet rocks to derive physically robust flow functions, i.e. relative permeability and capillary pressure. In this study, we use pore-network modelling, a reliable and physically based simulation tool, to predict the flow functions. We have developed a new pore-scale network model for rocks with variable wettability, from water- to oil-wet. It comprises a constrained set of parameters that mimic the wetting state of a reservoir. Unlike other models, it combines three main features: (1) A novel thermodynamic criterion for formation and collapse of oil layers. The new model hence captures wetting film and layer flow of oil adequately, which affects the oil relative permeability at low oil saturation and leads to accurate prediction of residual oil. (2) Multiple displacement chains, where injection of one phase at the inlet triggers a chain of interface displacements throughout the network. This allows for the accurate modelling of the mobilisation of many disconnected phase clusters that arise, in particular, during higher order WAG floods. (3) The model takes realistic 3D pore-networks extracted from pore-space reconstruction methods and CT images as input, preserving both topology and pore shape of the sample. For water-wet systems, we have validated our model with available experimental data from core floods. For oil-wet systems, we validated our network model by comparing 2D network simulations with published data from WAG floods in oil-wet micromodels. This demonstrates the importance of film and layer flow for the continuity of the various phases during subsequent WAG cycles and for the residual oil saturations. A sensitivity analysis has been carried out with the full 3D model to predict three-phase relative permeabilities and residual oil saturations for WAG cycles under various wetting conditions with different flood end-points.     

Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media

In this article, the numerical simulations for one-dimensional three-phase flows in fractured porous media are implemented. The simulation results show that oil displacement in matrix is dominated by oil–water capillary pressure only under certain conditions. When conditions are changed to decrease the amount of water entering into the fractured media from the boundary of the flow field, water in fracture may be vaporized to superheated steam. In these cases, the appearance of superheated steam in fracture rather than in matrix will decrease the fracture pressure and generate the pressure difference between matrix and fracture, which results in oil flowing from matrix to fracture. Assuming that oil is wetting to steam, the matrix steam–oil capillary pressure will decrease the matrix oil-phase pressure as the matrix steam saturation increases. After the steam–oil capillary pressure finally exceeds the pressure difference due to the appearance of superheated steam in fracture, the oil displacement in matrix will stop. It is also shown that variations of the water relative permeability curve in matrix do not result in different mechanisms for oil displacement in matrix. The simulation results suggest that the amount of liquid water supply from the boundary of flow field fundamentally influence the mechanisms for oil displacement in matrix.     

Efficiency of Capillary Imbibition Dominated Displacement of Nonwetting Phase by Wetting Phase in Fractured Porous Media

The critical and optimum injection rates as well as the critical fracture capillary number for an efficient displacement process are determined based on the experimental and numerical modeling of the displacement of nonwetting phase (oil) by wetting phase (water) in fractured porous media. The efficiency of the process is defined in terms of the nonwetting phase displaced from the system per amount of wetting phase injected and per time. Also, the effects of injection rate on capillary imbibition transfer dominated two-phase flow in fractured porous media are clarified by visualizing the experiments. The results reveal that as the injection rate is increased, fracture pattern begins to become an effective parameter on the matrix saturation distribution. As the rate is lowered, however, the system begins to behave like a homogeneous system showing a frontal displacement regardless the fracture configuration.     

Pore-Scale Imaging and Analysis of Wettability Order,Trapping and Displacement in Three-Phase Flow in Porous Media with Various Wettabilities

Three-phase flow in porous media is encountered in many applications including subsurface carbon dioxide storage, enhanced oil recovery, groundwater remediation and the design of microfluidic devices. However, the pore-scale physics that controls three-phase flow under capillary dominated conditions is still not fully understood. Recent advances in three-dimensional pore-scale imaging have provided new insights into three-phase flow. Based on these findings, this paper describes the key pore-scale processes that control flow and trapping in a three-phase system, namely wettability order, spreading and wetting layers, and double/multiple displacement events. We show that in a porous medium containing water, oil and gas, the behaviour is controlled by wettability, which can either be water-wet, weakly oil-wet or strongly oil-wet, and by gas–oil miscibility. We provide evidence that, for the same wettability state, the three-phase pore-scale events are different under near-miscible conditions—where the gas–oil interfacial tension is ≤?1 mN/m—compared to immiscible conditions. In a water-wet system, at immiscible conditions, water is the most-wetting phase residing in the corners of the pore space, gas is the most non-wetting phase occupying the centres, while oil is the intermediate-wet phase spreading in layers sandwiched between water and gas. This fluid configuration allows for double capillary trapping, which can result in more gas trapping than for two-phase flow. At near-miscible conditions, oil and gas appear to become neutrally wetting to each other, preventing oil from spreading in layers; instead, gas and oil compete to occupy the centre of the larger pores, while water remains connected in wetting layers in the corners. This allows for the rapid production of oil since it is no longer confined to movement in thin layers. In a weakly oil-wet system, at immiscible conditions, the wettability order is oil–water–gas, from most to least wetting, promoting capillary trapping of gas in the pore centres by oil and water during water-alternating-gas injection. This wettability order is altered under near-miscible conditions as gas becomes the intermediate-wet phase, spreading in layers between water in the centres and oil in the corners. This fluid configuration allows for a high oil recovery factor while restricting gas flow in the reservoir. Moreover, we show evidence of the predicted, but hitherto not reported, wettability order in strongly oil-wet systems at immiscible conditions, oil–gas–water, from most to least wetting. At these conditions, gas progresses through the pore space in disconnected clusters by double and multiple displacements; therefore, the injection of large amounts of water to disconnect the gas phase is unnecessary. We place the analysis in a practical context by discussing implications for carbon dioxide storage combined with enhanced oil recovery before suggesting topics for future work.

     

A New Mathematical Model for Force Gravity Drainage in Fractured Porous Media

In force gas/oil gravity drainage process in fractured porous media, gas is flowing in both matrix and fractures leading to produce a finite gas pressure gradient. Consequently, viscous force plays an important role for displacing matrix oil toward fractures in addition to gravity force that is required to be modeled appropriately. A new analytical model for estimation of steady state oil saturation distribution with assumption of fixed gas pressure gradient throughout the matrix is presented. Moreover, based on some results of this analytical model a different numerical formulation is developed to predict the performance of oil production process. Comparison of the results obtained from this numerical model with the results of a conventional simulator demonstrates that the newly developed model can be applied with satisfactory accuracy. Numerical simulations show that the viscous displacement in fractured porous media can reduce the capillary threshold height, and thus it suggests the force gravity drainage as a favorable production mechanism when the matrix length is close to the threshold height.     

Numerical Analysis of Imbibition Front Evolution in Fractured Sandstone under Capillary-Dominated Conditions

Fractures serve as primary conduits having a great impact on the migration of injected fluid into fractured permeable media. Appropriate transport properties such as relative permeability and capillary pressure are essential for successful simulation and prediction of multi-phase flow in such systems. However, the lack of a thorough understanding of the dynamics governing immiscible displacement in fractured media, limits our ability to properly represent their macroscopic transport properties. Previous experimental observations of imbibition front evolution in fractured rocks are examined in the present study using an automated history-matching approach to obtain representative relative permeability and capillary pressure curves. Predicted imbibition front evolution under different flow conditions resulted in an excellent agreement with experimental observations. Sensitivity analyses, in combination with direct experimental observation, allowed exploring the competing effects of relative permeability and capillary pressure on the development of saturation distribution and imbibing front evolution in fractured porous media. Results show that residual saturations are most sensitive to matrix relative permeability to oil, while the ratio of oil and water relative permeability, rock heterogeneity, boundary condition, and matrix–fracture capillary pressure contrast, affect displacement shape, speed, and geometry of the imbibing front.     

Modeling Water Leak-off Behavior in Hydraulically Fractured Gas Shale under Multi-mechanism Dominated Conditions

Fracturing-fluid leak-off in fractured gas shale is a complex process involving multiple pore/fluid transports and interactions. However, water leak-off behavior has not been modeled comprehensively by considering the multi-pores and multi-mechanisms in shale with existing simulators. In this paper, we present the development of a comprehensive multi-mechanistic, multi-porosity, and multi-permeability water/gas flow model that uses experimentally determined formation properties to simulate the fracturing-fluid leak-off of hydraulically fractured shale gas wells. The multi-mechanistic model takes into account water transport driven by hydraulic convection, capillary and osmosis, gas transport caused by hydraulic convection, and salt ion transport caused by advection and diffusion. The multi-porosity includes hydraulic fracture millipores, organic nanopores, clay nanopores, and other inorganic micropores. The multi-permeability model accounts for all the important processes in shale system, including gas adsorption on the organics’ surface, multi-mechanistic clay/other inorganic mineral mass transfer, inorganic mineral/hydraulic fracture mass transfer, and injection from a hydraulically fractured wellbore. The dynamic water saturation and pressure profiles within clay and other inorganic matrices are compared, revealing the leak-off behavior of water in rock media with different physicochemical properties. In sensitivity analyses, cases with different clay membrane efficiency, volume proportion of source rock, connate water salinity, and saturation are considered. The impacts of shale properties on water fluxes through wellbore, hydraulic fracture and matrix, and the total injection and leak-off volumes of the well during the treatment of hydraulic fracturing are investigated. Results show that physicochemical properties in both organic and inorganic matrices affect the water leak-off behavior.     

Study of Microscopic and Macroscopic Displacement Behaviors of Polymer Solution in Water-Wet and Oil-Wet Media

Performance of a polymer flood process requires the knowledge of rheological behavior of the polymer solution and reservoir properties such as rock wettability. To provide a better understanding of effects of polymer chemistry and wettability on the performance of a polymer flood process, a comprehensive experimental study was conducted using a two-dimensional glass micromodel. A series of water and polymer flood processes were carried out at different polymer molecular weights, degrees of polymer hydrolysis, and polymer concentrations in both water-wet and oil-wet systems. Image processing technique was applied to analyze and compare microscopic and macroscopic displacement behaviors of polymer solution in each experiment. From micro-scale observations, the configuration of connate water film, polymer solution trapping, flow of continuous and discontinuous strings of polymer solution, piston-type displacement of oil, snap-off of polymer solution, distorted flow of polymer solution, emulsion formation, and microscopic pore-to-pore sweep of oil phase were observed and analyzed in the strongly oil-wet and water-wet media. Rheological experiments showed that a higher polymer molecular weight, degree of hydrolysis, and concentration result in a higher apparent viscosity for polymer solution and lower oil–polymer viscosity ratio. It is also shown that these parameters have different impacts on the oil recovery in different wettabilities. Moreover, a water-wet medium generally had higher recovery in contrast with an oil-wet medium. This experimental study illustrates the successful application of glass micromodel techniques for studying enhanced oil recovery (EOR) processes in five-spot pattern and provides a useful reference for understanding the displacement behaviors in a typical polymer flood process.     

Miscible and Immiscible Foam Injection for Mobility Control and EOR in Fractured Oil-Wet Carbonate Rocks

Foam injection is a proven enhanced oil recovery (EOR) technique for heterogeneous reservoirs, but is less studied for EOR in fractured systems. We experimentally investigated tertiary \(\text {CO}_{2}\) injections, and \(\text {N}_{2}\) - and \(\text {CO}_{2}\) -foam injections for enhanced oil recovery in fractured, oil-wet limestone core plugs. Miscible \(\text {CO}_{2}\) and \(\text {CO}_{2}\) -foam was compared with immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foam as tertiary recovery techniques, subsequent to waterfloods, in fractured rocks with different wettability preferences. At water-wet conditions waterfloods produced approximately 40 % OOIP, by spontaneous imbibition. Waterflood oil recovery at oil-wet conditions was below 20 % OOIP, due to suppressed imbibition where water predominantly flowed through the fractures, unable to mobilize the oil trapped in the matrix. Tertiary, supercritical \(\text {CO}_{2}\) -mobilized oil trapped in the matrix, particularly at weakly oil-wet conditions, by diffusion. Recovery by diffusion was high due to small core samples, high initial oil saturation and a continuous oil phase at oil-wet conditions. Both immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foams and miscible, supercritical \(\text {CO}_{2}\) -foam demonstrated high ultimate oil recoveries, but immiscible foam was less efficient (30 pore volumes injected) compared to miscible foam (2 pore volumes injected) to reach ultimate recovery. This is explained by the capillary threshold pressure preventing the injected \(\text {N}_{2}\) gas from entering the matrix, verified by computed X-ray tomography, and the mobilized oil was displaced by the aqueous surfactant in the foam. At miscible conditions, there exists no capillary entry pressure between the oil-saturated matrix and the injected \(\text {CO}_{2}\) , allowing foam to invade the matrix for efficient oil recovery.     

A Parametric Model for Predicting Relative Permeability-Saturation-Capillary Pressure Relationships of Oil–Water Systems in Porous Media with Mixed Wettability

A parametric two-phase, oil–water relative permeability/capillary pressure model for petroleum engineering and environmental applications is developed for porous media in which the smaller pores are strongly water-wet and the larger pores tend to be intermediate- or oil-wet. A saturation index, which can vary from 0 to 1, is used to distinguish those pores that are strongly water-wet from those that have intermediate- or oil-wet characteristics. The capillary pressure submodel is capable of describing main-drainage and hysteretic saturation-path saturations for positive and negative oil–water capillary pressures. At high oil–water capillary pressures, an asymptote is approached as the water saturation approaches the residual water saturation. At low oil–water capillary pressures (i.e. negative), another asymptote is approached as the oil saturation approaches the residual oil saturation. Hysteresis in capillary pressure relations, including water entrapment, is modeled. Relative permeabilities are predicted using parameters that describe main-drainage capillary pressure relations and accounting for how water and oil are distributed throughout the pore spaces of a porous medium with mixed wettability. The capillary pressure submodel is tested against published experimental data, and an example of how to use the relative permeability/capillary pressure model for a hypothetical saturation-path scenario involving several imbibition and drainage paths is given. Features of the model are also explained. Results suggest that the proposed model is capable of predicting relative permeability/capillary pressure characteristics of porous media mixed wettability.     

Displacement of Cr(III)–Partially Hydrolyzed Polyacrylamide Gelling Solution in a Fracture in Porous Media

During waterflooding of a fractured formation, water may channel through the fracture or interconnected network of fractures, leaving a large portion of oil bearing rock unswept. One remedial practice is injection of a gelling solution into the fracture. Such placement of a gelling mixture is associated with leak-off from the fracture face into the adjoining matrix. Design of a gel treatment needs understanding of the flow of gelling mixture in and around the fracture. This flow is addressed here for Cr(III)–partially hydrolyzed polyacrylamide formulation through experiments and conceptual model. A fractured slab was used to develop a lab-model, where the flow along the fracture and simultaneous leak-off into the matrix can be controlled. Also, the fracture and matrix properties had to be evaluated individually for a meaningful analysis of the displacement of gelling solution. During this displacement, the gelling fluid leaked off from the fracture into the matrix as a front, resulting in a decreasing velocity (and pressure gradient) along the fracture. With pressure in the fracture held constant with time, the leak-off rate decreased as the viscous front progressed into the matrix. The drop in leak-off rate was rapid during the initial phase of displacement. A simple model, based on the injection of a viscous solution into the dual continua, could explain the displacement of Cr(III)–polyacrylamide gelling mixture through the fractured slab. This study rules out any major complication from the immature gelling fluid, e.g., build-up of cake layer on the fracture face. The model, due to its simplicity may become useful for quick sizing of gel treatment, and any regression-based evaluation of fluid properties in a fracture for other applications.     

Diffusion Mass Transfer in Miscible Oil Recovery: Visual Experiments and Simulation

Co-and counter-current type transfers due to diffusion and -free- convection caused by the buoyant forces between fracture and matrix were studied experimentally using 2-D glass-bead models. Mineral oil and kerosene were used as the displaced phase. The model saturated with oil was exposed to solvent phase (pentane) under static conditions (no flow in fracture) to mimic matrix-fracture interaction during gas or liquid solvent injection in naturally fractured reservoirs. Displacement fronts and patterns were analyzed and quantified using fractal techniques to obtain correlations between the fractal properties and displacement type. Displacements resulted in a mixture of bulk diffusion and -free- convection mainly depending on the interaction type (co- or counter-current), oil type, and displacement direction (horizontal and vertical). Conditions yielding different types of displacement patterns were identified. Finally, a stochastic model that was inspired from invasion percolation and diffusion limited aggregation algorithms was developed for the horizontal displacement cases. The experimental observations were matched to the displacement patterns obtained through the stochastic modeling.     

Unstable Immiscible Displacements in Oil-Wet Rocks

Displacement of a viscous fluid by a lower viscosity immiscible fluid (such as waterflood of a viscous oil) in a porous medium is unstable. The displacement front generates viscous fingers which lead to low oil recovery efficiency. These fingers are much smaller in width than typical reservoir simulation grid blocks, and capturing their effect in reservoir simulation is important. A dimensionless scaling group (viscous finger number) had been suggested in the past, which has a power-law relationship with the breakthrough recovery and cumulative recovery in unstable core floods. The relative permeability used in large grid block simulations had been modified to so-called pseudo-relative permeability on the basis of the dimensionless group, thus incorporating the effect of fingers in waterflood predictions. However, the previous proposed models were constructed from experiments in only water-wet rocks. This paper extends the recent viscous fingering models to oil-wet systems. Sandstone cores were treated to alter the wettability to oil-wet. Adverse viscosity water floods were performed in oil-wet cores. Viscosity ratio, velocity and diameter were varied. It is shown that the previously developed viscous finger number does not work for the oil-wet experiments. The correlating dimensionless number is modified for oil-wet systems; it is also different from the dimensionless group identified by Peters and Flock (Soc Petroleum Eng, 1981. doi: 10.2118/8371-PA) for oil-wet cores. A pseudo-relative permeability model has been developed for oil-wet cores. Corefloods have been matched by the new pseudo-relative permeability model to determine the model parameters. This pseudo-relative permeability model can be used in reservoir simulations of water and polymer floods in viscous oil-wet reservoirs.     

Dynamic Relative Permeability and Simulation of WAG Injection Processes

In this work, we investigate the impact of mobility changes due to flow reversals from co-current to counter-current flow on the displacement performance of water alternating gas (WAG) injection processes. In WAG processes, the injected gas will migrate toward the top of the formation while the injected water will migrate toward the bottom of the formation. The segregation of gas, oil and water phases will result in counter-current flow occurring in the vertical direction in some portions of the reservoir during the displacement process. Previous experimental and theoretical studies of counter-current flow have shown that the relative mobility of each of the phases in a porous medium is considerably less when counter-current flow prevails as compared to co-current flow settings. A reduction of the relative permeability in the vertical direction results in a dynamic anisotropy in phase mobilities. This effect has, to the best of our knowledge, not previously been considered in the modeling and simulation of WAG processes. A new flow model that accounts for flow reversals in the vertical direction has been implemented and tested in a three-phase compositional reservoir simulator. In order to investigate the impact of flow reversals, results from the new flow model are compared to cases where counter-current flow effects on the phase mobilities are ignored. A range of displacement settings, covering relevant slug sizes, have been investigated to gauge the impact of mobility reductions due to flow reversals. Significant differences, in terms of saturation distribution, producing GOR and oil recovery, are observed between the conventional flow model (ignoring mobility reductions due to counter-current flow) and the proposed new model that accounts for reductions in phase mobility during counter-current flow. Accordingly, we recommend that an explicit representation of flow transitions between co-current and counter-current flow (and the related impact on phase mobilities) should be considered to ensure accurate and optimal design of WAG injection processes.     

Wetting Phase Bridges Establish Capillary Continuity Across Open Fractures and Increase Oil Recovery in Mixed-Wet Fractured Chalk

The effect of fractures on oil recovery and in situ saturation development in fractured chalk has been determined at near neutral wettability conditions. Fluid saturation development was monitored both in the matrix and in the fractures and the mechanisms of fracture crossing were determined using high spatial resolution MRI. Capillary continuity across open oil-filled fractures was verified by imaging the water bridges established within the fracture. Despite an alternate escape fracture for the water, separate water bridges were shown to be stable for the entire duration of the experiments. The established capillary contact resulted in oil recovery exceeding the spontaneous imbibition potential in the outlet-isolated cores by ca. 10% PV. This is explained by viscous recovery provided by water bridges across open fractures. The size of the bridges seemed to be controlled by the wettability of the rock and not by the differential pressure applied across the open fracture.     

Leak-off During Placement of Cr(III)-Partially Hydrolyzed Polyacrylamide Gelling Solution in Fractured Porous Media

During waterflooding of a fractured formation, water may channel through the fracture or interconnected network of fractures, leaving a large portion of oil bearing rock unswept. One remedial practice is injection of a gelling solution into the fracture. Such placement of a gelling mixture (referred as gelant) is associated with leak-off from the face of the fracture into the adjoining matrix. As the gelant gets more crosslinked, the gelant encounters more resistance in flowing into the porous matrix. This article addresses the build-up of flow resistance as the Cr(III)-partially hydrolyzed polyacrylamide gelant, at various stages of crosslinking flows into the matrix. Flow experiments were conducted at constant injection pressure in unfractured Berea rocks that represent a matrix adjoining a fracture. Before entering the core, gelants underwent post-mixing delays, shorter than their gel time. On continued displacement, flow resistance developed that reduced the flow rate further. More delay, after mixing of gelant hastened, the build-up of resistance to flow and the resistance was contained nearer to the inlet face. Effect of flow over fracture face on the build-up of flow resistance in the matrix was also evaluated by conducting displacement of gelant in two fractured slabs. In one case, a part of the injected fluid came out of the fracture outlet with the rest leaking off into matrix. In the other case, all the fluid that entered into the fracture leaked off into the matrix. Build-up of flow resistances in the matrix for the two cases was compared. A simple conceptual model is presented that could explain the flow of gelant and build-up of resistance in porous rock at constant injection pressure.     

Immiscible Displacement of a Wetting Fluid by a Non-wetting One at High Capillary Number in a Micro-model Containing a Single Fracture

Most reservoirs in Iran are heterogeneous fractured carbonate reservoirs. Heterogeneity causes an earlier breakthrough and an unstable front which leads to a lower recovery. A series of experiments were conducted whereby the distilled water displaced n-Decane in strongly oil-wet glass micro-models containing a single fracture. Experimental data from image analysis of immiscible displacement processes are used to modify the Buckley?CLeverett and fractional flow equations by a heterogeneity factor. It is shown that the heterogeneity factor in the modified equations can be expressed as a function of fracture length and orientation.     

致密油藏动态裂缝扩展机理及应用

致密油藏采用注水吞吐补充地层能量取得了一定效果. 但多轮次注水吞吐后, 地层压力和产量降低快. 本文考虑了致密油藏复杂的裂缝形态, 根据艾尔文理论及弹性力学剖析I型裂缝尖端附近的应力场分布, 基于渗流力学、裂缝性致密油藏特征及动态裂缝渗流规律, 建立了多裂缝交叉裂缝扩展渗流模型, 结合注水诱导裂缝扩展机理及断裂力学能量守恒原理, 得到裂缝扩展长度. 依据致密油藏逆向渗吸原理, 提出将注水吞吐转为不稳定脉冲注水. 对比分析注水吞吐、脉冲注水2种能量补充发方式, 预测10年累计采油、压力及剩余油分布. 结果表明, 裂缝净内压随着注水量的增加而升高, 当应力场强度因子达到断裂韧度, 在裂缝尖端会发生扩展. 扩展及延伸的天然裂缝相互沟通, 呈现不规则复杂缝网, 在复杂缝网中主要发生逆向渗吸作用. 脉冲注水累计产油高、注水波及面积广、逆向渗吸作用强. 裂缝性致密油藏水平井注水吞吐转变为脉冲注水方式, 能够充分发挥动态缝网的逆向渗吸及线性驱替作用, 实现有效驱油的目的.