632 Reservoir Formation Damage Models Separating the Internal and External Filtration Processes For convenience in the modeling, Pang and Sharma (1994) divide the entire filtration process into two phases: (1) the initial internal cake filtration, and (2) the later external cake filtration. They separate these two filtration phases by a "transition time" after which the particle migration into porous formation becomes negligible and an external filter cake begins forming over the injection well formation face. Sharma and Pang (1997) assumed that transition from the internal to the external cake filtration occurs when the porosity, ((), of the formation face decreases to a minimum critical value, ({>*, by particle deposition, below which particle invasion into porous media is not possible. Their models apply for single phase water flow in the near wellbore region. Hence, the effect of the oil-water two-phase flow during the initial water injection period is neglected because this initial period is relatively short. Transition Time Wennberg and Sharma (1997) estimate the transition time based on the expressions given by Iwasaki (1937) for particle deposition rate and the filtration coefficient, respectively, as: do/dt = (19-14) and (19-15) where b is an empirical constant and X, 0 is the filtration coefficient without particle deposition. Although more sophisticated expressions are available in the literature (See Chapter 8), they used Eq. 19-15 for simplicity. Thus, invoking Eq. 19-15 into Eq. 19-14 yields the following expression, similar to the rate equation for particle deposition in pluggable pathways given by Gruesbeck and Collins (1982): (19-16) They obtain the analytic solution of Eq. 19-16 for constant flow rate and suspension particle concentration for two cases as the following: 1 * 1 f\ / t f\ t T\ <5 = huct , b = 0 (19-1/) Injectivity of the Waterflooding Wells 633 (19-18) At the transition time, the porosity attains the minimum critical porosity, <|>* , and the maximum critical volume fraction of the deposited particles becomes a* = (|) 0 - (j)* . Under these conditions, Eqs. 19-17 and 18 can be used to obtain the following expressions, respectively, for the transition time: (19-19) (19-20) Internal Filtration Models Pang and Sharma (1994, 1995) and Sharma et al. (1997) have pursued their derivations in terms of the following variables: V = M/(]) (19-21) K' = Kv (19-22) n = o/<|) (19-23) Here, v represents the interstitial velocity of the fluid phase and n denotes the fraction of the pore space occupied by the particle deposits in porous media. X' is the product of the deposition rate constant, K, and the interstitial velocity of the flowing suspension. In the following, their formulations are presented in a manner consistent with the rest of the presentation of this chapter. The damage of a core plug by the injection of a dilute particulate suspension can be described by means of the volumetric balance equations of the suspended and deposited particles in porous media, given, respectively, by (Wennberg and Sharma, 1997): „ —+ w —+ —= 0 dt 3;c 3r and (19-24) (19-25) 634 Reservoir Formation Damage subject to the initial and boundary conditions given by (Pang and Sharma, 1994): C ~ C o\ X )-> Q~®o\ X ) i > X ^ V ' t — \J c = c f (t), x = 0 , t>0 The instantaneous porosity is given by: (19-26) (19-27) (19-28) The analytical solution of Eqs. 19-24 through 27 used by Pang and Sharma (1994) implies some simplifications. It applies for the injection of dilute suspension of particles. Therefore, the effect of small amount of particle deposition, a, compared to the initial porosity can be neglected, the deposition rate coefficient is assumed constant, and a constant rate injection is considered. Thus, the analytic solution for constant §~$ 0 , A,~X 0 , and u = u 0 can be adopted from Rhee et al. (1986) as: c(x,t) = c 0 \x-— exp - Kut ut x> — <D (19-29) , x< — (19-30) where the term inside the square brackets expresses that c f is a function of (t-tyx/u). Considering that c f (t) = c f is constant and c 0 (x) = Q and a 0 (jc) = 0 in the laboratory core flow tests, Pang and Sharma (1994) simplify Eqs. 19-29 and 30, respectively, as: c(x,t) = 0, x>utl§ (19-31) c(x,t) = c / exp(-Xjc), x<ut/$ (19-32) and using Eqs. 19-31 and 32, they obtain the solution of Eq. 19-25 as: o(jc,0 = 0, x>ut/ty (19-33) Injectivity of the Waterflooding Wells 635 G(x,t) = c f kutexp(-X*), x>ut/ty (19-34) Pang and Sharma (1995) assumed that the permeability reduction primarily occurs by pore throat plugging. Therefore, they estimated the permeability of the porous formation as a harmonic average permeability of the combined plugged and unplugged regions as: (19-35) where f p represents the volume fraction of the deposited particles con- tributing to pore throat plugging, K p denotes the permeability of the plugged region near the pore throats, and K m is the permeability of the formation matrix, assumed to remain constant, which is equal to the initial nondamaged permeability, K 0 (i.e., K m =K 0 }. Therefore, Pang and Sharma (1995) rearranged Eq. 19-35 for the relative or fractional retained permeability of the porous formation undergoing particle deposition from dilute suspensions as, inferred by Payatakes et al. (1974): K(x,t) = /L where P is an empirical damage factor given by: P = /,(*./*,- 1 ) (19-36) (19-37) Thus, they calculate the harmonic average permeability of the damaged portion of the core by: K(t\- Xf { dx K(t) -T 0 n^t (19-38) Substituting Eqs. 19-33, 34, and 36 into Eq. 19-38, and then inte- grating, they derive the following expression: (19-39) in which 636 Reservoir Formation Damage N = (19-40) Note that Eq. 19-39 was previously derived by Wojtanowicz et al. (1987, 1988) as described in Chapter 10. Applying Eqs. 19-1 through 3, 10, 12, and 13, the injectivity ratio for the linear cases, applicable to laboratory core plugs without external cake formation (/? c = 0), is given by: a(t) = where *(*/) (19-41) x f = ut ut <L (19-42) x f = , —>L (19-43) Pang and Sharma (1994) simplify Eq. 19-41 by considering that the injection front reaches the outlet end of the core rapidly. Therefore, neglecting the damage during the short period of time until the front reaches the core outlet, Eqs. 19-39 and 41 yield for x f - L the following equation indicating that the reciprocal injectivity ratio is a linear function of time: 1 (19-44) External Filtration Models Considering the formation of an incompressible external cake with- out any particle invasion into the core plug, Pang and Sharma (1994) expressed the harmonic average permeability of the cake and the core system (Figure 19-2) as: h h _ K(t) ~ K c (t} K 0 (19-45) Injectivity of the Waterflooding Wells 637 Sharma et al. (1997) determine the thickness of the external filter cake by means of a volumetric balance of the particles in the cake as: h c = 0, t < t * (19-46) h = , t>f (19-47) where A is the cross-sectional area of the core plug and c f is the volume fraction of the fine particles in the water injected at a flow rate of q. For constant c f and q, Eq. 19-47 simplifies as: (19-48) Thus, substituting Eq. 19-48 into 45 and considering that the filter cake thickness is much smaller than the length of the core plug (i.e., h c « L), they obtained the following expression indicating that the reciprocal injectivity ratio is a linear function of time: 1 <x(0 where (19-49) (j>L (19-50) Filtration Coefficient Wennberg and Sharma (1997) point out that the filtration coefficient varies by particle deposition according to: (19-51) in which K 0 is the filtration coefficient with no deposited particles and F(o) is a function of the volume fraction of particles deposited. Their 638 Reservoir Formation Damage review of the various expressions available for prediction of the filtration coefficient is summarized and presented in the following. Ives (1967) proposed a general expression as: (19-52) in which x, v, z, and p are some empirical parameters and a M is the maximum of the volume fraction of the deposited particles necessary to make the filtration coefficient of porous media zero. This equation indicates that the filtration coefficient is equal to one when there is no deposited particles in porous media, and the filtration coefficient becomes zero when the volume fraction of deposited particles reaches a certain characteristic value of maximum a M . Chiang and Tien (1985) developed an empirical correlation as: X 2 )>a<10~ 2 (19-53) where N r is the suspended particle to porous media grain diameter ratio: N r = D p /D g (19-54) Rajagopalan and Tien (1976) developed the following expression: A, = Q. + 2.4x10' <0.18 (19-55) in which the dimensionless groups are defined as following. A s is Happel's dimensionless geometric parameter. A^ is the London parameter given by (19-56) where H = 3.0x10 I3 erg is Hamaker's constant. N g is the gravity number given by: (19-57) Injectivity of the Waterflooding Wells 639 N pe is the peclet number given by: N pe =vD g /D (19-58) v is velocity, ILL is viscosity, and p p and p e denote the particle and fluid densities, respectively. (3 is a packing parameter given by ,1/3 (19-59) where 0 is the porosity in fraction. D is the coefficient of diffusion for the Brownian motion of particles. Wennberg and Sharma (1997) analyzed the measurements of the filtration coefficient reported by various investigators and determined that these data mostly indicate power law-type relationships to the volumetric flux, the suspended particle size, and the porous media grain size as: X~M- a ,Z)P,zy (19-60) where a, (3, and y are some empirical exponents. They determined that 0<a < 2; (3 >0 for D p > liim and (3 < 0 for D p < Ijim; 0.6 <|3 < 1.2 for Ottawa sand; and y = 0.9 and 2. Diagnostic-Type Curves for Water Injectivity Tests Pang and Sharma (1994, 1997) identified four distinct type curves that can be used for interpretation of the water-quality tests. They justified these type curves with experimental data obtained from the literature as shown in Figure 19-3. Type curve 1 is a straight line indicating the formation of an incompressible external filter cake or a thin internal cake near the injection face of the core plug according to Eq. 19-49. The slope remains constant. Type curve 2 is for the similar case, but applies for compressible cakes. In this case, the porosity and permeability of the cake decrease by increasing filtration pressures. As a result, the slope of the curve increases with the filtration time or pore volume injected. Type curve 3 refers to a deep particle invasion and pore filling in the core plug, leading to a slower gradual permeability decrease. As a result, the slope of the curve decreases with the filtration time. Type curve 4 may be an S-shaped or other types of curves indicating a shift of the dominance of the different damage mechanisms during the filtration process. 8 •5 •a 500 1000 1500 Pore Volumes Injected (a) 2000 4000 8000 Pore Volumes Injected (b) 12000 o O P 1-4 ,3 1-3 1.2 1.1 2000 4000 6000 Pore Volumes Injected (c) 8000 o £ g +3 B S •8 I 10 20 30 40 Pore Volumes Injected (d) Figure 19-3. Diagnostic-type curves: (a) Type curve 1 (data from Todd et al., 1979), (b) Type curve 2 (data from Todd et al., 1984), (c) Type curve 3 (data from Todd et al., 1979), and (d) Type curve 4 (data from Pautz et al., 1989) (after Pang and Sharma, ©1997 SPE; reprinted by permission of the Society of Petroleum Engineers). Injectivity of the Waterflooding Wells 641 Models for Field Applications Pang and Sharma (1994, 1995, 1997) have also applied the above methodology and derived the models for prediction of the injectivity ratio given in Table 19-1 for other cases. Specifically, the open-hole, perforated and fractured wells as depicted in Figures 19-1, 19-4 and 19-5 are con- sidered separately in Table 19-1 for the internal and external filtration phases. Figure 19-4. Details of perforation for analysis of internal filtration (after Pang and Sharma, ©1995 SPE; reprinted by permission of the Society of Petroleum Engineers). wkf Lf Figure 19-5. Details of vertical hydraulic fracture for analysis of internal filtration (after Pang and Sharma, ©1995 SPE; reprinted by permission of the Society of Petroleum Engineers). [...]... "Characterization and Prediction of Formation Damage in Two-Phase Flow Systems, SPE 25429 paper, proceedings of the SPE Production Operations Symposium, March 21- 23, 19 93, Oklahoma City, Oklahoma, pp 231 -248 Liu, X., & Civan, F, "Formation Damage and Skin Factors Due to Filter Cake Formation and Fines Migration in the Near-Wellbore Region," SPE 2 73 6 4 paper, Proceedings of the 1994 SPE Formation Damage Control Symposium,... 4-6, 19 87 Wennberg, K E., & Sharma, M M., "Determination of the Filtration Coefficient and the Transition Time for Water Injection Wells," SPE Paper 38 181, Proceedings of the 19 97 SPE European Formation Damage Conference held in the Hague, The Netherlands, June 2 -3, 19 97, pp 35 3 -36 4 Wojtanowicz, A K., Krilov, Z., & Langlinais, J P., "Study on the Effect of Pore Blocking Mechanisms on Formation Damage, "... Oil/Gas FALSE +0 +15 TRUE Injector Well FALSE +0 +6 ^ >30 % Old ^> 30 % Condition Young +0 ^ 30 % +10 Age TVD 15,000/t ^ 30 % 15% 0.6 pn/ft 0 .3 pn/ft < 0 .3 pn/ft O 2 •4^ O U OS -« c 1 » •* -S aa «*& e! zz OPQQ Q< zz Q< W Reservoir Sand Migration and Gravel-Pack Damage o . /,(*./*,- 1 ) (19 -36 ) (19 - 37 ) Thus, they calculate the harmonic average permeability of the damaged portion of the core by: K(t- Xf { dx K(t) -T 0 n^t (19 -38 ) Substituting Eqs. 19 -33 , 34 , and 36 . Study," SPE Paper 38 180, Proceedings of the 19 97 SPE European Formation Damage Conference held in the Hague, The Netherlands, June 2 -3, 19 97, pp. 34 1 -35 1. Todd, A. C., et . 36 into Eq. 19 -38 , and then inte- grating, they derive the following expression: (19 -39 ) in which 636 Reservoir Formation Damage N = (19-40) Note that Eq. 19 -39 was previously