944 resultados para ES-SAGD. pressure drop. heavy oil. reservoir modeling and simulation
Resumo:
Diplomityön tavoitteena oli selvittää tapoja hallita voimalaitoksen höyryjärjestelmän paineensäätöä. Työ on tehty Neste Oilin Porvoon jalostamon energialaitoksen höyryjärjestelmästä. Energialaitos käsittää kaksi höyrykattilaa ja kaasuturbiinivoimalaitoksen, sekä öljynjalostamolla sijaitsevat höyrykattilan ja kaasuturbiinivoimalaitoksen. Höyryverkon häiriötilanteet aiheuttavat paineenlaskua höyryverkossa, jos jonkun höyryntuottajan tuotanto keskeytyy, tai paineennousua höyryverkossa, jos joku höyrynkuluttajista ei kuluta enää höyryä. Nykyisessä tilanteessa vakavin paineen laskua aiheuttava häiriötilanne olisi energialaitoksen höyrykattilan tuotannon keskeytyminen. Vakavin paineennousua aiheuttava häiriötilanne olisi suurimman kuluttajan,eli energialaitoksen vastapaineturbiinin, irtaantuminen höyryverkosta. Paineen hallinta voidaan toteuttaa oikeanlaisilla automaation ratkaisuilla, joita on mahdollista kehittää käyttökokemukseen perustuen järjestelmäpäivitysten yhteydessä. Häiriötilanteita voidaan yrittää ennakoida kunnossapidon keinoin niin, että vikaantuminen havaitaan ennen kuin vakava häiriö tapahtuu. Erityisesti kunnossapitostrategian kehittäminen auttaisi kohdistamaan kunnossapitotoimenpiteet oleellisimpiin kohteisiin kustannustehokkaasti. Lisäksi käyttöhenkilökunnan oikeanlaisella toiminnalla häiriön sattuessa voidaan häiriöiden aiheuttamia seuraamuksia tuotantoon ja turvallisuuteen lieventää. Käyttöhenkilökunnan toimintaa voidaan parantaa koulutuksella.
Resumo:
Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the Alluvial architecture above and below the maximum erosion depth of the modem channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits. The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements. This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.
Resumo:
Kaksifaasivirtauksen kuvaamiseen käytettävät mallit, ja menetelmät kaksifaasivirtauksen painehäviön määrittämiseksi kehittyvät yhä monimutkaisimmiksi. Höyrystinputkissa tapahtuvien painehäviöiden arvioinnin vaatiman laskennan suorittamiseksi tietokoneohjelman kehittäminen on välttämätöntä. Tässä työssä on kehitetty itsenäinen PC-ohjelma painehäviöiden arvioimiseksi pakotetulle konvektiovirtaukselle pystysuorissa höyrykattilan höyrystinputkissa. Veden ja vesihöyryn aineominaisuuksien laskentaan käytetään IAPWS-IF97 –yhtälökokoelmaa sekä muita tarvittavia IAPWS:n suosittelemia yhtälöitä. Höyrystinputkessa kulloinkin vallitsevan virtausmuodon määrittämiseen käytetään sovelluskelpoisia virtausmuotojen välisiä rajoja kuvaavia yhtälöitä. Ohjelmassa käytetään painehäviön määritykseen kirjallisuudessa julkaistuja yhtälöitä, virtausmuodosta riippuen, alijäähtyneelle virtaukselle, kupla-, tulppa- ja rengasvirtaukselle sekä tulistetun höyryn virtaukselle. Ohjelman laskemia painehäviöarvioita verrattiin kirjallisuudesta valittuihin mittaustuloksiin. Laskettujen painehäviöiden virhe vaihteli välillä –19.5 ja +23.9 %. Virheiden itseisarvojen keskiarvo oli 12.8 %.
Resumo:
This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.
Resumo:
Natriumhypokloriittia voidaan valmistaa kloorista ja lipeästä jatkuvatoimisessa absorberissa. Tässä työssä tutkittiin kokeellisesti, miten kaasun ja nesteen virtausnopeudet vaikuttavat täytekappalekolonnin tulvimiseen ja painehäviöön, kuinka nopeasti kloori absorboituu lipeään ja kuinka suuri hypokloriittiliuoksen kierrätys tarvitaan, ettei hypokloriitti ala hajota. Lisäksi luotiin matemaattinen malli, jolla voidaan mitoittaa jatkuvatoiminen vastavirtaperiaatteella toimiva natriumhypokloriittireaktori. Kloori–lipeäsysteemin havaittiin tulvivan suuremmilla virtausnopeuksilla kuin ilma–vesisysteemin. Tosin osa kloorista absorboituu jo ennen täytekappalekerrosta, minkä vuoksi kaasun todellinen virtausnopeus täytekappalekerroksen alaosassa on pienempi kuin mitattu arvo. Kolonnin painehäviö nousee erittäin jyrkästi tulvimispisteen läheisyydessä. Koska kloori absorboituu lähes täydellisesti ja vain kolonnin alaosa tulvii, voidaan kolonnia painehäviön kannalta operoida lähellä tulvimispistettä. Sekä mallinnuksen että koetulosten perusteella yli 99,99 % kloorista absorboituu koeolosuhteissa kahden metrin täytekappalekerroksessa. Nopea absorptio johtuu erittäin nopeasta, irreversiibelistä kloorin reaktiosta ja prosessille tyypillisestä natriumhydroksidikonversion rajoittamisesta alle 94 %:iin. Jotta varmistetaan, ettei hypokloriitti ala hajota, valmista hypokloriittiliuosta täytyy kierrättää kolonniin vähintään noin 4-kertainen määrä tuoreen lipeän syöttömäärän nähden.
Resumo:
Kandidaatintyössä luotiin CFD-malli mallintamaan jäähdytevirtausta kuulakekoreaktorin sydämessä käyttämällä Ansys Fluent -ohjelmaa. Mallin avulla tarkasteltiin virtauksen käyttäymistä ja painehäviötä ja saatuja tuloksia verrattiin aiempiin tutkimuksiin. Kandidaatin työssä on myös kerrottu mallintamisen etenemisestä ja laskentateoriaa.
Resumo:
Venturi scrubbers are high efficiency gas cleaners in which suspended particles are removed from gas streams by droplets formed by liquid atomisation, usually in the venturi throat. The size of the droplets formed is of fundamental importance to the performance of the equipment, both in terms of pressure drop and collection efficiency. In this study, drop sizes in a cylindrical laboratory scale venturi scrubber were measured using a laser diffraction technique. Gas velocity and liquid to gas ratios varied from 50 to 90 m/s and 0.5 to 2.0 l/m3, respectively. Water was inserted as perpendicular jets at the beginning of the throat. Measurements were performed at three positions: two located along the throat, and the last one at the end of the diffuser. The data presented here are a typical example of pneumatic atomisation and can be relevant to other industrial applications such as combustion and engine technology. Finally, results are compared to available correlations and the validity of these equations is discussed.
Resumo:
Microreactors have proven to be versatile tools for process intensification. Over recent decades, they have increasingly been used for product and process development in chemical industries. Enhanced heat and mass transfer in the reactors due to the extremely high surfacearea- to-volume ratio and interfacial area allow chemical processes to be operated at extreme conditions. Safety is improved by the small holdup volume of the reactors and effective control of pressure and temperature. Hydrogen peroxide is a powerful green oxidant that is used in a wide range of industries. Reduction and auto-oxidation of anthraquinones is currently the main process for hydrogen peroxide production. Direct synthesis is a green alternative and has potential for on-site production. However, there are two limitations: safety concerns because of the explosive gas mixture produced and low selectivity of the process. The aim of this thesis was to develop a process for direct synthesis of hydrogen peroxide utilizing microreactor technology. Experimental and numerical approaches were applied for development of the microreactor. Development of a novel microreactor was commenced by studying the hydrodynamics and mass transfer in prototype microreactor plates. The prototypes were designed and fabricated with the assistance of CFD modeling to optimize the shape and size of the microstructure. Empirical correlations for the mass transfer coefficient were derived. The pressure drop in micro T-mixers was investigated experimentally and numerically. Correlations describing the friction factor for different flow regimes were developed and predicted values were in good agreement with experimental results. Experimental studies were conducted to develop a highly active and selective catalyst with a proper form for the microreactor. Pd catalysts supported on activated carbon cloths were prepared by different treatments during the catalyst preparation. A variety of characterization methods were used for catalyst investigation. The surface chemistry of the support and the oxidation state of the metallic phase in the catalyst play important roles in catalyst activity and selectivity for the direct synthesis. The direct synthesis of hydrogen peroxide was investigated in a bench-scale continuous process using the novel microreactor developed. The microreactor was fabricated based on the hydrodynamic and mass transfer studies and provided a high interfacial area and high mass transfer coefficient. The catalysts were prepared under optimum treatment conditions. The direct synthesis was conducted at various conditions. The thesis represents a step towards a commercially viable direct synthesis. The focus is on the two main challenges: mitigating the safety problem by utilization of microprocess technology and improving the selectivity by catalyst development.
Resumo:
Cardiac hypertrophy that accompanies hypertension seems to be a phenomenon of multifactorial origin whose development does not seem to depend on an increased pressure load alone, but also on local growth factors and cardioadrenergic activity. The aim of the present study was to determine if sympathetic renal denervation and its effects on arterial pressure level can prevent cardiac hypertrophy and if it can also delay the onset and attenuate the severity of deoxycorticosterone acetate (DOCA)-salt hypertension. DOCA-salt treatment was initiated in rats seven days after uninephrectomy and contralateral renal denervation or sham renal denervation. DOCA (15 mg/kg, sc) or vehicle (soybean oil, 0.25 ml per animal) was administered twice a week for two weeks. Rats treated with DOCA or vehicle (control) were provided drinking water containing 1% NaCl and 0.03% KCl. At the end of the treatment period, mean arterial pressure (MAP) and heart rate measurements were made in conscious animals. Under ether anesthesia, the heart was removed and the right and left ventricles (including the septum) were separated and weighed. DOCA-salt treatment produced a significant increase in left ventricular weight/body weight (LVW/BW) ratio (2.44 ± 0.09 mg/g) and right ventricular weight/body weight (RVW/BW) ratio (0.53 ± 0.01 mg/g) compared to control (1.92 ± 0.04 and 0.48 ± 0.01 mg/g, respectively) rats. MAP was significantly higher (39%) in DOCA-salt rats. Renal denervation prevented (P>0.05) the development of hypertension in DOCA-salt rats but did not prevent the increase in LVW/BW (2.27 ± 0.03 mg/g) and RVW/BW (0.52 ± 0.01 mg/g). We have shown that the increase in arterial pressure level is not responsible for cardiac hypertrophy, which may be more related to other events associated with DOCA-salt hypertension, such as an increase in cardiac sympathetic activity
Resumo:
The objective of this master’s thesis was to design and simulate a wind powered hydraulic heating system that can operate independently in remote places where the use of electricity is not possible. Components for the system were to be selected in such a way that the conditions for manufacture, use and economic viability are the as good as possible. Savonius rotor was chosen for wind turbine, due to its low cut in speed and robust design. Savonius rotor produces kinetic energy in wide wind speed range and it can withstand high wind gusts. Radial piston pump was chosen for the flow source of the hydraulic heater. Pump type was selected due to its characteristics in low rotation speeds and high efficiency. Volume flow from the pump is passed through the throttle orifice. Pressure drop over the orifice causes the hydraulic oil to heat up and, thus, creating thermal energy. Thermal energy in the oil is led to radiator where it conducts heat to the environment. The hydraulic heating system was simulated. For this purpose a mathematical models of chosen components were created. In simulation wind data gathered by Finnish meteorological institute for 167 hours was used as input. The highest produced power was achieved by changing the orifice diameter so that the rotor tip speed ratio follows the power curve. This is not possible to achieve without using electricity. Thus, for the orifice diameter only one, the optimal value was defined. Results from the simulation were compared with investment calculations. Different parameters effecting the investment profitability were altered in sensitivity analyses in order to define the points of investment profitability. Investment was found to be profitable only with high average wind speeds.
Resumo:
Koalesenssi on ilmiö, jossa dispergoidun faasin pisarat pyrkivät muodostamaan suurempia pisaroita kunnes erotettava faasi muodostuu. Koalesenssi tapahtuu kolmessa päävaiheessa, jotka ovat lähestyminen, kiinnittyminen ja irrotus. Lähestymiseen vaikuttavat mekanismit ovat muuan muassa sieppaus, diffuusio, törmäysvaikutus, sedimentaatio, sähköiset repul-siovoimat ja van der Waalsin voimat. Kiinnittymisvaiheessa dispergoidun faasin pisarat syrjäyttävät väliaineen nestekalvon samalla kostuttaen väliaineen pinnan. Irrotusvaiheessa pisaran hydrodynaaminen voima voittaa pisaran ja väliaineen välisen adheesiovoiman. Koalesenssin tehokkuuteen vaikuttavat useat eri parametrit kuten virtausnopeus, pedin ominaisuudet, väliaineen ominaisuudet sekä emulsion ominaisuudet. Nämä kaikki asiat tulee ottaa huomioon koalesenssisuodatuksen suunnittelussa. Koalesenssisuodatus lukeutuu syväsuodatusmenetelmiin, jotka on ollut käytössä jo yli 100 vuotta. Koalesenssisuodatusmenetelmä on tehokas menetelmä pienten pisaroiden erottami-seen. Menetelmää käytetään esimerkiksi öljyisten jätevesien puhdistuksessa. Teollisen öljyn syväsuodatuksen etuihin kuuluu muun muassa sen kompakti koko, alhaisemmat käyt-tökustannukset, korkea erotusaste, kyky erotella pienetkin pisarat sekä helppo operointi, automatisointi ja huolto. Suurin haittapuoli on kuitenkin väliaineen tukkeutuminen, joten prosessi vaatii puhdistuksen tai väliaineen uusimisen. Tämän kandidaatintyön tarkoituksena oli koota kirjallisuustyö öljyn koalesenssisuodatuk-sesta. Työssä kartoitettiin koalesenssisuodatuksen lähtökohdat, teoria, tärkeimmät teolli-suuden sovellukset sekä väliaineet.
Resumo:
Rats fed a high-fructose diet represent an animal model for insulin resistance and hypertension. We recently showed that a high-fructose diet containing vegetable oil but a normal sodium/potassium ratio induced mild insulin resistance with decreased insulin receptor substrate-1 tyrosine phosphorylation in the liver and muscle of normal rats. In the present study, we examined the mean blood pressure, serum lipid levels and insulin sensitivity by estimating in vivo insulin activity using the 15-min intravenous insulin tolerance test (ITT, 0.5 ml of 6 µg insulin, iv) followed by calculation of the rate constant for plasma glucose disappearance (Kitt) in male Wistar-Hannover rats (110-130 g) randomly divided into four diet groups: control, 1:3 sodium/potassium ratio (R Na:K) diet (C 1:3 R Na:K); control, 1:1 sodium/potassium ratio diet (CNa 1:1 R Na:K); high-fructose, 1:3 sodium/potassium ratio diet (F 1:3 R Na:K), and high-fructose, 1:1 sodium/potassium ratio diet (FNa 1:1 R Na:K) for 28 days. The change in R Na:K for the control and high-fructose diets had no effect on insulin sensitivity measured by ITT. In contrast, the 1:1 R Na:K increased blood pressure in rats receiving the control and high-fructose diets from 117 ± 3 and 118 ± 3 mmHg to 141 ± 4 and 132 ± 4 mmHg (P<0.05), respectively. Triacylglycerol levels were higher in both groups treated with a high-fructose diet when compared to controls (C 1:3 R Na:K: 1.2 ± 0.1 mmol/l vs F 1:3 R Na:K: 2.3 ± 0.4 mmol/l and CNa 1:1 R Na:K: 1.2 ± 0.2 mmol/l vs FNa 1:1 R Na:K: 2.6 ± 0.4 mmol/l, P<0.05). These data suggest that fructose alone does not induce hyperinsulinemia or hypertension in rats fed a normal R Na:K diet, whereas an elevation of sodium in the diet may contribute to the elevated blood pressure in this animal model.
Resumo:
En Colombia las actividades de exploración, explotación, transporte y procesamiento de hidrocarburos que se vienen realizando desde comienzos del siglo XX son responsables de grandes procesos de transformación del territorio y de degradación de los ecosistemas en los que se realizan. Estos procesos han impactado negativamente la seguridad de las comunidades indígenas poniendo en riesgo su cultura y en algunos casos su existencia misma. Aunque históricamente los derechos de estas poblaciones frente a la explotación petrolera, y minera en general, han cambiado y su autonomía e integridad es protegida por la Constitución de 1991, las comunidades siguen teniendo una alta vulnerabilidad frente a la intervención de los ecosistema que habitan.La degradación ambiental producida directamente por las actividades petroleras y por los procesos de colonización que estas impulsan se constituye en una amenaza a la seguridad de las comunidades, cuyos territorios y recursos de subsistencia se ven disminuidos. La colonización y la presión sobre los recursos naturales que esta produce son motivadas principalmente por la pobreza de poblaciones campesinas que buscan nuevas tierras para habitar, a su vez estos dos procesos son causa de degradación ambiental que empobrece a las comunidades étnicas debido a que afecta sus fuentes de sustento, situación que genera inseguridad para los indígenas. Adicionalmente, la degradación ambiental y la disminución de los territorios ponen en riesgo la cultura de estos grupos humanos, pues afecta sus valores, tradiciones, autoridades y, en general, su forma de vida lo que constituye una amenaza a su seguridad.-----In Colombia, the exploration, exploitation, transport, and processing of hydrocarbons since the beginning of the 20th century have caused great territory transformations and ecosystem degradation. These processes have impacted adversely the indigenous communities security, exposing their culture and, in some cases, their existence itself. Even though, facing oil and, in general, mineral exploitation, the rights of this population have changed historically and their autonomy and integrity is protected by the 1991 Constitution, the communities are still highly vulnerable to the intervention on the ecosystem they inhabit.The environmental degradation directly arisen from the oil exploitation activities and the colonization they have driven, has become a threat to the security of the communities whose territories and subsistence resources have been reduced. Colonization and the resulting natural resource pressure are mainly caused by the poverty of the country population that seek new lands to occupy and these two facts cause in turn the environmental degradation that impoverish the ethnic communities by affecting their living sources, thereby causing insecurity to the indigenous population. In addition, environmental degradation and territory reduction risk these human groups’ culture by impacting their values, tradition, authorities and, in general, their way of living, and therefore turn into a threat to their security.
Resumo:
Atmospheric downwelling longwave radiation is an important component of the terrestrial energy budget; since it is strongly related with the greenhouse effect, it remarkably affects the climate. In this study, I evaluate the estimation of the downwelling longwave irradiance at the terrestrial surface for cloudless and overcast conditions using a one-dimensional radiative transfer model (RTM), specifically the Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART). The calculations performed by using this model were compared with pyrgeometer measurements at three different European places: Girona (NE of the Iberian Peninsula), Payerne (in the East of Switzerland), and Heselbach (in the Black Forest, Germany). Several studies of sensitivity based on the radiative transfer model have shown that special attention on the input of temperature and water content profiles must be held for cloudless sky conditions; for overcast conditions, similar sensitivity studies have shown that, besides the atmospheric profiles, the cloud base height is very relevant, at least for optically thick clouds. Also, the estimation of DLR in places where radiosoundings are not available is explored, either by using the atmospheric profiles spatially interpolated from the gridded analysis data provided by European Centre of Medium-Range Weather Forecast (ECMWF), or by applying a real radiosounding of a nearby site. Calculations have been compared with measurements at all sites. During cloudless sky conditions, when radiosoundings were available, calculations show differences with measurements of -2.7 ± 3.4 Wm-2 (Payerne). While no in situ radiosoundings are available, differences between modeling and measurements were about 0.3 ± 9.4 Wm-2 (Girona). During overcast sky conditions, when in situ radiosoundings and cloud properties (derived from an algorithm that uses spectral infrared and microwave ground based measurements) were available (Black Forest), calculations show differences with measurements of -0.28 ± 2.52 Wm2. When using atmospheric profiles from the ECMWF and fixed values of liquid water path and droplet effective radius (Girona) calculations show differences with measurements of 4.0 ± 2.5 Wm2. For all analyzed sky conditions, it has been confirmed that estimations from radiative transfer modeling are remarkably better than those obtained by simple parameterizations of atmospheric emissivity.
Resumo:
We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a two-dimensional reservoir in an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting or extracting fluid. Numerical solution of this problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l. This is a situation which occurs frequently in the application to oil reservoir recovery. Under the assumption that epsilon=h/l<<1, we show that the pressure field varies only in the horizontal direction away from the wells (the outer region). We construct two-term asymptotic expansions in epsilon in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive analytical expressions for all significant process quantities. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the reservoir, epsilon, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighborhood of wells and away from wells.