Search Documents

By R. L. Perrine

This paper shows how stability theory can be used to optimize solvent recovery of oil. Application of the theory leads to definition of the limiting conditions required for stable displacement to occur. One of these conditions is that the size of a solvent "slug" must at least equal a prescribed minimum. The criteria to be satisfied are miscibil-ity and stability. Stability implies that no viscous fingering will occur and that mixing will be caused only by a dispersion process. Miscibility implies that complete recovery will be obtained from the swept regions. Thus, for any specified set of reservoir conditions, an optimum use of solvent is defined. INTRODUCTION The early outlook for solvent flooding as a means to increase oil recovery was very favorable. Some laboratory results indicated that a "slug" containing perhaps 2 to 3 per cent of a hydrocarbon pore volume could be successful. However, other data suggested as much as 30 per cent would be required. The difference is of considerable economic importance. A theoretical explanation for these divergent results has been advanced by the author.1,2 Stability theory defines conditions for two distinct flow regimes. In stable displacement, solvent and oil become mixed by a dispersion process, and the solvent requirement is small. On the other hand, unstable displacement can degenerate into viscous fingering. The practical result is a considerable increase in the extent of mixing and, thus, in the solvent required. The present paper shows how to use stability theory to optimize solvent recovery of oil. Oil is usually considered to be displaced by a small amount, or a slug, of solvent. Gas in turn follows solvent. Any fingering will permit contact of nearly solvent-free gas and oil, leading to immiscibility and reduced recovery. Thus, our optimum is defined by two restrictions — stable flow and miscibility. Although an immiscible solvent process could prove more profitable than any miscible process, such processes lie outside the selected definition of an optimum. The usefulness of computations based on the present work will depend on the validity of the general theory. Because the important idealizations required tend to minimize the predicted solvent requirement, the results should state correctly the limiting conditions required. Therefore, theoretical calculations should have practical value. This paper first reviews several principles that result from stability theory. The following section describes use of these ideas in relatively simple systems which symmetry makes appear one-dimensional. Of particular importance is the step-by-step design of 3 stable minimum-solvent slug. The use of stability theory in cases of multi-dimensional displacement is discussed and, finally, suggestions are given for simplified practical application of the theoretical results. PRINCIPLES OF STABILITY THEORY The use of perturbation methods to derive stability conditions has been shown by the author.l The discussion is presented in terms of the fractional concentration of solvent in the single hydrocarbon phase. Quite generally, any disturbance can be represented as a Fourier series, of which only the least stable term need be considered. Denoting this term at time t = 0 by cn,

By H. J. Beach

This paper describes the properties and uses of cementitious mixtures containing 10 to 20 per cent bentonite and 2 per cent or more sodium chloride. Salt additions improve bentonite cement by increasing early strength and stabilizing viscosity without changing thickening time. The volume yield of cement is increased 20 per cent or more without lowered early strength. The unique properties of salt in this type of system make savings of more than 10 per cent possible where light-weight, low-strength cementitious material is adequate to fill behind pipe. The density range of these mixtures is 12 to 13 1b/gal and the filtration rates are medium-low, aiding the placement of long columns in primary cementing. Where available, sea water is an ideal salt source in preparing these slurries. INTRODUCTION Accelerators increase the rate of cement hydration. In speeding reaction with water, they can affect cement performance by reducing thickening time and increasing early strength after set. Accelerators that cause both effects can be considered "total", and those which produce only one effect can be called "partial". Materials that act conversely to accelerators are retarders. and they also can be total or partial. Partial retarders. make high-temperature cementing possible by prolonging thickening time while impairing early strength only slightly. Most known accelerators are total. A useful partial accelerator would increase strength without altering thickening time. This report concerns material combinations in which a normally total accelerator exerts unique partial properties. Utilizing these limited effects improves bentonite cement. SEPARATE EFFECTS OF SODIUM CHLORIDE AND BENTONITE ON CEMENT Sodium chloride and calcium chloride are commonly used in cement, and at low concentrations their effects are similar. Ludwig1 studied sodium chloride in API Class A and Class E cements.' Figs. 1 and 2 sum- marize his data taken at 140°F on slurry pumpability and early strength. Normal water-to-cement ratios are 0.46 and 0.40, respectively, for Class A and Class E cements. At saturation the salt content is 12 to 14 per cent of the weight of cement solids. The illustrations, therefore, cover the range of salt solubility in the mixing water. The data show that sodium chloride is a total accelerator of the Class E cement studied over the range of solubility. It is total in Class A cement to 6 per cent. Above this concentration salt becomes a retarder— partial from 6 to 10 per cent and total above. These results are general for different brands of Class A cement, although the shift from acceleration to retardation may occur at different salt concentrations for different temperatures. The results are specific for the Class E brand studied. In discussing his data, the investigator emphasized that different Class E cements might respond erratically to salt, depending on the brand chosen. Morgan and Dumbauld3 escribe the properties of bentonite cement in a paper that is a standard reference in using these systems. Bentonite functions in cement as an enabling agent; each bentonite addition equal to 1 per cent of the weight of dry cement makes possible the further addition of 4.5 per cent water.' Thus, 10 per cent bentonite in Class A cement, normally slurried neat with 46 per cent water, enables the slurry to contain another 45 per cent before settling. Water is one of the lightest and cheapest cement extenders; it reduces density and cost. Bentonite cements produce uniform slurries having low susceptibility to solids separation. After set, they offer low resistance to perforation.

By A. Eilers

NO one who has visited our Western mining districts, and studied the economical part of the beneficiation of the ores occurring all over that vast extent of country, can underrate the high importance of the above subject. By far the larger number of the districts which contain smelting ores, i. e., argentiferous and auriferous lead, or copper ores, are situated in the Great Basin, that great plateau between the Rocky Mountains on the east, and the Sierra Nevada on the west, almost the whole of which is comprised at present in the boundaries of Nevada, Utah, and part of Arizona. This region is essentially a barren country. The extreme dryness of the atmosphere permits but a scanty vegetation in the plains; and even in the detached mountain-chains running through it-generally from N. to S., or from N. W. to S. E.-there are no trees found except dwarfed pines and mahogany, at the head of sheltered ravines, and a few cottonwoods and willows, which fringe the insignificant streams, before the water sinks in the arid plains. Nearly all the mountain-chains in this region are rich in silver ores. That class of these ores which is adapted to amalgamation, and rich in silver, has been worked with profit for more than ten years. But before the advent of the transcontinental railways, mining was restricted to these ores alone, and the consumption f fuel could be met with the scanty supply of forest trees in the immediate vicinity of the mining districts. Since, however, the Union and Central Pacific Railroads have brought the Great Basin nearer to the commercial centres of the East and the Pacific coast, thus reducing the expenses of freight and labor materially, other silver deposits, containing poorer ores in greater abundance, have been rapidly taken

By J. E. Warren, J. H. Hartsock

An asymptotic solution to the equation that describes the flow of a slightly compressible fluid in an infinite porous medium has been used to estimate the interaction between two adjacent wells producing from a common reservoir. A direct method for approximating the interference time defined by Stevens and Thodosl has been suggested. An alternative definition for the time of interference, based on the minimum pressure change in the interference region, has been proposed; also, a direct method of determination has been prescribed. Examples have been employed to illustrate the use of both methods.

By K. A. Slagle

Hydraulic analysis of the wellbore has become increasingly inzportant for designing cementing operations and selecting equipment, materials and techniques to complenzent modern well-c-ompletion practices. Non-Newtonian fluid technology has advanced beyond the point where former empirical methods of analysis adequately define the hydraulic system and fluid properties. In view of these factors, this paper describes a series of rheological calculations which have been found practical, through field usage, for assistance in selecting a cementing program. A relatively simple laboratory method using standard viscometric equipment is suggested for determination of the rheological properties of slurries, and clrrta are presented on some of the more common cementitrg conzposition.A. A criterion for divergence from laminar-flow characteristics has been proposed. Usefulness of the calculations is indicated by examples of cementing operations where they have been used. INTRODUCTION With the changing aspects of well-completion practices during the past few years, it has been increasingly important to have a relatively simple method of analyzing the flow conditions existing in the well during cementing operations. This is particularly true in view of the improved economics toward which most of the changes have been directed. Rheological characteristics of slurries used for cementing should be a major consideration in the trend toward smaller casing sizes, either single or multiple strings. Receiving increased attention is the practice advocated in 1948 by Howard and Clark' of attaining turbulent flow with the fluids circulated during a primary cementing operation. While there may still be a difference of opinion concerning this technique, most available information indicates that superior primary-cementing results are generally obtained when high displacement rates are employed. Fluid properties of the slurry to be used must be available, as well as calculation methods, to determine what flow rates should be attained and the probable consequences in terms of frictional pressure and horsepower utilization. It would certainly be inappropriate to attempt high displacement velocities if sufficient pressure might be developed to create lost circulation. Since cementing slurries are non-Newtonian fluids, it is not possible to define their rheological or fluid properties by the single factor of viscosity and then make calculations for the quantities just described. Because the shear stress-shear rate ratio is not constant: it becomes necessary to establish at least two parameters for adequate fluid-flow calculations. It is not the purpose of this paper to delve into the mathematical development of non-Newtonian technology, nor to discuss the arbitrary classification system under which a single fluid may resemble two or three different classes depending upon experimental conditions. Rather, the intention is to present a useful series of calculations based on a concept applicable to both Newtonian fluids and to the preponderance of non-Newtonian fluids encountered in the oil-producing industry. Development of this approach was begun some 32 years ago,' and has most recently been brought to fruition by Metzner and his co-workers at the U. of Deleware. Some non-Newtonian fluids encountered in the petroleum industry, other than cementing slurries, have also had the benefit of this method of analysis."' The two parameters required to define the fluid are usually denoted by the symbols n' and K' and, for the purposes of this discussion, are called "flow behavior index" and "consistency index", respectively. These two slurry properties permit calculation of the Reynolds' number and the "critical" velocity, or the velocity at which departure from laminar flow begins. EXPERIMENTAL DETERMINATIONS The two principal instruments used for rheological studies are the pipeline (capillary-tube) viscometer and the rotational viscometer. When conveniently possible, a capillary-tube viscometer (where the pressure drop and flow rate of the material can be measured) is the better method for rigorous determination of the flow behavior index and consistency index for non-Newtonian fluids. With pressure-drop data at various flow rates, it is then possible to prepare a logarithmic plot of shear rate as the abscissa-shear stress as the ordinate. For fluids which do not exhibit time-dependency, these data will usually produce a straight line. The flow behavior index n' represents the slope of this line, while the consistency index K' becomes the intercept of this line at unity shear rate in accordance with the mathematical derivations associated with this concept of rheology. Due to the difficulties anticipated in maintaining a uniform, pumpable cement slurry for the time interval required to obtain measurements from the pipe viscometer, the n' and K' data reported herein were obtained using a direct-indicating rotational viscometer (Fig. 2). The

By M. G. Corson

The iron-silicon alloy series has always been one of the most puzzling among the binary alloys. Examining the well-known mechanical properties of the iron-rich alloys only we meet the following situation: The ultimate strength and yield point of alloys worked and annealed under the usual conditions represent nearly a straight line function of the silicon content up to 1.5 per cent. silicon (Fig. I). Starting at the ultimate strength of 35,000 lb. and at a yield point of 16,000 lb. for pure electrolytic iron, one reaches a maximurn strength of 93,000 lb. and a maximum yield point of 74,000 lb. at 4.5 per cent. silicon. This corresponds to an increase in ultimate strength by 13,000 lb, per each per cent. of silicon added and nearly the same rate applies to the yield point. The data mentioned are taken from the investigation by T. D. Yensen,1 and while they do not exactly corroborate the results obtained by earlier investigators, they seem to be the most correct. Plastic properties of the alloys are controlled by a much more complicated relationship. Up to 2..5 per cent. Si, the elongation stays close to 50 per cent. (in 2 in.) and the area reduction amounts to 85 per cent., a situation that forms the usual feature of all alloys representing strictly solid solutions. The plasticity of the latter is always nearly equal to that of their base metal. Beyond 2.5 per cent. silicon, the situation becomes much more complicated. Yensen observed an enormous drop in elongation and area reduction at about 2.65 per cent. silicon, where the first hardly reaches 10 per cent. and the second is still lower; he found, however, a considerable recovery in plasticity at about 2.8 per cent. Si and up to 1.5 per cent. Si it is represented by 20 per cent. for the elongation and 30 per cent. for area reduction. Yensen's graphs, as shown in Fig. 1, indicate also a considerable drop in strength and yield point at the same point of 2.65 per cent. Si,

By W. C. Hardy, B. W. McArthur

The purpose of this work is to show application of laboratory data in calculating solution gas-drive performance of the Cisco K-I reservoir. Included herein is a diagram showing the graphical relationship between six variables common to all solution gas-drive mechanisms. With an iterative process one may assume existence of any of these six variables at a particular pressure and calculate the production performance likely to occur in a reservoir produced by solution gas drive. Four unique pressure values are described at which the material balance equation alone is sufficient to predict the cumulative volumes of gas produced. Information derived from these pressure points establishes a criterion for determining compatibility of the pressure-volume-temperature (PVT) analysis of the produced crude with the gas production history. Two averaging procedures are used to arrive at a relative permeability-saturation relationship which describes the field-wide performance. INTRODUCTION The Cisco K-1 reservoir is a lenticular-type sand structure lying at approximately 6,300-ft subsurface. It is located in the H. and T. C. Railroad Co. Block 97, SW 1/4 S 392 and NW 1/4 S 385, Scurry County, Tex. Correlation of data from logs indicated that the Cisco K-1 reservoir occupied approximately 3,218 acre-ft of space with a closure of from 15 to 20 ft. Analysis of the production data indicates that the initial formation pressure and temperature were 2,450 psig and 121°F, respectively. The bubble point occurred at 1,715 psig after production of 12,256 STB oil and 11,520 Mcf gas. Only traces of water have been produced since initiation of production. Successive material balance calculations from June, 1952, to June, 1957, showed that no initial gas cap or water influx has contributed to production of oil and that the Cisco K-1 reservoir originally contained approximately 1.44 X 10" STB of oil. THEORY Methods for calculating solution gas-drive perform-

By J. W. Harden

THE observations here following are the result of an examination made as to the feasibility of extending the South Mountain Railroad from Pine Grove to Waynesborough, and of the commercial value of the ore through which such extension would be made to pass. The road at present built and in operation has a junction with the Cumberland Valley Railroad, half a mile east of Carlisle, and extends southward to Pine Grove, a distance of 176 miles. The extension proposed will add to it near 27 miles, making its total length 44$, more than 23 of which will run through, or in close proximity to, beds of ore. Along the route taken there are a number of openings at which the ore is gotten, of which the Craighead bank, Huntzinger's bank, Lowrey's bank, Wenkoop's bank, the Henry Clay bank, the Fuller

By Andrew S. McCraeath

The treatment which a steel receives, and the uses to which it may be applied, are frequently determined by the percentage of carbon which it contains; and especially is this the case in the different steel works in this country. The part which this element plays in the composition of steel is such an important one that a process for its accurate determination becomes an absolute necessity.

By E. A. Ernst, N. F. Carpenter

The benefits derived from an acidizing treatment are a function of the penetration achieved by the acid before complete spending. Additional penetration may be achieved by (1) controlling acid leak-08 into formation pores in the channel faces, and (2) retarding the reaction rate of the acid. A recently developed chemical additive consists of a synthetic polymeric material which absorbs hydrochloric-acid solutions, when suspended therein, swelling up to 40 times its original volume. These swollen particles have the ability to deform and seal-08 formation pores, providing fluid-loss control. In addition, they provide a diffusion barrier between the fracture face and the acid solution, prolonging the spending time of the acid. Field applications of this new technique have shown promising results. A method of conducting acid fluid-loss tests, using carbonate cores, is believed to provide fluid-loss data that are more representative of formation conditions than the conventional filter-paper determinations. INTRODUCTION The concept of oilwell acidizing has changed since its first commercial application, 30 years ago. Originally, it was visualized that the acid penetrated thousands of tiny pores and flow channels in the matrix rock, enlarging them by dissolving the carbonate walls. The resultant permeability increase was assumed to be the responsible factor in increasing production from the well. Recent laboratory studies,' however, have shown that this does not provide the complete picture. Although this type of individual pore penetration by the acid does take place during acid "soaks", designed to overcome "skin effect" due to mud invasion in the immediate vicinity of the wellbore, many years of experience have shown that considerable pressure is required to attain any appreciable injection rate into the fine capillary pores of the rock. During most acidizing treatments, the bottom-hole pressure build-up due to the restriction of flow into the formation exceeds the "breakdown" pressure of the rock so that a fracture is induced. In most cases, such fractures open up along natural, incipient fissures and zones of weakness in the rock and, therefore, tend to follow the natural stress pattern of the rock—whether it be horizontal, vertical or inclined. Because of the comparatively greater permeability of the channel in relation to that of the matrix, the bulk of the acid volume is diverted into the newly opened fracture. Here it quickly penetrates the formation, opening and ex- tending the fracture in much the same manner as a conventional fracturing fluid. Unlike the fracturing fluid, however, most acidizing solutions contain no propping agent; thus, the open fracture will again close when the injection pressure is relieved. Laboratory studies2 have shown that in many cases the etching of the fracture faces, resulting from the reaction between the acidizing solution and the carbonate rock, is nonuniform due to the heterogeneity of the rock structure. As a result, the two fracture faces no longer match when pressure is released, and support pillars and intermediate voids remain, forming a high-conductivity channel for well fluids. Unfortunately, this is not true over the entire area of the fracture, but only over that portion of the fracture where the rock has been partially dissolved by the acid. The acid solution spends as its travels away from the wellbore; once it has completely spent, even though it may provide additional mechanical fracture extension, no additional benefit due to etching of fracture faces can be expected. Studies of acid reaction rates under formation conditions,3 observing the effect of different variables upon spending time, have shown that the reaction was often so rapid that very little penetration of the formation occurred before the acid was spent. Study was undertaken to devise methods of increasing the penetration of the acid before spending, so as to provide greater benefit from the acidizing treatment by etching a greater portion of the fracture faces. Several techniques were devised to accomplish this purpose. First, chemical additives were developed which were designed to retard the reaction rate of the acid, causing it to penetrate a greater distance from the wellbore before finally becoming spent. Another method was to increase the injection rate of the acid. However, it was found that the resultant increased shear tended to accelerate the reaction rate of the acid, partially offsetting the benefits of the higher injection rate insofar as achieving increased penetration before spending was concerned.' Another approach to the problem of achieving increased penetration was the development of fluid-loss additives for acid solutions, which would minimize the volume of acid lost into formation pores in the fracture faces and provide maximum fracture extension for the volume of acid injected during the treatment. The use of fluid-loss additives is now considered the most effective method of providing maximum fracturing-fluid efficiency.~ Unfortunately, this latter technique does not solve the problem of rapid reaction rate, with consequent limitation of the fracture area benefited by reaction with unspent acid. A newly developed acid additive overcomes many of these limitations by providing the dual benefits of fluid-loss control and mechanical retardation of acid reaction

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 G. W. Binkley, R. E. Collins, G. K. Dumbauld

A mathematical analysis has shown that the primary factors affecting the deposition of cement in unfractured perforations during squeeze-cementing operations are: the properties of the cement slurry, the geometry of the perforations, and the filtration time. The validity of the analysis has been substantiated by laboratory investigations. The information presented in this paper makes possible the prediction of the rate of deposition of cement in and about unfractured perforations under given squeeze-cementing conditions, the proper filtration time to control the amount of deposition to any desired value, and the deposition properties needed in cement slurries to restrict to any desired degree the formation of nodes of cement .Solids inside the casing. INTRODUCTION The deposition of solids against a permeable medium by the process of filtration is an important property of a cement slurry. General recognition of the role of filtration in the use of cement in wells is evident in view of the fact that several patents on low-water-loss cements have been issued during the past few years. In the conventional, high-pressure squeeze-cementing technique, a large volume of cement slurry is usually displaced into the area surrounding the perforated interval after the formation has been fractured hydraulically.&apos; Near the end of the job, the pumping rate is slowed or stopped intermittently in an effort to obtain a high pressure in the wellbore. The way in which a high pressure is obtained is not understood thoroughly. but it is be- lieved to be caused by the deposition of a rigid mass of compacted cement particles either in the fractures of the formations or inside the casing opposite the perforated interval. In a method of squeeze-cementing introduced several years ago, a small amount of cement slurry is squeezed into the perforations at a low pressure to avoid fracturing the formation. The filling of the unfractured perforations with cement particles by the process of filtration is important to the success of this technique. Early in the development of this low-pressure squeeze-cementing technique, laboratory and field-size experiments showed that conventional, neat cement slurries were unsatisfactory if the tubing was extended into or through the perforated section of the casing. The extremely high filtration rate of neat cement slurries caused the deposition of a compacted mass of cement inside the casing and made it impossible to remove excess cement by reverse circulation and difficult to withdraw the tubing. Because the need for cements with a reduced filtration rate was obvioils and urgent, modified cements containing 12 per cent bentonite were used until research was completed on modified cements containing 25 per cent bentonite.< These modified cements were useful in the solution of the immediate problem, but the following questions remained unanswered: (I) what is the optimum filtration rate of cement slurries for this application, (2) how long should the squeeze pressure be applied. and (3) how great are the effects of squeeze pressure and well temperature on the rate of deposition? Furthermore, the relationship of the filtration properties of cement slurries to the dimensions of the perforations and the wellbore was unknown. In order to obtain information relative to these basic considerations, a systematic evaluation of all the factors involved in the deposition of Cement during low-pressure squeeze-cementing operations was initiated. Although the work was directed toward an evaluation of the factors influencing the deposition of Cement in the low-pressure squeeze-

By B. G. Matson, M. A. Whitfield, G. R. Dysart

This paper describes the development of u computer program to calculate changes that occur in the length of tubular goods due to temperature and pressure changes during stimulation operations. Due to the numerous variables involved and the uncertainty of all static and dynamic conditions that could exist, it becomes a staggering task for individuals charged with completions to perform the necessary mathematical calculations. The computer program permits advance calculations for several sets of conditions. INTRODUCTION In the Delaware basin of West Texas alone, 50 wells were contracted or drilled to 15,000 ft or deeper in 1965. Deep well activity is continuing in this and other areas on an expanding scale. Many of these deep wells require extensive stimulation for successful commercial production, and during these operations, pressures and temperatures are encountered that have a pronounced effect on the length of tubular goods. This length change during a large-volume, high-pressure stimulation treatment utilizing fluids considerably cooler than bottom-hole temperature can be of such a magnitude that permanent damage to casing and tubing will result unless mechanical design, pressures and fluid temperatures are evaluated and controlled. These pressure and temperature effects can be calculated. However, the process lends itself well to computer solutions because of the mathematical nature of the problem and the calculating hours involved in arriving at an answer. The engineering-hour demand becomes more severe as tapered strings are involved. On initial treatments on a given well, surface pressure and injection rate conditions are unknown, and offset well conditions have not proven to be a reliable method for making predictions. For these reasons, it has become rather standard procedure for operators to compensate for these uncertainties by placing unnecessary pressure and fluid temperature restrictions on stimulation design. On a number of occasions treating fluids have been preheated to as much as 160F as a means of compensating for thermal contmction resulting from pumping cool fluids. The maintenance of packer seals has been treated by Lubinski, Althouse and Logan&apos;,&apos; and the problem of therma1 effects on pipe has been explored by Ramey." These works were expanded and the results made applicable to everyday oilfield terminology before submitting them to computer programming. The pressure and temperature effects on tubing movement previously mentioned occur simultaneously as fluid moves through the pipe. The pressure changes, for purposes of explanation, are categorized here as to the various effects these pressures have on a tubing string. These divisions are (1) the piston-like results of forces acting on horizontal surfaces exposed to pressure, (2) swelling or ballooning of the tubing along its entire length due to the forces of pressure acting against the tubing walls, (3) the elongation of tubing due to frictional drag and (4) corkscrewing of the pipe due to internal pressure. Thermal changes are also of great importance, as their results may be more significant than any of the pressure effects. Steel is an excellent conductor of heat and the earth is a relatively poor conductor. It has been calculated that pipe temperatures at depths of more than 20,000 ft approach within as little as 25" the temperfature of the surface fluid after pumping for 2 hours, or a drop in temperature in some treatments of more than 220F. The equations presented in this paper were developed for computer programming and simplicity of input information; therefore, numerical constants such as Young&apos;s modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson&apos;s ratio for steel (0.3) are included with unit conversion factors. The moment of inertia of tubing cross-sectional area with respect to its diameter was changed to a constant times (D&apos; — d&apos;) where D is outer diameter and d is inner diameter. Units in the equations are length in feet, diameter in inches, density in pounds per gallon, pressure in psi, rate in barrels per minute and time in hours. PISTON-LIKE REACTIONS A change in tubing internal dimensions and the exposure of other horizontal surfaces to different pressures on the inside and outside of the tubing result in a reaction much like a piston under pressure. Such is the case when the internal diameter changes in a combination string of pipe, when seals of a slick joint assembly are subject to pressure and in the end effects of a tubing string. The change in tubing length due to the piston effects of a slick joint packer is affected by the various diameters involved, the tubing pressure Ap,, the casing pressure ,Ap,, length of pipe L, densities of fluid in the tubing before and during pump-

By T. M. Williams

DURING the year just ended we have had three great fires in the mines in the Wilkes-Barre district. One at the Empire Colliery, one at the Prospect shaft, and the other at the Baltimore old mine. It is proposed in this paper to give a brief history of the fire at the Empire Colliery, and the methods employed in extinguishing it. The fire was first discovered breaking out in the air-stack, alongside what is known as the " Kidder Plane," about one o&apos;clock on the morning of December 31st, 1873. Subsequently it was ascertained that the fire originated in the return airway somewhere near the steam boiler, located inside the mine, near the head of No. 5 slope. The coal from this slope was hoisted to the shaft gangway, at a point nearly four thousand feet distant from the foot of the shaft. As there were from four to five hundred mine cars hoisted per day from this slope the steam boilers were located at this place. The mines being of considerable extent required a large volume of pure air for ventilation, and no fan having been erected at this slope at the time, the heat from the boiler was made use of to assist ventilation, by conveying the return air from No. 5 slope into the large airway made from these boilers to the surface through the old "Kidder slope" to the stack where the fire was first discovered. The length of this slope or airway was about one thousand feet, with a sectional area of about sixty feet. The quantity of air passing through it must have been over sixty thousand cubic feet per minute