Search Documents

By R. H. Jacoby, M. J. Rzasa

Experimental equilibrium vaporization ratios were obtained for nitrogen, methane, carbon dioxide, ethane, and hydrogen sulfide in natural gas-condensate systems. Two different overall mixture compositions were investigated over the pressure range 500 to 4,000 psia at 100°, 150°, and 200 °F. Partial phase diagrams of the mixtures were obtained; one of these showed a critical point at 204°F and 3,875 psia. These and earlier data' were compared with the Kellogg correlation' at! 100°F and 1,000 psia. INTRODUCTION The authors' have presented equilibrium vaporization ratios for nitrogen, methane, carbon dioxide, ethane, and hydrogen sulfide in complex mixtures, such as natural gas-crude oil and natural gas-absorber oil. Since mixtures of the natural gas-condensate type also occur widely in petroleum reservoirs and often contain the non-hydrocarbons mentioned, it would be useful to have appropriate K data for calculating equilibrium vaporization conditions. Roland, Smith, and Kave-ler' obtained K's for methane and the light hydrocarbons in gas-condensate mixtures; and Poettmann and Katz3 obtained K's for carbon dioxide in a gas-condensate system. The data presented in this paper extend this coverage to include K's for nitrogen and hydrogen sulfide and phase diagrams for the experimental mixtures. EXPERIMENTAL WORK The equipment and experimental procedure was essentially as described previously, except for analysis of the vapor samples. The vapor sample bombs were filled at the pressure and temperature of the run. Immediately after sampling, a portion of this sample was flushed from the bomb into a 650 cc glass bulb, filling it to near atmospheric pressure. The sample in the glass bulb was then analyzed by the mass spectrometer. Some checks were made on this analytical procedure to find the accuracy of measurement of the C + content of the vapor samples. The modifications tried were complete fractionation of the vapor sample and the use of several liquids (depropanized condensate and retrograde liquid from several vapor samples) as calibrating materials for the mass spectrometer. The results obtained from these various checks differed by about 10 per cent relative and this is believed to be the order of accuracy of the vapor C+ concentrations.

Jan 1, 1953

By P. A. Witherspoon, S. P. Neuman

Evaluating the permeability of a caprock overlying a potential gas storage reservoir is a very critical problem. Pumping water from the reservoir can be used as an evaluation tool in analyzing this problem. Fluid level changes that occur in the aquifer us well as in the caprock can be measured with appropriately placed wells. If the leakage of water from the caprock into the aquifer is considerable, the effects will be apparent in the aquifer. If the leakage is slight, however, it will not be possible to detect it with certainty from observations in the aquifer alone. Fluid level measurements in the caprock must be relied upon. and improved methods of analyzing such effects have been developed which are based on a theoretical analysis of fluid flow through a caprock of infinite thickness. An example applying these methods to field data is discussed. INTRODUCTION One of the most critical problems in evaluating an aquifer gas storage project is determining the tightness of the caprock overlying the formation to be used as the storage reservoir. A formation that has previously held oil or gas obviously has a suitable caprock, but an aquifer that contains only water gives no such assurance. A number of aquifer projects in the United States have been troubled by gas leaking out of the intended storage zone, and the ensuing difficulties have led to the development of new evaluation methods. One of these new methods is pump testing wherein water is removed from the aquifer at some controlled rate prior to injection of gas. This fluid withdrawal causes a pressure drop to move out through the aquifer for considerable distances in a matter of days or weeks. Depending on the properties of the caprock, a pressure transient can also pass upward (as well as downward) through the caprock layers adjacent to the aquifer. Thus, if the operator has placed observation wells at appropriate distances from the pumping well, the rapidity with which the pressure transients reach different points in the system can be used to investigate the fluid transport properties of both the aquifer and its caprock. The usefulness of pump testing has been recognized by groundwater hydrologists for many years as a means of determining the potential yield and properties of aquifers used in water supply They have introduced the term "leaky aquifer" for a system in which an aquifer is overlain (or underlain) by semipermeable caprock layers. The ease with which water leaks into the aquifer during pumping can, of course, be very beneficial in bringing additional water to the pumped well. Hydrologists have therefore devoted considerable attention to this prob-lem. From the gas storage standpoint, however, the tighter the caprock layers that overlie the intended storage reservoir, the better are the conditions for minimizing or eliminating any vertical migration of gas. Thus, after a suitable geologic structure has been found, the emphasis in aquifer storage projects is in determining that the caprock is tight. Attention has recently been focused on the use of pump testing as one approach to solving this problem.23,24 This paper presents a further development on evaluating the permeability of a slightly leaky caprock when the caprock is of infinite thickness. From the practical standpoint, this means that the caprock layers are thick enough that pressure transients do not reach the outer boundaries of the system during the pumping test. In a subsequent paper, an analysis of the case where the caprock is of finite thickness will be presented. PREVIOUS WORK ON LEAKY AQUIFERS Jacob" developed a partial differential equation describing the flow of water in an aquifer of permeability k that is overlain by a leaky caprock of permeability k'. Fig. 1 shows a schematic cross-section of the system under consideration. One of his principle assumptions was that if k > > k', the direction of flow is essentially vertical in the caprock and horizontal in the aquifer. Neuman19 confirmed the validity of Jacob's assumption using a mathematical model. Another assumption was that a permeable source layer overlies the caprock (Fig. 1) and is able to maintain a constant hydraulic head at the upper boundary of the caprock. By neglecting the effects of compressibility within the caprock, Jacob1* developed a solution for a bounded circular aquifer. Later, Hantush and Jacobx5 used the same assumptions to solve the case of an infinite radial aquifer that is pumped at a constant rate. Their solution may be expressed in di-mensionless parameters by

By P. A. Witherspoon, R. W. Donovan, T. D. Mueller

The use of petroleum-barren aquifers for underground storage has become extremely important to the natural-gas industry. A critical problem in assessing the feasibility of a specific aquifer for such use is the permeability determination of the caprock over the proposed storage project. The approach used here is to conduct both static and dynamic field tests on the aquifer being analyzed. Valuable information on the possibility of communication between the storage aquifer and any other aquifers above can be obtained by measuring hydrostatic water levels and water analyses. Significant differences in such data give evidence of the lack of communication between the intended storage reservoir and other horizons. The dynamic approach requires that one well be pumped in the storage aquifer, and changes in fluid levels recorded in both the aquifer and its caprock. The interpretation of the data from such pumping tests involves the solution of nonsteady radial flow in an infinite aquifer and the influence on such flow of a leaky caprock. A finite-difference model has been used to investigate this problem, and the transient behavior has been solved numerically with a digital computer. It has been found that the pressure transients in the storage aquifer are not affected significantly by moderate caprock leakage. The pressure behavior of the caprock is a much better indicator of the degree of leakage, and generalized solutions for this behavior are included. Field data are presented to demonstrate both the static and dynamic approach. If is concluded that appropriate investigation of the groundwater hydrology in an aquifer-type gas-storage project can provide much valuable information for determining the effectiveness of the caprock to hold gas. INTRODUCTION Underground storage of natural gas in the United States has been developing at a rapid rate over the past few years. In 1955, the total gas-storage capacity was about 1.6 trillion cu ft; by 1961, this figure was almost 3.2 trillion cu ft, an increase of 100 per cent in six years.' This trend un- doubtedly will continue because the economics favor the development of gas storage, as opposed to the construction of new pipelines, to meet the inherent cyclic demand for fuel in the metropolitan areas of this country.' About 15 per cent of the current underground gas storage has been developed in petroleum-barren aquifers, i.e., geological domes or anticlines in which no commercial quantities of oil or gas had been produced prior to the storage operations. The necessity for using barren aquifers outside many metropolitan areas of this country has been due to the lack of depleted oil or gas fields that were near enough and large enough to meet the demands of such consuming areas. Pipeline companies have developed aquifer storage along their transmission lines to meet the fluctuating needs of their complex systems. Considerable thought has also been given to the problem of storing gas in a structureless aquifer, both in this country' and in the Soviet Union outside the city of Leningrad.'," Conditions such as these have led to the development of aquifer gas-storage projects in many parts of the U. S. Most of these developments have centered in the Mid-Continent area, and the greatest amount of activity has been concentrated in Illinois.6 Thus, the use of petroleum-barren aquifers for gas-storage purposes has become extremely important to the natural-gas industry. There are three basic problems in developing aquifer-type storage: (1) finding an adequate geologic structure, (2) finding a suitable storage reservoir within the structure and (3) determining the tightness of the caprock over the intended storage zone. The first two problems can be solved by applying conventional methods of exploration geology, but once these problems are solved, the question arises as to why no oil or gas is present in an otherwise favorable setting. Two situations are possible: (1) an adequate source bed was never present, or (2) a source bed was present but the petroleum seeped away because of a leaky caprock. Determining the tightness of the caprock is one of the most critical problems in assessing the feasibility of a specific aquifer for storage purposes. In attacking this problem, one usually takes cores of the caprock and subjects them to a rigorous investigation. Such core data are desirable, but they only detail the matrix properties and cannot be expected to reveal the gross characteristics of the caprock. Several gas-storage projects in the U. S. have had considerable leakage where

By Donald L. Katz, Leslie J. Noaker

Experimental data are presented for hydrate formation conditions for mixtures of hydrogen sulfide and methane. Vapor-solid equilibrium constants for hydrogen sulfide are also presented. Equilibrium constants, defined here as the ratio of the mole fraction of hydrogen sulfide in the vapor phase (dry basis) to the mole fraction of hydrogen sulfide in the solid phase (dry basis), were calculated using the vapor-solid equilibrium constants for methane presented in an earlier paper. The conditions at which a natural gas will form gas hydrates with water may be computed from its composition, when hydrogen sulfide is present in addition to carbon dioxide and the paraffin hydrocarbons. INTRODUCTION Many gases are known to form solid hydrates in the presence of liquid water at temperatures above 32°F. This hydrate formation is of importance to the natural gas industry inasmuch as most of the constituents found in natural gases will form hydrates. The attention of the industry was called to hydrates by Hammer-schmidtl0 in 1934. A number of papers have been written giving the temperatures and pressures at which these hydrates will form3,4,5,6,7,8,9,10,13,14,15,18,19 and a bibliography on the subject has been compiled13. Deaton and Frost4,5,6,7,8 report work for various pure gases and for some natural gases. Katz and co-workers3,4,13,14,18,19 have reported data for various gases at high pressure. They developed a method of predicting conditions for hydrate formation from equilibrium constants3,19 and from gab gravity." The effects of anti-freeze agents on the formation of hydrogen sulfide hydrate were reported by Bond and Nelson.' The hydrogen sulfide-water system, including hydrate, has been the subject of recent investigation by Sage and co-workers." Mixtures of hydrates behave as solid solutions and an equilibrium exists between the vapor phase and the solid phase that appears to be comparable to the equilibrium between vapor and liquid phases.3,10 By means of the vapor-solid equilibrium constants which have been published for methane, ethane, propane, iso-butane, and carbon dioxidel6 it is possible to determine the temperature and pressure at which hydrates will form for a specific gas mixture when water is present. Hydrogen sulfide which is frequently found in natural gases has a very strong tendency to form hydrate and must enter the solid solution when present in a natural gas forming a hydrate. To determine the effect of hydrogen sulfide when computing the conditions for gas hydrate formation, it is necessary to have the equilibrium constant giving the ratio of the concentration of hydrogen sulfide in the gas phase to the concentration in the solid phase. This paper presents data on hydrate formation for mixtures of hydrogen sulfide and methane and presents the equilibrium constants which were calculated from these data. EXPERIMENTAL PROCEDURE The apparatus used to determine the conditions of hydrate formation for the methane-hydrogen sulfide system was similar to that employed by previous investigators.3,18,19 A charge gas cylinder was used to supply the hydrogen sulfide-methane mixtures, called the charge gas, and to supply the pressure up to 1.000 psia. This

Jan 1, 1955

By John Riegle

To date. utility company underground storage of gas has generally been restricted to depleted dry gas fields. The Playa del Rey project is probably the first to successfully store gas in a partially depleted oil reservoir with the recovery of large volumes of pas at high rates as the main objective. Oil recovery has been a secondary consideration. Problems encountered unusual to those of storage in a dry gas zone were: (1) the removal and retardation of formation of emulsion. (2) the upstructure movement of fluid during withdrawal periods which formed fluid blocks, and (3) reservoir shrinkage resulting from encroachment of edgewater. The solution of these problems as outlined in the paper has resulted in increasing the withdrawal rate from the original design capacity of 4,000 Mcf per hour to the current 10.000 Mcf per hour. Additional increase is anticipated as the operations continue. Migration. reservoir performance. operational procedure and a historical record are included to make this a resume of the project from its war-time inception in 1942 to the first of 1952. INTRODUCTION The Playa del Rey underground gas storage project is probably the first project in which gas has been stored in a partially depleted oil reservoir, with the recovery of large volumes of gas at high rates as the main objective and the recovery of oil a secondary consideration. Problems unusual to those of storage in dry gas zones have made this a pioneer endeavor. The total storage capacity of approximately 1,500,000 Mcf of this reservoir is small in comparison with other underground gas storage projects. However. by study and experiment the deliverability into transmission lines. operating at pressures in excess of 150 lbs., has been nearly tripled in the past three years. to the current 10.000 Mcf per hour late, without the use of compression. It is anticipated this rate will gradually increase as additional fluid is removed from the formation. LOCATION AND HISTORY Playa del Key is located on Santa Monica Bay about two miles south of the beach commuity of Venice and 15 miles southwest of the renter of the City city Los Angeles. Fig. 1. The discovery well was completed at 6,194 ft in 1929. An upper zone at 3,905 ft was found on June 18. 1930. Townlot development was very rapid and by the end of 1930 there were 141 producing wells in the field. The Union Oil Co. completed a lower zone well south of the Townlot Field in May, 1931, which was believed to be an extension, until additional drilling of about 50 wells during 1934 and 1935 indicated this area to be a separate accumulation as shown on Fig. 2. The lower oil horizon, which is now the storage zone, had no original gas cap and is believed to be a continental detrital deposit of conglomerate that accumulated in valleys and depressions of a weathered schist high. which later subsided and was overlain by an impervious nodular shale. This type of formation has rapid variation of porosity. permeability and thickness due to lack of sorting action common to marine deposition. and the uneven surface of the original schist high. Production was limited on the upstructure sides by the pinching out of the conglomerate, and downstructure by edgewater. Unrestricted production, together with high permeabilities, led to large wastage of gas from this new area as facilities were not available for transport of gas to market. The decline was very rapid and by 1942 the portion of the new area south of Ballona Creek was nearly depleted. ACQUISITION OF PROPERTY The rapid increase in industrialization of Southern California during World war 11 forced upon the local gas com-

Jan 1, 1953

By R. J. Burgess, A. R. Ramey, A. R. Adams

During interpretation of pressure buildup tests on gas wells in a tight dolomite gas reservoir, peculiar behavior was noticed. Two straight lines were apparent. Effective permeability to gas taken from either straight line was about the same, and the Miller-Dyes-Hutchinson dimensionless time check for the straight line was proper for both straight lines. Geological data indicated the likelihood of scattered trending fractures in the reservoirs. Since the first straight Iine yielded permeability values close to the geometric mean permeability from core analyses, it was postulated that the reservoir model was that of an acidized well completed in the tight dolomite, but that widely scattered hairline fractures caused the mean permeability of the reservoir distant from the well to be higher than the matrix permeability. Because all other studies of fractured reservoirs to the authors' knowledge assumed that the fracture matrix was dense enough to communicate directly with the well, no interpretative methods were available. The Hurst line-source solution for a radial change in permeability for interference between oil reservoirs was adapted to pressure buildup testing. The result indicated that the first straight line should yield the proper matrix permeability and wellbore skin effect. The second straight line may be extrapolated to obtain static pressure. The time of bend between the straight lines was used to estimate distance to a fracture. Application to field test data is shown. It is believed that the methods developed and the case history presented will add to present tools available for pressure buildup interpretation. Introduction Since the pioneer studies by Miller, Dyes, and Hutchin-son1 and Horner' in 1950 and 1951, well test analysis has become recognized as one of the most powerful tools available to both production and reservoir engineers. Well test analysis serves as a logical basis for well stimulation and completion analysis, and for long-term reservoir engineering. Since the early 19501s, much effort has been placed on the development of well-test analytical methods. Reservoir and well conditions of increasing complexity have been considered systematically to provide the analyst with a catalog of causes and effects. Matthews and Russella state that some 200 papers dealing with this subject have been published in the last 35 years. Developments in well test analysis appear to have originated in one of two ways. Either a physically realistic field condition was anticipated and analytical solutions for the condition achieved, or anomalous field test behavior was recognized and interpretative methods sought for the anomaly. In recent years, it has appeared that the latter has inspired an increasing number of studies. The analyst today finds an increasing number of known cause and effect studies available for well test analysis, the classic of which is that of finding the specific flow problem that generated the answer — the well behavior. Although it may be impossible to achieve this goal uniquely, the analyst often is able to select a useful interpretation that combines all known performance and geologic data — or to show that various logical alternatives would not significantly affect the interpretation. During a recent reservoir study, we observed gas well test behavior that did not appear to fit behavior described previously. Although it cannot be said that we have found a unique interpretation, we shall present in this paper the peculiar behavior observed, and describe the reservoir and interpretative methods developed. Reservoir Description The subject gas reservoir is a 9-mile-long, narrow dolomite reservoir lying within a limestone of Ordovician age. (See Fig. 1.) The dolomitized rock in the field consists of dark brown to buff, dense to coarsely crystalline, vugular dolomite containing numerous hairline fractures, many of which may have been closed in the reservoir and parted when cores were brought to the surface. Larger fractures are also apparent in core, but usually are filled and sealed with euhedral dolomite crystals. Portions of the north flank of the reservoir are known to be cut by a sealing fault downthrown to the north. Gas wells located near the fault have higher open flow potentials than those more distant from the fault. This is believed to be a result of higher permeability near the fault due to more extensive and open fractures. Detailed coring and core analysis have been performed on several of the wells in this reservoir. Fig. 2 presents permeability variation' plots for both horizontal and vertical

Jan 1, 1969

By R. A. Wattenbarger, H. J. Ramey

A systematic study has been made of the application of the real gas pseudo-pressure m(p) to short-time gas well testing. The m(p) function can be used in real gas flow problems to account for the variation of viscosity and z-factor with pressure. A mathematical model was solved numerically to generate solutions of various real gas flow problems. The effects of turbulence, formation damage and wellbore storage were included in the model. The analysis of simulated well tests showed that the interpretation methods normally used for liquid flow are generally accurate when the m(p) is used. For practical rates without turbulence, the solutions are not rate-sensifive, and the flow capacity.. kh and skin effect can be determined accurately from either a buildup or a drawdown test. However, the kh calculated from a drawdown test can be significantly low when turbulence is present. A case was simulated in which this error wa.s 36 percent. Turbulence does not affect the determination of kh from buildup tests. The proper determination of kh from p and pa buildup and drawdown plots is developed by analyzing their relationship to the m(p) method. A simple equation is given that can be used for long-range gas well performance forecasting. This expression is compared with a method pre-sended by Russell et al." Introduction The formation flow capacity and wellbore damage condition can be determined for liquid-producing tests by means of buildup and drawdown tests. These tests make use of short-time pressure transient data rather than stabilized flow tests that are often used for gas well testing. Tracy' presented a method of gas well testing that was based on ideal gas equations and utilized p2 buildup plots. This method was shown to be good for low pressure wells. MatthewsZ suggested plotting p and using an average slope and average gas properties. This method was more successful on high pressure wells. Al-Hussainy, Ramey and Crawfordh howed that variation in gas properties could be simplified by using the real gas pseu do-pressure m(p). Some cases showed that the use of the m(p) function Original manuscript received in Society of Petroleum Engineers office Aug. 8. 1967. Revised manuscript received June 11, 1968. Paper (SPE 1835) was presented at SPE 42nd Annual Fall Meeting held in Houston Tex., Oct. 1-4. 1967, and at SPE 4th Annual Eastern Regional ~eetin$ held in Pittsburgh. Pa., Nov. 2-3. 1967. 0 Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. 'References given at end of paper. provided an accurate method of interpreting buildup and drawdown tests. Because the work of Al-Hussainy et al. included only a few actual cases, there was a need to explore the use of the m(p) method for a greater variety of flow conditions. The occurrence of turbulent flow around the wellbore often is an important factor in gas well testing. Swift and Kiel,' and Carter et a1.' treated turbulence in gas well testing, but did not consider the variation of viscosity and z-factor with pressure. This paper presents the results of an investigation of the application of the m(p) method to buildup and drawdown testing. Vhe investigation included the effects of turbulence, formation damage and wellbore storage. Gas Flow Equations To formulate the mathematical model, many of the assumptions usually used in well testing theory are applied. The system has radial geometry with a closed outer boundary and is composed of a horizontal porous formation that has uniform and isotropic rock properties and uniform thickness. Allowance is made, however, for a radial region of reduced permeability near the wellbore. This region represents formation damage. The geometry of the system is shown in Fig. 1. Darcy's law does not always apply to gas flow. A more general expression is needed for non-Darcy, or "turbulent" REGION OF DAMAGED PERMEABILITY SEALED BOUNDARIES ClRCULAR WELL BORE Fig. 1—Radial flow model. AUGUST, BT7

Jan 1, 1969

By Jens P. Nielsen, W. H. Justice

This paper describes an improved acidizing technique which has been applied in acidizing gas wells in the La Gloria Field. Wells acidized in this manner exhibited a much greater increase in deliverahility than wells acidized in the conventional manner. The acid is injected into the well in small slugs separated by a small volume of high pressure gas. The acid and gas are displaced from the tuhing by high pressure gas to reduce the hydrostatic head on the well when the acid is recovered from the well. The well is opened for production as soon as the acid and gas are displaced. This procedure permits a very rapid recovery of the acid which results in large increases in the deliverability of the gas wells. It has been found advantageous in the acidization of gas injection wells to inject the high pressure gas directly after the acid without backflowing the acid out of the well. This practice has made it possible to inject the gas in the injection wells of this field with a much lower pressure differential. This reduces the horsepower required to inject the gas and also decrease:, the number of injection wells required per reservoir. INTRODUCTION The results obtained by acidizing gas wells ill the La Gloria Field have varied considerably in the past. For this reason. a planned study of the acidizing procedure was undertaken in order to obtain more consistent results. The Frio Sands of the La Gloria Field are rather permeable and it was not necessary to acidize many of the wells when they were originally completed. At the completion of an eight-year cycling program in several of the gas reservoirs, it became necessary to recomplete many wells in other reservoirs. It was found that the opan Flow potential of the wells recompleted in different sands was not as large as wells that were originally completted in that reservoir. A study of the drilling records of several of these wells indicated that the sands had originally been subjected to considerable hydrostatic drilling fluid overload while the well was being drilled to the deeper sands of the field. It is indicated that during this eight-year period, the drilling mud and water that were lost into the formation when the well was being drilled had formed a relatively impermeable Eone around the well bore. This impermeable zone around the well bore made it very difficult to obtain high deliverability wells when the original wells are recompleted in other zones. Muskat' has shown that in uniformly permeable and porous formations a zone of partial plugging one ft in radius around the well would drastically decrease the productivity of the well. The analysis for this problem shows that the production capacity of a well is very sensitive to the value of the permeability of the zone immediately surrounding the well bore. This is to be expected in view of the highly localized character of the pressure drop in a radial flow system about the well center. The reduction of the deliverability is not proportional to the thickness of the mud-affected zone as the major portion of the reduction in deliverahility IS caused hy the

Jan 1, 1952

By R. V. Smith, E. J. Dewees, R. H. Williams

Before this paper was presented, the terms "f" and "v1/f in flow equations had been called "friction coefficient" and "friction factor," respectively. However, many authors have used the names for the terms interchangeably, with the result that considerable confusion has existed in the literature on flow of fluids. In an effort to establish a clear and concise terminology for these terms, the authors have designated f as "resistance coefficient" and dm as "transmission factor." Transmission factors for the flow strings of gas wells were determined by tests on wells in six fields of Oklahoma, Texas, and Louisiana. Data were collected on 14 gas wells while gas was flowing through 1 1/4, 2, 21/2, and 3-in. tubing and 51/2 and 7-in. casing. The transmission factors developed from these data were found to follow the resistance equation published by Nikuradse1 for turbulent flow in "rough" pipes with an average absolute roughness of 0.000,65 in. This average absolute roughness is believed to be representative generally of the flow strings of gas wells with the exception of those that have been exposed to corrosion or other conditions causing undue roughness. Also, use of an absolute roughness of 0.000,65 in. probably will result in transmission factors that are more representative of gas-well conditions than those in common use today, provided that the transmission factors are calculated by the methods outlined in this report. The authors believe that use of these transmission factors will permit calculation of more accurate subsurface pressures in gas wells than those in use today. Engineers use transmission factors in flow equations to compute subsurface pressures in flowing gas wells from measured wellhead pressures. These subsurface pressures are necessary for determining the productivity of gas wells from "back-pressure" tests. The methods of calculating pressures in flowing gas wells, as outlined by Rawlins and Schellhardt2 in Bureau of Mines Monograph 7, are based on the transmission factors in the well-known Weymouth formula" Natural gas engineers have realized that errors introduced by using transmission factors as given by the Weymouth formula are relatively unimportant in testing low-capacity gas wells; but they also know that such factors are important considerations in testing large-capacity wells. Accordingly, the research described in this report was undertaken to fulfill the need for more accurate means of calculating subsurface pressures in gas wells. This report is part of the Bureau's study of the flow of natural gas through pipe and commercial pipelines. Several of the conclusions regarding flow in the flow strings of gas wells were influenced by the results of the study of flow in horizontal pipelines. These instances are mentioned in the text of this report. METHODS OF COMPUTING TRANSMISSION FACTORS Methods and procedures for computing transmission factors for turbulent flow in circular pipes from test data taken on gas wells have been derived and explained in detail in a paper by Smith4. As these methods were used to calculate the transmission factors presented in this report, a summary is included. Starting with

Jan 1, 1955

By F. W. Jessen, R. W. Durie

The development in the laboratory of mall salt cavities illustrated that the rate at which salt dissolved was affected by the rate at which fresh water was injected into the cavity. Experiments carried out with small salt cores indicated, however, that the rate of salt removal with forced convection was not significantly different from that with natural convection, provided the flow was in the laminar range. The rate of salt removal was apparently controlled by the water salinity, which was indirectly determined by the rate at which fresh water was injected. With turbulent flow at the salt surface, the rate of salt removal could be increased several-fold with the same salinity in the cavity. A quantitative evaluation of the rate of salt removal under natural convection was attempted by using techniques developed for the analogous heat convection system. The values of the rate of salt removal obtained experimentally were considerably greater than the calculated values due to surface irregularities. To determine the effect of these irregularities, the rate of salt removal was found for a smooth surface by experimentally determining the rate of salt removal as an irregular surface developed and then extrapolating this trend back to initial time. Excellent agreement between the experimental and calculated values then were obtained. INTRODUCTION Underground cavities leached from massive salt formations have been used extensively for many years for storing LPG products. The development techniques for storing methane gases in the liquid state at a low temperature leads to speculation that underground salt cavities of controlled geometry may be quite important in this field. In any case, the increasing number of salt cavities used for storage indicates an impending need for a better knowledge of the washing process to efficiently develop maximum storage volumes within specified regions. It is the purpose of this paper to discuss some aspects of -the mechanism of the dissolution of salt in the formation of cavities of controlled geometry. Although no attempt is made to resolve the problems that may be associated with the storage of methane in the liquid state in salt formations, the fact that salt cavities can be formed at reasonable cost would appear to justify investigation in this direction. Further, underground cavities have the advantage of allowing storage at higher pressures (and hence higher temperatures) than surface installations. Although the thermal conductivity of salt is somewhat higher than frozen earths, 0.017 cal/(sec) (cm)(c) for salt at 32F vs 0.005 cal/(sec) (cm) () for ice at this temperature, this disadvantage is offset by the higher storage temperature in underground cavities. For example, at atmospheric pressure, the temperature of liquid methane would be —258F. In a storage cavity at 2,000 ft, the pressure would be approximately 1,000 psia, which would allow a storage temperature of -9OF. The formation temperature would be roughly 110F at this depth, hence the total temperature differential would be 200F. Under surface storage conditions the corresponding temperature difference would be some 333 F (75 F + 258 F). The results of this difference can be compared more readily by substituting the appropriate values of thermal conductivity, temperature difference, heat capacity and density into the solution of the diffusion equation for the unsteady linear conduction of heat. For the underground salt storage, the heat flux is 15cal/(sec) (sq cm), where t is the time in seconds from the initiation of the cooling. For the surface system, and using the thermal properties of ice at the average temperature, the corresponding figure is 21cal/(sec) (sq cm). Hence, the heat losses through ice under surface conditions would be some 30 per cent higher than the underground salt storage. This ratio would vary, considerablv. Of course, for other conditions, but it indicates that an advantage over surface facilities may be obtained by storage in underground salt cavities.

Jan 1, 1965

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By M. R. Tek, F. H. Poettmann, M. L. Grove

The general problem of single-phase natural gas flow through porous media has been considered both by mathematical analysis and field experiments. Theoretical studies on the fluid dynamics of natural gas wells resulted in derivations of expressions which permit a ber-ter understanding of the flow mechanism for conditions under which Darcy's law is no longer valid. In applying these results to natural gas production and testing technology, special emphasis war placed on the study of low permeability, slowly stabilizing formations. These applications resulted in a method of predicting the performance coefficients of tight natural gas wells. The calculation of the performance coefficient requires the determination of the "test-index", pkh/µ, the drain- age radius, R, and the effective wellbore radius, rw. Methods for estimating these parameters from field data were also developed. Field data from 15 wells were collected and analyzed. The measured values of the performance coefficients were compared with the predicted values. This statistical evaluation gave a bias of 7.14 per cent and a standard deviation of 11.6 per cent. INTRODUCTION The performance of a natural gas well is determined by the physical properties of the reservoir rock, the extent and the geometry of the drainage area, the properties of the flowing fluid, and the conditions of pressure distribution within the drainage matrix. The relation between the quantity of gas delivered and the pressure drop within the reservoir is characteristic of the behavior of each well and is often referred to as the back-pressure behavior of the well. It is well established that for normal, single-phase gas wells the relationship between the gas delivery rates and the bottom-hole pressure take, in general, the form, Q = c(P3 8 -P2 1)n. ,.......(1) where C is the performance coefficient, and n is the exponent corresponding to the slope of the straight-line relationship between Q and (P2 s, — P2 l) plotted on logarithmic coordinates. The determination of a well's performance by the back-pressure behavior not only permits the production engineers to analyze many operating problems, but also provides necessary information useful and essential to predict the future development of the field.' Important engineering and production problems such as calculation of gas deliverability into a pipeline at a predetermined line pressure, design and analysis of gas-gathering lines, determinations of spacing and number of wells to be drilled during field development to meet future market or contract obligations, etc., all depend, at least in part, on the availability and use of reliable back-pressure curves.' In the past the behavior of gas wells was evaluated

Jan 1, 1958

By R. E. Schilson, F. H. Poettmann

The direct determination of the stabilized performance behavior of low capacity, slowly stabilizing gas wells is extremely time-consuming and wasteful of gas. From both field experience and theoretical considerations, a test procedure has been evolved by which the stabilized hack pressure behavior of such gas wells can he predicted without having to revert to long time flow tests. The method consists of using the isochrona1 test procedure to establish the slope of the back pressure curve, "n", and the short time variation of the performance coeficient, "C", with time. From this short time transient flow data and theoretical considerations, the value of C at large times can he established. By assuming the radius of drainage of a well to be half the distance between wells, one can calculate the stabilization time for various well spacing patterns. Once the stabilization time for a given .spacing has been determined, the value of C can be calculated and the stabilized back-pressure curve can he. establbrhed. The calculated perfortnance coefficient as a function of tinze was cotnpared to the experintentally measured values for a number of gas wel1.e. The deviation of the calc~*lated from the cxperimental res1tlt.s vary depending on the set of short time experimental points used to evabrate the parcrttzeters of the equation. The longer the tirne far the flow test data user1 in the calculations, the better was the agreement with the experimental results. The time necessary to obtain this data from well tests varies consitlerably, depending on the physical natrcre oj the rt3scrvoir under consideratiot~. INTRODUCTION For many years, the U. S. Bureau of Mines Monograph 7' has served as a guide for testing and evaluating the performance of gas wells by means of the back-pressure method. The back-pressure performance of a gas well is expressed by the following equation: Q~ CW-W.........(I) where the characteristics of the back-pressure equation are determined by C, the performance coefficient, and n, the exponent which corresponds to the slope of the straight line when Q and {P* - PN2) are plotted on logarithmic paper. Q is gas flow rate at standard conditions, and P, and P, are equalized and flowing bottom-hole pressures, respectively. Prior to the development of the back-pressure tcst, the "open flow" capacity method of testing a well was common. By this method, the wells were flowed wide open to the air and the flow rate measured. Such procedure was wasteful of gas and did not provide information on the deliverability of the gas to the pipe line. Monograph 7 Procedure The back-pressure method of testing wells was developed to overcome these shortcomings. Although much has been learned regarding the laws of the flow of gas through porous formations, the original development of the back-pressure relationship was based entirely on empirical methods. The back-pressure behavior provides the engineer with information essential in predicting the future development of a field. It permits him to calculate the deliverability of gas into a pipe line at predetermined line pressures, to design and analyze gas gathering lines, to determine the spacing and number of wells to he drilled during the development of a field to meet gas purchasers' requirements, and to solve many other technical and economic problems. As described in Monograph 7, the flow-after-flow method of back-pressure testing, when applied to fast stabilizing and usually high capacity wells, correctly characterized the behavior of the well. However, as the value of the gas at the wellhead increased, small capacity gas wells having slow rates of stabilization became economically operable. The flow-after-flow method of testing could not be used to describe the behavior of these slowly stabilizing wells. The procedure of Rawlins and Shellhardtl for establishing the back-pressure behavior of a gas well was based on the rcquirement that the data be obtained under stabilized flow conditions; that is, that C is constant and does not vary with time. C depends on the physical properties of the reservoir, the location, extent and geometry of the drainage radius, and the properties of the flowing fluid. In a highly permeable formation, only a very short period of time is required for the well to reach a stabilized condition, and, consequently, the requirements for the test procedure described in Monograph 7 are met. For a given well, n is also constant

By H. O. McLeod, J. M. Campbell

This paper presents the results of the data obtained in the first stage of a long-range study at high pressures of the system, vapor-hydrate-water rich liquid-hydrocarbon rich liquid. The data presented are for the three-phase systems in which no hydrocarbon liquid exists. Tests were performed on 10 gases at pressures from 1,000 to 10,000 psia. One of these was substantially pure methane, and the remainder were binary mixtures of methane with ethane, propane, iso-butane and normal butane. Several conclusions may be drawn from the data. 1. Contrary to previous extrapolations, the hydrocarbon mixtures tested form straight lines in the range of 6,000 to 10,000 psia which are parallel to the curves for pure methane, when the log of pressure is plotted vs hydrate formation temperature. 2. The hydrate formation temperature may be predicted accurately at pressures from 6,000 to 10,000 psia by using a modified form of the Clapeyron equation. The total hydrate curve may be predicted by using the vapor-solid equilibrium constants of Carson and Katz' to 4,000 psia and joining the two segments with a smooth continuous curve between 4,000 and 6,000 psia. 3. The use of gas specific gravity as a parameter in hydrate correlations is unsatisfactory at elevated pressures. 4. The hydrate crystal lattice is pressure sensitive at elevated pressures. INTRODUCTION Prior to 1950 many studies had been made of the hydrate forming conditions for typical natural gases to pressures of 4,000 psia.""'"'"" Most of these attempted to correlate the log of system pressure vs hydrate formation temperature, with gas specific gravity as a parameter. One of the more promising correlations was made by Katz, et al, which utilized vapor-solid equilibrium constants. The only published data above 4,000 psia are those of Kobayashi and Katz7 for pure methane to a pressure of 11,240 psia. In the intervening years, most published charts for the high-pressure range have represented nothing more than extrapolations of the low-pressure data, with the methane line serving as a general guide. The reliability of these charts has become increasingly doubtful (and critical) in our present technology as we handle more high-pressure systems. The portion of our high-pressure hydrate research program reported here was designed to: (1) investigate the reliability of existing charts; (2) obtain actual data on gas mixtures to 10,000 psia; and (3.) develop a simple hydrate correlation that was more reliable than those which simply used specific gravity as a parameter. Binary mixtures of methane and ethane, propane normal butane, or iso-butane were injected into a high-pressure visual cell containing an excess of distilled water. Hydrates were formed and then melted to observe the decomposition temperature of the hydrates at pressures from 1,000 to 10,000 psia. EQUIPMENT The equipment consisted of a Jerguson 10,000-lb high-pressure visual cell, a 10,000-1b high-pressure blind cell and a Ruska 25,000-1b pressure mercury pump. The visual cell was placed in a constant-temperature water bath controlled by a refrigeration unit and an electric filament heater. A Beckman GC-2 gas chromatograph was used in analyzing the gas mixtures after each run was completed. EXPERIMENTAL PROCEDURE After evacuating the gas system, the heavier hydrocarbon was injected into the high-pressure mixing cell to that pressure necessary to give the desired composition. This cell then was pressured to 1,100 to 1,200 psia by methane from a high-pressure cylinder. The mixing cell holding the gas contained a steel flapper plate and was shaken intermittently over a period of 15 minutes. After mixing, the valve to the high-pressure visual cell containing excess distilled water was opened, and the gas mixture was allowed to flow into the cell. The temperature in the water bath was lowered 10" to 15'F below the estimated hydrate decomposition point. As a first check, the temperature was increased at a rate of 1°F every six minutes to find the approximate point of decomposition. It was again lowered 1.5° to 5°F to form hydrates. The temperature was raised to within l° of the estimated decomposition point and then increased 0.2F every 10 to 15 minutes until the hydrates decomposed. This procedure was repeated at various pressures to obtain 7 to 13 points for each mixture between 1,000 and 10,000 psia. After completion of the hydrate decomposition tests, the gas mixture composition was analyzed with a calibrated gas chromatograph. These gas analyses have an estimated error of ± .1 per cent.

By Anthony Folger, Howard R. Lowe, J. F. Dougherty, E. G. Trostel, Joffre Meyer

The status of the natural gas industry in Alberta, Canada, is described with particular reference to the current extent of natural gas reserves and possibilities for additional development. Certain major fields are discussed as representative of the types of gas accumulations within various broad geographical divisions of the provinces. The estimation of the future performance of a gas reservoir is discussed as the basis for an estimate of the future availability of natural gas from presently known reserves. The future availability of the provincial reserves and the estimated future demand for gas within the province are described. INTRODUCTION This paper is a condensation of a detailed study of the natural gas reserves and the future availability of pipeline gas in the province of Alberta as prepared for Trans-Canada Pipe Lines Limited and presented in public hearings at Edmonton and Calgary before the Petroleum and Natural Gas Conservation Board of the Province of Alberta. The study began in October, 1950, and resulted in the preparation of five volumes of data at the expenditure of more than 20,000 man-hours. This paper is designed to afford a broad picture with detail limited to a type example of both natural ga-accumulation and future availability of gas supply. The province of Alberta has an area of 255.285 square miles which, for sake of comparison. is 8,359 square miles smaller than the state of Texas. Natural gas was first discovered in 1883 at the now abandoned Alderson well located some 20 miles southeast of the town of Princess. The first gas in commercial quantities was found at Medicine Hat in 1890. During the last 68 years. 187 localities have been discovered which are capable of gas production, and 73 per cent of these have been found since 1945. Sixty fields currently are producing gas, including 22 non-associated gas accumulations. Location of the oil and gas fields and prospects in the province of Alberta is shown in Fig. 1, together with a list of the localities of measurable gas, Fig. 1A. The estimated total provincial gas reserve as of Jan. 1, 1952, for the 104 fields estimated, is 11.7 trillion cu ft proved and probable, and for the 187 known areas capable of production is in excess of 16 trillion cu ft on a total proved, probable, and possible basis. Five fields (Pincher Creek, Leduc-Woodbend, Viking-Kinsella, Cessford and Medicine Hat) are estimated to have recoverable reserves of more than one trillion cu ft each. A sixth field, Turner Valley, was initially in this group, but its reserves have been depleted substantially. Twenty-four fields, of which the above mentioned six fields are a part, have estimated gas reserves (including the possible category) in excess of 100 billion cu ft. These 24 fields account for 79 per cent of the total provincial reserves. Comparative systemic analyses reflect that 56 per cant of Alberta's gas reserves are contained in sediments of Cretaceous age; 30 per cent in sediments of Mississippian age; 10 per cent in sediments of Devonian age; three per cent in sediments of Triassic age; and the remaining one per cent in Jurassic and Permo-Pennsylvanian sediments. It is expected that the relative importance of the Mississippian, Devonian, and Triassic sediments will be increased as exploration progresses. GENERAL GEOLOGY The western Canada basin comprises the sedimentary area between the Canadian Shield and the Rocky Mountain Front Range, the Arctic Ocean, and the International Boundary. The basin has an approximate area of 765,000 square miles and embraces parts of the provinces of Manitoba. Saskatchewan. Alberta, British Columbia and the Northwest Territories. Generalized columnar sections for the Province of Alberta are shown in Fig. 2. Substantially all of Alberta, except its northeast corner, has sedimentary cover. Sediments of all the geologic systems are represented. except those of Silurian age.

Jan 1, 1952

By H. J. Ramey

The wellbore acts as a storage tank during drawdown and build-up testing and causes the sand-face flow rate to approach the constant surface flow rate as a function of time. This effect is compounded if non-Darcy flow (turbulent flow) exists near a gas wellbore. Non-Darcy flow can be interpreted as a flow-rate dependent skin effect. A method for determining the non-Darcy flow constant using this concept and the usual skin effect equation is described. Field tests of this method have identified several cases where non-Darcy flow was severe enough that gas wells in a fractured region appeared to be moderately damaged. The combination of wellbore storage and non-Darcy flow can result in erroneous estimates of formation flow capacity for short-time gas well tests. Fortunately, the presence of the wellbore storage eflect permits a new analysis which can provide a reasonable estimate of formation flow capacity and the non-Darcy flow constant from a single short-time test. The basis of the Gladfelter, Tracy and Wilsey correction for wellbore storage in pressure build-up was investigated. Results led to extension of the method to drawdown testing. If non-Darcy flow is not important, the method can be used to correct short-time gas well drawdown or build-up data. A method for estimation of the duration of wellbore storage effects was developed. INTRODUCTION In 1953, van Everdingen and Hurst generalized results published in their previous paper3 concerning wellbore storage effects to include a "skin effect", or a region of altered permeability adjacent to the wellbore. Later, Gladfelter. Tracy and Wilsey4 presented a method for correcting observed oilwell pressure build-up data for wellbore storage in the presence of a skin effect. The method depended upon measuring the change in the fluid storage in the wellbore by measuring the rise in liquid level. To the author's knowledge, application of the Gladfelter, Tracy and Wilsey storage correction to gas-well build-up has not been discussed in the literature. It is, however, a rather obvious application. Gas storage in the wellbore is a conlpressibility effect and can be estimated easily from the measured wellbore pressure as a function of time. Several approaches to the wellbore storage problem have been suggested. As summarized by Matthews, it is possible to minimize annulus storage volume by using a packer, and to obtain a near sand-face shut-in by use of down-hole tubing plug devices. Matthews and Perrine have suggested criteiia for determining the time when storage effects become negligible. In 1962, Swift and Kiel' presented a method for determination of the effect of non-Darcy flow (often called turbulent flow) upon gas-well behavior. This paper provided a theoretical basis for peculiar gas-well behavior described previously by Smith. Recently, Carter, Miller and Riley observed disagreement among flow capacity k,,h data determined from gas-well drawdown tests conducted at different flow rates for short periods of time (less than six hours flowing time). In the original preprint of their paper, Carter et al. proposed that the discrepancy in flow capacity was possibly a result of wellbore storage effects. Results of an analytical study of unloading of the wellbore and non-Darcy flow were recorded by carter.14 In the final text of their paper, Carter et al.!' stated that they no longer believed wellbore storage was the reason for discrepancy in their kgh estimates. In view of the preceding, this study was performed to establish the importance of non-Darcy flow and well-bore storage for gas-well testing. In the course of the study. a reinspection of the previous work by van Everdingen' and Hurst' was made, and the basis for the Gladfelter, Tracy and Wilsey' wellbore storage correction was investigated and extended to flow testing. WELLBORE STORAGE THEORY As has been shown by Aronofsky and Jenkins,11-12 Matthews," and others, flow of gas can often be approximated by an equivalent liquid flow system. The following developnlent will use liquid flow nomenclature to simplify the presentation. Application to gas-well cases will be illustrated later. First, we will use the van Everdingen-HursP treatment of wellbore storage in transient flow to establish (1) the duration of wellbore storage effects, and (2) a method to correct flow data for wellbore storage. DURATION OF WELLHORE STORAGE EFFECTS When an oil well is opened to flow. the bottom-hole pressure drops and causes a resulting drop in the liquid level in the annulus. If V. represents the annular volume in cu ft/ft of depth, and p represents the average density of the fluid in the wellbore, the volume of fluid at reservoir conditions produced from the annulus per unit bottom-hole pressure drop is approximately: res bbl-- (V, cu ft/ft) (144 sqin./sq ft) psi -(5.615 cu ft/bbl)(pIb/cuft) ........(I)

Jan 1, 1966

By A. Satter, J. M. Campbell

Reported herein are the results of a careful and detailed study of the non-ideal behavior of pure gases and their mixtures. Included are: (1) new data on five ternary systems composed of methane, ethane and H2 S; (2) a simple compressibility factor correlation that is inherently superior to present correlations, particularly for gases containing H2S and CO2; and(3) a detailed study of combination rules and the effect of system composition on the choice thereof. This study makes use of the rather large mass of data already available in the literature. A complete re-examination of the data and ideas presented in the last 25 years was considered desirable as a prelude to our basic concern — the effect of diluents on gas behavior. A consideration of both the macroscopic and microscopic properties of gases provides a better insight which, in turn, gives a firmer basis for improved correlation techniques. Such a study has shown that expressing the compressibility factor Z as a function of acentric factor w, as well as reduced temperature and pressure, yields a correlation that is broader in scope. The study of various combination rules has shown that better results are obtained by "tailoring" the rule used to the system composition. To do so improves the basic reality of results by overcoming some of the anomalies often found when using Kay's rule alone. Tentative recommendations are made regarding the most reliable combination rule for use with a given class of gas. The data presented are useful for estimating the direction and magnitude of the expected deviation when using a given rule. Although more work is needed, particularly around the critical region and with CO2 mixtures, the advantage of the classification scheme proposed is apparent. INTRODUCTION When one attempts to write a PVT equation to fit the data for actual gases, greater precision is obtained by the use of a multiple number of empirical constants. This has lead to multiple-constant equations such as Benedict-Webb-Rubin, Beattie-Bridgman, Keyes, etc., which are capable of yielding very precise results for pure gases in a range for which data to get the constants are available. As a matter of practicality, though, the use of such equations for gas mixtures is limited. Because of the infinite number of gas analyses available, any attempt to compile the constants needed requires a prohibitive amount of experimental data. This could be overcome by the use of a combination rule, but there is no real advantage in doing so because the end result offers no practical impovement over the Z factor correlation. The most widely used method of predicting the volumetric properties of pure gases is based upon the "theorem of corresponding states". According to this theorem, "all pure substances have corresponding molal volume at corresponding temperature and pressure if the reference point of correspondence is the critical point". Generalized compressibility charts for gases were prepared first by Cope and associates1 in 1931 and later by Brown and co-workers2 in 1932. However, the most commonly used charts are those of Dodge,3 Nelson and Obert,4 Hougen and atsson: and Standing and Katz.6 The work of Katz and co-workers has provided us with basic data for the hydrocarbons most widely used today. Their original chart6 was compared with a relatively large amount of multi-component data for gases consisting almost entirely of normal paraffin hydrocarbons. A deviation of only + 1.2 per cent was obtained.39 In the 20 years following publication of this work it has been found that the behavior of most mixtures of paraffin hydrocarbons could be predicted by this correlation within at least 5 per cent. Where difficulty has been encountered it has largely involved one or more of the following circumstances: pressures above 4,000 psig, mixtures containing large amounts of heavy ends and/or aromatics, systems in the critical region and mixtures containing polar compounds and/or CO2. The abnormal error sometimes found with such gases, not too unexpected for this method, is

By R. D. Carter

Methods are presented for calculating the performance of multiwell gas reservoirs. These methods account for two-dimensional, unsteady-state flow of a non-ideal gas through a heterogeneous formation. The performance can be obtained for gas reservoirs with known physical properties and configuration subjected to specified development and production schedules. These methods may prove helpful in optimizing the exploitation of gas reservoirs. The numerical solrrtion procedures given for this non-linear problem are extensions of known techniques for the linear diffusion equation. Two of these procedures are based on a new difference equation given by Saul'ev. Limitations and advantages of these procedures are discussed. An example solution is presented and is used to compare the calculation methods. INTRODUCTION Li planning the operation of a multiwell gas reservoir, it is desirable to predict its performance for each alternate development and producing schedule under consideration. This paper presents three numerical procedures that can be used to calculate such gas reservoir performance using a digital computer. These methods account for two-dimensional transient flow of a non-ideal gas within the reservoir. Extension of these methods to three-dimensional flow is possible. The paper is divided into three parts. The first describes the problem for which a solution is wanted and three numerical procedures that can be used to obtain solutions. The second discusses the information necessary to set up a problem. This part will be of greatest interest to the field engineer. The third gives data and results for an example problem. PROBLEM AND SOLUTION METHODS In the following development, a gas reservoir is regarded as a sealed, porous, permeable, heterogeneous body in which a single, non-ideal gas phase is flowing. Flow is permitted in the x- and y-horizontal directions, but is considered to be negligible in the vertical or z-direction. The reservoir can be tapped by one or more wells. Methods to be described provide a basis for calculating the complete performance of the reservoir, including rates and cumulative production from the wells, and pressure distribution and decline in the reservoir. If the isothermal flow of a non-ideal gas obeying Darcy's Law occurs in two dimensions in a region R, the governing equation in this region is:

Jan 1, 1967

By F. Kurata, P. C. Davis, T. L. Gore

An experimental method and apparatus for the study of the low-temperature phase and volumetric behavior of volatile mixtures are described. The phase diagrams, including the critical points, and the volumetric behavior of a natural gas and a mixture of the natural gas with methane at temperatures from -40°F to -200°F and pressures to 1.200 11si are givrn. INTRODUCTION Knowledge of the phase and volumetric behavior of mixtures of the light hydrocarbons at low temperatures is becoming increasingly important as interest in low-temperature separation processes increases. Distillation at low-lemperatures has been used for some time to recover ethane and ethylene from natural gas and refinery gas for petrochemical uses. Removal of excess nitrogen from natural gas by distillation has been proposed as a means of increasing the capacity of gas transmission lines and of improving the marketability of the gas.1,2 Evaluation of the latter process has been hampered, however, by the scarcity of data from which the phase behavior of natural gases at low temperatures might be predicted. While existing correlations of vaporization equilibrium constants may be used to predict the behavior of such mixtures at lower pressures, high-pressure computations often require knowledge of the location of the critical point of the system involved. Although correlations have been developed for the prediction of the critical temperature and pressure of certain types of hydrocar-bon systems, they do not apply to very volatile systems such as nitrogen-bear-ing natural gases. This paper describes the initial stages of an investigation of the low-temperature behavior of a numher of natural gases and other gases of commercial interest. The experimental method and apparatus are described. and the phase and volumetric behavior of a natural gas and a mixture of the natural gas with methane are presented. PREVIOUS WORK The published low-temperature vapor-liquid equilibrium data for hydrocarbons and mixtures of hydrocarbons with nitrogen and helium are listed in Table I. Most of these data are for binary systems, and the data of Stutz-man and Brown" at 100 psi are the only data on a natural gas. The data of Bloomer and Parent3 on the methafie-nitrogen system, which have been phblished since this study was completed, are a valuable addition to the knowledge in this field. Prior to the publication of these data the only published critical point data on systems having critical temperatures below the ice point were those of Ruhemann.1 on the methane-

Jan 1, 1952

By K. H. Kilgren

At high formation pressures the distillate produced from a gas-condensate reservoir may be black in color. In this event the dense gas phase existing above the dew point is correspondingly dark. Volumetric phase data and an analysis of a reservoir fluid exhibiting these characteristics, together with a description of the visual equilibrium cell in which these observations were made, are presented in this paper. INTRODUCTION Previously the author, like many others in the oil and gas industry perhaps, tacitly assumed that the expressions black or dark oil system, crude oil system and bubble-point system were synonymous. Crude oil reservoir fluids are bubble-point systems and yield a black or dark stock-tank oil of relatively low API gravity. Conversely, a clear or amber colored trap distillate of high API gravity is assumed indicative of a dew point system or a gas-condensate reservoir fluid. This broad classification appears satisfactory for shallow reservoirs, but as the following study demonstrates, may be misleading when applied to deep reservoirs. Theoretically, there is no reason to exclude the possibility of producing a dark, low-gravity distillate from a gas-condensate reservoir. At sufficiently high values of pressure and temperature, heavy, dark-colored hydrocarbons may exist in the vapor state of a multi-component system. If enough dark-colored components are present in the reservoir vapor phase, the resulting condensate will be dark. The reservoir fluid investigated in the present study supports this contention. Stock-tank production was black in color and measured 29" APT gravity. From outward appearances, the liquid closely resembled a medium gravity crude oil. Experimental measurements proved the reservoir fluid was in reality a gas-condensate system. Volumetric phase data for this high-pressure system and a description of the visual cell in which the study was conducted successfully are presented in this paper. THEORY Phase behavior of a reservoir fluid can be predicted accurately with reference to a pressure-temperature phase diagram. If the reservoir temperature is lower than the critical temperature of the hydrocarbon fluid in place, bubble-point behavior will be observed. If the reservoir temperature lies between the critical and cricondentherm temperature, dew point behavior and retrograde condensation will occur. For reservoir temperatures above the cri- condentherm, only a single gas phase can exist in the reservoir regardless of pressure. Providing the composition of the reservoir fluid were known, it would be possible to predict the critical temperature and estimate the phase behavior from equilibrium relationships. However, the usual practice is to obtain a sample of reservoir fluid, subject it to varying pressures at the reservoir temperature and observe the phase behavior experimentally. The latter method was used to obtain the data reported here. WELL AND TRAPPING INFORMATION A summary of pertinent data for the well from which the reservoir fluid was sampled is presented in Table 1. This well is located offshore Louisiana. Except for the pressure which substantially exceeds hydrostatic pressure, the information does not appear unusual. Prior to the sampling program, the well was produced for 22 hours at a stock tank oil rate of 139 B/D. Average trapping conditions and gauging data for the six-hour test period that followed are summarized in Table 2. Samples of the first-stage trap gas and liquid were obtained during the latter portion of the test period. Ambient temperature remained 5 to 10F below trap temperature and presented no problem for sampling. Surface wind and moderate foaming of stock-tank oil presented some difficulty in obtaining accurate stock-tank gauges. SAMPLE ANALYSIS Compositions of the gas and liquid samples are shown in Table 3. The trap gas was analyzed by isothermal chromatography which revealed only a trace of heptane in the stream. The trap liquid was initially analyzed by low-temperature fractional distillation, yielding a bottom product of heptane and heavier components. Specific gravity of this fraction was measured and the mol weight was determined by freezing point depression. The hexane and lighter overhead gas collected during distillation was

Jan 1, 1967