Technical Notes - Development of a Generalized Darcy Equation

General equations relating the pressure drop necessary to sustain the flow of a fluid through a porous matrix at a given rate have been developed. The results indicate that at high values of flow rate the pressure-flow behavior may not necessarily satisfy the usual Darcy equation. The mathematical analysis, carried through the micro-pore geometry and extended through the macro-reservoir scale, indicate that Darcy's law, of limited applicability to certain ranges of Reynolds numbers, can be generalized through the inclusion of some additional parameters. The "generalized Darcy equation" has also been formulated in dimen-sionless form permitting the evaluation of its predictive accuracy with regard to literature data. A comparison between predicted and experimental values indicates that the generalized Darcy equation predicts the pressure drops with good agreement over all possible ranges of Reynolds numbers. INTRODUCTION The limits and the nature of validity of Darcy's law' has been a subject of every-day interest to the industry for many years. It is well known that as the Reynolds number, characteristic of the fluid flow through porous media, becomes large, Darcy's law gradually loses its predictive accuracy and ultimately becomes completely void. For the last 20 years much has been said and written on this subject. Unfortunately little has been accomplished to bring about a satisfactory agreement, at least on the nature of the threshold of validity of Darcy.'s law. Fluid dynamists, geo-physicists, and engineers all had their individual views, explanations, interpretations and concepts on the subject. To some, a mechanistic analogy with pipe-flow proved a satisfactory explanation.' To others,' turbulence, in its random character, was incompatible with the geometric structure of consolidated porous systems. To some,4 turbulence merely represented a factor influencing the permeability measurements and again to others5,6,7 em-pirical or semi-empirical correlations proved satisfactory from an engineering viewpoint. Deviations from Darcy's law at high flow rates have been studied by systematic experiments by Fancher, Lewis, and Barnes.' In an article on the flow of gases through porous metals, Green and Duwezs conclude that the onset of turbulence within the pores appears unsatisfactory to explain deviations from Darcy's law. This view is held by many others. While the subject remained controversial for many years, the development of vast natural gas reserves throughout recent years further justified considerable interest on this problem from the standpoint of gas reservoir behavior. As large amounts of field data became available from the operation of many gas fields, it became evident that the steady-state behavior of gas wells was not, in general, in agreement or compatible with Darcy's law. This suggested a careful reconsideration of all mechanisms which may account for pressure drops in addition to viscous shear. In a series of articles9,10 . Hou-peurt indicated that deviations from Darcy's law may be explained on the basis of kinetic energy variations and jetting effects without resorting to assumptions on turbulent flow conditions. Another article by Schneebeli11 indicates that special experiments by Lindquist clearly demonstrated that the onset of turbulence does not necessarily coincide with conditions of deviation from Darcy's law. This view is also held by M. King Hubbert.12 Starting with the basic pressure-flow relations suggested by Houpeurt, the derivation, development and extension of analytical expressions to -supplement and generalize Darcy's law has been the objective of this work. MATHEMATICAL ANALYSIS Derivation of Dimensionless Pressure-drop, Flow-rate Relations In considering the flow of a fluid through a porous matrix geometrically represented by a succession of capillary passages in the shape of truncated cones,810 an approximate expression may be derived relating viscous and inertial, i.e., total pressure drop to the physical properties of the fluid, geometric properties of the rock matrix and the rate of flow: ?P/?r = µ/k V [ 1 + c(m4 - 1) p V/16n" mµ w] ..........(1) Let us formally set: c (m4 - 1) / 16n" m = a d ......(2) Such a representation is equivalent to assert that the term [c(m4 — 1)/ 16n"m], variable with various porous media and probably highly variable within a given porous medium, may be macroscopically defined as equal to a lithology factor times the aver-age grain diameter d. In view of the usual grain and pore size distribu-