The majority of the formation damage models were developed for single phase fluid systems. This assumption is valid only for very specific cases such as the production of particles with oil flow and for special core tests. Nevertheless, it is instructive to understand these models before looking into the multi-phase effects. Therefore, the various processes involving single-phase formation damage are discussed and the selected models available are presented along with some modifications and critical evaluation as to their practical applicability and limitations. The methodology for determination of the model parameters are presented. The parameters that can be measured directly are identified. The rest of the parameters are determined by means of a history matching technique. The applications of the models and the parameter estimation method are demonstrated using several examples.

An evaluation and comparison of six selected models bearing direct relevance to formation damage prediction for petroleum reservoirs are carried out. The modeling approaches and assumptions are identified, interpreted, and compared. These models are applicable for special cases involving single-phase fluid systems in laboratory core tests.

Porous media is considered in two parts:

1 the flowing phase, denoted by the subscript /, consists of a suspension of fine particles flowing through and

2 the stationary phase, denoted by the subscript s, consists of the porous matrix and the particles retrained.

The Thin Slice Algebraic Model

Model Formulation

Wojtanowicz et al. (1987, 1988) considered a thin slice of a porous material and analyzed the various formation damage mechanisms assuming one distinct mechanism dominates at a certain condition. Porous medium is visualized as having tortuous pathways represented by Nh tubes of the same mean hydraulic equivalent diameter, Dh, located between the inlet and outlet ports of the core as depicted. The crosssectional area of the core is A and the length is L. The tortuosity factor for the tubes is defined as the ratio of the actual tube length to the length of the core.

in which Cl is an empirical shape factor that incorporates the effect of deviation of the actual perimeter from a circular perimeter. As a suspension of fine particles flows through the porous media, tubes having narrow constrictions are plugged and put out of service. If the number of nonplugged tubes at any given time is denoted by Nnp and the plugged tubes by Np, then the total number of tubes is given by:

The Darcy and Hagen-Poiseuille equations given respectively by

are considered as two alternative forms of the porous media momentum equations, q is the flowrate of the flowing phase and Ap is the pressure differential across the thin core slice. Thus, equating Eqs. 10-5 and 10-6 and using Eqs. 10-1 and 10-2 the relationship between permeability, K, and open flow area, A is obtained as:

The permeability damage in porous media is assumed to occur by three basic mechanisms:

1 gradual pore reduction (pore narrowing, pore lining) by surface deposition,

2 single pore blocking by screening (pore throat plugging) and

3 pore volume filling by straining (internal filter cake formation by the snowball effect).

Gradual pore reduction is assumed to occur by deposition of particles smaller than pore throats on the pore surface to reduce the cross-sectional area, A, of the flow tubes gradually as depicted in this article. the number of tubes open for flow, Nnp, at any time remains the same as the total number of tubes, Nh, available. Hence,

Then, using Eq. 10-9 and eliminating A between Eqs. 10-4 and 10-7 leads to the following equation for the permeability to open flow area relationship during the surface deposition of particles:

Single pore blocking is assumed to occur by elimination of flow tubes from service by plugging of a pore throat or constriction, that may exist somewhere along the tube, by a single particle to stop the flow through that particular tube. Therefore, the cross-sectional areas of the individual tubes, Ah, do not change. But, the number of tubes, Nnp, open for the flow is reduced as depicted in Figure 10-3. The area of the tubes eliminated from service is given by:

The number of tubes plugged is estimated by the ratio of the total volume of pore throat blocking particles to the volume of a single particle of the critical size.

The critical particle size, d, is defined as the average size of the critical pore constrictions in the core. fd is the volume fraction of particles in the flowing phase, having sizes comparable or greater than d. pp is the particle grain density. pp/ is the mass concentration of particles in the flowing suspension of particles. Because Ah is a constant, Eq. 10-7 leads to the following permeability to open flow area relationship:

Pore filling occurs near the inlet face of the core when a suspension of high concentration of particles in sizes larger than the size of the pore throats is injected into the core as depicted in this article. The permeability, Kc, of the particle invaded region decreases by accumulation of particles. But, in the uninvaded core region near the outlet, the permeability of the matrix, Km, remains unchanged. The harmonic mean permeability, K, of a core section (neglecting the cake at the inlet face) can be expressed in terms of the permeability, Kc, of the Lc long pore filling region and the permeability, Km, of the Lm long uninvaded region as

Rc(t) and Rm are the resistances of the pore filling and uninvaded regions defined by

The rate of increase of the filtration resistance of the pore filling particles is assumed proportional to the particle mass flux of the flowing phase according to:

kc is the pore filling particles resistance rate constant.The instantaneous porosity of a given cross-sectional area is given by:

(|)0 and (() denote the initial and instantaneous porosity values, e is the fractional bulk volume of porous media occupied by the deposited particles, given by

mp is the mass of particles retained per unit volume of porous media and pp is the particle grain density. For convenience, these quantities can be expressed in terms of initial and instantaneous open flow areas, Afo and Ay> and the area covered by the particle deposits, A , as

Substituting Eqs. 10-25 through 10-27, Eqs. 10-23 and 24 become, respectively

at themlet and outlet of the core. Wojtanowicz et al. (1987, 1988) omitted the accumulation of particles in the thin core slice and simplified Eq. 10-30 to express the concentration of particles leaving a thin section by:

The rate of particle retention on the pore surface is assumed proportional to the particle mass concentrations in the flowing phase according to:

The rate of entrainment of the surface deposited particles by the flowing phase is assumed proportional to the mass of particles available on the pore surface according to:

Then, the net rate of deposition is given as the difference between the retention and entrainment rates as:


Cernansky, A., & Siroky, R. "Deep-bed Filtration on Filament Layers on Particle Polydispersed in Liquids," Int. Chem. Eng., Vol. 25, No. 2, 1985, pp. 364-375.

Cernansky, A., & Siroky, R., "Hlbkova Filtracia Polydisperznych Castic z Kvapalin na Vrstvach z Vlakien," Chemicky Prumysl, Vol. 32 (57), No. 8, 1982, pp. 397-405.

Civan, F. "A Generalized Model for Formation Damage by Rock-Fluid Interactions and Particulate Processes," SPE Paper 21183, Proceedings of the SPE 1990 Latin American Petroleum Engineering Conference, October 14-19, 1990, Rio de Janeiro, Brazil, 11 p.

Civan, F. "Evaluation and Comparison of the Formation Damage Models," SPE 23787 paper, Proceedings of the SPE International Symposium on Formation Damage Control, February 26-27, 1992, Lafayette, Louisiana, pp. 219-236.

Civan, F., & Knapp, R. M. "Effect of Clay Swelling and Fines Migration on Formation Permeability," SPE Paper No. 16235, Proceedings of the SPE Production Operations Symposium, Oklahoma City, Oklahoma, 1987, pp. 475-483.

Civan, F. "A Multi-Phase Mud Filtrate Invasion and WellBore Filter Cake Formation Model," SPE Paper No. 28709, Proceedings of the SPE International Petroleum Conference & Exhibition of Mexico, October 10-13, 1994, Veracruz, Mexico, pp. 399-412.

Civan, F., Knapp, R. M., & Ohen, H. A. "Alteration of Permeability by Fine Particle Processes," J. Petroleum Science and Engineering, Vol. 3, Nos. 1/2, October 1989, pp. 65-79.

Civan, F., Predictability of Formation Damage: An Assessment Study and Generalized Models, Final Report, U.S. DOE Contract No. DE-AC22- 90BC14658, April 1994.

Civan, F. "Modeling and Simulation of Formation Damage by Organic Deposition," Proceedings of the First International Symposium on Colloid Chemistry in Oil Production: Asphaltenes and Wax Deposition, ISCOP'95, Rio de Janeiro, Brazil, November 26-29, 1995, pp. 102-107.

Civan, F. "A Multi-Purpose Formation Damage Model," SPE 31101, Proceedings of the SPE Formation Damage Symposium, Lafayette, Louisiana, February 14-15, 1996, pp. 311-326.

Civan, F. "Interactions of the Horizontal Wellbore Hydraulics and Formation Damage," SPE 35213, Proceedings of the SPE Permian Basin Oil & Gas Recovery Conf., Midland, Texas, March 27-29, 1996, pp. 561-569.

Gruesbeck, C, & Collins, R. E. "Particle Transport Through Perforations," SPEJ, December 1982b, pp. 857-865.

Gruesbeck, C., & Collins, R. E. "Entrainment and Deposition of Fine Particles in Porous Media," SPEJ, December 1982a, pp. 847-856.

Khilar, K. C., & Fogler, H. S. "Colloidally Induced Fines Migration in Porous Media," in Amundson, N. R. & Luss, D. (Eds.), Reviews in Chemical Engineering, Freund Publishing House LTD., London, England, January-June 1987, Vol. 4, Nos. 1 and 2, pp. 41-108.

Khilar, K. C., & Fogler, H. S. "Water Sensitivity of Sandstones," SPEJ, February 1983, pp. 55-64.

Liu, X., Civan, F, & Evans, R. D. "Correlation of the Non-Darcy Flow Coefficient, J. of Canadian Petroleum Technology, Vol. 34, No. 10, 1995, pp. 50-54.

Metzner, A. B., & Reed, J. C. "Flow of Non-Newtonian Fluids—Correlation of the Laminar, Transition, and Turbulent Flow Regions," AIChE J., Vol. 1, No. 4, 1955, pp. 434-440.

Nayak, N. V, & Christensen, R. W. "Swelling Characteristics of Compacted Expansive Soils," Clay and Clay Mineral, Vol. 19, No. 4, December 1970, pp. 251-261.

Ohen, H. A., & Civan, F. "Predicting Fines Generation, Migration and Deposition Near Injection and Production Wells," Proceedings of the First Regional Meeting, American Filtration Society, Houston, Texas, October 30-November 1, 1989, pp. 161-164.

Ohen, H. A., & Civan, F. "Simulation of Formation Damage in Petroleum Reservoirs," SPE Advanced Technology Series, Vol. 1, No. 1, April 1993, pp. 27-35.

Ohen, H. A., & Civan, F. "Simulation of Formation Damage in Petroleum Reservoirs," SPE 19420 paper, Proceedings of the 1990 SPE Symposium on Formation Damage Control, Lafayette, Louisiana, February 22-23, 1990, pp. 185-200.

Schechter, R. S., Oil Well Stimulation, Prentice Hall, Englewood Cliffs, New Jersey, 1992, 602 p.

Seed, H. B., Woodward, Jr., R. J., & Lundgren, R. "Prediction of Swelling Potential for Compacted Clays," /. Soil Mech. Found. Div., Proc. Am. Soc. Civ. Eng., 88(SM3), June 1962, pp. 53-87.

Wojtanowicz, A. K., Krilov, Z., & Langlinais, J. P. "Study on the Effect of Pore Blocking Mechanisms on Formation Damage," SPE 16233 paper, presented at Society of Petroleum Engineers Production Operations Symposium, Oklahoma City, Oklahoma, March 8-10, 1987, pp. 449-463.

Wojtanowicz, A. K., Krilov, Z., & Langlinais, J. P. "Experimental Determination of Formation Damage Pore Blocking Mechanisms," Trans, of the ASME, Journal of Energy Resources Technology, Vol. 110, 1988, pp. 34-42.