Search Documents

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',' 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's modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson'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' — d') 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 H. A. Lesser, H. L. Stone, G. H. Bruce

ESSO PRODUCTION RESEARCH CO., HOUSTON, TEX. A mathematical model that represents the conduction heating of a rock formation of limited permeability is formulated. Heat is introduced by the injection of a hot condensing gas into horizontal fractures through the formation. The fractures are assumed to be equally spaced; the flow of fluid is linear. The model consists of three coupled partial differential equations and two functional relationships. Temperature histories for both fractures and formation are obtained by numerical procedures. Application of the model in predicting in situ heating of oil shale is illustrated with four sample cases. Influences of shale thermal diffusivity, fluid pressure in the fracture and fracture spacing upon heating rate are demonstrated. In all cases considered, several years were required to heat all the formation to 700F with 1,000F steam. Less than 32 per cent of the injected heat was utilized for heating the shale. To facilitate extensive calculations of this type, a set of seven dimension-less quantities are derived. With this dimensionless system, the effect of controllable parameters on economics can be explored. INTRODUCTION Vast deposits of hydrocarbons are contained in oil shale formations in the western United States and in other countries. These hydrocarbons, known as kerogen, are difficult to recover because they are immobile solids and because the rock formations in which they are found have very limited permeability. The kerogen can be converted through pyrolysis to recoverable liquids and gases if the formation temperature is elevated to the range 600 to 800F. To accomplish this heating, the formation might be artificially fractured and a hot fluid injected into the fractures. If the injection and production wells are vertical and the fractures horizontal, interwell communication will be established and the hot fluid will flow from an injection well to a production well. In this process, some of the heat in the fluid will be given up and conducted into the formation. Products of pyrolysis will flow from the heated rock matrix into the fractures and be produced, commingled with the injected fluid. The practicality of such a process is strongly influenced by the rate at which heat can be transferred from the hot fluids to the formation. An analysis of this heat conduction problem is the subject of this paper. Linear flow of a condensing gas through a set of horizontal, equally spaced fractures is considered. Temperature histories for both fractures and formation are obtained from numerical procedures. The analysis also applies to formations containing streaks of natural permeability instead of fractures. Heat flow problems similar to the one treated here have been studied by other workers. Thomas1 has solved a radial problem in which he considers the injection of a noncondensing gas into a single, horizontal fracture and approximates the temperature profile in the fracture by a step function. The solutions to other problems, which are less similar to the present problem than the one treated by Thomas, have been reviewed by Spillette.2 In addition to these mathematical studies, a field test of an in situ shale oil recovery process using hot natural gas as the injected fluid has been reported.3 The first part of this paper is concerned with the mathematical development of the problem. The differential equations, boundary conditions and functional relationships which describe the physical model are formulated and the equations are cast into finite difference form. The method of solution of the difference equations is discussed (Appendix B). The second part of this paper presents results computed for a typical oil shale formation for which steam at 1,000F is the injected fluid. Fracture spacing, horizontal and vertical thermal diffusivities

Jan 1, 1967

By G. W. White

Over the past several years, the steady progress of gas lift from a mysterious art into a developing science has been traceable through the many fine papers that have appeared in Journal of Petroleum Technology. Consequently, this magazine has rightly become a standard for the serious students of gas lift all over the world. The following discussion is offered in a sincere attempt to accurately present certain topics covered in the subject paper in the light of acknowledged state of the art procedures. INSTALLATION DESIGN The first point to be questioned is the author's installation design, beginning with Eq. 6. The use of 1,200 psi in the spacing equation instead of the 1,275 psi noted as available seems, at first glance, to be a simple arithmetic error. Upon further examination, it appears evident that both this 75 psi plus the 100-psi safety factor must be subtracted from available kickoff pressure because of conditions imposed by the characteristics of the ASC valve (see Fig. B-1). Pressure-charged valve design is incorrectly penalized by spacing the first valve using only 1,100 psi. The use of such a high back-pressure makes it appear that a typical ASC installation would require no more valves than a comparable conventional valve installation. Using a value of 100-psi tree pressure, the well could be valved to 11,000 ft very conservatively with seven pressure-operated valves, but would still require 11 ASC valves. The use of 400-B/D spacing gradient for a 140-B/D rate is not an accepted practice using precharged valves. The author would have done well to note the advantages to be gained by reducing the 500-psi tree pressure, which is clearly explained by Beadle, Harlan and Brown.' In Step 1 of Part III of the continuous flow discussion, the author incorrectly applies the Poettmann and Carpenter curves.2 Item 2 is an incorrect statement. An excellent reference on this point is the paper by Baxen-dell and Thomas.' VALVE MECHANICS The next topic relates to the mechanics of precharged valves. Eq. 2 for the ASC valve can be used to describe the condition of equilibrium, wherein the valve is closed with opening pending. When this equation is used to determine relative values of casing, dome and tubing pressure with conditions approaching the transition value, some strange results are obtained. The following uses the author's symbols and values given on page 689. Fig. B-1 shows the condition that exists at the moment an ASC valve has been uncovered during the unloading process (Condition I). The obvious result is a locked condition. The basis for the curve is shown on the figure which includes 1,000 psig casing pressure to uncover the valve. If injection gas is not bled off the casing until the pressure is decreased by an amount exceeding the stabilizer discharging valve pressure, upward adjustment of the casing pressure is mandatory for the valve to open and unloading operations to continue. In Condition I the dome pressure exceeds the casing pressure by 20 psi, and the tubing pressure required to open the valve is 1,212 psig. This tubing pressure would not exist under the assumed conditions. In Condition I1 the casing pressure is increased to the dome pressure, and the tubing pressure required to open the valve is 1,102 psig. It is apparent that the casing pressure must be increased to a value exceeding the dome pressure before unloading operations could continue. In Condition I11 a maximum casing pressure of 100 psi (stabilizer charging valve set pressure) above the dome pressure is reached, and the tubing pressure required to open the valve is 549 psig. The author's statement after Eq. 5 that "the minimum tubing pressure required to open an ASC valve is equal to the casing pressure at depth less the transition pressure of the valve" is incorrect except for Condition I11 only. Otherwise, the equation is invalid.

Jan 1, 1965

By M. Kaitz

A continuing discussion in both the petroleum engineering and economic literature is directed to the difficulties encountered in the use of discounted cash flow rate of return (DCF) as a measure of investment worth. Although useful in most Instances, DCF has been criticized because it is time-consuming in its trial-and-error solution, theoretically invalid, not adaptable to cash flow streams yielding multiple return solutions and not entirely reliable in selecting between mutually exclusive investments. The theoretical invalidity of DCF stems from the reinvestment assumption implicit in the calculation that earnings are reinvested at the DCF rate. The emerging consensus of the economic literature is that the net present value or the present worth of the net cash flow stream (discounting at the average opportunity or cost of capital rate) is more correct and reliable. Other criteria proposed have been ratios of net present value divided by initial investment or by present value of all investments in a project. All of these criteria are simple to determine and explicitly assume a reinvestment rate for [he income generated by a project. This paper develops and discusses a theoretically valid profitability criterion which is simple to compute and retains the appeal of a percent return on investment. It is called "percentage gain on investment" or PGI. It measures the gain an investment is expected to realize over like capital invested in the average opportunity and explicitly considers reinvestment potential. Why add another Concept to the large array of investment criteria now available, any one of which, or perhaps a combination of several, appears to embrace the form's (or individuals) objectives? The answer is that not one of the existing criteria provides both a readily comprehensible and theoretically valid measure of risk coverage that has general application. The proposed PGI does fulfill these requirements. INTRODUCTION An ancient expression warns that "one must yield to the times" — there are better ways of doing things. A review of the petroleum engineering and economic literature on one topic alone, measurement of investment worth, certainly is witness to this truth. In use for a number of years, the DCF has recently received attention, directed mainly to its theoretical invalidity. Several alternatives to DCF have been proposed to provide a valid, simply determined criterion to describe investment worth and to overcome the criticisms previously mentioned. This paper introduces another method called percentage gain on investment (PGI) and is proposed as but one of several yardsticks that should be used in making investment decisions. MEASURES OF INVESTMENT WORTH This paper will consider only those criteria which give weight to the time pattern of future earnings. These criteria are usually compared with an average opportunity rate or cost of capital of the firm to judge the relative worth of the investment. For purposes of demonstration, a 9 percent average opportunity rate will be used throughout this paper. Implicit in the DCF calculation is the assumption that earnings are reinvested at the DCF rate. Some argue, though, that there is no reinvestment assumption,' that the DCF rate is simply that maximum rate of interest one can pay on the investment over the life of the project and break even. The determination of DCF is accomplished by discounting the net cash flow stream at that rate (DCF) which will yield zero. The question is: why should reinvestment potential be explicitly considered in calculating return on investment or other criteria measuring economic worth? Perhaps the answer lies in consistent or equal treatment of future cash flow. It appears entirely illogical to give different present worth value to $1 received, say, 10 years from now, which is the circumstance resulting in comparing projects with different DCF returns. In fact, $1 received in 10 years has the same value regardless of which project generated the income. DCF return thus favors investment projects which are expected to provide early income as compared to those providing long-term income. Not surprisingly, the controversy on reinvestment assumption is an old one. Hoskold' discussed this same problem in 1877. He considered the future income from mineral properties as an annuity or a series of fixed future payments. (As will be demonstrated, his equation can be modified for variable income.) Prior to Hoskold, the value or what one could pay for the mine, was determined with the use of standard, single-interest tables. Here is what Hoskold said with regard to these tables. "This table, and others of its kind, to be found in most works on annuities, is constructed correctly according to

By J. L. Huitt, B. B. McGlothlin

The introduction of new fracture propping agents that are brittle but much stronger than sand created the problem of what loading strength is required for a propping agent to be effective in a given formation. It is shown that the load at which the propping agent crushes should exceed the load at which total embedment in the fracture faces is possible. Simple laboratory tests to determine loading strength of the propping agent and embedment in the fracture faces, and use of these data in selecting a propping agent for a given formation, are discussed. INTRODUCTION One of the most important factors in the design of hydraulic fracturing treatments is the selection of a propping agent that can effectively provide the fracture flow capacity needed for stimulation of a well. Sand, once generally accepted as being synonymous with propping agent in hydraulic fracturing, is now recognized as having limited effectiveness in many formations because of its low resistance to crushing. Sand particles are brittle and have relatively low strength. Because of this property, sand particles are crushed in rocks that offer high resistance to the penetration of fracture faces by the proppant particles when the fracture attempts to close under the action of the overburden load. For rocks that offer a high resistance to penetration, deform able particles are more effective propping agents than sand. However, for this same type of rock, a propping agent that does not deform, yet does not crush, is often more effective. Thus, a rigid propping agent with sufficient strength to prevent crushing is desirable. A method for determining the strength required for a rigid propping agent to function effectively in given formations is discussed. BEHAVIOR OF RIGID PROPPANTS AND FRACTURE FACES RELATED STUDIES An early qualitative description of the reaction of propping sand in fractures was given by Hassebroek et al.' In discussing fracturing in deep wells, the authors mentioned that even though propping sand entered the fractures, a high flow capacity did not result due to crushing or embedding of the propping sand. Dehlinger et al.2 in discussing the reaction of propping sand surmised that, because of the hardness of sand particles, deformation occurred in the fracture faces contacting the propping sand. In later studies,3,4 methods of determining the embedment of propping sand in fracture faces of soft rock and the critical load at which propping sand is crushed by the fracture faces in hard rock were discussed. In working with de-formable proppants, Kern et al. considered proppant particles to be deformed into cylindrical disks by action of the overburden and then pressed slightly into the fracture faces by further action of the overburden. Rixie et al.'0 reported on embedment pressure and presented a method of selecting a propping agent for use in given formations. The propping agents included sand, walnut shells and aluminum pellets. All these studies have contributed materially to a better understanding of propping agent behavior; however, the strength of brittle proppants (sand, glass and ceramics) required to result in embedment rather than crushing has not been discussed. This topic will be covered in the ensuing discussion. PROPPANT PARTICLE CRUSHING—-EMBEDMENT For this discussion, a rigid propping agent is considered to be one that is brittle and fails under tensile stress when loaded to a critical value. In an earlier study4 it was shown that the Hertzian4 loading theory could be applied to a spherical brittle propping agent if the propping agent and fracture faces behaved elastically. At the failure of the proppant, the ratio of the load to the square of the diameter of the particle should be constant for a given material combination, or: Lc/dp2=C ............(1) A partial derivation of this equation from proppant and formation properties is included in the Appendix. Should a rigid particle not be crushed as a load is applied, it embeds in the fracture faces. A study3 of particle embedment in fracture surfaces has been published. The embedment can be described by an equation based on Meyer's metal penetration hardness relationships: d1/dp=B 1/2[L/dp2]m/2..........(2) In Eq. 2, B and m are constants that are characteristic of the rock; the significance of the other terms is shown in Fig. 1. A STANDARD DEFINITION FOR PROPPANT LOADING STRENGTH Eq. 1 is useful in appraising propping agent strength," but it is strictly applicable only when the area of contact between a particle and a fracture face (or loading plate)

Jan 1, 1967

By L. R. Raymond, J. L. Hudson

This paper proposes a comparatively short-term (8 to 10 hours) well test for detecting and characterizing well-bore damage and for measuring mean formation permeability. The proposed test is made by injecting fluid at constant pressure, recording injection rate as a function of injection time. After one to four hours of injection, the well is shut in and fall-off of bottom-hole pressure is obtained as a function of shut-in time. Formation permeability is estimated by an iterative technique. First, a value of formation permeability is assumed. Then, a plot of the recorded injection rate as a function of dimensionless time is made, using the assumed pertneability value. From the slope of the injection-rate curve. a new value of formation permeability is calculated. If the new value agrees with the original assumed value, the assumption was the correct formation permeability. If the values do not agree, the process is repeated using the new permeability value in the calculation. Convergence is rapid, and a reliable permeability value results. Pressure fall-off data are used to check the result. Graphs of pressure and injection rate us functions of time given in the paper show that changes in permeability of the formation in the neighborhood of the wellbore are disclosed by this technique. Thus, the short-term test can he used to detect formation damage. Also, a rough measure of the radial extent of damage can be inferred, which is helpful in designing stimulation treatments. The mathematical model used for this work was a single-zone, horizontal reservoir with a damaged zone in which permeability decreased continuously as radial distance to the wellbore decreased. This model is more realistic than the usual two-zone, discontinuous permeability model used in published works; calculations indicate the realistic model is valid. Vertical variations in horizontal permeability were studied with this model, and results indicate that the permeability measured by the short-term test is the mean horizontal permeability for the vertical interval tested. The proposed short-term test thus should be useful in detecting and characterizing formation damage and in measuring formation permeability needed in calculating reservoir transmissibility. INTRODUCTION To plan the most efficient production or injection schedule for a well and to design or evaluate the optimal stimu- lation treatment, it is necessary to know the properties of the reservoir adjacent to the well, particularly the reservoir transmissibility and characteristics of a damaged zone, if one exists. Several techniques for determining reservoir transmissi-bility from well tests have been presented in the literature. 1,2,3,4 All these techniques rely on conducting constant-rate well tests that often are difficult to execute. A constant-pressure well test is generally easier to carry Out. and this paper contains the first available method for the analysis of constant-pressure well tests. Determination of wellbore damage from transient well tests has been the subject of several papers."" From these studies it is apparent that information necessary for determination of the characteristics of a damaged zone is available shortly after the transient well test is initiated. Consequently, it may not be necessary to carry out an extensive well test (for example, a pressure build-up test) if the primary purpose of the test is to detect the existence of wellbore damage. All previous studies of well testing to determine wellbore damage have been based on the two-zone perrneability model. In this model the damaged zone has a permeability k,, extending to a radius r,,, and the formation permeability k obtains from r, to the drainage radius r,.. Consequently, there is a discontinuity in permeability at r = r,,. This discontinuity can be eliminated by assuming a continuous variation in permeability through the damaged zone. The effect of this assumption on transient well tests is discussed in following sections of this paper. In addition, all formations have within them vertical permeability variations associated with lithology changes throughout the zone of interest. This paper also considers the effect of these variations on transient well tests. ANALYSIS OF CONSTANT-PRESSURE WELL TESTS The mathematical analysis associated with the injection of fluid at constant wellbore pressure into a single-zone, horizontal reservoir completely filled with a fluid of small and constant compressibility and constant viscosity is given in Appendix A. In this analysis it is assumed that the well is located at the center of an undamaged, circular drainage area. From this analysis, the formation permeability can be obtained as follows. 1, Estimate a value for the formation permeability k. 2. Prepare a plot of injection rate q vs

Jan 1, 1967

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner' makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,' is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors'' derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel' terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones'" initiates steady-state

Jan 1, 1967

By C. T. Philipp

Recycling of waste is much preferred over disposal. This paper discusses the Enviroscience, Inc. (ESI) process, now being commercialized by MetWool, Inc., to recycle hazardous wastes using an extractive metallurgical process which is followed by vitrification of non-reducible metal oxides. In this process, volatile and reducible metals are reclaimed as their respective oxides and elements.

Jan 1, 1995

By James K. Gowans

Polaris Mine, a Cominco Ltd. property located on Little Cornwallis Island, Northwest Territories, Canada, (Lat. 75°30'), is the world's most northerly base metal mine. Polaris Mine, which commenced operations in late 1981, is characterized by its unique process barge. The lead/zinc concentrator, powerhouse, shops, offices, dry and warehouse have all been incorporated on a 400 ft. long by 100 ft wide by 65 ft high barge. This process barge was constructed in Southern Canada, at Trois Rivieres, Quebec, then towed 3500 miles into its final resting place on the high arctic. This paper describes the unique design of the barge concentrator and deals with the start-up and continuing operation of the Polaris Concentrator.

Jan 1, 1984

By E. H. Bentzen

The purpose of this paper is to explain the production of sand for glass manufacturing and the reasoning behind the sand specifications set by the glass industry. In order to discuss the production of glass sand, there is a need to understand what glass is and how it is made. Therefore this paper discusses glass, as well as sand beneficiation.

Jan 1, 1977

By J. M. Bray, R. R. H. Cohen, P. Kiousis, S. R. Al-Abed

"INTRODUCTION Conventional tailing storage and management practices typically utilize engineered facilities, such as ring dikes and valley impoundments. To reduce the environmental costs associated with disposing of mine wastes in tailing storage facilities (TSFs), several beneficial re-use strategies have been developed over the years. For example, thickened tailings have been used for paste backfill. It is common for tails to comprise up to 60% of the material’s total weight and it has been estimated that TSF volumes can be reduced by 60% when using this strategy, (Benzaazoua, et al. 2004). Other beneficial re-use strategies have produced unpaved temporary access roads, hot and cold asphalt, chip seals, self-compacting concrete, and golf course sand traps. This project sampled mine tailings from two different TSFs in Colorado and produced structural concrete that could potentially be used for commercial foundations, sidewalks, and paved roads. The tailings were substituted for the fine aggregates typically found in conventional concrete. The two main objectives of this study were to investigate the structural integrity of the produced concrete samples and the leaching of heavy metals from the resulting concrete. The ratio of tailings to fine aggregate (T/FA) was carefully studied to find if there was a correlation between the substitution rate and the structural integrity. The potential maximum substitution rate was also investigated for the produced concrete specimens. The potential of leaching mining impacted water (MIW) was evaluated by comparing the leachability of the raw mine tailings to the encapsulated tailings by using standard shake extraction prediction tests. MATERIAL AND METHODS Materials The four components that make up the concrete mix design for this project are: cement, coarse aggregate (gravel), fine aggregate (sand and mine tailings), and water. The type I/II ordinary Portland cement (OPC), gravel, and sand came from local cement and concrete companies in Denver, Colorado. Mine tailings were collected from two different TSFs in the Rocky Mountain Region, where Pride of the West TSF in Silverton, CO will be denoted as Site A. The primary metals of interest at the mill site are gold, silver, lead, copper, and zinc. The location of Site B will not be mentioned in this report due to confidentiality agreements with the mining company"

Jan 1, 2016

By Roger Potvin

Abstract The demand for pure zirconium and titanium metals in Germany increased considerably during the war. Pure zirconium was used in flash bulbs, in vacuum tubes, and in time fuses for bombs, where it has a decided advantage because it gives no gas in burning. Under the form of alloy, it was used with other rare elements for making substitutes for platinum; it was also used as a scavenging agent in high-grade steels. Metallic titanium was mainly used in connection with the manufacturing of special alloys. The authors were sent to Germany by the Department of Reconstruction and have investigated the German process developed for producing these rare metals. They describe the procedure and special equipment used for their commercial production in Germany during the war. General Metallic zirconium and titanium were .produced on a commercial scale in Germany by the thermal reduction, with metallic sodium, of potassium fluozirconate and titanium chloride, respectively, using a mixed potassium and sodium chloride flux. According to the information gathered, the sole German producer of these rare metals in pure form was the Deutsche Gold und Silber Scheideanstalt (commonly called Degussa), at their Werk II, 215 Gutleutstrasse, Frankfurt/Main, which was visited by the authors on September 18th, 1945. As far back as 1934, Degussa was carrying on active research on the manufacture of zirconium, and nearly one ton of the pure metal had been produced before September, 1939. Concurrently, a quite similar process was worked out for the recovery of titanium and it was in commercial operation in the year 1938, during which the production totalled about half a ton (metric). The demand, especially for pure zirconium, increased considerably during the war, attaining a total of about 48 tons during the period October, 1939, to the end of 1944, and reaching a peak of more than 15 tons in 1943. In 1944, the plant was slowed down due to the frequent and disastrous bombing by the Allies, especially in March and September.

Jan 1, 1946

By R. L. Whiting, H. J. Ramey

The basic study from which this paper was prepared was started as the result of the growing need throughout the world for increasing quantities of energy in all forms. Quite obviously, natural forms of energy that are readily available at low development cost are those in greatest demand. The underdeveloped countries — and particularly those having little or no petroleum resources — are the countries in which the most interest is being shown in the newer energy sources. One of the least expensive energy sources is natural geothermal steam.'-* Although this form of energy has been recognized for centuries, it has been only during the past 20 years that serious efforts have been made to harness it. Natural geothermal steam energy in Italy, New Zealand, Mexico, Japan, and California is now being produced through wells to drive turbines and generate electricity. Further, active exploration for natural geothermal steam is being conducted in Hawaii, Fiji, Taiwan, Chile, Russia, Greece, and Katanga. It is surprising to find that most geothermal steam exploration is in the "steam seep" stage. That is, surface studies are made and exploratory wells are drilled in the general area of steam seeps. However, in the larger geothermal steam areas, there have been efforts to apply the most modern geological and reservoir engineering principles in order to define the reservoir parameters, particularly those relating to estimates of reserves and future productivity. This paper is concerned with the development of appropriate equations and techniques to facilitate these estimates. Production of natural steam or hot water presents problems different from those experienced in the production of oil and gas. For example, steam or hot water systems may be essentially single-component systems, while hydrocarbon systems are most frequently multicomponent fluid systems. Heat effects are much larger for water than for hydrocarbon systems; and the natural steam production may or may not be isothermal, while production of petroleum reservoirs is considered normally to be isothermal. Petroleum reservoir engineering principles may be applied to natural steam or hot water reservoirs if the inherent differences in the systems are considered. The basic considerations involved in geothermal steam reservoir engineering are: thermodynamics, physical and thermal properties of water, materials and energy balances, fluid influx, and performance matching and predicting. Thermodynamics Fig. 1 is a pressure-temperature diagram for the liquid-vapor region for pure water, showing the critical point and five other points representing possible initial conditions for a geothermal steam reservoir. To generate a physical appreciation for the mechanics of natural steam production, it is helpful to consider

Jan 1, 1970

By J. E. Strassner

Oilfield emulsions are stabilized primarily by film-forming asphaltenes and resins containing organic acids and bases. Adding inorganic acids and bases radically changes the physical properties of these films and their emulsion-stabilizing properties. Interfacial films show one of three different physical characteristics when retracted in pendent drop tests. These compressed films have rigid, soap-type or minimum film properties depending on the quantity of asphaltenes and resins in the crude oil and on the pH and composition of the brine. There exists for each crude oil-brine system an optimum pH range over which the adsorbed film exhibits minimum contracted filtn properties. In this range, intmfacial tension is high — frequently near its maximum value — indicating the absence of highly surface-active specier; crude oil-brine emulsions generally show minimum stability; and surfactant requirements for breaking these emulsions are significantly reduced — sometimes no surfactant is needed. About two-thirds of the world's crude oil is produced as emulsions that must be treated (demulsified) to be marketed.' Approximately three-fourths of all domestic crude oil produced must be treated. The percentage of crude requiring treatment increases as fields become older and more water is produced. Demulsification most frequently is attained through use of surfactants, heat and electric treaters. It is estimated that over 75 percent of all oilfield emulsions are broken by chemicals generally containing alkylene oxide condensates. Although acids and bases are known to affect emulsion stability:-' they rarely are used in the oil field except in batch treatment of the most severe emulsions. Nevertheless, pH adjustment by using acids or bases in demulsification has some advantages over other methods. Mainly, it is cheaper. Acids and bases frequently are cheaper to use than surfactants, or can appreciably reduce the amount of surfactant required for demulsification. Secondary advantages are that less oil is lost by emulsi-fication into the aqueous phase, and that less corrosion is encountered when the optimum pH is basic. Limitations to the use of pH include (1) use in brines containing high concentrations of reactive ions, (2) closer supervision for pH control than is required for surfactants and (3) need for corrosion protection when the optimum pH is sufficiently acidic. These factors indicate that the economics of pH con trol in demulsification will be maximized when relatively large treating facilities are needed, when relatively large surfactant concentrations are required to break an emulsion, when a relatively low concentration of the aqueous phase is emulsified in the oil and when the emulsified aqueous phase contains relatively low concentrations of reactive ions. Because there are no broad guide lines for using pH control in breaking oillield emulsions (and because the limited guides available are highly specific and frequently contradictory), this paper will attempt to clarify the principles involved in using acids and bases to break emulsions. These principles also should be useful in refinery separations and in other oilfield applications where stable or unstable petroleum emulsions are involved. Principles of Emulsion Stability Oilfield emulsions generally occur as water droplets dispersed in oil (w/o) and occasionally as oil droplets in water (o/w). These emulsion droplets are usually in the 0.1- to 10-micron range. However, more coarsely dispersed water droplets frequently are present in the emulsion. Pure hydrocarbons and water will not form stable emulsii except in extreme dilutions where the large distance brtween droplets inhibits coalescence. An emulsion breaks because interfacial tension acts to minimize the surface area of the dispersed drops by coalescence. However, emulsions can be stablized for long periods of time by the presence of a third phase that adsorbs at the oil-water interface and inhibits coalescence of the dispersed phase."' Emulsions generally are stabilized by (1) repulsive charges on the surfaces of the dispersed phase and (2) adsorbed films that, being preferentially wetted by the continuous phase, act as a physical barrier to inhibit contact between the dispersed droplets. Both these factors generally are important in stabilizing and in breaking o/w emulsions. However, repulsive charges are less effective in stabilizing w/o emulsions which generally are stabilized by adsorbed films.V or this reason the nature and importance of these films at oil-water interfaces were examined closely; as will be shown, pH changes were found to affect the films significantly, and this work led to laboratory tests of pH- adjustment in emulsion breaking. This paper will present evidence to show that emulsions are stabilized primarily by films that form around droplets

Jan 1, 1969

By D. A. Flanagan, B. L. Landrum, P. B. Crawford, B. C. Norwood

Herern is described a heat conduction model which has been developed to predict the pressure distribution within a reservoir being subjected to rrnsteady-state depletion. Although potentiometric models have been widely used to study reservoir pressure distribution, their application is limited to the study of steady-state cases. Electrical analyzers may be used to study unsteady-state phenomena, but their application and simulation have certain limitations. Much of the now accepted theory of unsteady-state fluid flow in porous media has evolved from previous heat conduction knowledge. An analogy between the flow of fluid in an unsymmetrical reservoir and con-duct ion of heat in a similarly shaped unsymmetrical metal plate is made to show adaptability of the hear conduction principle to transient reservoir cases. A description of the model and necessary instrumenrarion is included. Data show influence of well location on pressure distribution during unsteady-state deple[ion. Application of results is made to a hypothetical reservoir and to a reservoir located in an irregular aquifer. INTRODUCTION Development of methods for predicting behavior of petroleum reservoirs has resulted partially from adaptation of electrical and thermal flow concepts to the flow of fluids through porous media. Extensive investigation of electrical and thermal flow problems has resulted in mathematical solutions applicable to many reservoir problems. The relative ease of measuring electrical flow has brought wide acceptance of the use of the potentiometric model to study certain field problems. Potentiometric models are restricted to the study of steady-state conditions because of the instantaneous response of electrical flow to changes in the applied potential. The electrical reservoir analyzer has been one of our most useful instruments for studying certain phases of unsteady-state behavior1. Other than approximation techniques, no method exists for determining pressure distribution in an irregular reservoir during primary depletion. Heat flow, like fluid flow, is inherently transient. The model described here was developed to use this transient nature of heat flow, measurable in the laboratory, to simulate transient flow of fluids in the porous medium of the reservoir above the bubble point. A metal plate shaped like the reservoir to be studied was heated by a "well" source. Temperature increase in the plate is analogous to pressure decline in the field if the field is being depleted, or pressure increase if fluid is being injected. This study is intended to determine the feasibility of using the transient nature of the heat flow to study the transient behavior of fluid flow. It is hoped that this new technique will aid in the understanding of one of the more complex phases of reservoir engineering. THE FLUID FLOW-HEAT CONDUCTION ANALOGY The similarity of the equations which describe the flow of fluids through a porous medium and the conduction of heat through a solid is well known2. Rate of flow of fluids is given by Darcy's law3. while the rate of flow of heat is given by Fourier', (I,, - KA 4T........(2) Combining Eq. 1 with the equation of continuity for a compressible liquid", i.e.,results in Expressing y in terms of the compressibility, i.e., pr = peey ,..........(5) and by expanding the exponential and neglecting higher orders of cp one can obtain

By J. V. Vogel

Jan 1, 1969

By T. Seldenrath

In the development of a fluid -operated hammer drill' for accelerated penetration of hard rock formations in oil wells, a research investigation was conducted to evaluate the percussion effects obtained with different design characteristics and to determine the possibilities of percussion in the low frequency range. Tests were conducted on granite blocks and comparable impact forces were measured with a load cell under selected test conditions. These full-scale laboratory tests provided an evaluation of the effectiveness of percussion which was further supported by actual downhole field performance. While much pertinent data have been presented by other investigators the method of evaluation described in this article resulted in good correlation between laboratory and field performance and may be applicable in other low frequency percussion developments. EXPERIMENTAL METHOD For these tests the bit was held stationary and the test block was rotated and forced upward against the bit by a hydraulically-lifted rotary table. "Weight on bit" was calculated from pressure readings on the fluid system of the lifting cylinders—corrected for tare and friction. Percussion was obtained with several tubular hammers of different weights up to 320 lb. These were lifted mechanically and allowed to drop by gravity from various heights onto a slidably-mounted bit sub or "floating anvil" for calibration of the load cell response with known kinetic energy of the hammers. During evaluation of percussion effects, the hammers were hydraulically operated to provide a range of percussion frequencies up to 1,020 blows per minute delivered to the anvil. To obtain comparable test readings, the load cell of the strain gauge type, was attached to the bottom of the anvil or bit sub in place of the bit. The load cell rested on a typical granite test block mounted on the drilling table. The latter was hydraulically supported as in the actual drilling operation. The effect of support resistance was the same for all tools when calibrating. Impact forces so determined werethusdirectly comparable with one another for the purposes of this investigation. When a roller bit was substituted for the load cell for actual drilling, the impact force would probably be reduced due to the resiliency of the bit body, legs. cones and teeth and the penetration of the teeth into the rock. However, since the bits and rock specimens were uniform in these tests, an evaluation of actual tool performance, in terms of parameters previously determined with the load cell at the selected standard conditions of calibration, proved satisfactory. Throughout these tests, 8¾-in. standard tricone, hard rock toothed roller bits and carbide-insert roller bits were used. Rotary speeds were 50 and 100 rpm; static weight on bit was applied up to 40,000 Ib, and percussion frequencies up to 1,020 blows per minute were available. During this investigation it was found that toothed bits on which the teeth had developed flat ends about 3/16 in. wide did not change appreciably for the duration of a test. Hence, all toothed bits used in this study were in this condition. Conventional reference curves obtained in these tests show that increasing weight on hit by 2:1 yielded an increased penetration rate of approximately 2.7:1, while an increase in rotary speed of 2:1 netted an increased penetration rate of roughly 1.7:l. Results are thus of the same order as those reported for hard rock drilling by other investigators." CALIBRATION Under the selected standard test conditions previously described, with the lifting table hydraulically supported, and the anvil interposed between hammer and load cell, the hammers were dropped from known heights. Thus, with the kinetic energy of the hammer known, the corresponding maximum impact force shown by the cell was determined. (It may be of passing interest to note here that special drop tests with the hammers showed that no appreciable reduction in impact force resulted from introduction of the anvil or hit sub.) From these calibrations, it was possible to determine the theoretical kinetic energy available in the hammer from the maximum impact force shown by the load cell while operat-

By C. D. Hall, F. E. Dollarhide

Fluid Ioss agent.s for crude oil and for water have been studied in dynamic tests. A treatment using a spearhead with a fluid loss agent followed by plain fluid appears feas ible in crude oil, but not in water. An equation for spearhead depletion shows that spurt loss relative to fracture width must be low, if the portion of spearhead fluid in the treatment is to be small. The presence of colloidal matter in crude oils aids the fluid Ioss agent. Unlike in kerosene, where flow limited the agent deposition, in crude oils the filter cake continually formed and leak-off declined. The volume-time relation varied somewhat for different crudes, but was best described by a square root of time function. Spurt loss was inversely proportional to agent concentration. After the fluid loss agent initiated the filter cake, the crude oil colloids built on it effectively. A 2-minute or a 5-minute spearhead with double the normal agent concentration gave the same fluid Ioss curve as the same concentration did for a 30-minute test. The agents tested in water gave fluid Ioss plots on which, for the first few minutes, volume was proportional to the square root of time, but later became proportional to time. For fracture area calculation the customary square root of time function is a satisfactory approximation. Leak-off rates and spurt losses were higher in water systems than in oils. The spurt Ioss tended to be inversely proportional to concentration. In spearhead tests, the filter cakes were not eroded by water flow. However, the rather high spurt loss values make spearhead treatments impractical for water-based fluids. Introduction The effects of dynamic testing conditions on the performance of fluid loss agents in kerosene have been studied previously.' We have extended the work to include crude-oil- and water-based fracturing fluids. An understanding has been gained of the mechanisms of formation and functioning of the filter cakes of fluid loss agents. The practical aspects of evaluating performance of agents in relation to fracture area calculations also are considered. The feasibility of using the fluid loss agent in a spearhead stage of the treatment is examined further for both types of fluids. Experimental Procedure The dynamic fluid loss tests were performed in an apparatus similar to the high-pressure apparatus described in a previous publication.' A fracturing fluid was circulated over a rock surface located in a closed pressurized loop. The fluid flowed axially over the cylindrical surface of a core 2 in. in diameter X 3.5 in. long, mounted (with the flat ends sealed off) in a pipe, with 0.117 in. annular clearance. The filtrate was collected in a central hole in the core and led through valves to graduated cylinders. Provision was made for changing quickly the circulating fluid during the test (spearhead runs) without interrupting the filtration pressure. The only modifications were to add heating tapes and water jackets for the tests with crude oils, all conducted at ISOF, and to change all parts exposed to the test fluid to stainless steel for the tests with water-based fluids. The latter tests were made at room temperature, 80F. Three crude oils were tested. A mixed crude, obtained from a local refinery, contained a considerable amount of light ends. For safety reasons, it was stripped to 250F vapor temperature before use in the fluid loss tests. The other two oils were used as obtained from lease tanks. One was a greenish-brown, 37" API paraffinic crude, and the other was a black, 32" API asphaltic crude. The fluid loss agent for oil, here designated for brevity as Agent A, was Adomite@ Mark II*, a granular solid commercial agent, the same as previously tested in kerosene.' Three different compositions of fluid loss agents were tested in Tulsa tap water. Agent B was adomit& Aqua*, a solid commercial fluid loss agent, comprising clays and hydrophilic gums principally derived from starch. Agent C was a mixture of three parts of Agent B with two parts of silica flour. Agent D was Dowel1 J137, a mixture of guar gum and silica flour. The test cores were cut from contiguous blocks of Berea or Bandera sandstones. For the oil tests, the cores were oven dried, evacuated, saturated with kerosene, and the kerosene permeability was measured. The cores used with the water-based fluids were pretreated by saturating with 3 percent calcium &loride solution to minimize pemeability damage by the fresh water due to clay migration. The

Jan 1, 1969

By H. D. Hodges, J. K. Godbey

In order to better understand the fracturing process, bottom-hole pressures were measured during a number of typical fracturing operations. A recently developed system was used that allows simultaneous surface recording of both the bottom-hole and wellhead pressures on the same chart. The results from six fracruring treatments are summarized on the basis of the pressure data obtained. Al-though no complete analysis is attempted, the value of accurate pressure measurements is emphasized. Important characteristics of the bottom-hole pressure record do not appear at the wellhead because of the damping effect of the fluid-filled column. In four of the six treatments described, the formations apparently fractured during the initial surge of pressure with only crude oil in the well. The properties of the fluids used during the treatments are given and the fluid friction losses are obtained directly from the pressure records. This technique is also shown to be adequate for determining when various fluids, used during the process, enter the formation. INTRODUCTION Hydraulic fracturing for the purpose of increasing well productivity is now accepted in many areas as a regular completion and workover practice. Numerous articles have appeared in the literature discussing the various techniques and theories of hydraulic fracturing'. In general, three basic types of formation fractures are recognized today. These are the horizontal fracture, the vertical fracture, and fractures along natural planes of weakness in the formation'. Any one or all three of these fracture types may be present in a fracturing operation. However, with only the wellhead pressure record as a guide, it is difficult at best to determine if the formation actually fractured, and is almost impossible to determine the type of fracture induced. These difficulties arise in part because the wellhead pressure record, especially when fracturing through tubing, does not accurately reflect the pressure variations occurring at the formation. Several factors contribute to this effect and preclude the possibility of using the wellhead pressure as a basis for accurately calculating the bottom-hole pressure. These factors are: 1. The compressibilities of the fluids which damp the pressure variations. 2. The changes in the densities of the fluids or apparent densities of the sand-laden fluids. 3. The flowing friction of the various fluids and mixtures, which is dependent on the flow rates and the condition of the tubing, casing, or wellbore. 4. The non-Newtonian characteristics of a sand-oil mixture and its dependence upon the fluid properties, the concentration of sand, and the mesh size used. 5. The unknown and variable temperatures throughout the fluid column. Because of these reasons it was determined that in order to obtain a more accurate knowledge of the nature of fracturing, the bottom-hole pressure must be measured along with the pressure at the surface during a fracturing treatment. Even with accurate pressure data, a reliable estimate of the nature of fracturing is still dependent upon knowledge of the tectonic conditions. However, the hydraulic pressure on the formation is basic to any approach to a complete analysis. In order to accomplish this objective a system was developed to record the wellhead and bottom-hole pressures simultaneously at the surface. By recording both pressures on a dual pen strip-chart recorder, it was possible to greatly expand the time scale so that rapid pressure variations would be faithfully recorded. By such simultaneous recording, time discrepancies inherent in separate records are eliminated, thus overcoming one of the most difficult problems associated with bottom-hole recording systems. This paper illustrates the results obtained by using this system during six typical fracturing operations. All of these tests were taken in wells that were treated through tubing. By a direct comparison of the wellhead and bottom-hole pressures, the importance of obtaining complete pressure information during a fracturing treatment is emphasized. THE INSTRUMENTATION AND PROCEDURES The bottom-hole pressure measuring instrument consisted of a pressure-sensing element, a telemetering section, and a lead-filled weight or sinker bar. The pressure-sensing element used was an isoelastic Amerada pressure-gauge element. By using an isoelastic element, no temperature compensation was necessary in the tests described, since the temperature was believed to be well below the maximum temperature limit of 270°F. The rotary output shaft of this helical Bourdon tube element was coupled to a precision miniature potentiometer. The rotation of the pressure-gauge shaft thus changed the resistance presented by the potentiometer

By A. G. Weber, R. L. Ledbetter, A. L. Hicks

A procedure which makes use of both analog and digital computers has been developed for predicting the pressure-production behavior. of water-drive reservoirs. The electric analyzer is used for matching the past pressure and production history of the reservoir and inferring the physical characteristies of the aquifer. When the history has been matched, a .set of base curves called influence functionce" are derived and the future behavior of the reservoir is calculated from these curves. INTRODUCTION Electric analog methods are particularly well-suited for complex water-drive problems. Like all other methods of analyzing water-drive reservoirs, however, analog methods are subject to certain limitations. One of the more obvious and serious of these results from inaccurate or insufficient knowledge about physical characteristics of the reservoir and aquifer. However, it is not insurmountable because the unknown characteristics can usually be inferred from the reservoir's previous behavior. Simons and Spain have shown that even with scant reservoir data and limited production history, electric analyzer studies can provide useful and reliable results. Aside from this limitation on analyzer predictions, the main disadvantage has been lack of a technique for keeping predictions abreast of latest production forecasts. Pressure behavior of a reservoir depends on the rate at which fluid is withdrawn; it is therefore necessary to specify the future withdrawal schedule before the future pressure behavior can be calculated. Conventionally, the electric analyzer predicts pressure behavior for a specific production schedule. If this schedule is abandoned in favor of some new one, as is frequently the case, either the study must be rerun or results adjusted in some manner. The usual procedure is to predict pressures for several different forecasts covering a wide range of production rates. Pressures for specific rates other than those studied can then be approximated by interpolating between the analyzer predicted curves. However, reliability of this type of prediction is doubtful. This paper describes a technique for predicting future behavior of water-drive reservoirs that retains all advantages of electric analyzer analysis and eliminates the main disadvantage. That is, the analyzer is used conventionally to infer unknown reservoir parameters from the pressure history, but not to predict the future behavior directly. Instead, a derived set of base curves, called influence l'unctions, can be used to calculate future behavior for any production schedule.