Search Documents

By Owen R. Rice

Changes in coke hardness affect the working of the blast furnace, for soft coke is an obstacle to proper furnace operation. Soft coke is due to a low hydrogen-oxygen ratio in the coal charged; increasing this ratio increases the hardness of the coke. The best means of determining the hardness of coke is the combined tumbler and ball-mill. What is Good Coke? It is not always appreciated that the correct answer may be such a reliable forecast to the physical condition of the blast furnace, and incidentally its entire operation, that among the ordinary variations in other raw materials and also in operating conditions, it assumes a position of primary significance. Inasmuch as coke forms two-thirds of the total bulk of material charged, this assertion will not seem exaggerated, and in the subsequent discussion its applicability will be presented. Coke is the ideal blast-furnace fuel because it is strong and hard enough to resist the crushing and grinding effects of the furnace burden, and its porosity allows the ready passage of gas upward through the charge, and rapid combustion before the tuyeres. Coke deprived of its strength and with a cell structure physically frail and chemically susceptible to solution by carbon dioxide has its chief asset eliminated. When the melting zones are approached the fluid cinder trickling toward the hearth will envelope the small pieces of coke and form a sticky mass that will retard the free passage of gas. Blast pressure will then rise, driving will slow up and fragments of coke will lie inert before the tuyeres, rendered practically incombustible by their slaggy envelop. It is believed also that coke messes at cast are attributable to quantities of such coke floating unburned in the hearth. It has been correctly claimed that coke rendered small outside the oven possesses greater combustibility than that so developed in the

Jan 1, 1922

By Oliver Bowles

WHILE vast quantities of limestone and dolomite. are used in metallurgy, the estimated production in 1926 being 23,860,000 tons, there are many problems connected with their use which have not received adequate study. The literature of metallurgy is notably lacking in comprehensive discussions of fluxing or furnace stone. Approximately 900 lb. of limestone is used for every long ton of pig iron produced in the blast furnace, but this important constituent of the charge receives little attention compared with the intensive study that has been applied to ores and fuels, the other important constituents of the charge. In an effort to elucidate some of the doubtful points, the writer has made a wide review of metallurgical literature, and has had the disappointing experience of finding not more than a paragraph or two on metallurgical stone' for each two or three thousand pages of literature reviewed. Furthermore, the limited information supplied is in some instances contradictory and even absurd. For example, a correspondent in a reputable metallurgical magazine classifies magnesia and silica together as impurities in flux "that must be melted and run off as slag." From the furnace operator's viewpoint it is desirable that as complete knowledge as possible be obtained concerning the effects of chemical composition and physical character of the stone on its action in the furnace and on the quality of the metal produced. A better knowledge of these features is advantageous also to the quarry operator in that it enables him to work more intelligently in producing stone that will best satisfy the metallurgists' requirements. The following pages comprise a brief compilation of data now available. Much of the information is well established and probably acceptable to most metallurgists. On some points there is room for debate, and the chief purpose in presenting a paper at this time is to arouse discussion which will throw more light on controversial points and thus assist the author in preparing a more comprehensive and accurate discussion of metallurgical limestone problems.

Jan 1, 1928

By Oliver Bowles

While vast quantities of limestone and dolomite are used in metallurgy, the estimated production in 1926 being 23,860,000 tons, there are many problems connected with their use which have not received adequate study. The literature of metallurgy is notably lacking in comprehensive discussions of fluxing or furnace stone. Approximately 900 lb. of limestone is used for every long ton of pig iron produced in the blast furnace, but this important constituent of the charge receives little atten. tion compared with the intensive study that has been applied to ores and fuels, the other important constituents of the charge. In an effort to elucidate some of the doubtful points, the writer has made a wide review of metallurgical literature, and has had the disappointing experience of finding not more than a paragraph or two on metallurgical stone for each two or three thousand pages of literature reviewed. Furthermore, the limited information supplied is in some instances contradictory and even absurd. For example, a correspondent in a reputable metallurgical magazine classifies magnesia and silica together as impurities in flux " that must be melted and run off as slag." From the furnace operator's viewpoint it is desirable that as complete knowledge as possible be obtained concerning the effects of chemical composition and physical character of the stone on its action in the furnace and on the quality of the metal produced. A better knowledge of these features is advantageous also to the quarry operator in that it enables him to work more intelligently in producing stone that will best satisfy the metallurgists' requirements. The following pages comprise a brief compilation of data now available. Much of the information is well established and probably acceptable to most metallurgists. On some points there is room for debate, and the chief purpose in presenting a paper at this time is to arouse discussion which will throw more light on controversial points and thus assist the author in preparing a more comprehensive and accurate discussion of metallurgical limestone problems.

By John. C. Reed

THIS paper presents briefly some of the principal results of 3 ½ months field work during the summer of 1938 in the Chichagof mining district, southeastern Alaska. The report is preliminary and presumably will be followed later by a more complete report to be published by the Geological Survey. The paper is concerned principally with the structural geology of a block of steeply dipping graywacke as revealed underground in the Hirst Chichagof and the Chichagoff mines, the two principal mines of the district. The two mines have produced about $16,000,000 worth of gold and silver. The attitude of the graywacke, which is locally cut by dikes and sills of aplite, is remarkably uniform and trends about N. 62° W. The average dip is about 67° SW. The graywacke is cut by large faults, the best known of which are the Hirst and the Chichagof, whose average trends are about N. 30° W. and average dips about 75° SW. The faults at some places split into the bedding of the graywacke and at others are deflected from their normal trends for considerable distances along the bedding. The relative movements along the faults was northwest and up at about 30° for the south-west or hanging-wall sides. The ore bodies are elongate quartz lenses in the faults. Most of the quartz is mineralized but locally large quartz masses are too low grade to be mined. The emplacement of the ore bodies appears to have been structurally controlled in large part by the intersections of the faults with the bedding surfaces. Places where the faults split or are deflected appear to have been particularly favorable places for ore deposition. There has been some post-quartz movements along the faults but it is thought that such movements have not been extensive enough to greatly modify the original positions and shapes of the ore bodies.

Jan 1, 1939

By D. G. F. Hedley, D. F. Coates, K. Barron

INTRODUCTION Preceding papers presented in this symposium have considered the influence of various parameters, such as geologic structure, groundwater, etc., on the stability of rock slopes and have given some analyses and guide lines for slope design. To carry out such analyses it is necessary to put numbers to these parameters; where possible, this is best done on the basis of in- situ observation and measurement. Since none of the design analyses can be considered exact and since conditions vary during excavation, it is also desirable to monitor the behaviour of the slopes as the excavation proceeds. In this manner a feedback of information on the actual excavation conditions can allow a continuous reassessment of the design parameters to be made and can give early warning of incipient failures. The purpose of instrumenting rock slopes is thus to provide information which fits into one of two categories: (a) To give basic data to assist on-site assessment and design of the rock excavation. (b) To provide data for the control and assessment of safety during excavation. The purpose of this paper is to review briefly some of the instruments and measuring techniques available for gathering this data and to indicate where such data may prove useful. This is not a comprehensive review of all available instrumentation but, rather, is a survey of some of the instruments and techniques with which the authors have had some experience and which have proven useful. STRUCTURAL ASSESSMENT Whilst most structural assessment of a site is done on the basis of geologic mapping and core examination, there exists some instrumentation which can, in particular circumstances, assist in this assessment. Seismic instruments can be used to determine the geometry of formations and to give an indication of the relative competency of rocks. These instruments operate using seismic reflection or refraction techniques and, although mainly confined to shallow sub-surface exploration, can be useful in distinguishing between soil/ bedrock contacts and, in some cases, the extent of weathered surface rock.

Jan 1, 1971

THIS paper presents briefly some of the principal results of 3 ½ months field work during the summer of 1938 in the Chichagof mining district, southeastern Alaska. The report is preliminary and presumably will be followed later by a more complete report to be published by the Geological Survey. The paper is concerned principally with the structural geology of a block of steeply dipping graywacke as revealed underground in the Hirst Chichagof and the Chichagoff mines, the two principal mines of the district. The two mines have produced about $16,000,000 worth of gold and silver. The attitude of the graywacke, which is locally cut by dikes and sills of aplite, is remarkably uniform and trends about N. 62º W. The average dip is about 67º SW. The graywacke is cut by large faults, the best known of which are the Hirst and the Chichagof, whose average trends are about N. 30ºW. and average dips about 75º SW. The faults at some places split into the bedding of the graywacke and at others are deflected from their normal trends for considerable distances along the bedding. The relative movements along the faults was northwest and up at about 30º for the south- west or hanging-wall sides. The ore bodies are elongate quartz lenses in the faults. Most of the quartz is mineralized but locally large quartz masses are too low grade to be mined. The emplacement of the ore bodies appears to have been structurally controlled in large part by the intersections of the faults with the bedding surfaces. Places where the faults split or are deflected appear to have been particularly favorable places for ore deposition. There has been some post-quartz movements along the faults but it is thought that such movements have not been extensive enough to greatly modify the original positions and shapes of the ore bodies.

Jan 1, 1939

By C. D. Stahl, S. M. Farouq Ali, W. E. Culham

One- and two-dimensional mathematical models have been developed that simulate transient, two-phase flow of hydrocarbon mixtures in porous media in a manner that accounts for interphase mass transfer. Numerical simulations of one-dimensional depletion-drive experiments using a two-component hydrucarbon fluid were used to establish the validity of the mathematical models. In addition, the experimental and numerical data were used to demonstrate that production rate had a relatively insignificant effect on the recovery of individual hydrocarbon components from the experimental system, and that attainment of equilibrium between phases is possible for a wide range of liquid and vapor velocities in reservoirs containing light hydrocarbon fluids. Results of some two-dimensional numerical simulations are also presented. INTRODUCTION This study was undertaken to develop a mathematical model that would simulate transient, two-phase flow of hydrocarbon mixtures in porous media under conditions that result in interphase mass transfer and to test the validity of the assumptions used to set up the model. In addition, the study was designed to determine if production rate influences the recovery of individual hydrocarbon components from reservoirs producing by depletion drive. Two-phase flow in porous media, with interchange of components between the two phases, is important in many petroleum recovery processes. Studiesl-5 conducted within the last 3 years have outlined methods of solving multiphase flow problems incorporating mass transfer. Some of these studies have also indicated the importance of accounting for mass transfer under various producing conditions. An earlier work6 first demonstrated the importance of combining relative permeability data with equilibrium ratios in compositional balance methods. The mathematical model presented in this paper is formulated so that a phase behavior package, as described in previous papers,5,17 is not required as an integral part of the routine employed to solve for the primary dependent variables. The finite difference formulation is designed so that all the primary variables can be solved for simultaneously. This is accomplished by utilizing one basic set of equations. These innovations, which are in contrast to other models2,3,5 but are similar in some respects to the approach used by Taylor,l render the total problem computationally simpler than any of the previously referenced formulations. The mathematical model was developed by combining Darcy's law with a continuity equation for each hydrocarbon component. The principal assumptions invoked in the formulation were (1) that capillary forces and diffusional effects are negligible, and (2) that thermodynamic equilibrium exists in the reservoir at all times. No assumption as to the type of vaporization process was made in formulating this model. Experimental data were required to complete this study. These were generated by conducting several depletion drive experiments. The experimental apparatus consisted of a sandstone core enclosed in a pressurized casing. The apparatus was designed in such a manner that the core could be charged with a liquid hydrocarbon mixture and depleted at different production rates. The experimental tests were designed to determine the effect of production rate on component recovery. In addition, direct comparison of experimental and mathematically predicted

Jan 1, 1970

By B. A. Eaton

The subject of many discussions and technical papers in the last 20 years has been the prediction of the well-bore pressure gradients that are required to induce or extend fractures in subsurface formations. The subject merits this attention because of the frequently recurring problems that arise from an inability to predict fracture pressure gradients. Encountered in several common types of operations in the oil industry are problems associated with the prediction of formation fracture pressure gradients. When wells are being drilled in both new and old fields, lost circulation is often a very troublesome and expensive problem. Complete loss of circulation has been disastrous in some cases. Many times, such disasters could have been avoided if techniques for calculating fracture pressure gradient had been employed in the well plans, and if casing strings had been set, and mud weight plans had been followed accordingly. In areas of abnormally pressured formations, the prediction of fracture gradients during the well-planning stage is extremely important. In fact, it is as important as the prediction of formation pressure gradients, which has received a great deal of attention in recent years. There are several published methods used to determine fracture pressure gradients. However, none of these methods appears to be general enough to be used with much reliability in all areas. In 1957, Hub-bert and Willis published a classical paper that included the development of an equation used to predict the fracture extension pressure gradient in areas of incipient normal fau1ting.l Overburden stress gradient, formation pore pressure gradient and Poisson's ratio of rocks were the independent variables that were shown to control fracture pressure gradient, the dependent variable. In 1967, Matthews and Kelly published another fracture pressure gradient equation that is different from that of Hubbert and Willis in that a variable "matrix stress coefficient" concept was utilized.3 Later the same year, Costley wrote about a similar idea.5 Goldsmith and Wilson used a least-squares curve-fitting technique and field data from the Gulf Coast area to correlate fracture pressure gradient with formation pore pressure gradient and formation depth.' They noted that the fracture pressure gradient increased with increasing depth while the pore pressure gradient remained constant. In each of these cases, the problem for which a solution was sought was to determine the bottom-hole pressure gradient required to initiate or extend a fracture. Results of the previous work show that fracture pressure gradient is a function primarily of overburden stress gradient, pore pressure gradient, and the ratio of horizontal to vertical stress. There is argument for a fourth variable in that in many cases breakdown fracture pressure gradient is greater than the fracture extension pressure gradient. However, if the fracturing fluid is able to penetrate the formation through the pores or existing cracks, there is very little

Jan 1, 1970

By J. S. Aronofsky

The first step in a quantitative analysis ot the mechanism of oil displacement by water in a fractured reservoir is usually conceded to be the solution of the differential equation describing the saturation distribution of two immiscible fluids flowing in a porous medium, where the capillary pressure is taken into account. In such a system the production mechanism may consist of displacement of oil both by the flow of water due to natural or artificially imposed pressure gradients and by imbibition, which implies a flow of water not due to external pressure gradients. Owing to the presence of the two oil displacement mechanisms, the mathematical model given by the differential equation intended to describe the system may not properly represent the behavior of the physical system. In fact, in the reservoir the rate of water advance may be very slow, and in the case of a fractured reservoir with a great number of large fractures, the pressure difference determining the flow of water through the matrix may be much less than 1 lb/'psi over lengths of a few feet. In such a case, imbi-bition (the exchange between oil in the matrix and water in the fractures resulting from capillary forces) may become, with time. a significant element of the production mechanism. It occurred Lo the authors, however, that without going into a physical analysis of the process of production, it might be possible by means of simple abstract reasoning to throw some light on the variation of recovery with time under conditions occurring in a highly fractured oil reservoir with rising water table. The object of this paper is to present both the reasoning and its application to a reservoir of the highly fractured type. Specifically, the analysis given here was undertaken to try to explain the increase of recovery (as defined later) with time as observed in this reservoir, without having to assume unlikely variations in the reservoir parameters with depth. This attempt has been successful as will become clear upon comparison of the computed recoveries with the actual field data. ABSTRACT MODEL Let us consider a small volume. of porous matrix saturated with oil at time, t = 0. Let the process of oil displacement by water start at time, t = t,,. At some time, t, the process will have terminated. Then a volume of oil equal to or smaller than the original oil contained in the matrix will have been produced. The first basic assumption that describes the model and guides the forthcoming reasoning is that the oil production from the small volume. dv, is a continuous monotonic function of time and that it converges to a finite limit. Such an assumption is not inconsistent with the results of laboratory waterflood tests as well as with results of imbibition tests where this is, in fact, observed. The second basic assumption is that none of the properties which determine the rate of convergence change sufficiently during the process to affect this rate or the limit. Let it be assumed that the form of the function of time relative to production from the matrix volume. dv, is given by where V,(t) is the volume of oil produced up to that time t. R is the limit toward which the recovery converges, A is a constant giving the rate of convergence, and V,(t) is the volume of oil originally in place in the volume, dv. It should be noted that recovery at time I will be understood here to he It follows that r. tends to R as t tends to infinity. CONSTRUCTION OF RESERVOIR FROM ABSTRACT MODEL Let the reservoir consist of a series of identical blocks of porous matrix stacked vertically and separated by fractures. Let each of these blocks satisfy the conditions of our abstract model. These conditions are: (1) that recovery is a continuous mono-tonic function of time converging to finite limit, and (2) none of the properties that determine the rate of convergence change sufficiently to affect the rate or the limit. Let water be rising in fractures so that oil production from any part of the block starts when the water comes in conLact with it. TOTAL RECOVERY FROM THE RESERVOIR COMPARED TO RECOVERY FROM A SINGLE MODEL ELEMENT As stipulated in Eq. 1, in the case of the abstract model,

By MORTON WEBB

Discussion of the paper of Max Boehmer, presented at the Pittsburg meeting, March, 1910, and printed in Bulletin. No. 38, February, 1910, pp. 119 to 122. W. MORTON WEBB, Germiston, Transvaal, South Africa (communication to the Secretary*):-The experience of Mr. Boehmer in the Leadville district and his reputation as an engineer assure the interest of anything he may write on the ore-bodies of that district; and all engineers who have operated in Leadville will agree that the theory of the origin of the Leadville deposits, as given by Mr. Emmons in his famous monograph, does not adequately explain the genesis of all the deposits. The solution of the problem does not, however, appear to be quite so simple as Mr. Boehmer would make it; and the advance pamphlet issued by the U. S. Geological Survey, pending the completion of the main report, shows that there is still much ground for controversy. No doubt the later eruptives had a greater influence on the deposition of the ore-bodies than was thought at the time of the first investigation by the Survey; but it has always appeared to me that the influence exerted by these dikes was more in the nature of a physical than a chemical one. The weight of evidence in recent years is mostly in the direction of ascension rather than lateral secretion; and one of the strongest evidences is the fact that most of the ore-bodies lie under one of the comparatively impervious layers which help to make up the structure of these hills. On Iron hill I have found a very persistent white porphyry sheet, split off from the overlying porphyry and lying from 30 to 60 ft. below the main white porphyry and blue limestone contact. This sheet varies from a few to 50 ft. in thickness, and where found adds another "contact" to that region ; some of the largest ore-bodies on Iron hill lie immediately under this sheet. Above this is the main white porphyry and * Received Mar. 3, 1910.

Feb 1, 1911

By C. R. McEwen

The displacement equations of Buckley and Lever-ett' have been successfully applied to the prediction of oil recovery in frontal drives for a number of years. Commonly, the capillary pressure term is omitted from Leverett's fractional flow formula. The resulting satura-tion-vs-distance curves become multi-valued, however. Buckley and Leverett, recognizing that the neglect of capillarity was most likely the cause, proposed that this part of the solution he ignored. An ingenious method of propagating a discontinuity in place of the multi-valued portlon of the curve was devised by Cardwell using "non-capillary displacement theory". However, it would seem that a completc description of the displacement pracess would require the inclusion of capillary forces. Three publications on displacement theory have appeared in which capillary pressure was included. Ter-williger. et al', calculated saturation profiles which matched those which they observed when gas displaced water vertically downward. Based on their study, Jones-Parra and Calhoun5 roposed a method of computing saturation distributions sufficiently removed in time from the initial conditions so that all saturations in the flood front region were assumed to move with the same velocity. This latter assumption corresponds to the tangent construction on the flowing fraction-vs-satura-tion curve proposed by Welge." A numerical solution to the problem of water displacing oil from a water-wet porous medium was recently presented by Blair, ef at'. The initial conditions included a uniform saturation slightly above interstitial water saturation. that is, the permeability to water ahead of the front was not zero. This appeared to result in a region of slowly rising water saturation ahead of the front, at least for the particular set of relative permeability and capillary pressure data used to illustrate the computation. The authors also postulated an "end-effect" at the outflow face of the system, whereby no water could pass out of the system until the water saturation had built up to a high value. The saturation profiles computed under these assumptions show water accumulating at the downstream end of the system while the flood front proper is still traveling toward the outlet. The purpose of this paper is to present the results of a numerical solution to the one-dimensional water-flooding equation in which no outlet end-effect is considered and where there is no flow of water ahead of the flood front. It should be noted that the contemporary paper of Fayers and Sheldon8 deals with the problem under discussion. For a thorough treatment including both capillarity and gravity the reader is referred to their paper.

By W. B. Gogarty, H. R. Bailey

Methods are presented for determining molecular diffusion coefficients by using data from, capillary flow experiments. These methods are based on a numerical solution (presented in a previous paper) of the partial differential equation describing the combined mechanisms of flow and diffusion. Results from this numerical solution are given and compared with the approximate analytical solution of G. I. Taylor. The numerical solution is valid over a much larger time range. These methods are applied to experimental results for the fluid pairs water-potassium permanganate solution and amyl acetate-orthoxylene. Both of these fluid pairs have approximately equal densities and viscosities. Graphical and numerical techniques are presented for deter mining diffusion coefficients from the flow data. Values obtained by these techniques are compared with values obtained by other methods. INTRODUCTION The molecular diffusion coefficient is known to be a variable in determining the amount of mixing in a miscible displacement process. The effect of molecular diffusion on dispersion in longitudinal flow through porous media -has been examined by different investigators. These investigators concluded that at low velocities of flow, the amount of dispersion is approximate1y proportional to the molecular diffusion coefficient. The influence of diffusion on fingering, channeling, and overriding has been mentioned by other investigators.4 Recent studies5,6 have been made on the effects of molecular diffusion in connection with the problem of gravity segregation. Many different methods have been developed for the experimental determination of molecular diffusion coefficients? These methods differ mainly according to boundary conditions selected and analytical procedures used. Nevertheless, all of these methods have the condition in common that the bulk fluids in which diffusion is occurring are stationary with respect to each other. In connection with a series of papers on mixing in capillary flow,8-10 Taylor suggested the use of a flow method for determining molecular diffusion coefficients. Additional studies 11-13 have been conducted on miscible displacements in capillary tubes, but the data from these studies were not used for the specific purpose of determining diffusion coefficients. The flow method proposed by Taylor results in a single value of the diffusion coefficient for the fluid pair used in the displacement experiments. This single value represents the true value for the fluid pair when the diffusion coefficient is independent of concentration. If the diffusion coefficient is a function of concentration, the single value obtained by the flow method gives an average value for the coefficient of the fluid pair. These average values are based on diffusion taking place over the entire range of concentration, i.e., from 0 per cent of one fluid to 100 per cent of that same fluid. In field applications of the miscible displacement process, gradients occur over the same range of concentration as are found in the displacements in capillary tubes. Molecular diffusion coefficients obtained from the capillary flow method should, therefore, be especially relevant to field operations. This investigation was undertaken to evaluate the feasibility of obtaining molecular diffusion coefficients from capillary -flow experiments. In making this evaluation, diffusion coefficients were first determined for two systems from data obtained in capillary flow experiments. These values of the diffusion coefficient were then compared to values obtained by other methods. MIXING IN CAPILLARY FLOW—THEORETICAL The theoretical basis for determining molecular diffusion coefficients from capillary flow experiments is the partial differential equation relating the mechanisms of flow and diffusion. For the case of laminar flow and with a fluid pair having equal densities and viscosities, the equation is given by

By G. Thodos, R. B. Grieves

A method has been developed For predicting the critical temperatures and critical pressures of binary mixtures of carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen with the normal paraffin hydrocarbons. For carbon-dioxide and hydrogen-sulfide systems, relations are presented that take into account the peculiar behavior of mixtures with closely boiling components, such as carbon dioxide-ethane and hydrogen sulfide-propane mixtures which exhibit minimum critical temperature points. For hydrogen, nitrogen and carbon-monoxide systems, the extreme critical behavior caused by wide differences in pure component properties is established In addition, those fixed gas-paraffin systems which resemble paraffin-paraffin systems are also accounted for. For a mixture of known composition, the pure component critical temperatures, critical pressures and normal boiling points are all that are required to determine its critical point. Graphical relations are presented relating Tc and PC , of the mixture to the pure component properties. From the treatment of 12 carbon-dioxide and hydrogen-sulfide systems reported in the literature (74 mixtures), the expected error lor the critical temperature is approximately 1.5 per cent, and for the critical pressure, approximately 2 per cent. From the treatment of six hydrogen, nitrogen and carbon-monoxide systems reported in the literature (30 mixtures), the expected error for both the critical temperature and critical pressure is approximately 2.5 per cent. The relationships, which have been developed with only normal paraffins as the hydrocarbon components, may be extended to those isoparaffins and olefins which fall within the allowable volatility ranges. INTRODUCTION Many of the fixed gases — carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen — occur in natural mixtures with hydrocarbons. Carbon dioxide and hydrogen sulfide are frequent components of the fluids produced from underground petroleum reservoirs. Nitrogen, carbon dioxide and hydrogen sulfide are present in varying quantities in most natural gases and gas-condensate well effluents. Hydrogen mixtures are of considerable interest in many phases of refining processes of petroleum. The determination of the critical temperatures and critical pressures of such mixtures is of value in vapor-liquid equilibrium studies, for the prediction of the characteristics of underground reservoirs, and for reduced-state correlations of PVT, transport and thermodynamic properties. The accurate estimation of the critical point for binary mixtures is an important initial step toward a complete analysis for the establishment of the critical temperatures and pressures of multi-component mixtures. Methods for predicting the critical temperatures and critical pressures of binary hydrocarbon systems have already been presented in the literature.8 It is possible to apply these existing methods to fixed gas-paraffin mixtures, but due to their unusual critical behavior, values calculated deviate considerably from experimental values. For systems containing trace quantities of the fixed gases, these methods are acceptable; however, for systems containing more than 5 mol per cent of the fixed gases, these utterly fail to produce reasonable critical values. Consequently, in this study a method has been developed for handling such binary mixtures over the entire composition range. CARBON-DIOXIDE AND HYDROGEN-SULFIDE SYSTEMS The critical behavior and the vapor pressure behavior of mixtures of carbon dioxide and hydrogen sulfide with paraffinic hydrocarbons may be quite similar or quite dissimilar to that of paraffin-paraffin mixtures, depending on the volatilities of the components involved. The critical temperature and normal boiling point of carbon dioxide are very close to the corresponding values for ethane, while its critical pressure is considerably higher than that of ethane. The normal boiling point of

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge's method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey's theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey's method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge's6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

By C. Chu

A theoretical investigation has been made of the forward combustion process using a one-dimensional linear mathematical model, taking into consideration the effect of the vaporization-condensation which occurs on the leading edge of the heat wave. This work involves the solution of five coupled partial differential equations. Besides the vaporization-condensation phenomenon, these equations account for conduction, convection, combustion, heat loss, diffusion and bulk fluid flow. For the one-dimensional linear model studied, the vaporization-condensation phenomenon does not induce appreciable change in the temperature at the combustion front; and its primary effect is to create a steam plateau and to increase the length of the heated zone ahead of the combustion front. This effect becomes more pronounced at lower pressures, higher porosities or reduced gas saturations. The peak temperature and the temperature profile on the leading edge of the heat wave stabilize after a certain period. The length of the steam bank remains practically constant, although the length of the water bank increases as the heat wave advances. INTRODUCTION The existence of the vaporization-condensation phenomenon in the heat-wave process and the important role played by the phenomenon have been recognized by several investigators. Kuhn and Kochl stated that steam plateaus were frequently observed on the temperature records of their experiments. The steam plateaus were attributed primarily to the vaporization and subsequent condensation of the interstitial water existing in the oil sand. Szasz2 suggested that both lighter hydrocarbons and water are vaporized on the leading edge of the heat wave, carried forward in the gas stream, and then condensed to create banks of oil and water. Martin et al suggested that the vaporization-condensation phenomenon is one of the main mechanisms of the heat-wave process, along with thermal expansion and viscosity reduction. Wilson et al. reported the existence of a steam plateau several inches in length in their small-scale tube-run experiments. However, this important phenomenon has never been taken into consideration in the numerous theoretical analyses by various authors. 5 The purpose of this work was to study the thermal aspects of a linear heat wave, taking into consideration the vaporization-condensation on the leading edge of the wave, to determine the effect of this phenomenon on the temperature profile of the reservoir, and to investigate how this effect varies when other process variables are changed. THEORY We consider a reservoir of porous medium of cross-sectional area A, extending from x = 0 to x = L. This reservoir contains, aside from the solid matrix itself, a gas phase and a "combined liquid phase" which is a combination of two immiscible liquid phases — namely, an oil phase and a water phase. The oil present in the reservoir is assumed to consist of three fractions, a noncondensable gas, a nondistillable residuum, and a vaporizable oil fraction which may be present in both the gas phase and liquid phase. Before the heat-wave process begins, preheating has taken place and has imparted an initial temperature distribution T,(x) to the reservoir. At the start of the process, a stream of oxygen-bearing gas is introduced through the face at x = 0. This gas supports the combustion of the residual fuel and supplies the heat throughout the process. In order to develop a tractable system of equations, the model properties are postulated as follows: (1) the system is isotropic and homogeneous; (2) there is a uniform fuel content throughout the reservoir; (3) combustion takes place in a moving vertical planar region of infinitesimal thickness; (4) the temperatures of the fluids and the solid matrix are identical; (5) radiation is negligible; (6) the gas and liquid saturations are uniform throughout the reservoir and remain unchanged

Jan 1, 1965

By F. W. DeVries

The subject article is an excellent comparison between two potential oxidants, molecular oxygen and hydrogen peroxide, as used in the Palangana field at the time of the trials. The tests showed no advantage to using hydrogen peroxide; however, it would be unwise to generalize that the experience in these trials was universally true. During the time of the work cited, extremely high levels of an unusual impurity were recycling in the lixiviant, which substance was a catastrophically destructive catalyst for peroxide. Since this material was not removed before recycle, it was not surprising that much of the peroxide fed was decomposing in the feed piping and never made contact with the uranium mineralization. Under these conditions, indeed, peroxide showed no improvement over oxygen. This is, however, perhaps the most adverse environment encountered to date in US in-situ uranium leaching. In other operations, peroxide has proved advantageous where two factors have been exploited: 1.Despite the tendency for Hz02 to decompose to water and oxygen in the in-situ leach environment, its higher redox potential before decomposition has provided higher leach rates. In some formulations, overloading beyond the O2 saturation equivalent with peroxide has allowed higher oxidant levels, and higher leaching rates, when overloading with molecular O2, led to vapor-locking.

Jan 1, 1983

By D. A. Flanagan, P. B. Crawford

A study has been made of the feasibility of storing liquid meihane at low pressures in undergrohd caverns. Methane liquefies at — 258°F at atmospheric pressure. It is shown that the methane evaporation rates will rapidly decrease and cool the surrounding rock so that at the end of one month they would be between 50 Mcf/hr for a cavern of 25-ft radius and 700 Mcf/hr for a cavern of 100-ft radius. At the end of 10 years, the evaporation rates would be 18 and 100 Mcf/hr, respectively, for caverns of the same radii. The evaporation rates may be reduced by a factor of 2 to 10 by the application of insulation. The cost for mining the caverns is estimated to be $.75 to $1.25/Mcf of storage. This is substantially less than surface storage; it is believed to be safer and to result in lower maintenance, savings in space and savings in strategic materials. INTRODUCTION During the past few years, there has been an increasing interest in the economic feasibility of liquefying methane. Methane liquefaction is being considered for tanker transport; in addition, liquefaction is being reconsidered for shaving peak gas demands. Several articles have described natural gas liquefaction and the progress of the tanker in making trial runs from the United States to Great Britain to determine the feasibility of tanker transport. At the present time, pipelines are not designed to supply peak gas loads during extremely cold periods such as are oiten encountered in the North and Northeast. Gas is being stored in underground reservoirs en-route to its destination, but in many instances satisfactory storage in porous reservoirs has not been practical, especially along the Eastern seaboard where few petroleum reservoirs have been found. In England and other foreign countries, it is unlikely that satisfactory porous structures could be found, and it may be desirable to mine or excavate the rock to obtain storage. By storing the methane near the consumption point, product availability can be increased during periods of need. Methane liquefies at - 258OF at atmospheric pressure, and 1 cu ft of liquid methane will make about 600 cu ft of gaseous methane. Methane liquefaction for shaving peak gas demands was conducted in Cleveland, Ohio in the early 1940's. The Cleveland plant was designed to liquefy 4 MMcf/D of natural gas.' Development of several low-temperature processes has resulted in some improvements in the design of low-temperature storage vessels. Liquefied-methane storage vessels can be designed for various evaporation rates. The desired rate will depend on the reason for liquefaction. For shaving peak gas demands, the desired evaporation rates would be low, possibly between 0.2 and 1.5 per cent/day. For methane tanker transport and supply, the evaporation rates must be high and must equal the total daily rate of supply. It is apparent that the cost of storage facilities will depend on the objective and will vary appreciably. Surface storage of methane may cost between $3,000 and $20,000/MMscf, depending on the size and design. The failure of the surface storage vessels at Cleveland, with the resulting injury and loss of many lives, has caused considerable emphasis to be placed upon developing new, safer and improved methods of storing liquid methane. Storage of liquefied propane and butane in underground salt domes has been very satisfactory. The per-barrel cost of underground storage for any sizeable capacity is very small compared to the cost of surface storage for LPG mixtures. The larger the capacity of the underground storage, the less is the per-barrel cost. Numerous advantages appear to exist if liquefied methane could be satisfactorily stored underground, the principal advantages being increased safety, lower initial cost, lower maintenance cost, savings in space and savings in strategic materials. The purpose of this paper is to determine the feasibility of storing liquid methane underground in mined caverns. If a large spherical cavern is dug underground and is filled instantaneously with liquid methane, the surface of the sphere may approach the liquid methane temperature almost instantaneously. The temperature distribution for a sphere of radius r. in an infinite medium, initially at one temperature with spherical surface kept at T, from time t = 0, is given by Eq. 1.'

By D. R. Wieland, H. T. Kennedy

The object of the work reported here was to determine the content of elemental sulfur in gaseous methane, carbon dioxide, hydrogen sulfide, and in mixtures of these gases, at pressures and temperatwes encountered in natural gas reservoirs. Sulfur content at equilibrium is reported for pure methane, carbon dioxide, hydrogen sulfide, and on three binary mixtures of each of the three pairs of gases at pressures of 1,000, 2,000, 3,000, 4,000, 5,000 and 6,000 psia and at temperatures of 150°, 200° and 250°F. In addition, the sulfur content of three ternary mixtures at the same temperatures and pressures are reported. The results indicate that the sulfur content is higher in the gases at higher temperatures and pressures. The content is highest in hydrogen sulfide, intermediate in carbon dioxide and lowest in methane. INTRODUCTION At ordinary pressures and temperatures, the concentration of a nonvolatile material, such as sulfur in a gas at equilibrium, is a function of the vapor pressure of the material and is independent of the nature of the gas. As the pressure and temperature increase, the gas assumes some of the properties of liquids, including the power to dissolve other liquids and solids, to an extent dependent on the nature of both the gas and the material dissolved. At equilibrium, the content of sulfur in the several gases considered here may thus be considered as solubilities which are fixed for a given composition of gas, temperature and pressure. The study of the solubility of elemental sulfur in gases is of interest because sulfur is sometimes present in reservoirs producing natural gas and must be present in the vapor phase. Upon reduction of pressure and temperature, the sulfur precipitates from solution in the reservoir and in the tubing and fittings. Due to the fact that the greatest pressure drop in the reservoir is around the wellbore, the volume of sulfur precipitated will be greatest in this locality and can cause a substantial reduction in the permeability of the formation in this area. As the natural gas flows up the tubing string, the pressure and temperature are further decreased, causing further sulfur precipitation. If the volume of free sulfur is large and remedial measures are not taken, complete plugging of the tubing can occur. Knowledge of the content of sulfur as a function of composition, temperature and pressure will enable the operator to determine what changes in pressure and/or temperature are necessary to keep the sulfur in solution. If the sulfur solubility is large and can be controlled, the operator may desire to produce the well at a high wellhead pressure and temperature and reclaim the sulfur at the surface. The high solubility of sulfur in hydrogen sulfide (5.6 per cent by weight at 5,000 psi and 200°F) suggests the feasibility of recovering sulfur by injecting this gas into sulfur-bearing reservoirs and expanding the produced gas to precipitate out the sulfur. The present project was designed to measure the solubility of sulfur in carbon dioxide, methane, hydrogen sulfide and mixtures of these three gases at various temperatures and pressures. These gases were the major constituents of the gas well where plugging by sulfur was encountered. A study of the published literature shows no data on the solubility of sulfur in any gas, although Hannay and Hogarth,' in 1880 reported that it was soluble in carbon disdlide above the critical temperature of carbon disulfide. EQUIPMENT AND PROCEDURE Fig. 1 shows the layout of equipment employed to bring gases to equilibrium with sulfur and to isolate the quantity of sulfur contained in a known volume of gas. Pure gases were measured in the charging bomb at known temperatures and pressures and displaced into the reservoir bomb to make up the desired mixtures. When equilibrium was achieved

By R. J. Sandrea, S. M. Farouq Ali

The results of an experimental and theoretical study of the effects of rectilinear impermeable barriers and highly permeable channels on the sweep efficiency and conductivity of a five-spot network are presented. The study was carried out using Hele-Shaw, and conductive sheet analogs for both. normal and inverted flood patterns. A functional relationship was developed which provides quantitative prediction of the above parameters; i.e., sweep efficiency and conductivity, for the nonhomogeneous system. To evaluate this, it is necessary only to measure the interference modulus of the obstacle using a simple expression. The results indicate that impermeable barriers, with an interference modulus greater than zero, always cause a decrease in the conductivity and the sweep efficiency of the pattern. On the other hand, permeability channels do not have a significant effect on these parameters, unless they are located along the streamlines. The study also discusses a method of formulating the difference equations pertaining to the two-dimensional Laplacian in a system containing rectilinear barriers. In addition, a technique is described for tracing the progress of a free surface, using marked particles, in a stream containing discontinuities SuCh as those occurring at the extremities of barriers. The above computational scheme was used in a digital computer program to simulate results of the Hele-Shaw analog. The computed values of sweep efficiency mere in good agreement with the experimental values. INTRODUCTION At the present time waterflooding remains the most widely used technique for the secondary recovery of oil. Among the numerous factors that affect the performance of a waterflood, the areal sweep efficiency constitutes one of the most critical parameters. It is a function only of the geometric distribution of the wells for a homogeneous medium and unit mobility ratio. Geological considerations, however, such as the presence of sealing faults or solution channels in the rock matrix, would generally exert an adverse effect on the areal sweep efficiency of a waterflood network. Surprisingly little information is presently available on the quantitative implications of such permeability interferences. It was the object of this investigation to consider the presence of rectilinear, impermeable and highly permeable interferences in an otherwise homogeneous five-spot pattern, and to analyze their effects on the resulting sweep efficiency and pattern conductivity. The experimental phase of the study was conducted using two different models: the Hele-Shaw flow cell and the field plotter. For experimental convenience, the geometric boundaries of the five-spot were assumed to constitute streamlines even in the case where the interference was asymmetrically placed with respect to a repeated five-spot pattern. A second object of the study was to develop a computational scheme for numerically calculating the pressure distribution in a five-spot flow pattern containing rectilinear, impermeable barriers, and to use this distribution to compute the sweep efficiency. EXPERIMENTAL PROCEDURE HELE-SHAW ANALOG The Hele-Shaw analog was used both as a visual model and for determining the sweep efficiency for systems involving impermeable barriers. The physical limitations of this model as a means of experimentally studying the motion of a fluid in porous medium have been adequately discussed elsewhere.2,7 Likewise, it can be easily verified that the analogy remains valid when a portion of the region between the plates contains an obstacle of the same thickness as the distance between the plates. The model used in this study was composed of

THE BULOLO gold-dredging operations were described at some length by Frank W. Griffin at the Oct. 11 meeting of the San Francisco Section. This enterprise is on the island of New Guinea, about 1800 miles north of Sydney, Australia. The area in which the present placers exist originally was German territory but later the location was taken over by Australians. Gold was discovered in 1923 by prospectors who required 8 to 9 days to pack into the country over severe mountain trails. In 1926 a rich discovery was made on Edie Creek and in 1927 the deposit on the Bulolo River was discovered.

Jan 1, 1932