Freddy Humberto Escobar 1, Yu Long Zhao 2 and Mashhad Fahes 3 1.

Size: px
Start display at page:

Download "Freddy Humberto Escobar 1, Yu Long Zhao 2 and Mashhad Fahes 3 1."

Transcription

1 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. CHARACTERIZATION OF THE NATURALLY FRACTURED RESERVOIR PARAMETERS IN INFINITE-CONDUCTIVITY HYDRAULICALLY- FRACTURED VERTICAL WELLS BY TRANSIENT PRESSURE ANALYSIS Freddy Humbero Eobar 1, Yu Long Zhao and Mashhad Fahes 3 1 Universidad Surcolombiana/CENIGAA, Avenida Pasrana - Cra 1, Neiva, Huila, Colombia Sae Key Laboraory o Oil and Gas Reservoir Geology and Eploiaion, Souhwes Peroleum Universiy, Xindu Sree, Xindu disric, Chendu, Sichuan, P. R. China 3 The Universiy o Oklahoma, Mewbourne School o Peroleum and Geological Engineering, E. Boyd S. SEC Rm, Norman, Ok, USA eobar@uo.edu.co ABSTRACT I has become common o hydraulically racure a naurally racured ormaion o increase he well s producion poenial. Since he mass ranser beween racured nework and hydraulic racure is much higher han ha rom mari o racures, he hydraulic racure-racure nework inerporosiy low parameer is much higher ha ha o mari-racure nework. As a resul, he ransiion period behavior rom naurally racured o homogeneo akes place beore radial low regime during he early bilinear, linear or ellipical (birradial) low regimes. The purpose o his paper is o provide epressions by boh convenional analysis and TDS echnique or characerizing he naurally racured parameers when he ransiion period inerrups he response o an ininie-conduciviy racure. The developed epressions or boh mehodologies were saisacorily esed wih simulaed eamples. Keywords: pressure ransien analysis, inerporosiy low parameer, TDS echnique, ellipical low, naurally racured reservoirs. 1. INTRODUCTION Even wells in naurally racured ormaions are being subjeced o hydraulic racuring since he bulk permeabiliy is low and well produciviy is no as high as epeced. Typical such cases have been seen in some wells o he Orinoco basin oohill in Colombia where he wells were racured j aer being drilled. Tiab (1994) was he irs o ideniy he ellipical low regime in pressure es daa o ininie-conduciviy hydraulically racured verical wells. He called i birradial low, provided he governing pressure and pressure derivaive equaions and implemened he TDS echnique or is characerizaion. This low regime can also be seen in horizonal wells beore he lae radial low regime shows up and is due mainly o horizonal anisoropy. The las model or he ellipical low in horizonal wells was presened by Marinez, Eobar and Bonilla (01). The original pressure governing model or ellipical low inroduced by Tiab (1994) involves he drainage area, which may aec he calculaions i he pressure es is oo shor or he developmen o he lae pseudoseady-sae period. An eperienced er o he TDS echnique may deal wih ha problem ing he poin o inersecion beween he radial low and he birradial low regime lines, Tiab (1994). However, o overcome his problem, Eobar, Bonilla and Ghisays-Ruiz (014) ormulaed a new ellipical/birradial low model or verical wells in eiher racured or unracured ormaions. Alhough Tiab s model, Tiab (1994), works perecly, i requires he knowledge o well-drainage area or esimaing he hal-racure lengh. However, someimes his condiion canno be me becae he pressure es needs o be run long enough or he developmen o he lae ime pseudoseady-sae period or he deeraion o he drainage area. The new model presened by Eobar e al. (014) ecludes he reservoir drainage area and slighly modiies Tiab s model, Tiab (1994), o accoun or naurally-racured double-porosiy sysems. Tiab and Beam (007) presened a pracical inerpreaion o he pressure behavior o a inieconduciviy hydraulically racured verical well locaed in a naurally racured reservoir. The inerpreaion is based on analyical equaions derived o deere permeabiliy, racure sorage capaciy raio, inerporosiy low coeicien, skin and wellbore sorage rom he pressure derivaive plo wihou ing ype-curve maching echnique. In oher words, hey implemened he TDS echnique, Tiab (1993), or such sysems. Par o he work presened by Eobar, Marinez and Monealegre (009) was oced on he implemenaion o convenional analysis or he work presened by Tiab and Beam (007). In heir work, Tiab and Beam (007) developed epressions or bilinear and linear low regimes bu ecluded he presence o he ellipical or birradial low regime. Thereaer, his work consiues an eension o 535

2 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. he research o Tiab and Beam (007) or ininieconduciviy racures where early birradial low is seen. This implemenaion was saisacorily esed wih synheic eamples.. MATHEMATICAL FORMULATION.1. Governing equaion The dimensionless pressure behavior in he Laplace domain or a hydraulically racured well in double-porosiy reservoirs was inroduced by Ozkan and Raghavan (1988): 1 s s ( s)[1 0.73) PD () s K 0 0( ) d (1) s s( s) khp PD (6) 141.qB D kh(* P ') * PD' (7) 141.qB Finally, he dimensionless racure conduciviy inroduced by Cinco-Ley, Samaniego and Doguez (1976) is deined as: k w C D (8) k For a pseudoseady-sae mari low model, he uncion o sorage capaciy and inerporosiy low parameer is given by: (1 ) () s () (1 ) s 9/5 D *P D ' * Eobar e al. (014) provided he ollowing pressure and pressure derivaive governing equaions or ellipical low regime: 5 D PD 9 6 D D* PD' 6 When = 1, Equaions (3) and (4) accoun or homogeneo reservoirs. For he case o naurallyracured ormaions, =, which is he dimensionless soraiviy coeicien. Consequenly, or he purpose o his work, will be replaced by he soraiviy raio... Dimensionless quaniies The dimensionless ime, based on hal-racure lengh, is given below: D (3) (4) k (5) c The dimensionless pressure and pressure derivaive parameers or oil reservoirs are given by: 1.E-0 [( 9/5 * PD ')* ] 0.03 D 1.E-06 1.E-05 1.E-04 1.E-03 1.E-0 1.E+0 1.E+03 1.E+04 1.E+05 9/5 * D Figure-1. Eec o he inerporosiy low parameer on he pressure and pressure derivaive during he ellipical low regime, = TDS echnique By plugging he dimensionless parameers in Equaions (3) and (4) and solving or he hal-racure lengh, we obain: qb hpbr1 c k qb h (* P') BR1 c k (9) (10) Figures-1 hrough 4 were generaed ing he analyical soluion provided by Ozkan and Raghavan (1988), Equaion (1), which considers low beween hydraulic racure and racure nework. As seen in Figure-1, once he ransiion period vanishes, he reservoir behaves as i i were homogeneo. In his case, Equaions (9) and (10) become: 5353

3 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. qb hpbr1 c k (11) qb.03 (* P') h ( c ) k (14) qb h (* P') BR1 c k (1) Where BR1 corresponds o he sraigh line drawn on he second ellipical low regime which occurs aer he ransiion period. However, beore he ransiion period eiher ellipical low or linear low may ake place. I linear low occurs a early ime, hen Equaions (9) and (10) m be replaced by hese wo epressions developed by Tiab and Beam (007); Fracure conduciviy can be ound ing an epression presened by Tiab (003); kw k e r s w (15) I is observed in Figure-1 ha he imum poin o he pressure derivaive imes he inerporosiy low parameer o he power is a consan: qb Ph ( c ) k (13) 9/5 ( D * D ') 0.03 P (16) The ollowing epression o deere he inerporosiy low parameer is obained once Equaion (7) is subsiued in Equaion (16): 9/5 D *P D ' * 9/5 [P & *P ' ]* D D D 1.E-0 [( 9/5 * PD ')* ] ma 0. D ( * / ) 0.475hr D ma E-03 1.E-0 1.E+0 1.E+03 1.E+04 1.E+05 1.E+06 * D Figure-. Correlaion o he maimum and imum pressure derivaive coordinaes during ellipical low period, = E+0 1.E-0 Heerogeneo behavior E-04 1.E-03 1.E-0 1.E+0 1.E+03 1.E+04 * D Homogeneo behavior Uni-slope line during ransiion period Figure-3. Eec o he soraiviy raio on he pressure and pressure derivaive during he ellipical low regime, = 100. qb kh(* P ') 5/9 (17) Figure- conains a log-log plo o pressure and pressure derivaive boh muliplied by he inerporosiy low parameer raised o he power 9/5 (or ) vers dimensionless ime muliplied by inerporosiy low parameer and divided by he soraiviy raio. During he early heerogeneo behavior, all curves uniy no maer he values o soraiviy raio. A maimum poin is displayed beore he ransiion period which coordinaes are given by: [( * P ') ] 0. (18) 9/5 D D ma D ma (19) Once Equaion (7) is subsiued ino Equaion (18), we obain: qb kh(* P ') ma 5/9 (0) Subsiuing Equaion (5) ino Equaion (19) yields wo epressions o ind eiher inerporosiy low parameer or soraiviy raio one as a uncion o he oher so; 5354

4 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. k ma (1) c is inerruped by he ransiion period, ollowed by he ellipical or birradial low. In such cases, dividing Equaion (14) by Equaion (1) leads o: c () k ma h * P' BR1 k qb * P ' (8) Figure-3 was prepared o uniy he uni-slope line showing up during he ransiion period. Such sraigh line was ied o: ( * P ') (3) 16/5 D D D Aer subsiuing he dimensionless parameer and solving or he inerporosiy low parameer, i yields: c h (* P') qb 5/16 (4) Equaion (4) works ing he coordinaes o any poin along he uni-slope line. During radial low regime, he value o he dimensionless pressure derivaive is one hal. Equaing his value o Equaion (3) and solving or he inerporosiy low parameer, we ge: The poin o inersecion beween he early birradial line, Equaion (4), and he uni-slope line during he ransiion period, Equaion (3), leads o he ollowing epression: 1 k BRi c 9/16 (9) Tiab and Beam (007) inroduced he governing pressure derivaive equaion or early linear low: D o: 1 D * PD' (30) Which, in combinaion wih Equaion (3), leads c k i 5/16 (5) c (31) k 3/5 Li Where i is he inersecion poin o he uni-slope line and he radial low regime. From observaion o he imum poin in Figure-, an epression wih a correlaion coeicien o was obained: 9/5 P D * P 9/5 Dpss ma ma (6) Dividing Equaion (10) by Equaion (1) and solving or he soraiviy raio, we ge:.7778 * P' BR1 * P' BR1 (7) Equaion (7) is valid whenever birradial/ellipical low eiss a boh side o he ransiion period. However, in ininie-conduciviy racure, i is possible ha he linear low occurs a early ime, and hen 1.E-0 1.E+0 1.E+03 1.E+04 1.E+05 1.E+06 * D Figure-4. Dimensionless pressure imes he inerporosiy low parameer o he power vers dimensionless ime muliplied by /, = Convenional analysis As seen in Figure-4, here is a poin in which he dimensionless pressure converges o a value during he ransiion period. From here we can read ha: 9/5 P (3) Dpss 5355

5 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. Subsiuing Equaion (3) in and solving or he inerporosiy low parameer yields: qb khppss 5/9 (33) As seen in Figure-4, once he ransiion period vanishes a line wih slope o 0.4 is developed. Draw such line and draw a parallel line hroughou he poins j beore he developmen o he ransiion period. Read wo poins in boh lines a he same ime value and ake is raio, which is ed in he ollowing iing equaion ha has a correlaion coeicien o The soraiviy raio is ound rom such separaion: Phigh Plow 5e (34) 3. SYNTHETIC EXAMPLES Table-1. Inpu daa or eamples. Eamples Parameer 1 k, md 5 1 h, 100 0, % 5 10 B, rb/stb , cp c, 1/psi q, BPD 30 50, /9 (50)(1.5)(1.) (5)(100)(3.) 3.1. Synheic eample 1 Figure-5 repors simulaed pressure drop vers ime daa obained wih he inormaion given on he second column o Table-1. I is required o ind he naurally-racured reservoir parameers. Soluion by Convenional Analysis In Figure-5 a horizonal line was drawn going hrough he ransiion poins. The ollowing value o pressure drop was read: P pss = 3. psi Using he above value in Equaion (33), i yields: Then, a 0.4-slope line is drawn aer he ransiion period. Anoher parallel line goes j aer he iniiaion o he ransiion period so he ampliude can be esimaed. Two pressure drop values are read a he same ime rom boh sraigh lines, so: P high = psi P low = 6 psi Equaion (34) will lead o he esimaion o he soraiviy raio: e E+0 P, psi P pss 3. psi ( P) psi high ( P) 6psi low The hal-racure lengh is ound rom he slope o a Caresian plo (m L= psi/hr 0.5 ) o pressure vers he square roo o ime, Figure-7, ing Equaion (A.7): 30(1.) (100)(360.8) (0.05)(0.0001)(5)(0.01) Soluion by TDS Technique 1.E-06 1.E-05 1.E-04 1.E-03 1.E-0, hr Figure-5. Pressure drop vers ime log-log plo or eample

6 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. P & *P ', psi 1.E+03 1.E+0 ( * P') 0.75 psi ma ( * P') psi ma hr 1.E-0 1.E-06 1.E-05 1.E-04 1.E-03 1.E-0, hr ( * P') 190 psi i hr Li 0.06 hr 1. hr BRi hr (* P') 8.5psi ( * P') psi hr Figure-6. Pressure and pressure derivaive vers ime log-log plo or eample 1. 5 BR1 The inerporosiy low parameer is esimaed wih Equaions (17), (0), (), (4) and (5): 5/9 30(1.5)(1.) (100)(0.085) 5/9 30(1.5)(1.) (100)(0.75) (0.1)(0.05)(1.5)(0.0001)(10 ) ( ) 5/16 (0.05)(0.0001)(10 )(100)(0.305) (30)(1.)(0.003) P, psi m psi/ hr L 5/16 (0.05)(0.0001)(10 ) (0.06) The soraiviy raio is esimaed wih Equaions (6), (8), (9) and (31): , hr Figure-7. Caresian plo o pressure vers he square roo o ime or eample 1. The ollowing inormaion was read rom he pressure and pressure derivaive log-log plo provided in Figure-6. ma = hr (*P ) ma= 0.75 psi (*P ) = psi (*P ) = psi BRi = hr (*P ) BR1= 8.5 psi Using Equaion (14); (*P ) BR1= 8.5 psi (*P ) = 190 psi = hr i = 0.06 hr Li = 1. hr = hr 30(1.) (190)(100) (0.0001)(5) The hal-racure lengh is also ound wih Equaion (1); 30(1.) (8.5) 0.05(0.0001) (8.5) (1.) (190) (50)(0.015) (0.05)(1.5)(0.0001)(10 ) (0.05)(1.5)(0.0001)(10 ) 3/5 9/ (5)(50 )(1.) Synheic eample Figure-8 presens synheic pressure drop vers ime daa obained wih he inormaion given on he hird column o Table-1. I is required o ind he naurallyracured reservoir parameers. Soluion by Convenional Analysis In Figure-8 a horizonal line was drawn going hrough he ransiion poins. The ollowing value o pressure drop was read: P pss = 184 psi P high = 1600 psi P low = 373 psi 5357

7 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. m L = psi/hr 0.5 read rom Figure-10. Equaion (33) leads o an esimaion o a value or he inerporosiy low parameer and Equaion (34) provides a soraiviy raio o A hal-racure lengh o 35.1 is ound wih Equaion (A.7) 1.E+04 The ollowing inormaion was read rom he pressure and pressure derivaive log-log plo provided in Figure-9. ma = hr (*P ) ma= 44 psi (*P ) =.4 psi BRi = hr (*P ) = 6.53 psi (*P ) BR1=105 psi (*P ) = 1400 psi = 0.09 hr Li = 340 hr = 0.04 hr P, psi 1.E+03 P pss 184 psi ( P) high 1600 psi ( P) 373psi low The calculaions were perormed in he same way as in eample 1 and a summary o he resuls is given in Table-3. 1.E+0 Table-. Resuls rom convenional analysis. 1.E-04 1.E-03 1.E-0 1.E+0 1.E+03 Figure-8. Pressure drop vers ime log-log plo or eample 1., hr Parameer Eample 1 Inpu Eample Inpu Equaion Number, A E+05 Li 340 hr Table-3. Resuls rom TDS echnique. P & *P ', psi 1.E+04 1.E+03 (* P') 44psi 1.E+0 ma ( * P') 1400 psi ( * P') BR1 105 psi i 7hr (* P').4psi ma hr 0.09 hr ( * P') 6.53 psi 0.04 hr 1.E-04 1.E-03 1.E-0 1.E+0 1.E+03, hr BRi hr Figure-9. Pressure and pressure derivaive vers ime log-log plo or eample. P, psi m psi/ hr L Parameer Eample 1 Eample, Equaion Number Inpu Inpu Inpu , hr Figure-10. Caresian plo o pressure vers he square roo o ime or eample. Soluion by TDS Technique 4. COMMENTS ON THE RESULTS The provided eamples in his work demonsrae he applicabiliy o he developed mahemaical epressions. I is worh o red ha he esimaion o he naurally racured parameers is very sensiive, and as a resul, deviaion errors o one order o magniude are allowed. The mehod provides saisacory resuls especially or he inerporosiy low parameer. I urher 5358

8 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. accuracy is required, hen, e o compuer program is highly recommended. These eamples were perormed by hand and he values orm he plos were read by eye, consequenly, deviaion rom inpu values is epeced. Equaions or gas low are given in Appendi B and he equaions or he esimaion o hal-racure lengh are given in Appendi A. CONCLUSIONS a) New mahemaical epressions or boh TDS echnique and convenional analysis were inroduced and successully esed wih synheic cases or he ull characerizaion o pressure ess conduced in naurally-racured reservoirs whose wells have been subjeced o hydraulic racuring. b) For he case o he TDs echnique, ive epressions or he esimaion o he inerporosiy low parameer were developed. Four equaions or he esimaion o he soraiviy raio were also inroduced in his work. ACKNOWLEDGEMENTS The auhors hank Universidad Surcolombiana, he Universiy o Oklahoma and he Naional Science Fund or Disinguished Young Scholars o China (Gran No ) and he Naional Program on Key Basic Research Projec (973 Program, Gran No. 011CB01005), and he Naural Science Foundaion o China (Gran No ). Nomenclaure A Draining area, B b C D c h k k w m P P w q q g r w s Oil volume acor, rb/stb Shores disance rom a laeral bounday o a well, Dimensionless racure conduciviy Compressibiliy, 1/psi Formaion hickness, Formaion compressibiliy, md Fracure conduciviy, md- Slope Pressure, psi Well-lowing pressure, psi Oil low rae, STB/D Gas low rae, MSCF/D Wellbore radi, Hal-racure lengh, Skin acor. Laplace parameer Tes ime, hr T (* P ) Reservoir emperaure, R Pressure derivaive, psi ( D*P D ) Dimensionless pressure derivaive Greek Change Porosiy, raccion Mari-racure nework inerporosiy low parameer Fracure nework-hydraulic racure inerporosiy low parameer μ Viosiy, cp Suies BR BR BR BR1 BRi D D ell Li ma w i Variable o ideniy homogeneo (=1) or heerogeneo (=) reservoirs Dimensionless soraiviy coeicien Birradial low in he homogenoe región (second birradial) Second birradial low a he ime o 1 hr Birradial Birradial a 1 hr Inersec o birradial and uni-slope lines Dimensionless Dimensionless based on hal-racured lengh Ellipical (same as birradial) Linear a 1 hr Inersec o linear and uni-slope lines Maimum poin Minimum poin Well Time REFERENCES Uni-slope during heerogeneo o homogeneo ransiion period Inersec o radial al uni-slope lines Agarwal G Real Gas Pseudo-ime a New Funcion or Pressure Buildup Analysis o MHF Gas Wells. Paper SPE 879 presened a he 54 h echnical conerence and 5359

9 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. ehibiion o he Sociey o Peroleum Engineers o AIME held in Las Vegas, NV, Sepember 3-6. Cinco-Ley H., Samaniego F. and Doguez N Transien Pressure Behavior or a Well wih a Finiy- Conduciviy Verical Fracure. Paper SPE 6014 presened a he SPE-AIME 51 s Annual Fall Technical Conerence and Ehibiion, held in New Orleans, LA, Ocober 3-6. Eobar F.H., Monealegre-M. M. and Marínez J.A Convenional Pressure Analysis or Naurally Fracured Reservoirs wih Transiion Period beore and Aer he Radial Flow Regime. CT and F - Ciencia, Tecnología y Fuuro. 3(5): Eobar F.H., Ghisays-Ruiz A. and Bonilla L.F New Model or Ellipical Flow Regime in Hydraulically- Fracured Verical Wells in Homogeneo and Naurally- Fracured Sysems. Journal o Engineering and Applied Sciences. ISSN Vol. 9. Nro. 9: Marinez J.A., Eobar F.H. and Bonilla L.F. 01. Reormulaion o he Ellipical Flow Governing Equaion or Horizonal Wells. Journal o Engineering and Applied Sciences. 7(3): Ozkan E. and Raghavan R. 1988, Aug 17. Some New Soluions o Solve Problems in Well Tes Analysis: Par - Compuaional Consideraions and Applicaions. Sociey o Peroleum Engineers. doi:null Tiab D Analysis o Pressure and Pressure Derivaive wihou Type-Curve Maching: 1- Skin and Wellbore Sorage. Journal o Peroleum Science and Engineering. 1: Tiab D Analysis o Pressure and Pressure Derivaive wihou Type Curve Maching: Verically Fracured Wells in Closed Sysems. Journal o Peroleum Science and Engineering. 11: Tiab D Advances in pressure ransien analysis - TDS echnique. Lecure Noes ManualThe Universiy o Oklahoma, Norman, Oklahoma, USA. p Tiab D. and Beam Y Pracical Inerpreaion o Pressure Tess o Hydraulically Fracured Wells in a Naurally Fracured Reservoir. Sociey o Peroleum Engineers. doi:10.118/ ms. January 1. Appendi-A. Esimaing hal-racure lengh by Convenional Analysis Eobar e al. (014) provided he ollowing equaion; qb k P9.486 kh c Or, ell (A.1) P m (A.) Which implies ha a Caresian plo o P vs. (or drawdown) or ΔP vs. [( p+δ) - Δ ] (or buildup) provides a sraigh line which slope, m ell, provides he halracure lengh, qb k khmell c (A.3) During he homogeneo par, Equaion (A.3) becomes: qb k khmell c I he early low regime is linear, hen: qb P h ck P m (A.6) L (A.4) (A.5) So rom he Caresian plo o eiher pressure or pressure drop vers he square roo o ime: qb hml c k (A.7) Appendi-B. Gas low Equaions (1), () and (3) are applied o gas wells i he viosiy and oal sysem compressibiliies are given a iniial condiions, i means (c ) i. However, i hose epressions are epressed ing he pseudoime concep, Agarwal (1979), i yields: 5360

10 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved k ( ) Da a P (B.1) For gas wells, Agarwal (1979) also included he pseudopressure deiniion, mp ( ) D hk mp ( i ) mp ( ) (B.) 14.5qT which dimensionless pseudopressure derivaive is given by: hk * m( P)' D* m( P) D' (B.3) 14.5qT Aer replacing Equaions (B.1), (B.) and (B.3) ino Equaions (3) and (4) leads o obain: a( ) BR BR qt P h( m( P) k (B.4) qt k khmell For gas wells, Equaions (13) and (14) become: qt 1 h( m( P)) k qt 1 h* m( P)' k 0.5 (B.11) (B.1) (B.13) The analogo orm o Equaion (A.7) or gas low is hen: qt 1 hml k The gas version o Equaions (0) o (3) is: 0.5 (B.14) qt ( P) h (* mp ( )' k a BR BR qt h (* mp ( )' BR1 k (B.5) (B.6) qt kh[* m( P ') ma ] (B.15) ka( P) ma (B.16) 5/9 For he homogeneo zone: a( ) BR BR qt P h( m( P) k a( ) BR BR qt P h (* mp ( )' k qt h (* mp ( )' BR1 k (B.7) (B.8) (B.9) In convenional analysis or gas low case, Equaion (3) becomes, (B.17) k ( P) a ma Equaion (4) or gas low is: h [* mp ( ')] qba ( P) 5/16 Epressions (6) and (7) or gas wells: a( P) a( P) a( P) ma a( P) ma (B.18) (B.19) Dividing Equaion (B.6) by Equaion (B.9) and solving or he soraiviy raio: qt k P ( ) a P kh 0 Then, he slope o he Caresian plo will give, (B.10).7778 * m( P)' BR1 * m( P)' BR1 (B.0) 5361

11 VOL. 10, NO. 1, JULY 015 ISSN ARPN Journal o Engineering and Applied Sciences Asian Research Publishing Nework (ARPN). All righs reserved. Dividing Equaion (B.13) by Equaion (B.9) and solving or he soraiviy raio: [* m( P)'] BR1 h k qb [* m( P)'] (B.1) Equaions (9) and (31) rewrien or gas low are: 1 k a ( P) BRi /16 (B.) (B.3) k P 3/5 a( ) Li are: Finally, he gas versions or Equaion (33) and (34) qt khm( P) pss 5/9 (B.4) m( P) high m( P) low 5e (B.5) 536

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions M. Onur 03.04.008 PET467E-Analysis of Well Pressure Tess/008 Spring Semeser/İTÜ Miderm Examinaion (Duraion 3:00 hours) Soluions Name of he Suden: Insrucions: Before saring he exam, wrie your name clearly

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

DETERMINATION OF RESERVOIR DRAINAGE AREA FOR CONSTANT- PRESSURE SYSTEMS USING WELL TEST DATA

DETERMINATION OF RESERVOIR DRAINAGE AREA FOR CONSTANT- PRESSURE SYSTEMS USING WELL TEST DATA TRMINATION OF RSRVOIR RAINAG ARA FOR CONSTANT- PRSSUR SSTMS USING WLL TST ATA Freddy-Humbero scobar *, uly-andrea Hernández and jebbar Tiab Universidad Surcolombiana, Av. Pasrana Cra., Neiva, Huila, Colombia

More information

SPE ω =... (5)

SPE ω =... (5) SPE 9003 Use of Producion Daa Inversion o Evaluae Performance of Naurally Fracured Reservoirs T. Marhaendrajana, Insiu Teknologi Bandung; T.A. Blasingame, Texas A&M Universiy; and J.A. Rushing, Anadarko

More information

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri.

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri. UNIVERSIY OF RENO DEPARMEN OF INFORMAION AND COMMUNICAION ECHNOLOGY 385 Povo reno Ialy Via Sommarive 4 hp://www.di.unin.i MEASUREMENS OF RANSIEN PHENOMENA WIH DIGIAL OSCILLOSCOPES Anonio Moschia Fabrizio

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

A New Correlating Parameter to Quantify Productivity of Extended Hydraulic Fractured Wells

A New Correlating Parameter to Quantify Productivity of Extended Hydraulic Fractured Wells Inernaional Journal o Applied Engineering Research ISSN 97-56 Volume, Number 5 (6) pp 979-985 A New Correlaing Parameer o Quani Producivi o Eended Hdraulic Fracured Wells T. Marhaendrajana* Associae Proessor,

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Effect of Heat Loss on Pressure Falloff Tests

Effect of Heat Loss on Pressure Falloff Tests Inernaional Journal of Engineering Trends and Technology (IJETT) Volume 39 Number 1- Sepember 016 Effec of Hea Loss on Pressure Falloff Tess Ashkan Jahanbani G. #1, Tom A. Jelmer # #1, eparmen of Peroleum

More information

Lecture #6: Continuous-Time Signals

Lecture #6: Continuous-Time Signals EEL5: Discree-Time Signals and Sysems Lecure #6: Coninuous-Time Signals Lecure #6: Coninuous-Time Signals. Inroducion In his lecure, we discussed he ollowing opics:. Mahemaical represenaion and ransormaions

More information

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS Scieniic Research o he Insiue o Mahemaics and Compuer Science DIERENT VARIANTS O THE BOUNDARY ELEMENT METHOD OR PARABOLIC EQUATIONS Ewa Majchrzak,, Ewa Ładyga Jerzy Mendakiewicz, Alicja Piasecka Belkhaya

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS

A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS THERMAL SCIENCE: Year 7, Vol., No. A, pp. 33-4 33 A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS by Xiao-Jun YANG a and Feng GAO a,b * a School of Mechanics and Civil Engineering, China Universiy

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

SPE Copyright 2001, Society of Petroleum Engineers Inc.

SPE Copyright 2001, Society of Petroleum Engineers Inc. SPE 757 Decline Curve Analysis Using Type Curves Evaluaion of Well Performance Behavior in a Muliwell Reservoir Sysem T. Marhaendrajana, Schlumberger, and T. A. Blasingame, Texas A&M Universiy Copyrigh,

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

MINI FRAC & AFTER CLOSURE ANALYSIS REPORT

MINI FRAC & AFTER CLOSURE ANALYSIS REPORT MINI FRAC & AFTER CLOSURE ANALYSIS REPORT Company Ld. Name Field Name Formaion February 6 18, 013 DISTRIBUTION: Clien Name/Company Ld. DEFINITIVE PROJECT COORDINATOR: PREPARED BY: Definiive Opimizaion

More information

Hybrid Forecasting Methods for Multi-Fractured Horizontal Wells: EUR Sensitivities

Hybrid Forecasting Methods for Multi-Fractured Horizontal Wells: EUR Sensitivities Advances in Peroleum Exploraion and Developmen Vol. 3, No., 0, pp. -0 DOI:0.3968/j.aped.9554380030.5 ISSN 95-54X [Prin] ISSN 95-5438 [Online] www.cscanada.ne www.cscanada.org Hybrid Forecasing Mehods for

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams Combined Bending wih Induced or Applied Torsion of FRP I-Secion Beams MOJTABA B. SIRJANI School of Science and Technology Norfolk Sae Universiy Norfolk, Virginia 34504 USA sirjani@nsu.edu STEA B. BONDI

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves Rapid Terminaion Evaluaion for Recursive Subdivision of Bezier Curves Thomas F. Hain School of Compuer and Informaion Sciences, Universiy of Souh Alabama, Mobile, AL, U.S.A. Absrac Bézier curve flaening

More information

Math 221: Mathematical Notation

Math 221: Mathematical Notation Mah 221: Mahemaical Noaion Purpose: One goal in any course is o properly use he language o ha subjec. These noaions summarize some o he major conceps and more diicul opics o he uni. Typing hem helps you

More information

Economics 8105 Macroeconomic Theory Recitation 6

Economics 8105 Macroeconomic Theory Recitation 6 Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

A Limit Symmetry of Modified KdV Equation and Its Applications

A Limit Symmetry of Modified KdV Equation and Its Applications Commun. Theor. Phys. 55 011 960 964 Vol. 55 No. 6 June 15 011 A Limi Symmery o Modiied KdV Equaion and Is Applicaions ZHANG Jian-Bing Ï 1 JI Jie SHEN Qing ã 3 and ZHANG Da-Jun 3 1 School o Mahemaical Sciences

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

INDEX. Transient analysis 1 Initial Conditions 1

INDEX. Transient analysis 1 Initial Conditions 1 INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera

More information

I. OBJECTIVE OF THE EXPERIMENT.

I. OBJECTIVE OF THE EXPERIMENT. I. OBJECTIVE OF THE EXPERIMENT. Swissmero raels a high speeds hrough a unnel a low pressure. I will hereore undergo ricion ha can be due o: ) Viscosiy o gas (c. "Viscosiy o gas" eperimen) ) The air in

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18 A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly

More information

Orthogonal optimization of horizontal well fracturing method

Orthogonal optimization of horizontal well fracturing method Vol.2, No.9, 1030-1034 (2010) hp://dx.doi.org/10.4236/ns.2010.29126 Naural Science Orhogonal opimizaion of horizonal ell fracuring mehod Chi Dong*, Kao-Ping Song, Xiao-Qi Chen Key Laboraory of Minisry

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Advanced Organic Chemistry

Advanced Organic Chemistry Lalic, G. Chem 53A Chemisry 53A Advanced Organic Chemisry Lecure noes 1 Kineics: A racical Approach Simple Kineics Scenarios Fiing Experimenal Daa Using Kineics o Deermine he Mechanism Doughery, D. A.,

More information

10.6 Parametric Equations

10.6 Parametric Equations 0_006.qd /8/05 9:05 AM Page 77 Secion 0.6 77 Parameric Equaions 0.6 Parameric Equaions Wha ou should learn Evaluae ses of parameric equaions for given values of he parameer. Skech curves ha are represened

More information

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations Naural Gas Engineering Producion Analysis Reserves Esimaion in Unconvenional Reservoirs New Rae-Time Relaions T.A. Blasingame, Texas A&M U. Deparmen of Peroleum Engineering Texas A&M Universiy College

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder# .#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,

More information

Estimation of complex motions in video signals. Erhardt Barth. Institute for Neuro- and Bioinformatics University of Lübeck, Germany

Estimation of complex motions in video signals. Erhardt Barth. Institute for Neuro- and Bioinformatics University of Lübeck, Germany Esimaion o complex moions in video signals Erhard Barh Insiue or Neuro- and Bioinormaics Universiy o Lübec, Germany Nonlinear analysis o mulidimensional signals: LOCOMOTOR LOcal adapive esimaion o COmplex

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Fitting Parameters in a Differential Equation With a Non-Analytical Solution

Fitting Parameters in a Differential Equation With a Non-Analytical Solution NDSolve u.nb 1 Fiing Parameers in a Dierenial Equaion Wih a Non-Analyical Soluion Mahew Smih 9-Jan-005 In[]:= Remove "Global` " O General::"spell" O General::"spell1" Graphics` Oen dierenial equaions do

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

Solutionbank Edexcel AS and A Level Modular Mathematics

Solutionbank Edexcel AS and A Level Modular Mathematics Page of 4 Soluionbank Edexcel AS and A Level Modular Mahemaics Exercise A, Quesion Quesion: Skech he graphs of (a) y = e x + (b) y = 4e x (c) y = e x 3 (d) y = 4 e x (e) y = 6 + 0e x (f) y = 00e x + 0

More information

CHAPTER 2: Mathematics for Microeconomics

CHAPTER 2: Mathematics for Microeconomics CHAPTER : Mahemaics for Microeconomics The problems in his chaper are primarily mahemaical. They are inended o give sudens some pracice wih he conceps inroduced in Chaper, bu he problems in hemselves offer

More information

y h h y

y h h y Porland Communiy College MTH 51 Lab Manual Limis and Coninuiy Aciviy 4 While working problem 3.6 you compleed Table 4.1 (ormerly Table 3.1). In he cone o ha problem he dierence quoien being evaluaed reurned

More information

LAB # 2 - Equilibrium (static)

LAB # 2 - Equilibrium (static) AB # - Equilibrium (saic) Inroducion Isaac Newon's conribuion o physics was o recognize ha despie he seeming compleiy of he Unierse, he moion of is pars is guided by surprisingly simple aws. Newon's inspiraion

More information

Lecture 10: Wave equation, solution by spherical means

Lecture 10: Wave equation, solution by spherical means Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in

More information

a. Show that these lines intersect by finding the point of intersection. b. Find an equation for the plane containing these lines.

a. Show that these lines intersect by finding the point of intersection. b. Find an equation for the plane containing these lines. Mah A Final Eam Problems for onsideraion. Show all work for credi. Be sure o show wha you know. Given poins A(,,, B(,,, (,, 4 and (,,, find he volume of he parallelepiped wih adjacen edges AB, A, and A.

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

Single and Double Pendulum Models

Single and Double Pendulum Models Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double

More information

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN

More information

REFORMULATION OF THE ELLIPTICAL FLOW GOVERNING EQUATION FOR A MORE COMPLETE WELL TEST DATA INTERPRETATION IN HORIZONTAL WELLS

REFORMULATION OF THE ELLIPTICAL FLOW GOVERNING EQUATION FOR A MORE COMPLETE WELL TEST DATA INTERPRETATION IN HORIZONTAL WELLS VO. 7, NO. 3, MARCH 0 ISSN 89-6608 006-0 Aian Reearch Publihing Neor (ARPN. All righ reerved..arpnjournal.com REFORMUATION OF THE EIPTICA FOW GOVERNING EQUATION FOR A MORE COMPETE WE TEST DATA INTERPRETATION

More information

04. Kinetics of a second order reaction

04. Kinetics of a second order reaction 4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC EE 435 Lecure 3 Absolue and Relaive Accuracy DAC Design The Sring DAC . Review from las lecure. DFT Simulaion from Malab Quanizaion Noise DACs and ADCs generally quanize boh ampliude and ime If convering

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in Circui Variables 1 Assessmen Problems AP 1.1 Use a produc of raios o conver wo-hirds he speed of ligh from meers per second o miles per second: ( ) 2 3 1 8 m 3 1 s 1 cm 1 m 1 in 2.54 cm 1 f 12 in 1 mile

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

4.5 Constant Acceleration

4.5 Constant Acceleration 4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic Model for Cancer Cell Growth through Single Forward Mutation Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

10.1 EXERCISES. y 2 t 2. y 1 t y t 3. y e

10.1 EXERCISES. y 2 t 2. y 1 t y t 3. y e 66 CHAPTER PARAMETRIC EQUATINS AND PLAR CRDINATES SLUTIN We use a graphing device o produce he graphs for he cases a,,.5,.,,.5,, and shown in Figure 7. Noice ha all of hese curves (ecep he case a ) have

More information

RC, RL and RLC circuits

RC, RL and RLC circuits Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.

More information

The field of mathematics has made tremendous impact on the study of

The field of mathematics has made tremendous impact on the study of A Populaion Firing Rae Model of Reverberaory Aciviy in Neuronal Neworks Zofia Koscielniak Carnegie Mellon Universiy Menor: Dr. G. Bard Ermenrou Universiy of Pisburgh Inroducion: The field of mahemaics

More information

2.9 Modeling: Electric Circuits

2.9 Modeling: Electric Circuits SE. 2.9 Modeling: Elecric ircuis 93 2.9 Modeling: Elecric ircuis Designing good models is a ask he compuer canno do. Hence seing up models has become an imporan ask in modern applied mahemaics. The bes

More information

Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture

Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture Scienific Herald of he Voronezh Sae Universiy of Archiecure and Civil Engineering. Consrucion and Archiecure UDC 625.863.6:551.328 Voronezh Sae Universiy of Archiecure and Civil Engineering Ph. D. applican

More information

Age (x) nx lx. Age (x) nx lx dx qx

Age (x) nx lx. Age (x) nx lx dx qx Life Tables Dynamic (horizonal) cohor= cohor followed hrough ime unil all members have died Saic (verical or curren) = one census period (day, season, ec.); only equivalen o dynamic if populaion does no

More information

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling 2.39 Tuorial Shee #2 discree vs. coninuous uncions, periodiciy, sampling We will encouner wo classes o signals in his class, coninuous-signals and discree-signals. The disinc mahemaical properies o each,

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

PRESSURE AND PRESSURE DERIVATIVE ANALYSIS FOR ASYMMETRY FINITE-CONDUCTIVITY FRACTURED VERTICAL WELLS

PRESSURE AND PRESSURE DERIVATIVE ANALYSIS FOR ASYMMETRY FINITE-CONDUCTIVITY FRACTURED VERTICAL WELLS Revista Fuentes: El Reventón Energético Vol. 15 Nº 2 de 2017 - Jul/Dic - pp 71/78 PRESSURE AND PRESSURE DERIVATIVE ANALYSIS FOR ASYMMETRY FINITE-CONDUCTIVITY FRACTURED VERTICAL WELLS Freddy Humberto Escobar

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

Calculus Tricks #1. So if you understand derivatives, you ll understand the course material much better. a few preliminaries exponents

Calculus Tricks #1. So if you understand derivatives, you ll understand the course material much better. a few preliminaries exponents Calculus Tricks # Eric Doviak Calculus is no a pre-requisie or his course. However, he oundaions o economics are based on calculus, so wha we ll be discussing over he course o he semeser is he inuiion

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information