Estimation of parameters for simulation of steady state foam flow in porous media

Transcription

1 Estimation o parameters or simulation o steady state oam lo in porous media Kun Ma, Jose L. Lopez-alinas, ibani Lisa Bisal and Geore J. Hirasai Department o Chemical & Biomolecular Enineerin Rice University, Houston, TX 07/02/202

2 Outline. Foam simulators have many parameters. Ho do e determine them? 2. Compare the experimental results ith the oam models in a commercially available reservoir simulator. 3. Develop methodoloy to describe oam mobility rom common oam experiments.

3 ystem A ystem B

4 -D oam experiments ystem A andpac: silica sand 20/40 Lenth: 27.5 cm Inner diameter: 2.58 cm uractant: IO 58 ith.0% t NaCl R-CHOH-CH 2 -CHO 3 --R ~75% R-CH=CH-CHO 3 --R ~25%, here R+R = C2-5 Permeability: 58.0 darcy Porosity: 36.0%

5 -D oam experiments Total supericial velocity: 20 t/day oam, app p u u

6 -D oam experiments Total supericial velocity: 20 t/day

7 as mobility reduction /FM. Ashoori E, Heijden TLM, Rossen WR 200 Fractional-Flo Theory o Foam Displacements With Oil. PE Journal 5:pp Computer Modelin Group 2007 TAR TM User's Guide. Calary, Alberta, Canada Foam model u r p D r n r FM FM mmob F ater F sur Gas mobility is a unction o both ater saturation and suractant concentration.

8 FM TAR Foam model old mmob F ater F sur mmob: the reerence oam mobility reduction actor; F ater 0.5 arctan[ epdry mdry] mdry: the critical ater saturation volume raction above hich the maximum oam strenth is reached; Cs epsur Fsur or Cs msur msur or C s msur msur: the critical suractant concentration above hich as mobility is independent o suractant concentration.. Rossen, W. R. and Renema, W. J uccess o Foam AG Processes in Heteroeneous Reservoirs. PE Annual Technical Conerence and Exhibition. Anaheim, Caliornia, U..A., ociety o Petroleum Enineers.

9 Hih and lo quality reime r / / / r r r / / r mdry? n r FM r. Chen, L., Reme, A. B., et al imulatin Foam Processes at Hih and Lo Foam Qualities. PE/DOE Improved Oil Recovery ymposium. Tulsa, Olahoma. 2. Alvarez, J. M., Rivas, H. J., et al Uniied Model or teady-tate Foam Behavior at Hih and Lo Foam Qualities. PE Journal 63.

10 and mdry An example usin mmob = 0000, epdry=000 and mdry = 0.:. is close but not equal to mdry; 2. can be calculated throuh max oam, app oam, app mdry=0.000 =0.038

11 and mdry An example usin mmob = 0000, epdry=000 and mdry = 0.: - curve is very steep near and precise calculation o is needed. transition oam quality mdry=0.000 =0.038

12 The problem to solve r r oam app mdry mmob measured,,,,, mdry mmob measured r r max,, oam app app oam r r oam app, olve mmob, mdry and throuh the olloin equations:

13 Usin Equations a and b to determine a contour plot o as a unction o mmob and mdry Usin Equations c and d to determine a contour plot 2 o μ oam,app as a unction o mmob and mdry max oam, app oam, app Eqn a, measured r r Eqn c r r Eqn b oam, app r r Eqn d Perorm superposition o contour plots and 2 and indentiy the point mmob, mdry here =,measured in contour plot and μ oam,app = μ oam, measured in contour plot 2 cross over

14 Match experimental data =0.5 Computed rom: r r Computed rom: r r oam app,

15 Match experimental data mmob = mdry = 0.072

16 Match experimental data Total supericial velocity: 20 t/day mmob=44200 mdry=0.072 epdry=000

17 Eect o epdry mmob = 44200, and mdry = 0.072

18 Dependence on suractant concentration

19 Revised Foam model ne FM mmob F ater F sur F ater 0.5 arctan[ epdry mdry C s msur epmdry ] instead o mdry in the old model F sur C s msur epsur or C or C s s msur msur

20 urace tension msur hypothesized

21 Match experimental data

22 Comparison o

23 ystem A ystem B

24 Foam apparatus or hih temperature ystem B

25 Experimental data uractants: ZA blend Temperature: 94

26 Match experimental data mmob mdry 0.044

27 Match experimental data uractants: ZA blend Temperature: 94 Total supericial velocity: 3.0 t/day

28 Experimental data shear thinnin

29 hear thinnin eect in oam model r n r FM FM mmob F ater F shear F sur F ater 0.5 arctan[ epdry mdry] u F shear u, re epv F sur deault value

30 hear thinnin eect r r oam app u,,, r r u n c r c r r 0 FM n c r c r r 0 next slide

31 hear thinnin eect ] arctan[ 0.5, 0, n c r c oam app r epv re u u mdry epdry mmob n oam app r c r c, 0 next slide

32 Lo-lo it to experimental data lo 0 u x } { arctan 0.5 lo, 0, 0, 0 0 mdry epdry y n n oam app r c r c n oam app r oam app r lo, 0 epv re u mmob x epv y

33 Lo-lo it to experimental data =0.7

34 Match experimental data y y epv x lo 0 mmob epv u, re x epv mmob lo0 u epv i u, re 2.7 t / day mmob 3364, re 0.7 or shear-thinnin modelin

35 ,, Parameters obtained Dry-out modelin u total 3.0 t / day mmob mdry epv 0 u, re arbitrary lide 25 hear-thinnin modelin mmob 3364 mdry epv u re, t / assumed 0.7 assumed day Used in next slide

36 Match experimental data

37 Conclusions. A ne method o ittin the parameters in the TAR oam model is presented and a unique roup o parameters is ound or modelin the oam property in silica sandpac ith the suractant 0.02%-0.2% IO 58 in.0% NaCl solution. 2. A revised dry-out model or eect o suractant concentration is proposed. 3.The critical suractant concentration msur in the oam model is at least one order o manitude above the CMC. 4. The shear thinnin eect is modeled usin the TAR oam model.

38 Than you!

39 Parameters or oam simulation

40 Match experimental data ystem B mmob mdry 0.044

41 Eect o epdry An example usin mmob = 0000, and mdry = 0.: