<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJFD</journal-id><journal-title-group><journal-title>Open Journal of Fluid Dynamics</journal-title></journal-title-group><issn pub-type="epub">2165-3852</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojfd.2012.24015</article-id><article-id pub-id-type="publisher-id">OJFD-25719</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  External Cake Build Up at Surface of Porous Medim
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ooman</surname><given-names>Fallah</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>Ahmadi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohammad</surname><given-names>Afkhami Karaee</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hasan</surname><given-names>Rabbani</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Petroleum Engineering, Firoozabad Branch, Islamic Azad University, Firoozabad, Iran</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>hooman.fallah2@gmail.com(OF)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>13</day><month>12</month><year>2012</year></pub-date><volume>02</volume><issue>04</issue><fpage>145</fpage><lpage>148</lpage><history><date date-type="received"><day>August</day>	<month>4,</month>	<year>2012</year></date><date date-type="rev-recd"><day>September</day>	<month>15,</month>	<year>2012</year>	</date><date date-type="accepted"><day>September</day>	<month>23,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Deep bed filtration of carrier water with particles occurs in numerous processes of environmental, chemical, petroleum and civil engineering. For example in the petroleum industry, deep bed filtration of drilling fluid happens during well drilling; it also takes place near to injection wells during seawater injection or produced water re-injection causing in- jectivity reduction. In this work, the filtration theory and its mathematical model are given and the filtration coefficient will be discussed. Then the transition time, the time at which the initial layer of external filter cake is completely formed has been derived mathematically. The transition time can be obtained by calculating the trapping efficiency as a function of the number of previously deposited particles where very few particles can invade the rock and the external filtration begins to build. The prediction of particle capturing and particle retention by mathematical model is an essen- tial stage during planning and design of above-mentioned industrial processes.
 
</p></abstract><kwd-group><kwd>Porous Media; Trapping Efficiency; External Filtration; Transition Time; Filtration Coefficient</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Most oil and gas reservoirs have a natural water layer called formation water beneath the hydrocarbon layer. Also to enhanced oil recovery, seawater is usually injected into the reservoir, which may be associated with oil droplets and particles in the injected water causing injectivity reduction. Damage due to flow of dilute oil-inwater emulsions in porous media depends on several factors. The droplets concentration, the permeability of the core, the rate of the injection water and the properties of the crude oil all affect the permeability impairment.</p><p>Generally, during the particulate suspension flow and colloids through porous medium, the solid and the liquid phases in the suspension can be separated. The retention process of particles when flowing through porous medium is called Deep Bed Filtration. The deep filtration occurs because of several mechanisms: the contacting of particles with the retention site, the fixing of particles sites, and the breaking away of previously retained particles.</p><p>The problem of formation damage is split into two parts which includes external cake build up and internal filtration. In this study, deep filtration, its mathematical modelling and the effect of the different parameters are explained. Then the transition time, the time at which the initial layer of external filter cake is completely formed will be explained. Finally, the transition time calculation has been done, mathematically.</p></sec><sec id="s2"><title>2. Filtration Mechanism</title><sec id="s2_1"><title>2.1. Interception</title><p>It happens when the particle following a streamline hits the surface of a grain and attaches to it. Particles with a density equal to the carrier fluid density follow the streamline in the porous medium at the low velocity. The particles will be brought into contact with grain surfaces, because of the finite size of the particles.</p><p>When the particle is retained by a previously deposited particle we refer to it as bridging and the simplest bridging process that exists is blocking or sieving of the mobile particles.</p></sec><sec id="s2_2"><title>2.2. Impaction</title><p>When the density of particle is larger than that of the fluid, inertia forces will deviate the particle from the stream line and it may attach to the surface of a grain. The impaction is more significant for air filtration and not for water filtration.</p></sec><sec id="s2_3"><title>2.3. Diffusion</title><p>Small particles do random Brownian motion that increases the number of collisions between particles and grains. Diffusion is usually neglected for larger particles.</p></sec><sec id="s2_4"><title>2.4. Size Exclusion</title><p>When a particle enters a throat that is too small for it to pass, it gets stuck there. So the large particles are captured in small pores and pass through large pores. Several attempts to correlate the formation damage with size of particles and pores were unsuccessful. The phenomenological model for average concentration is not general/ universal enough.</p></sec><sec id="s2_5"><title>2.5. Surface Forces</title><p>Double layer repulsion forces and Van Der Walls are important to capture the particles. These forces will play role in the particle in contact with the grain surface. The dominance of the force is depending on filtration medium properties, particle properties and flow [<xref ref-type="bibr" rid="scirp.25719-ref1">1</xref>].</p></sec></sec><sec id="s3"><title>3. Filtration Theory</title><p>The filtration theory is comprised of a system of two equations: the mass conservation equation and a kinetic equation.</p><sec id="s3_1"><title>3.1. Mass Conservation</title><p>The volumetric balance conducted on either the liquid phase or the solid particles yields the following one dimensional convection-dispersion reaction equation [<xref ref-type="bibr" rid="scirp.25719-ref2">2</xref>].</p><disp-formula id="scirp.25719-formula108515"><label>(1)</label><graphic position="anchor" xlink:href="5-2320028\f492005f-48d7-42a8-8b19-8ff6f661644d.jpg"  xlink:type="simple"/></disp-formula><p>where t is time; x is the spatial coordinate in the direction of flow; <img src="5-2320028\34cbfbe0-0e7e-4c54-b157-e0ae4ac97b57.jpg" />is the effective porosity; c is the volumetric concentration of suspended particles with respect to the pore volume; <img src="5-2320028\43cbb78d-0a1f-401f-b1cd-2de90b141f4c.jpg" />is the volumetric concentration of retained particles with respect to the bulk volume.</p><disp-formula id="scirp.25719-formula108516"><label>(2)</label><graphic position="anchor" xlink:href="5-2320028\6578f703-b617-4756-960a-e491532d12c7.jpg"  xlink:type="simple"/></disp-formula><p>Here <img src="5-2320028\e8ddd8ad-c8b3-4867-9fc3-b3cbd0231bad.jpg" /> is the permeate superficial velocity.</p></sec><sec id="s3_2"><title>3.2. Kinetic Equation</title><p>The kinetic equation describes the dynamic interaction between the suspended particles, the deposited particles and the porous medium. Iwasaki (1937) introduced the filtration coefficient, λ, in his formulation of the kinetic equation [<xref ref-type="bibr" rid="scirp.25719-ref3">3</xref>].</p><disp-formula id="scirp.25719-formula108517"><label>(3)</label><graphic position="anchor" xlink:href="5-2320028\3f76078a-f461-4bf0-8652-84534f89bd80.jpg"  xlink:type="simple"/></disp-formula><p>It is generally accepted het for deep bed filtration phenomena covering the water injection field the filtration function is dependent on previous retention <img src="5-2320028\d814c102-f826-4b5e-a934-9a8ea18e2294.jpg" /> only. Thus,</p><disp-formula id="scirp.25719-formula108518"><label>(4)</label><graphic position="anchor" xlink:href="5-2320028\f6cec401-4961-4e2f-af5e-efc6d8cac54b.jpg"  xlink:type="simple"/></disp-formula><p>So, the filtration theory in the field of petroleum engineering, is defined as following equations:</p><disp-formula id="scirp.25719-formula108519"><label>(1)</label><graphic position="anchor" xlink:href="5-2320028\bbb6b5e3-c19a-48dd-a3d5-9201992526f2.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25719-formula108520"><label>(2)</label><graphic position="anchor" xlink:href="5-2320028\eb0d4457-d12e-43e5-8562-ad1e2aad1c54.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25719-formula108521"><label>(4)</label><graphic position="anchor" xlink:href="5-2320028\e6a70b94-206b-4dfb-9879-c69e4c8f7c24.jpg"  xlink:type="simple"/></disp-formula><p>The net increase of particles within an infinitesimal filter volume equals the particles entering the volume on the upstream side minus the particles leaving the volume on the downstream side minus the particles leaving the volume on the downstream side minus the loss of particles to the formation through deposition.</p></sec></sec><sec id="s4"><title>4. Filtration Coefficiant</title><p>The filtration coefficient changes as particle deposition proceeds. This can be written in as&#160;</p><disp-formula id="scirp.25719-formula108522"><label>(5)</label><graphic position="anchor" xlink:href="5-2320028\e550752e-2b41-4c10-a16c-e75ddf745542.jpg"  xlink:type="simple"/></disp-formula><p>The filtration coefficient is a dynamic property that changes with the specific deposit,<img src="5-2320028\b054b38f-9c0d-4c28-8634-f2a4a37ba150.jpg" />. The value of λ increases during the initial stages of filtration, this process called filter ripening. Iwasaki described the ripening period with an equation</p><disp-formula id="scirp.25719-formula108523"><label>(6)</label><graphic position="anchor" xlink:href="5-2320028\ffccdd4e-3e2b-435e-8b0a-17219b0969e3.jpg"  xlink:type="simple"/></disp-formula><p>The deposited particles can themselves act as collectors for subsequent retention of flowing particles. If, however, the surface forces between the deposited and suspended particles are strongly repulsive, λ is likely to decrease with time as the collectors are progressively coated with deposited particles. If particles continue to deposit, they will finally begin bridging the pore throats and if the conditions allow a majority of the pore throats to be bridged, the external filter cakes starts forming.</p></sec><sec id="s5"><title>5. Transition Time <img src="5-2320028\6b33a4af-7442-412c-9b2d-1eff619228d8.jpg" /> Concept</title><p>In developing filtration models, both internal and external filter cakes need to be accounted for because both are generally present in the filtration process. We postulate that during some initial period, an internal filter cake is formed. As more particles are trapped on the surface of the rock, a point is reached where very few particles can invade the rock and an external filter cake begins to build (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p><p>We refer to the time at which no more particles invade the rock, i.e. the time at which the initial layer of external</p><p>filter cake is completely formed, as the transition time<img src="5-2320028\9dc6ee6b-f556-4c70-afad-d325d4cc3d6a.jpg" />. The transition time can be obtained by calculating the trapping efficiency as a function of the number of previously deposited particles. The trapping efficiency is defined as the fraction of suspended particles trapped in the porous media [<xref ref-type="bibr" rid="scirp.25719-ref4">4</xref>].</p><p>If we can determine the conditions under which particles will form the internal and external filter cakes and the time required to form the initial layer of external filter cake<img src="5-2320028\fa77eb9c-12cf-4260-bc37-f227a641d44e.jpg" />, then the entire filtration process can be approximated by applying an internal cake-filtration model before the transition time and an external cakefiltration model after the transition time. Pure external filtration can be obtained in the limit <img src="5-2320028\b6288973-3d58-4d5e-af55-0d572edcd481.jpg" /> and pure internal filtration can be obtained in the limit<img src="5-2320028\b5b163d2-1739-45d7-b593-8c457bd139df.jpg" />.</p><p>Some authors declared that the effect of penetrated particles on formation permeability may be dominant only at the early stages of injection [<xref ref-type="bibr" rid="scirp.25719-ref5">5</xref>]. Once the cake is established and reaches a characteristic thickness, these effects are generally negligible compared to those of filter cake properties.</p>A Brief Review of Injectivity Decline Models<p>If the impairment is due to only to a filter cake buildup, the total pressure drop, <img src="5-2320028\5c8d96f6-3e4e-4cd8-9f47-7c2bb0186223.jpg" />(constant), can be equated to the sum of pressure drop across the filter cake, <img src="5-2320028\3b9ddedc-6d62-41f8-8c75-56b3fe4133be.jpg" />, and across the formation boundary,<img src="5-2320028\241bd941-66e3-4238-958d-03553922b729.jpg" />. With continuous injection, depending on the filter cake permeability and thickness, only the pressure drop across the filter cake increases. This reduces the formation (matrix) pressure drop, resulting in significant losses in injectivity [<xref ref-type="bibr" rid="scirp.25719-ref6">6</xref>].</p><p>Barkman and Davidson developed an equation which expressed the half-life of an injector as a function of formation characteristics, injection rate, Particles concentration and filter cake properties [<xref ref-type="bibr" rid="scirp.25719-ref7">7</xref>].</p><p>More assumption has been used in modeling the external filter cake and concluded since the thickness of the filter cake is small for waterflooding, due to lower solids concentration, the filter cake resistance can be assumed linear rather than radial.</p><p>The limitation of Barkman and Davidson’s and Eylander’s methods results from the assumptions that the matrix of the filter cake is Incompressible and the cake thickness varies linearly with time. However, must particles encountered downhole could be highly compressible (scale products and corrosion).</p><p>Khatib tried to find the relationship between permeability and porosity in order to calculate the pressure drop across the filter cake and thus determine the injectivity decline curve [<xref ref-type="bibr" rid="scirp.25719-ref5">5</xref>]. He the used improved model and concluded that the cake porosity is not constant and varies with applied pressure, particularly for compressible particles such as iron hydroxide and the filter cake characteristics of the solids can be defined by developing empirical permeability/porosity correlations.</p></sec><sec id="s6"><title>6. Simplified Approach for Water Re-Enjection Projects</title><p>To summarize, three possible scenarios exist for the time evolution of the filtration coefficient [<xref ref-type="bibr" rid="scirp.25719-ref8">8</xref>]:</p><p>1) λ increases monotonically. At some point in the time, the initial layers of grains no longer permit particles to pass. The filtration coefficient at this position diverges to infinity while the collection efficiency, η, goes to one, and external cake starts forming which we refer this time as the transition time.</p><p>2) after an initial ripening period, λ starts to decrease and approach zero and no external cake starts builing.</p><p>3) λ is constant or decrease slowly, the pore throats start to bridge, leading to external cake formation. The transition time can, however be orders of magnitude larger than for case 1.</p>The Transition Time<p>Here, we operate in terms of the filtration coefficient, rather than the collection efficiency and postulate that the time development of any filtration theory can be approximately described by the simplest stated equation:</p><disp-formula id="scirp.25719-formula108524"><label>(6)</label><graphic position="anchor" xlink:href="5-2320028\647a2ea2-84ba-48a5-ac15-79eacc85f469.jpg"  xlink:type="simple"/></disp-formula><p>To derive the transition time, <img src="5-2320028\8cfa13c7-4f03-42dc-8abd-a75accc10f4c.jpg" />, from Equation (6) we solve the differential equation</p><disp-formula id="scirp.25719-formula108525"><label>(7)</label><graphic position="anchor" xlink:href="5-2320028\888d510e-1104-422d-a3e9-2fd2606a2509.jpg"  xlink:type="simple"/></disp-formula><p>And find the solution for constant flow rate is</p><disp-formula id="scirp.25719-formula108526"><label>(8)</label><graphic position="anchor" xlink:href="5-2320028\5e42e9a9-0657-4024-9de0-80cbcf438bcf.jpg"  xlink:type="simple"/></disp-formula><p>Then, we state the only a certain percentage of the porous medium can be filled by particles. When the critical porosity, <img src="5-2320028\76fa1fb8-0b43-4bb9-99f6-64640c37b5cb.jpg" />, is achieved, no particles can go through. We denote this critical porosity and note that the transition time corresponds to <img src="5-2320028\0a92cc2a-385e-4229-995e-e391e5527e10.jpg" /> and insert it to Equation (8) will give us the transition time.</p><disp-formula id="scirp.25719-formula108527"><label>(9)</label><graphic position="anchor" xlink:href="5-2320028\99a5cf63-6beb-423d-8396-08b814a3e22a.jpg"  xlink:type="simple"/></disp-formula><p>Little experimental data are available on the transition time and the corresponding specific deposit just below the surface. The critical porosity <img src="5-2320028\f20f4a6c-7422-4677-ae9d-1070c9e46f63.jpg" /> therefore needs to be estimated from theoretical arguments.</p><p>If the initial pore space of a formation is completely filled with particles, the remaining porosity will be the product of formation porosity and the filter cake porosity. A reasonable guess could be that about 50% of the pore space is being filled before the cake starts forming <img src="5-2320028\79fb705c-e845-4dbe-b25e-ac8603fe1548.jpg" />.</p></sec><sec id="s7"><title>7. Conclusions</title><p>Both internal and external filter cakes are considered in the new model by introducing the term, transition time <img src="5-2320028\a37943d6-43ea-4819-baa0-db08c44cbb73.jpg" /> that has been derived mathematically.</p><p>The external filtration build up in front of the surface of the porous medium can be described by a filtration coefficient <img src="5-2320028\bb56a453-3a3d-48f3-949c-e86fc1212145.jpg" /> that may vary with time. The initial filtration coefficient can be determined from correlations derived for conditions under which attractive forces act between formation grains and injection particles.</p><p>Where the repulsive forces exist, the filtration coefficient must be determined from experiments. The cake porosity is not constant and varies with applied pressure, particularly for compressible particles.</p></sec><sec id="s8"><title>REFERENCES</title></sec><sec id="s9"><title>Nomenclator</title><p><img src="5-2320028\3f100546-c555-419f-a8af-4031131e80df.jpg" />: Time</p><p><img src="5-2320028\07220f67-d174-4e95-a019-379ecfa4380d.jpg" />: Spatial coordinate in the direction of flow</p><p><img src="5-2320028\0a6dd8da-27c1-4fb4-a3d8-acffa84333c9.jpg" />: Effective porosity</p><p><img src="5-2320028\679514cc-57b8-4b7d-94c6-7cae8c763057.jpg" />: Volumetric concentration of suspended particles with respect to the pore volume</p><p><img src="5-2320028\1393e3b2-56c0-48ae-a096-1ae3e92c4c16.jpg" />: Volumetric concentration of retained particles with respect to the bulk volume</p><p><img src="5-2320028\4b2310a5-41b1-4a03-93b2-9aa1f03823ed.jpg" />: Permeate superficial velocity</p><p><img src="5-2320028\bd4eadbe-d4f7-4121-812d-b1521579fc2c.jpg" />: Filtration coefficient</p><p><img src="5-2320028\7c042e46-48ec-4b1f-a455-19bf3c6a638b.jpg" />: Transition Time</p><p><img src="5-2320028\ce637ec2-984a-4fc1-be59-beb4f4603de4.jpg" />: Critical porosity</p></sec></body><back><ref-list><title>References</title><ref id="scirp.25719-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">R. Farajzadeh, (PWRI), “An Experimental Investigation into Internal Filtration and External Cake Build Up,” MSc. Thesis, Delft University of Technology, Delft, 2004.</mixed-citation></ref><ref id="scirp.25719-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">J. P. Herzig, D. M. Leclerc and P. Le Goff, “Flow of Suspension through Porous Media—Application to Deep Filtration,” Journal of Industrial and Engineering Chemistry, Vol. 65, No. 5, 1970, pp. 8-35.  
doi:10.1021/ie50725a003</mixed-citation></ref><ref id="scirp.25719-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">T. Iwasaki, “Some Notes on Sand Filtration.” Journal of the American Water Works Association, Vol. 29, 1937, pp. 1591-1602.</mixed-citation></ref><ref id="scirp.25719-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">P. Bedrikovetsky and F. D. Siqeira, “Modified Particle Detachment Model for Colloidal Transport in Porous Media,” Transport in Porous Media, Vol. 86, No. 2, 2011, pp.383-413.</mixed-citation></ref><ref id="scirp.25719-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Z. I. Khatib, “Prediction of Formation Damage Due to Suspended Solids: Modeling Approach of Filter Cake Build up in Injectors,” Paper SPE 28488, 1994.</mixed-citation></ref><ref id="scirp.25719-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">S. Pang and M. M. Sharma, “Supplement to SPE 28489, A Model for Predicting Injectivity Decline in Water Injection Wells,” SPE Formation Evaluation, Vol. 12, No. 3, 1997, pp. 194-201.</mixed-citation></ref><ref id="scirp.25719-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">J. H. Barkman and D. H. Davidson, “Measuring Water Quality and Prediction Well Impairment,” Journal of Petroleum Technology, 1972, pp. 865-868.</mixed-citation></ref><ref id="scirp.25719-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">F. Al-Abduwani, W. M. G. van den Broek and P. K. Currie, “Laboratory and Mathematical Study of Produced Water Re-Injection in Oil Reservoirs,” European Water Management Symposium, Aberdeen, 21-22 May 2004.</mixed-citation></ref></ref-list></back></article>