Search Documents

By R. A. Greenkorn

This project originated in a practical problem—we needed five tracers that could be used together to locate flow paths in a pilot flood. While tracers for subsurface liquids have been used since the turn of the century,1-15 none of those reported in the literature appeared to be either consistent or quantitative enough for our purposes. Most were used in field systems without controlled experiments to determine the accuracy and precision of analysis, and many were tracers that could not be used collectively. The ideal tracer, of course, would follow the fiuid of interest exactly, traveling at the same velocity as the fluid front. But the ideal is impractical to attain because adsorption-desorption effects cause the tracer to lag behind the front; these effects, plus diffusion-dispersion effects, cause the tracer front to spread more than the fluid front. Thus, our objective was not to locate a tracer that would be ideal for all circumstances but rather, to find one that would approximately follow the fluid, or one that under controlled conditions could be corrected to calculate the movement of the fluid front. We considered tracers satisfactory—(1) if they were easy to analyze; (2) if their breakthrough-elution curves were not too different from those for the chloride ion, a tracer believed to follow the fluid front closely; and (3) if we could calculate from the curves a material balance, at 1.25-PV (pore volume) injected, within 5 per cent of that calculated from chloride curves. Of a possible 35 materials, we selected 13 tracers that could be quickly and easily identified and whose analysis was claimed to be accurate within 5 per cent. Only one of these, tritiated water, was a radioactive tracer. Radioactive tracers are easy to detect even in small quantities, but they require special handling and special equipment. Also, those that can be used together are limited because special equipment is required to separate emissions from the various tracers. Two of the original 13 tracers were eliminated in static tests to determine how accurately and precisely they could be analyzed, and to check on gross adsorption. The remaining 11 were flowed through a 9-ft linear sandstone model, and breakthrough-elution curves were obtained. Finally, three tracers were field-tested as breakthrough tracers. These tests are described in the following sections. The 13 tracers considered in these experiments were EDTA (ethylene diamine tetra acetic acid), fluo-rescein, picric acid, salicylic acid and ammonium, boron (as borate), bromide, dichromate, iodide, nitrate and thiocyanate ions, plus chloride ion and tritiated water. All but the chloride ion and tritiated water were subjected to static sand tests to eliminate the tracers that could not be analyzed quickly and accurately and to eliminate those that showed excessive adsorption. All but two of the tracers, EDTA and salicylic acid, qualified for flow tests on this basis. The procedures used in the static tests (see Appendix A for analytical details) were as follows. Each tracer (in the amount shown in Table 1) was dissolved in 1,000 ml of water. Then, 800 ml of this solution was added to 500 gm of sand, and the remaining 200 ml was reserved as a control standard. The sand-tracer mixture was agitated once a day over a 10-week period, and during this period we analyzed five duplicate samples from this mixture, along with five duplicate samples from the control solution. From the control-sample analyses, we constructed a control chart on which the sand-tracer analysis results were plotted. The control chart, and its use to determine accuracy and precision of analytical results, as well as the amount of adsorption, is described in Appendix B. The results of the static tests are summarized in Table 1, which lists the concentration c (concentration at time t) divided by c, (initial concentration). If there was no adsorption, c/c, = 1. The acceptable deviation from this value varied from tracer to tracer, depending upon the limits of analytical precision established for the different tracers. The best results were obtained with boron, bromide and dichromate, with all values falling within limits of accuracy and precision. Ammonium, iodide, nitrate and picric acid also were satisfactory, although the ammonium results were erratic. Flourescein and thiocyanate were unsatisfactory ini-

By Frank W. Harris

The strength of hard drawn copper wire is a question of considerable importance to both manufacturer and consumer. Unlike steel and alloy wires, in which strength is governed by both chcniical and physical considerations, t,he high conductivity required for copper wire in the clcc-trical industry generally precludes the possibility of obt,aining tensile strength through the use of alloys. Consequently, this strenglh is largely a rnatter of the ~hysical structure, and all factors known to affect this structure, such as rolling temperature, die contours, lubricants, drawing speeds, etc., arc matters of prime importance. The experiments described here have becn carried out dl~ring the past year at the Baltimore plant of the American Smelting & Refining Co., in connection with a general investigat,ion of sorne of the questions mentioned above, with the object of studying thc physical ~t~ructure of hard drawn copper wire and the influence of ccrtain factors on t,hat st'ruct.ure. The work has been of an exploratory nature, rather than con~prehensivc., but it was thought that the result,s might prove of sr~fficient general interest tjo warrant their publicatlion. The bulk of the hard drawn coppcr wire used in places where strength is an important factor is rnanufacturcd to strict specification, based USU-ally in the United States on thab of t,he Arrlerican Society for Tcst'ing Materials (B-1-23). This specification calls for almost as high a, valuc of t,ensilc strength, consistent with elongation, as it is possible pract'icnlly to obtain, and only physically perfect mlaterial can be relied on to lliect such a specification. Obvious defects, such as laminat'ions protluc:etl through fins in bhe original rods or brittleness duo lo ovt~rdrawing, arc. usually recognized and need only passing refercncc. Laminatlions can be detected by twisting and untwisting the wircx several times; they will show up on the surface in the form of fins. Pig. 1 shows a very bad example of this defect. Thc wire may bc perfectly smooth before twisting hut, when such larrlinations arc3 prcstmt they will readily show up under this treatrncnt. This trouble is caliscd by sirriilar defects in the original rod produced by overfills in rolling. Fig. 2 shows the cross-section of the rod from which thc wire in Fig. 1 was drawn.

By C. E. Meissner

The need for a practical method of deciding upon new alloy steel analyses to widen its markets was the problem facing the Chrome Steel Works at the beginning of 1927. In addition to determining the physical properties of the alloys under consideration, it was necessary to test their value in service. Small melts of steel in a crucible, from which test specimens could be made, were insufficient to determine more than physical properties, therefore actual castings and ingots had to be made before the products could be placed in service. Furthermore, to save further experimental work, it was desirable that these heats should be made on a "pilot" basis; that is, actually duplicating regular open-hearth furnace practice. Therefore a one-ton (2000-lb.) acid open-hearth furnace was built. Furnace Design There are no special features in the design of this furnace. Standard acid open-hearth practice is followed, with the reproportioning of the dimensions necessitated by its small size; notably, a considerably increased length in proportion to width in order to get proper flame propagation. After the first heat was made it was found that the check-ers were insufficient. The checker chambers were increased in size and the checker openings were made much smaller. The furnace was started late in March, 1927, and made 65 heats. Experimental work was then stopped unti1, the data accumulated could be worked up and summarized. The run was resumed in November and the furnace is still operating, having made at this writing 60 additional heats. For the last series of runs, there was some rebuilding of both the furnace itself and the checkers. Slag pockets were added. Fig. 1 shows details of construction; Fig. 2, the general appearance. All designing was done by the operating and engineering staff. Construction materials were those available around the plant. The reversing valve is a standard butterfly type taken from an old furnace. The slag pockets and checkers are cleaned by digging down to arched openings which have been bricked up.

By N. T. Cotman, C. H. Pickering, P. B. Crawford

A beat-conduction model has been developed to study the flow of fluids in a stratified oil reservoir which is being subjected to unsteady-state depletion. To simulate stratification, plates of different metals were joined together so that cross flow could occur between strata. A description of the model, plate construction and necessary instruments is included. The data indicate that stratified oil reservoirs experiencing cross flow may be treated as a single, homogeneous, producing sand having properties intermediate between those of the layers making up the system. A technique for averaging core-analysis data to arrive at the proper average value of rock permeability and porosity is presented INTRODUCTION Although oil has been produced from reservoirs for 100 years and reservoir injection processes have been undertaken for 40 years, the industry is still inadequately informed on the problem of reservoir inhomogeneities. Stratification is perhaps the best-known type of reservoir heterogeneity. Urenl discussed the possible significance of this nonuniformity in his 1927 paper on the theoretical aspects of water flooding. Since this time the problem of vertical permeability variations has been the subject of much discussion and numerous papers. Levorsen2 simply describes the origin of layered beds in his Geology of Petroleum. Some beds are field-wide and may be traced from well to well, while others apparently pinch-out and become discontinuous. In addition, the physical characteristics of a layer may vary from place to place. For these reasons, a stratified system is not easily defined. In calculations involving oil reservoirs in which permeabi1ities vary, methods presently used contain many simplifying assumptions. The reservoir frequently is considered as consisting of permeable layers, each of uniform vertical and lateral permeability and each behaving as if separated from adjacent layers by thin impermeable layers. Other assumptions frequently made in reservoir calculations are steady-state flow conditions, unit mobility ratio, linear geometry, average or homogeneous porosity, and negligible capillary-pressure and gravity effects. Muskat3 has. developed solutions to the stratification problem for mobility ratios other than unity for both linear and exponential permeability distributions. He assumes lateral uniformity and continuity of all productive strata and does not consider porosity. This study will take into account transient effects, cross flow and varying porosity, as well as the variations in permeability, by means of a heat analogy. Muskat5 has shown the similarity of the equations which describe the flow of fluids in a porous medium to the flow of heat through a solid. Landrum, et al,6 made an analogy between the flow of fluid in an unsymmetrical reservoir and the conduction of heat in a similarly shaped metal plate. This analogy is summarized briefly as follows. If fluid is produced from an oil reservoir above the bubble point, the resulting pressure distribution will be identical to the temperature which would result in a metal or heat-conducting plate if heat were similarly removed. This means that a study of heat flow in stratified thermal models will behave like fluid flow in stratified oil reservoirs. For mathematical correspondencies, set- Ref. 6. The fluid mobility k/p is proportional to the thermal conductivity, and the porosity-times-compressibility product is proportional to the density times specific heat of the thermal plate. EXPERIMENTAL EQUIPMENT AND PROCEDURE The stratified reservoir model used in this study is illustrated schematically in Fig. 1. A Leeds and Northrup Speedomax Type G temperature recorder

By J. H. Henderson, L. A. Wilson, R. L. Gergins, R. J. Wygal, D. W. Reed

This paper presents a method of predicting the production history of an underground combustion recovery process. A rigorous solution of the thermodynamics and hydrodynamics involved is beyond the scope of this report. However, a practical scheme to appraise combustion recovery performance has been Worked out. By prudent assumptions, regarding the attainment of steady-state conditions, a trial-and-error solution, which satisfies both material balance and three-phase relative permeability requirements has been evolved and tested. As a combustion zone moves through a formation the interstitial water, the water of combustion and a portion of the oil wil1 be vaporized. These vapors will move downstream to a cooler region of the formation where they will condense. Part of the remaining oil will be displaced by the imposed gas drive and what remains behind will be consumed as fuel. The fluids will distribute themselves downstrean? in a manner which will satisfy the three-phase relative permeability characteristics of the formation. Three distinct zones will exist ahead of the combustion zone: (I) a zone termed a water bank which contains three mobile phases, oil, water and gas; (2) a region containing two mobile phases, oil and gas, called the oil bank; and (3) a zone having the fluid saturution distribution which existed prior to combustion. The mobile water will impose a water flood on the oil in the region containing three flowing phases. The process can be visualized as a simultaneous water .flood and combustion-supported gas drive. It is possible to estimate the saturation distributions in the water bank and oil bank as well as the production history for any combustion temperuture if the porosity, three-phase relative permeability characteristics, oil viscosity and initial saturations are known. This has been done for injection pressures ranging from I to 100 atm, oil viscosities of 10 to 1,000 cp and porosities from 20 to 40 per cent. Several of these calculations have been checked in the laboratory, with the result indicating rather good agreement between the predicted and observed production histories. The analytical and experimental techniques are described and the comparison of predicted and observed performance presented. INTRODUCTION Considerable research by the oil industry during the past 10 years has been directed toward the thermal recovery of crude oil. Of all schemes considered, the one receiving the most attention is in situ combustion. This is a process in which part of the oil is actually burned within the reservoir rock. The uniqueness of the process lies in its potential. Theoretically, this potential generates from the following: (1) it may open the door to vast reserves previously not economically producible, (2) it may induce increased production rates where low rates are now realized and (3) it may give high over-all recovery efficiency. The general idea of the process has received considerable attention in the literature '-' and the following brief description is intended for orientation only. In situ combustion can be visualized by considering a linear system, one end of which represents the injection well and the other the producing well. The formation at the inlet end is raised to ignition temperature by placing a gas or electric heater at the face of the formation or by the ignition of thermite, charcoal or similar material. While heating, the water and part of the oil in the heated region is vaporized and moved out into the reservoir by the injection of gas containing oxygen. In addition, part of the oil is displaced as liquid by the imposed gas drive. The oil not displaced is modified to an immobile, coke-like material (hereafter referred to as coke) by the high temperatures existing immediately ahead of the combustion zone. This coke is used as fuel, and the combustion front advances through the formation at a temperature which is probably in excess of 800°F. As the combustion front advancek it will drive ahead all interstitial oil not rendered immobile and all interstitial water plus the water formed by the combustion of hydrogen. As a result the sand behind the combustion zone is freed of all oil and water. Although the general description of the process has been outlined in the literature, there has been no consideration of the behavior of the fluids moving downstream from the combustion front. It is the purpose of this paper to present a technique

By R. E. Gilchrist, R. F. Nielsen

When a gas is displaced by mother in a porous medium, and there is a relative immobile liquid present. there is a transition zone in which the gas composition varies from essentially that of the original gas to that of the injected gas. The change in extent of this zone with time may be governed by a mixing process in the direction of flow and by the rate of establishment of gas-liquid equilibrium, if the gases are appreciably soluble in the liquid. The relative importance of these two factors in deterri~ining the growth of the transition zone is discussed, with appropriate equations and experimental data. INTRODUCTION The various solution methods of oil recovery involve a transfer of the components of a mixture, whether within a single phase or between a liquid and gas. The similarity of a natural porous medium to systems used in separation and chromatographic processes makes it seem reasonable that certain concepts used by chemical engineers in "mass transfer" might be applicable to petroleum reservoirs. One approach, the HETP (height equivalent to a theoretical plate), involves a distance or length of time over which equilibrium might be con. sidered as established. Another approach involves direct use of a rate of mixing or rate of establishment of equilibrium. The latter approach has been preferred in relation to production problems. DISCUSSION AND PROCEDURE When a gas or liquid in a porous medium is displaced by one of different composition but miscible with it there will be set up a transition zone in which the composition varies from cssentially that of the original fluid to essentially that of the injected fluid. If the fluid being displaced is a gas, and there is a relatively immobile liquid present, the transition zone may be

By W. E. Showalter

This paper discusses some of the results of combstion-drive tests which were made in a test cell using a sand bed 10 in. in diameter x 10-ft long. The test method is illustrated and described. The relationship between the API gravity of the in situ oil and the amount of air required for combustion drive is discussed in detail. Other things constant, the air requirement for combustion drive increases as the APl gravity of the in situ oil decreases. If the test results apply to actual reservoirs, the lowest-priced oils may cost the most to recover by this method. Information is shown which indicates that the effect of pressure on the amount of hydrocarbon burned is not large. A method of predicting air requirements from the API gravity of the in situ oil is presented. INTRODUCTION Combustion drive is the term used to identify the process of interstitial or in situ burning as an oil recovery method. Part of the in situ oil is burned to generate the energy needed to produce the remainder of the oil. Combustion drive as an oil recovery mechanism remains an economic uncertainty in spite of all the work that has been done by the industry in both laboratory and field?-6 This paper will show some of the results of tests which were made in a test cell for the purpose of studying the nature of the combustion-drive process. It will present data which indicate that the API gravity of the in situ oil is a significant indicator of the amount of air required to drive a burning front through oil sand. Air requirement varies inversely as the API gravity of the in situ oil. EXPERIMENTAL The tests were performed in a cell which utilized a cylindrical sand section 10 in. in diameter x 10-ft long. The thin-walled metal pipe which held the sand was wound with twenty external electrical resistance heaters which, by means of an automatic controller, maintained adjacent sections of the wall of the pipe at temperatures equal to the temperatures of the contained sand. Each heater covered 6 linear in. of the pipe. By this means lateral heat loss from the sand section was minimized, thereby causing the sand section to simulate more closely a horizontal increment of a combustion-drive reservoir. Fig. 1 shows a schematic diagram of the test assembly. Thermocouples to measure the temperature in the sand were located every 6 in. along the length of the sand section. The pipe containing the sand was enclosed in a cell designed for an operating pressure of 500 psig. The inlet air pressure was controlled at the inlet, and the gas flow rate was controlled and measured at the outlet of the cell. The oil sand used for the tests was prepared by mixing first water and then oil with the noncon-solidated sand using a closed mixer similar to a cement mixer. Table 1 shows a screen analysis of the sand. Ninety percent of the sand was 100 mesh or finer. This sand was a mixture of 80 per cent No. 120 Nevada White Sand and 20 per cent Tennessee Hi-Fusion Moulding Sand No. 3. The Nevada sand was a clean silica sand. The Tennessee sand con-

By William Hurst

Muskat's depletion performance equation is here derived considering the expansion behavior of the reservoir hydrocarbon system and a simple fractional-flow equation. This nietkod of derivation leads logically to the two extensions that follow. The first of these is concerned with gravity segregation in a depletion-drive reservoir. The second is concerned with including an empirically determined tern? for water influx in the performance equation. The more general equation for gravity segregation when there is a primary gas cap and empirically-determined water influx is stated for completeness. These equations have been found useful in reservoir performance calculations in Eastern Venezuela. A discussion on the methods of solving these equations follows, and considers firstly the effect of taking finite intervals in the numerical integration, and sec-ondly, methods of incorporating the time functions involved in segregation in with the expansion behavior. The paper concludes with a brief general discussion on further extensions to the depletion performance equation. INTRODUCTION The two fundamental sources of energy by which oil is produced from a reservoir result from pressure depletion inside the boundaries of the reservoir and fluid encroachment across the boundaries of the reservoir. The wells in either case form low pressure outlets through which oil and gas may be produced by the expansive force of the reservoir fluids and associated encroaching fluids. When the reservoir pressure is higher than the bubble-point pressure of the oil, so that there is no free gas in the reservoir, these expansive forces are the only ones available for the production of oil. However, when the reservoir pressure is less than the bubble-point pressure of the oil, free gas is vaporized as the pressure falls. With both oil and free-gas phases present, the additional forces of gravity and capillarity may operate on the gas-oil system, as they have previously operated on the oil-interstitial water system. Gravity tends to segregate the free gas from the oil due to their density difference. Capillarity opposes and eventually balances gravity as the more extreme free gas and oil saturations are reached, preventing the independent move- ment of free gas until it is above a certain saturation, and the independent movement of oil when it is below a certain saturation. The type of depletion performance equation chosen for predicting the future performance of a reservoir depends on the amount of past history available. When the reservoir is somewhere past the halfway mark in depletion, some form of decline curve is often used. With less past history, material balance equations which incorporate empirical factors based on the past performance are often used. When, however, the amount of past history is small, the Muskat depletion performance equation will usually be used. The distinguishing feature of this type of equation is that empirical factors based on the over-all or macroscopic reservoir behavior are almost or entirely absent. Each parameter affecting the reservoir performance is ascribed an independent set of values based on measurements made on laboratory samples; that is, incorporating microscopic empirical factors. In establishing Muskat-type depletion performance equations, it is necessary to consider the reservoir as consisting of a number of associated blocks, in each of which the saturations and pressures may be considered uniform, and in each of which all substances have uniform pressure-volume characteristics. Thus, a primary gas cap can usually be considered as one block and an aquifer as another. Gravity segregation may be negligible for practical purposes when the rock and oil properties are adverse and/or the dip or thickness of the reservoir is too small. In this case the whole oil leg may be considered as one block, except in very large reservoirs. In very large reservoirs the fluid and rock properties may vary enough, particularly in the dip direction, for it to be necessary to divide the oil leg into a number of blocks, in each of which the relevant quantities may be considered uniform. When gravity segregation of the oil and free gas is not negligible, it is necessary to consider the space occupied by the initial oil leg as divided into two blocks, a secondary gas cap and an effective oil leg, in each of which saturations may be considered to be uniform. The total volume of these two blocks is thus constant, but the secondary gas cap grows continuously at the expense of the oil leg. Muskat1 derived depletion performance equations for the basic case of an oil reservoir with closed boundaries and without segregation, and for the case of an oil

By W. B. Gogarty

With flow of non-Newtonian fluids in porous media, effective viscosities are needed for use in the Darcy equation. These viscosities depend on the rock parameters and flowing conditions. In this investigation, rheological and flow results were correlated to develop an expression for calculating effective viscosities. Surfactant-stabilized dispersions of water in hydrocarbon were used with consolidated sandstone cores of different permeabilities. The average shear rate in a core was related to the permeability and porosity of the core and the frontal velocity through the core. By using the correlation obtained for determining effective viscosity, experimental and theoretical values were compared for each fluid system by determining average and maximum errors. For the fluid system showing the largest average error, the correlation fit 38 experimental points with an average and maximum percentage error of 4.3 and 14.8, respectively. For the fluid system showing the least error, average and maximum percentage errors of 1.8 and 4.3, respectively, were calculated using 29 experimental points. INTRODUCTION Use of non-Newtonian fluids in the petroleum industry is not new. Fluids of this type have been used for many years as fracturing agents and drilling muds. Recently, attention has focused on the use of high molecular weight polymer solutions for secondary recovery. Results have been reported on studies relating to the flow behavior of polymer solutions in porous media.1-4 Each study gives the viscosity in the Darcy equation as a function of the rheological properties of the fluid, the characteristics of the porous medium and the pressure gradient. The functional forms of the reported expressions differ somewhat in each study. The first three studies used unconsolidated porous media while the fourth used a single consolidated core. Because of the limited information available on the flow of non-Newtonian fluids in porous media, further investigation was undertaken in this area. We selected non-Newtonian fluids that did not contain high molecular weight polymers in order to minimize core plugging. The over-all objective was to gain a better understanding of the flow behavior of non-Newtonian fluids in consolidated porous media. Specifically, our goal was to determine the relationship between the average shear rate in a core and the frontal velocity. Because of this goal, it was necessary to determine the effect of permeability on the shear rate-frontal velocity relationship. We wanted to replace the viscosity in the Darcy equation with an expression relating to the average shear rate in the core and the rheological properties of the non-Newtonian fluid. THEORY For a time-independent, non-Newtonian fluid, the apparent viscosity is a function of the shear rate. This relationship is given by µ = F(v')...........(1) With Eq. 1, a single value of shear rate determines the apparent viscosity at that shear rate. With flow in porous media, a spectrum of shear rates occurs at any average throughput. An average shear rate v' can be used to represent the spectrum of shear rates that exists within the core. By using v' and the rheological data given by Eq. 1, an effective viscosity is determined for the flowing fluid. If v' is the correct average, µe becomes the appropriate viscosity for use in the Darcy equation. Under these conditions, flow rate-pressure drop relationships can be determined for non-Newtonian fluids. HOW the correct average shear rate might be related to the frontal velocity through the core and the characteristics of the porous media will now be considered. Later, the experimental data are presented that indicate the validity of these

By Henry J. Miller, Marcella Linderman

During some experiments made by Henry J. Miller, partly in Germany and partly in the United States, in which it was found necessary to melt metals in quantities up to 60 kg. in a vacuum or under low pressures of hydrogen, molybdenum was tried as the resistor in the electric melting furnaces. Molybdenum ribbons 16 by 0.55 mm. and 18 by 0.45 mm. were used as resistors. The furnaces were of the type shown in Fig. 1. About the heating tube, Y, are wound parallel resistors of molybdenum ribbon from the water-cooled electrodes, El, Ez. The whole is within a copper container having outlets for the removal or introduction of gases and a window, W, for observation, the space about the heating tube and electrodes being filled with insulating material. The life of the molybdenum used as the resistor was found to depend very much on the manner of using it, its failure being due in some cases to thinning by evaporation or chemical reaction and in others to hardening, which prevented it from following the expansions and shrinking of the heating tubing without cracking. As such failures are often rather expensive in commercial work, it was thought worth while to determine the causes. The factors influencing the behavior of molybdenum resistors used in this manner may be briefly stated as follows:

By R. F. Krueger

Results are presented from tests of dynamic fluid-loss rates to cores from clay-gel water-base drilling fluids containing different commercial fluid-loss control agents (CMC, polyacrylate or smt,ch), organic viscosity reducers (quebracho and complex metal lignosulfonate) and oil at several different levels of concentration. In the dynamic system the most effective individual additives to the clay-gel drilling fluid, based on cost-equalized concentrutiom, were found to be starch and the viscosity reducers. These results do not conform with the rankings determined by API fluid-loss rests, which indicate CMC, polyacrylate and starch to be the most effective and comparable. Generally, minimum dynamic fluid-losr rates were attained at cost-equalized concentrations of additive (including thinner) of about $1.00/ bbl, or less. For chernically treated clay-gel drilling fluids, both the standard and the high-pressure API filter-loss tests were found to he inaccurate indicators of trends in dynamic fluid-loss rates under the test conditions used, particulurly for drilling muds containing viscosity reducers. From a practical field viewpoint, restrictions on the applicability of the API fluid-loss test are such that it is open to question whether or not results of this test can be used routinely with confidence as an indicator of control of down-hole fluid loss under field treating conditions. INTRODUCTION The petroleum industry spends large sums of money during drilling operations to control the fluid-loss properties of drilling fluids based on the standard API filter-loss test,' which is a static filtration system. Laboratory studies' ' of dynamic filtration have shown that in a given time period filtrate loss from a circulating mud stream is greater than from a static system and that it is a function of linear mud velocity, pressure and the properties of the drilling fluid. Ferguson and Klotz' and Horner, et al," observed that (I) the dynamic fluid-loss rates for the drilling fluids used were not related to the extrapolated API filter loss and (2) the drilling fluids with the lowest API filter losses did not have the lowest dynamic fluid-loss rates. However, there has been no published information on the relative effects on dynamic fluid-loss rate as a given drilling fluid is treated with increasing amounts of chemical additive to reduce the API filter loss. Such information is economically important because drilling-fluid costs rise rapidly as chemical requirements increase. This paper presents the results of a study of dynamic filtratioi rates to cores from a clay-gel water-base drilling fluid treated with various commercial viscosity reducers and chemical fluid-loss control agents. The dynamic fluid-. loss rates to cores are compared with the standard API filter-loss values at several different levels of additive concentration. Dynamic filtration rates were obtained in each experiment under two different simulated wellbore conditions: (1) filtration just above the bit through a new mud cake laid down dynamically on a freshly drilled formation and (2) filtration up-hole through a mud cake formed by deposition of a static filter cake on top of the initial dynamically formed cake. The latter case corresponds to the bottom-hole conditions existing above the bit when mud circulation is restarted after a stand of pipe has been added or a round trip has been made to change the bit. Except for the short-duration, high-rate filtration beneath the bit where no mud cake can form, these two conditions probably represent the two extremes of dynamic filtration. Because thickness of a dynamic mud cake formed on freshly exposed formation is limited by the shearing action of the mud stream, the filtration rate for this condition is high. On the other hand, once circulation is stopped and a static mud cake forms on top of the dynamic cake, re-starting circulation has only a small effect on the cake properties and filtration rate is much lower thereafter. A discussion of the mechanics of mud-cake deposition and dynamic filtration is outside the scope of this paper but may be found in more detail in publications by prior investigators. APPARATUS AND EXPERIMENTAL CONDITIONS The test equipment used to simulate the dynamic flow conditions existing during drilling was a modification of that described previously by Krueger and Vogel: A schematic flow diagram is shown in Fig. 1. In general, a power-driven, high-pressure mud pump capable of delivering up to 60 gallmin was used to circulate drilling fluid parallel to the faces of 1-in. diameter sandstone cores mounted in a 2 3/4-in. ID high-pressure test cell. Pump rates were controlled by means of a magnetic clutch to maintain an average axial fluid velocity of 110 ft/min in the annular space between the cell wall and a 1 1/2-in. rod positioned on the center line of the cell. The core specimens were Berea sandstone plugs sealed with plastic inside 1 1/8-in. OD tubes and were fluid-saturated prior to use. Burettes were used to accumulate fluid discharged from the cores. The mud sump shown was used for treatment and storage of the drilling-fluid samples during a particular test. The valve arrangement permitted either (1) circulating drilling fluid through the by-pass line while treating with

By Francis R. Conley, John A. Sievert, John N. Dew

A brief review is presented of the past performance of a number of large, thin, highly permeable reservoirs with low dips in the Bolivar Coastal fields of Venezuela. The performance of these reservoirs indicates that the fluids are segregated vertically within the sand section by gravity. With this assumption, equations are developed which describe the performance under pressure maintenance operations. Methods of solving these equations and results of simple example calculations are presented. Example calculations indicate that under pressure maintenance conditions injected gas tends to {tow preferentially along the top of the sands and that encroaching water has a tendency to flow preferentially along the bottom. The expected performance of segregated fluids is discussed and compared with that of fluids which are uniformly distributed in any sand section. INTRODUCTION The field performance of a number of reservoirs described herein indicates that the fluids are segregated vertically within the sand section by the force of gravity. Thus, it was felt that any method used to predict the future performance of these reservoirs should consider the effects of segregation in the sand sections. A study of the available information on some of the reservoirs suggests that increased recoveries may be expected if the reservoir pressure is maintained by either crestal gas injection or flank water injection. To predict future performance of reservoirs under pressure maintenance operations, a method of analysis was needed which would account for segregation of fluids in the sand sections. Several methods of analysis have been developed to take into account the segregation of fluids in the reservoir as a whole1,3 To our knowledge no method considers segregation of the fluids within sand sections in the manner indicated by the past performance of several reservoirs in the Mara-caibo Basin. This paper outlines part of the work done in studying some of these reservoirs and contains a description of their performance characteristics. The analysis presented is restricted to pressure maintenance conditions, since space limitations prevent a full discus- sion of the development of the mathematical relations and the various methods of solving them. Only a sketch of the development of the mathematical relations is given. The various methods of solving these relations are pointed out, but the actual determination of solutions to various problems are omitted except for one example. It is felt that these results may aid in the study of reservoirs outside the Maracaibo Basin. Some concepts on which the present analysis is based are outlined by D. N. Dietz. However, the analysis presented herein includes several factors not considered by Dietz: (1) variations in permeability and in the cross section, (2) various shapes of the cross section, and (3) the production of fluids. Also, the mathematical development presented by Dietz differs considerably from our corresponding analysis. RESERVOIR CHARACTERISTICS The reservoirs under consideration, large and thin, with low dip and high permeability, are large to the extent that they contain from 0.5 to over four billion bbl of oil initially in place. The section thicknesses vary from 100 to 400 ft and the lengths from 10,000 to 40,000 ft. Thus, a length-wise cross section of the producing formation appears to be long and thin. The dip angles vary from 0 to 6 degrees and the average permeabilities vary from 0.5 to over three darcies. The reservoirs contain from 20 to 50 per cent shale inter-laminated with the productive sands. Most shale breaks within the major sands do not correlate from one well to the next, and correlation is often difficult, even between wells drilled from the same location. However, it is often possible to correlate the shale breaks between major sands for some distance. Herein we are concerned with the performance of reservoirs containing oil of gravities over 20" API. Past reservoir performance indicates good pressure communication and, in cases where pressure sinks have developed, large amounts of fluid migration have occurred. Few of the reservoirs have had initial gas caps, but also, few have been found to be highly under-saturated. Due to gravity segregation, secondary gas caps usually form before the pressure has been reduced more than 25 per cent of its initial value. The effect of gravity has not only been apparent on a reservoir-wide basis, but has also caused segregation of the reservoir fluids in productive sections. Proof of this is found from the results of selective well tests, workovers, and electric and various other types of logs. Evidence that

By Wilbur H. Somerton

W. H. Somerton is to be commended for his application of the methods of dimensional analysis in combination with the results of laboratory drilling tests to obtain a significant formula for the rate of penetration for the process of rotary drilling. This formula is shown to be a consequence of a more general formulation developed by the discussor as a result of his investigations of the basic action of most methods of rock drilling by mechanical means—the momentary loading of a chisel-shaped tool on the surface of the rock at the bottom of the hole. Some of the author's conclusions are re-examined in the light of this general formulation. DRILLING RATE FORMULA The author selects a minimum number of the most significant variables associated with the rotary drilling process to set up his two dimensionless parameters. Only one variable is associated with the geometry of the bit —the bit diameter, D; only one variable is associated with the properties of the rock pertinent to drilling— the drilling strength, S; and only two independent variables are associated with the conditions of operation the weight on the bit. F, and the rate of rotation, N. The dependent variahle of interest is, of course, the rate of penetration, R. The results of the author's drilling tests show that the variable R that appears in the dimensionless parameter, RIND is approximately proportioned to the square of the variable, F, that appears in the other dimension-less parameter, F/D2S. This establishes the functional relationship between these two simple parameters that leads to the formulation of a third dimensionless parameter, RD1T2/NF1. This last parameter is a significant one for the rotary drilling process, in general, although the author properly limits the validity of his equations to his present test conditions. The author realizes that the numerical value of this parameter, C, is not a constant but may be a function of one or more additional dimensionless groupings of the many other variables associated with the rotary drilling process, e.g. some dimensionless measure of bit wear. The author's parameter is shown to be a special case of a more general dimensionless parameter, developed by the discussor. that is applicable to practically all methods of rock drilling by mechanical means, namely. RD2S/P. P is the rate at which mechanical work is done on the rock, and S is the drilling strength, a term originated by the discussor. S is defined as the ratio of energy input to the rock to the volume of rock broken

By F. D. Patchen

Changes in the properties of partland cement upon the addition of fine-ground silica are discussed. Data were collected from formulations cured for periods up to 60 days at temperatures varying from 180" to 350°F. Compressive strengths, permeabilities and thickening times are compared. It is demonstrated that addition of find-ground silica in proper proportion results in a hydrated cement with superior properties above 230°F; i.e., early compressive strengths are higher and increase with additional curing, permeabilities are significantly lnrver and thickening times are increased. X-ray analysis of hydrated cement samples have been used to determine more fully the complex chemistry of these systems. Dicalcium silicate alpha-hydrate was found in all neat cements cured above 260°F. A 10 weight per cent silica addition produced a maximum amount of this phase. Xonotlite and tobermorite have been identified as hydration products at 320°F in systems containing 35 or more weight per cent silica. Formation of these components appears to be primarily responsible for the gain in strength and reduction of permeability observed after the addition of silica. INTRODUCTION For several years it has been recognized that loss of compressive strength occurs when portland cement is cured at a temperature above 250°F. This phenomenon has been called "strength retrogression". More recently it has been recognized that increased permeability usually is associated with decreased compressive strength at these temperatures. At curing temperatures above 300°F, these factors reach serious proportions. A compressive strength loss of 50 per cent between curing times of 24 hours and 14 days has been reported,' and permeabilities as high as 10 md have been measured in neat cements cured at 320°F for seven days.' Earlier investigators,' although unsuccessful in eliminating strength retrogression, stated ". . . .it appears possible that by the proper selection of additives for cement, strength loss with age at high temperatures can be either eliminated or considerably lessened". One solution suggested was the use of a pozzolanic material, hydrated lime, and a retarder. Menze15 in 1934 discovered that portland cement containing finely divided silica possessed high compressive strengths when the curing temperature was 350°F. The use of such a composition is now an established practice in the manufacture of steam-cured cement blocks. It is reasonable to assume that the same beneficial effect can be used advantageously in high-temperature oilwell cementing operations. Ludwig first mentioned the use of sand for this purpose. Later, a portland cement containing moderate additions of silica flour and retarded with CMHEC, a cellulose derivative. was recommended for high-temperature application." These formulations have been successfully employed in the field."'" This paper presents the results of a research program to test the efficacy of finely ground silica in increasing the compressive strength and reducing the permeability of oilwell cement compositions containing portland cement. Included also are X-ray diffraction patterns of powdered samples which have been used to determine some of the important chemical reactions occurring in these systems. LABORATORY PROCEDURES AND MATERIALS All compressive-strength test specimens were prepared and cured from a slurry in accordance with API RP 10B, with two modifications. A final curing pressure of 5,000 psi was maintained for a one- or three-day period; but, for further curing, pressures were reduced to the vapor pressure of water at the curing temperature. Also, the water content was varied in samples containing silica to maintain the slurry viscosity within a range of 5 to 15 poise as measured on a consistometer. Test specimens at 180°T were cured in a water bath at atmospheric pressure. At higher temperatures the specimens were cured in a steam-jacketed high-pressure curing vessel at 5,000 psi for a one- or three-day period. cubes cured for longer periods were removed from the pressure-curing vessel after 24 hours, sealed in small bombs containing water and placed in an aging oven.

By M. J. Buerger

Possibly the most puzzling features observed during a single-crystal deformation test are the appearance of slip striae on the surface of the crystal and the strengthening of the specimen.' While the first phenomenon is mainly of theoretical interest, the latter is fundamental to all strain-hardening and cold-working processes in which the deformation occurs by translation. The plastic deformation of crystals has received considerable discussion2 and need not be reviewed here in any great detail. It has been observed, essentially, that after a certain amount of. initial stress is applied, the surface of the crystal becomes covered with sets of parallel striae or rings which appear to divide it into thin slabs or slip blocks.3 The striae are the traces of rational crystallographic planes along which shearing movements are localized, these movements providing the strain which accompanies stress applied beyond pure elasticity. As deformation under the applied stress proceeds the slip blocks shear past one another, although the essentially single-crystal character of the specimen is unimpaired. All of the shearing apparently takes place in the immediate vicinity of the original slip striae without the appearance of new planes of movement.4 Since each increment of strain is attended by an increased stress requirement to continue the deformation, the crystal strengthens during the process.

By C. W. Arnold, H. L. Stone, D. L. Luffel

The purpose of this research is to determine (I) the efficiency of small banks of enriched gar driven by methane in displacing oil from a porous medium and (2) the effects of variation in bank size and composition of that efficiency. Most of the experiments were conducted in a sand-packed tube 20-ft long and 1/2-in. in diameter. The hydrocarbon system generally used was methane, butane and decane at 2,500 psia and 160°F. The results of these experiments indicate that, in the regions contacted by the gas, a small bank of an oil-miscible gas driven by methane can displace all of the oil in a piston-like manner. If the enriched gas is of such composition as to remain immiscible with the oil, displacement of oil is less efficient than for the miscible case, and the gas bank travels through the sand with a velocity less than that of the driving gas. These data along with theories discussed imply that smaller banks and less total gas are required when the enriched gas and oil are miscible. INTRODUCTION Widespread application of enriched-gas drive to the recovery of oil rests upon a key factor — the use of limited quantities, or "banks", of enriched gas. At the present time, the value of liquefied petroleum gas or other enriching agents discourages their use in a continuous injection technique, or even in a large bank, except in a few isolated reservoirs. If small banks of enriched gas driven by methane were as effective in displacing oil as is continuous injection, the enriched-gas drive process might be applied to a larger number of reservoirs. Previous research on the mechanics of the enriched-gas drive process reported by Stone and Crurnpl and by Kehn, Pyndus and Gaskell has utilized continuous injection of enriched gas. This work has shown that two types of displacements occur. With gases containing sufficient intermediates. the oil is displaced misciblv and complete recovery is obtained from the regions swept. When gases are used which contain insufficient intermediate hydrocarbon for miscible displacement, oil is displaced immiscibly. In the latter type, selective solution of the intermediate hydrocarbons causes a swelling and reduction in viscosity of the oil and leads to an increased recovery over that obtained by dry-gas (methane) drive. The size of the enriched-gas bank necessary for efficient displacement of oil is determined by those factors which cause deterioration of the bank. A differentiation may be made between those factors which operate on a microscopic scale and those which act on a macroscopic scale. On the smaller scale, the enriched gas mixes in the direction of flow by diffusion and convection with the fluids immediately preceding and following it. On the larger scale, the gas may by-pass the oil by flowing through permeable streaks, by overriding the oil because of density difference, or by fingering because of unfavorable viscosity ratios. In such cases, the enriching material tends to mix with the oil both laterally and in the direction of flow. The increase in effective area available for diffusion and dispersion of the enriching components leads to a faster degradation of the bank and a need for a larger bank than is necessary for those cases in which no by-passing occurs. The effects of such macroscopic factors in the deterioration of enriched-gas banks have been reported in a separate paper by Blackwell, Terry and Rayne. The present study was confined to the factors which operate on the smaller scale, in particular to the behavior of banks of enriched gas in sands uniformly swept by the gas. Experiments were designed to answer the following questions. 1. Can small banks of enriched gas driven by methane be used to secure oil recoveries comparable to those obtained by continuous injection of enriched gas? 2. What is the optimum bank size (the minimum bank size necessary to obtain a recovery comparable to that obtained by continuous injection of the same enriched gas)? 3. How many total pore volumes of gas must be injected to obtain the maximum recovery when the optimum bank size is used? 4. What is the effect of varying the number of enriching components in a gas bank? This report describes the experimental investigations and discusses the results in terms of their significance to reservoir behavior.

By B. L. Landrum, P. B. Crawford

Theoretical and potentiometric model studies have been made of the effect of non-uniform lateral permeabilities on pattern sweep efficiency and production capacity in waterflood and gas-cycling programs. It is shown that a diflerence in directional permeability by a factor of three may result in a sweep efficiency of only 43 per cent for a five-spot pattern or a sweep of either 79 or 38 per cent for a direct line-drive square pattern, depending on the direction of the line-drive flood. Changes in the pattern conductivity varied from about 0.8 to 1.34 over this same permeability variation, depending on the pattern used. It is suggested that measurements be made to determine the possible magnitude and extent of the directional permeability phenomenon early in the field developinent and certainly prior to the initiation of any fluid-injection progmm. INTRODUCTION Irregularities in reservoir sand properties long have been a major difficulty to anyone attempting to explicitly describe the field characteristics of oil production. In particular, it is well known that vertical and lateral permeabilities often differ appreciably; however, the existence of large regions with lateral permeabflity variation is not widely recognized. A number of years ago, extensive studies were conducted by the Secondary Recovery Research Labcratory of the Pennsylvania Grade Crude Oil Association, primarily on the Bradford field. Johnson and Hughes' reported a permeability trend in the northeast-southwest direction. They indicated that flow in the preferred direction may be 25 to 30 per cent greater in that direction than in the northwest-southeast direction. They also reported that similar effects may be found in other nearby fields. The origin of the permeability variation has been discussed by Griffith Hutchinson' described the results of labora- tory tests on limestone cores, pointing out that per-ferential directional permeabilities were significant in one-half of 10 formations studied and that the average permeability ratio was 16:1. Barfield, Jordan and Moore analyzed the fluid-injection displacement performance in the Spraberry and reported that the effective ratio of permeability parallel to the major axis to the permeability perpendicular to this axis was 144:l and that it inclined along an angle of about 30". To prevent the reduction of sweep efficiency which would normally result from flooding or cycling in such reservoirs, well spacing could be varied from the customary square to a diamond-shaped pattern. Johnson and Hughes' presented a method for determining the desired shape. Another method will be discussed in detail in the next section. It is the purpose of this paper to investigate quantitatively the effect of directional permeability upon the sweep efficiency and production capacity of the five-spot and line-drive well patterns. Information currently in the literature will be applied to this particular problem and augmented by a potentiometric model study. It is stressed that the effects herein discussed may or may not be of great importance because field data on this phenomenon are far too scant to form such an opinion. It is suggested that, if conditions similar to those found by Johnson and Hughes frequently exist (as may well be the case), more consideration cdrtainly should be given to this problem.

By J. F. Gormley

Drill pipe, and the connections used to join the pipe together, have had a long history of development and improvement. With the growing use of high tensile-strength pipe, the rotary shouldered connections used to join the lengths of drill pipe have become more critical. Fatigue-type failures of rotary-shouldered connections, or tool joints, have become a cause of much conern to operators and manufacturers. Most fatigue tests on tool joints have been made on the basis of bending moment applied vs number of cycles, rather than on the basis of stress vs cycles. A desire to learn more about the nature of stress in a tool joint due to bending caused a rather intensive experimental study to determine what actually happens to cause fatigue failure. It was known that a tool joint with insufficient make-up torque was subject to early fatigue failure. Therefore, the establishment of the relationship between make-up torque and allowable bending moment would enable a testing engineer to better correlate the results of various tests. The problem of calculating the stress in a tool joint due to an applied bending moment was complicated by: the effect of unknown stress concentrations; the applicability of the well known formula S = Mc/I; and the question of whether to use the combined cross section of the pin and box in calculating the moment of inertia I. TEST EQUIPMENT The test equipment consisted of the following. 1. A lathe modified to act as a fatigue machine for full-size tool joint or drill-collar specimens. 2. An adapter plate bolted to the headstock of the lathe and into which a short adapter was screwed. 3. Tool-joint pins on sections of drill collars bucked into the adapter. 4. A hydraulic load cylinder which rested on the bed of the lathe and applied bending moments to the joints. 5. Paper-backed strain gauges (FAP-12-12, 1/8-in. gauge length and A-5-1, 1/2-in. gauge length). 6. Eastman 910 contact cement for bonding gauges to joint and Epoxylite 222 waterproofing. 7. A 12-position switchbox, SR-4 bridge balancing unit connected to an SR-4 strain indicator to read the gauge output. 8. A four-channel Sanborn recorder, Model 150, was used as the readout instrument for part of the testing. The fatigue machine was well suited as a loading device. Because of the long lever arm available, the change of moment with distance was small; consequently, the vertical shear stress was small, less than 300 psi. Some difficulty was experienced with the hydraulic load cylinder at high bending moments due to leakage. The strain gauge equipment was simple and quite reliable. A 41/2-in. extra hole (E. H.) joint with "as machined" threads and shoulders and a 2 7/8-in. internal flush (I. F.) joint with copper-plated threads and shoulders were selected for testing because of their common usage in the field. In addition, the 4 1/2-in. E. H. was suitable because it had an inside diameter large enough to permit the insertion and accurate positioning of strain gauges without using special jigs, while the 2 7/8-in. I. F. joint required much lower bending loads to cause separation of the shoulders of the joint when made up to the recommended torque. The joints were made to standard API specifications, except for very small bevels on the box and pin shoulders and cylindrical stress-relief grooves added to the pins. Larger bevels would have

By R. C. Craze

This paper discusses the characteristics of the velocity logs now available to the petroleum industry, and some of their advantages and limitations. The velocity log was designed as an aid to the geophysicist to determine velocity layering in the earth and to provide fast, accurate and economical vertical travel time data. After a number of logs had been made, studies were carried out on their use ia correlating strata in adjacent wells. It was found that adjacent strata usually show sufficient velocity contrast to make it possible to correlate these from well to well, even in areas where electric or radioactive logs fail.' In addition, velocity logs have been useful in the determination of porosity and thickness of reservoir rock.' These uses for velocity logs require an understanding of the characteristics of the log in response to the changes in velocity and in borehole conditions encountered normally in logging. GENERAL DISCUSSION OF INSTRUMENTATION The measurement made in velocity logging is more accurately described as a measurement of time rather than velocity. Ideally the log is obtained by measuring the length of time required for compressional sound waves to travel through a fixed distance in the formation. In practice, the time measurement may be made accurately but the distance through which the waves have traveled may vary slightly. Changes in borehole diameter are responsible largely for the difficulty in maintaining a fixed distance, especially if the diameter variations are so abrupt and frequent that the borehole instrument cannot follow contours of the hole. Several service companies have velocity logging equipment available or under development which have different transmitters, receivers, and measuring systems. While these differences are important, there are actually only two systems in use at the present time (Fig. I)."' The first system initiates a pulse at a trans-mitter and detects it at a single receiver. The transmitter and receiver are separated by a fixed distance, called the spacing. The time measured is that between the initiation of the pulse and the first arrival of acoustical energy at the receiver. The second system uses two receivers and a transmitter. The receivers are separated by a fixed distance and the time between arrivals of energy at the two receivers is measured. Here, the spacing is the distance between receivers. The two systems have much in common for it is possible to record both logs simultaneously. Usually the transmitter is an electro-mechanical transducer which produces discrete

By J. G. Savins

Analytical procedures are described for optimizing the selection of a rheological model when it is desired to express the functional relationship between true shearrate and shearing stress in analytical form. The procedures are extended to two principal model categories — "Generalized Newtonian" and "Viscoplastic" models. Graphical examples are presented to illustrate the utility of certain characteristic deriuative functions in distinguishing between these categories and in determining how well a particular model reflects the flow behavior of rheologically anomalous materials in simple shearing flows. INTRODUCTION Modern drilling fluids and fracturing fluids, being chemically complex, are frequently non-Newtonian — a peculiarity which complicates the solution of problems such as designing flowline and process equipment, and predicting laminar-turbulent behavior, solid-suspending ability and cutting-transport efficiency. Therefore, techniques of data analysis which accurately characterize the behavior of rheologically complex materials are needed for the accurate formulation of solutions of these and similar problems. It is the purpose of this paper to describe the application of a sensitive method of rheological data analysis — the dual differentiation-integration method of optimizing the selection of the functional relation between shear rate and shearing stress — to a variety of materials of interest in drilling and production practices. BASIC CONCEPTS THE DIFFERENTIATION METHOD To a first approximation, the viscosity function for a non-Newtonian system can be expressed in a form analogous to Newton's law of viscosity by the following generalization2.3 between the stress On introducing the restrictions that the system is incompressible and the motion is rectilinear, then l = 0 and l3 = 0 ; and with the further restriction that the flow is slip-free, the viscosity (an invariant) is then determined solely by 12. With this approach one can set up a suite of integral equations relating observed kinematical and dynamical parameters for certain types of flow experiments involving simple shearing motion, e.g., flow in a cylindrical tube, flow between rotating concentric cylinders, and flow between fixed parallel surfaces of large aspect ratio. Taken in the order cited, these integral equations are as follows.