Gruesbeck and Collins (1982a) developed a partial differential model based on the concept of parallel flow of a suspension of particles through plugging and nonplugging pathways, as depicted in this article. Relatively smooth and large diameter flowpaths mainly involve surface deposition and are considered nonplugging. Flowpaths that are highly tortuous and having significant variations in diameter are considered plugging. In the plugging pathways, retainment of particles occurs by jamming and blocking of pore throats when several particles approach narrow flow onstrictions.

Sticky and deformable deposits usually seal the flow constrictions (Civan, 1994, 1996). Therefore, conductivity of a flow path may diminish without filling the pore space completely. Thus, the fluid seeks alternative flow paths until all the flow paths are eliminated. Then the permeability diminish even though the porosity may be nonzero. Another important issue is the criteria for jamming of pore throats. As demonstrated by Gruesbeck and Collins (1982b) experimentally for perforations, the probability of jamming of flow constrictions strongly depends on the particle concentration of the flowing suspension and the flow constrictionto-particle diameter ratio. Gruesbeck and Collins (1982a) assumed that the liquid and particles have constant physical properties. The porous media is incompressible, homogeneous and isotropic.

There is hydraulic communication through the interconnectivity of the plugging and nonplugging pathways and therefore the pressure gradients and the particle concentrations of the suspension flowing through the plugging and nonplugging pathways are the same. The volume flux through the core is constant and only the external particle invasion is considered. The flow through porous media was assumed to obey the Darcy Law. In this section, the Gruesbeck and Collins (1982) model is presented with the modifications and improvements made by Civan (1995). (J)PO and §npg denote the initial pore volume fractions of the plugging and nonplugging pathways of the porous media (Civan, 1994, 1995).

These values can be determined experimentally for a given porous media and the particle size distribution. εp and εnp represent the fractions of the bulk volume occupied by the deposits. Thus, the instantaneous porosities are:

The fractions of the bulk volume containing the plugging and nonplugging pathways can be approximated, respectively, by:

However, Gruesbeck and Collins (1982) assume a constant value for ƒp (and therefore for ƒnp =l-ƒp), which is a characteristic of the porous medium and the particles. Total instantaneous and initial porosities are given, respectively, by:

The total deposit volume fraction and the instantaneous available porosity are given, respectively, by:

The rate of deposition in the plugging pathways is given by assuming the pore filling mechanism:

when the pore throat-to-particle diameter ratio is below its critical value. Civan (1990, 1994) recommended the following empirical correlation:

which is determined empirically as a function of the particle Reynolds number:

The rate of deposition in the nonplugging tubes is given as the difference between the rates of surface deposition and sweeping of particles (Civan, 1994):

where ðnp is the volume fraction of particles in the suspension of particles flowing through the nonplugging pathways. kd and ke are the surface deposition and mobilization rate constants. ke = 0 when tw<tcr. ɲe is the fraction of the uncovered deposits that can be mobilized from the pore surface, estimated by:

tcr is the minimum shear stress necessary to mobilize the surface deposits. inp is the wall shear-stress in the nonplugging tubes, given by the Rabinowitsch-Mooney equation (Metzner and Reed, 1955):

and the mean pore diameter is given by:

where C is an empirical shape factor. It can be shown that Eqs. 10-101 and 105 simplify to the deposition rate equations given by Gruesbeck and Collins (1982):

This requires that the effects of the permeability and porosity changes be negligible, the fraction of the uncovered deposits be unity, the suspension of particles be Newtonian, and the particle volume fractions of the suspensions flowing through the plugging and nonplugging pathways be the same, that is,

The permeabilities of the plugging and nonplugging pathways are given by the following empirical relationships (Civan, 1994):

Then, the average permeability of the porous medium is given by:

Note that Eq. 10-116 was derived independently by Civan (1992) and Schechter (1992) and is different than the expression given by Gruesbeck and Collins (1982). The superficial flows in the plugging and nonplugging pathways are given, respectively, by:

Considering that the physical properties of the particles and the carrier liquid are constant, the volumetric balance of particles in porous media is given by:

Substituting Eq. 10-100 into Eq. 10-119, and then rearranging, an alternative convenient form of Eq. 10-119 can be obtained as:

Following Gruesbeck and Collins (1982), Eq. 10-120 can be simplified for cases where e and a are small compared to ɸ0 and unity, respectively, and for constant injection rate, as:

The initial particle contents of the flowing solution and porous media are assumed zero:

where L is the length of porous medium. The particle content of the suspension of particles injected into the porous media is prescribed as:

Alternatively, the pressures of the inlet and outlet ends of the porous media instead of the flow rate can be specified. Then, the volumetric flux can be estimated by the Darcy law:

Then, the pressure obtained by solving Eqs. 10-126 through 10-128 is substituted into Eq. 10-124 to determine the volume flux. The preceding formulation of Eq. 10-119 or 120 applies to the overall system following Gruesbeck and Collins' (1982) assumption that the particle concentrations in the plugging and nonplugging pathways are the same according to Eq. 10-113. When different concentrations are considered, Eq. 10-120 should be applied separately for the plugging and nonplugging paths, respectively, as suggested by Civan (1995):

k is a particle exchange rate coefficient. A solution of Eqs. 10-129 through 132 along with the particle deposition rate equations, Eqs. 10-101 and 105, yields the particle volume fractions in the plugging and nonplugging flow paths.


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.