Naturally Fractured REservoirs

[spydaman]

Im trying to simulate my first model ever wid Eclipse. This is what ive tried doing, but its not working right. Error says "There are no Active Cells in the REservoir " and " Incomplete set of Corner points Specified". PLease Help !!


RUNSPEC

TITLE
Fracture Model Design

START
21 'MAR' 2000 /

METRIC

OIL

WATER

VISCD

DUALPERM

DUALPORO

MONITOR

RSSPEC

NOINSPEC

DIMENS
3 3 4 /

EQLDIMS
1 100 50 1 20 /

REGDIMS
1 1 0 0 /

TABDIMS
1 1 20 20 1 20 20 1 /

WELLDIMS
3 3 3 3 /


GRID

GRIDFILE
2 /

INIT

ZCORN
2500 2500 3000 3000 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 2800 2800
2500 2500 2500 2500 2800 2800 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500
3000 3000 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 3000 3000 2500 2500
2500 2500 2800 2800 2500 2500 2500 2500 2800 2800 2500 2500 2500 2500 3000 3000
2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500
3000 3000 2500 2500 2500 2500 2800 2800 2500 2500 2500 2500 2800 2800 2500 2500
2500 2500 3000 3000 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 3000 3000
2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 2800 2800 2500 2500 2500 2500
2800 2800 2500 2500 2500 2500 3000 3000 2500 2500 2500 2500 3000 3000 2500 2500
/
-- Grid Block Coordinate Lines
--
-- Grid Block Coordinate Lines
--
--
-- Grid Block Coordinate Lines
--
--
-- Grid Block Coordinate Lines
--
COORD
0 0 2500 0 0 2500
0.001 0 2500 0.001 0 2500
0.0011 0 2500 0.0011 0 2500
0.0021 0 2500 0.0021 0 2500
0 0.001 2500 0 0.001 2500
0.001 0.001 2500 0.001 0.001 2500
0.0011 0.001 2500 0.0011 0.001 2500
0.0021 0.001 2500 0.0021 0.001 2500
0 0.002 2500 0 0.002 2500
0.001 0.002 2500 0.001 0.002 2500
0.0011 0.002 2500 0.0011 0.002 2500
0.0021 0.002 2500 0.0021 0.002 2500
0 0.003 2500 0 0.003 2500
0.001 0.003 2500 0.001 0.003 2500
0.0011 0.003 2500 0.0011 0.003 2500
0.0021 0.003 2500 0.0021 0.003 2500 /

[DAH7542]

spydaman,

Upload the whole *.data file (and include files if used) to take a look of it.

[spydaman]

Thank you for replying. Ive built my model not using the coding, but used the User-Interphase of Eclipse Software. SO ive basically copy pasted the data files from notepad.REgret for any errors caused.



RUNSPEC

TITLE
Fracture Model Design

START
21 'MAR' 2000 /

METRIC

OIL

WATER

VISCD

DUALPERM

DUALPORO

MONITOR

RSSPEC

NOINSPEC

DIMENS
3 3 4 /

EQLDIMS
1 100 50 1 20 /

REGDIMS
1 1 0 0 /

TABDIMS
1 1 20 20 1 20 20 1 /

WELLDIMS
3 3 3 3 /


GRID
--
--Force NOECHO for GRID and EDIT Sections
NOECHO
--

GRIDFILE
2 /

INIT



INCLUDE
GRIDUNIT
--
-- Grid data units
--
'METRES' /
MAPAXES
--
-- Grid Axes wrt Map Coordinates
--
0 0 0 0 0 0 /
ZCORN
2500 2510 2510 2520 2520 2530 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500
2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500 2500
2500 2500 2500 2500 2510 2510 2510 2520 2520 2530 2510 2500 2500 2500 2500 2500
2510 2500 2500 2500 2500 2500 2510 2500 2500 2500 2500 2500 2510 2500 2500 2500
2500 2500 2510 2500 2500 2500 2500 2500 2510 2520 2520 2530 2530 2540 2510 2500
2500 2500 2500 2500 2510 2500 2500 2500 2500 2500 2510 2500 2500 2500 2500 2500
2510 2500 2500 2500 2500 2500 2510 2500 2500 2500 2500 2500 2520 2520 2520 2530
2530 2540 2520 2500 2500 2500 2500 2500 2520 2500 2500 2500 2500 2500 2520 2500
2500 2500 2500 2500 2520 2500 2500 2500 2500 2500 2520 2500 2500 2500 2500 2500
/
-- Grid Block Coordinate Lines
-- Grid Block Coordinate Lines
--
-- Grid Block Coordinate Lines
--
--
-- Grid Block Coordinate Lines
--
--
-- Grid Block Coordinate Lines
--
COORD
0 0 0 0 0 0
1000 0 0 1000 0 0
2000 0 0 2000 0 0
3000 0 0 3000 0 0
0 1000 0 0 1000 0
1000 1000 0 1000 1000 0
2000 1000 0 2000 1000 0
3000 1000 0 3000 1000 0
0 2000 0 0 2000 0
1000 2000 0 1000 2000 0
2000 2000 0 2000 2000 0
3000 2000 0 3000 2000 0
0 3000 0 0 3000 0
1000 3000 0 1000 3000 0
2000 3000 0 2000 3000 0
3000 3000 0 3000 3000 0 /

MAPUNITS
--
-- Units for MAPAXES Data
--
'm' /

INCLUDE
PORO
0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1
0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8 0.8 0.1 0.8
/
BOX
1 3 1 3 1 2 /
PERMX
50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 1 3 4 /
PERMX
50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 3 1 2 /
PERMY
50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 1 3 4 /
PERMY
50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 3 1 2 /
PERMZ
50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 1 3 4 /
PERMZ
50 1000 2*50 1000 50
/
ENDBOX
BOX
1 3 1 3 1 2 /
ACTNUM
--
-- Active Grid Blocks
--
18*1
/
ENDBOX

--
MULTY
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
/
MULTZ
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
/
MULTX
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
/
OLDTRAN
--
-- Block Centre Transmissibilities
--
SIGMAV
0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3
/


--
--Force ECHO back ON after GRID and EDIT Sections
ECHO
--

PROPS

END


EQUIL
--
-- Equilibration Data Specification
--
2500 172.368925 2500.9 1* 2500.2 1* 1* 1* 1*
/

RTEMP
--
-- Initial Reservoir Temperature
--
85 /

RPTSOL
--
-- Initialisation Print Output
--
'PRES' 'SOIL' 'SWAT' 'RESTART=2' 'FIP=2' /

RPTRST
--
-- Restart File Output Control
--
'BASIC=2' /


ROCK
--
-- Rock Properties
--
137.89514 5e-005
/

PVTW
--
-- Water PVT Properties
--
137.89514 1 4.35113231691849e-005 0.6 1*
/

PVCDO
--
-- Dead Oil PVT Properties
--
137.89514 1.1 0.000145037743897283 3.2 1*
/

RSCONSTT
--
-- Constant Rs per PVT Region
--
10 1
/

DENSITY
--
-- Fluid Densities at Surface Conditions
--
932.274372 1008.6824262 1
/



-- OFFICE-SCAL-HEADER-DATA
--
-- Water/Oil Saturation Functions
--
SWOF
--
-- Water/Oil Saturation Functions
--
0.19 0 0.8 0
0.5 0.3 0.4 0
0.8 0.6 0 0
1 0.8 0 0
/

END

[spydaman]

Also, is there anyway we can activate Natural Fractures in Eclipse other than through Dual Porosity Dual Permeability [DPDP]?? According to their manual, on activating DPDP model, each grid block is divided into 2 simulation cells, representing matrix and fracture volumes of the cell.

This is a line from Eclipse manual

"To model such systems[i.e DPDP], two simulation cells are associated with each block in the geometric grid, representing the matrix and fracture volumes of the cell."

I want to individually add properties to both matrix and fractures, so that i could vary my fracture properties like width, orientation and see the difference in performance when it is flooded wid Water, Polymer and ASP. Hope i got my message through.
Thank you for ur time !!

[DAH7542]

spydaman,

Error says "There are no Active Cells in the REservoir " and " Incomplete set of Corner points Specified":

Reason: your grid is 3x3x4 (36 cells total). You have defined only 144 values in ZCORN, but you must have 36*8 values in ZCORN since each block has 8 corners. That's the problem, ECLIPSE has insufficient data to calculate pore volumes and they are set to zero (as warned in PRT file)

[spydaman]

Dear DAH7542
Thank your for ur replies

i dont actually know the issue, because i input all the values and when i run it, the values disappear. Anyways ive sorted it out by running the model in using Block Center rather than Corner Point

But i have this issue simulating the fracture.According to the Eclipse Manual

"In a dual porosity run of ECLIPSE the number of layers in the Z-direction should be doubled.
ECLIPSE associates the first half of the grid (the first NDIVIZ/2 layers) with the matrix blocks,
and the second half with the fractures. In such runs NDIVIZ must therefore be even; ECLIPSE
checks that this is the case."

I have input the matrix properties for the first 4 Zlayers of my 7x7x8 model and fracture properties for the remaining 4 Zlayer. But running the simualtion in the Grid Section, the Permeabilities and Porosities entered appear the same for all the region.. Moreover, no matter how many layers i add in the Zlayer, it is considered as a single layer when simulated.
What do you think d possible reason is ?? When i remove the Fractured REservoir model from the Case Definition, it works perfect, but d moment i click on the Fractured REservoir, the model completely changes.
Could u please help me how to tackle modelling a Fractured REservoir ??

Thank you

[spydaman]

....

[DAH7542]

I suggest you do this:

1. Check the *.data file directly (i.e. read it with any text editor) before running the simulation to ensure that all the keywords and the information you defined are there, you're probably constructing your model by using Elipse-Office (Am I right?); from my experience, Eclipse-Office sometimes has problems, and I wonder if you are lacking some information just like it happened before with the ZCORN values.

2. Take a look at the example data files provided with Eclipse within the installation folder, there are many example runs including DualPoro-SinglePerm and DualPoro-DualPerm simulation.

[paolomaldini]

Thanks

[spydaman]

Thank you DAH7542

I checked out the manuals n found a Dual Poro example, all thanks to you..
I tried to model a polymer flood with it, but some how gave me the following error
" NON-LINEAR EQUATION CONVERGENCE FAILURE
ITERATION LIMIT REACHED - TIME STEP CHOPPED FROM 0.53"

Everythin else seems to be ok.. Ill post my data file here..

Hope you could help me out wid this one too
Thank you again


--
-- DUAL POROSITY TEST RUN NUMBER 1
--
-- THIS IS AN 8 BY 8 QUARTER FIVE SPOT, DIMENSIONED 8 BY 8 BY 2
-- TO ALLOW FOR THE MATRIX AND FRACTURE BLOCKS
--
--
--

RUNSPEC ================================================== ===========
POLMYER FLOODING IN FRACTURED RESERVOIR TEST1
= NDIVIX NDIVIY NDIVIZ QRDIAL NUMRES QNNCON MXNAQN MXNAQC QDPORO QDPERM
8 8 2 F 1 T 0 0 T F /
= OIL WATER GAS DISGAS VAPOIL
T T F F F /
= UNIT CONVENTION
'FIELD ' /
= NRPVT NPPVT NTPVT
5 13 1 /
= NSSFUN NTSFUN QDIRKR QREVKR
13 2 F T /
= NDRXVD NTEQUL NDPRVD
2 1 100 /
= NTFIP
1 /
= NWMAXZ NCWMAX NGMAXZ NWGMAX
2 13 1 2 /
= QIMCOL NWFRIC NUPCOL
F 0 4 /
= MXMFLO MXMTHP MXMWFR MXMGFR MXMALQ NMMVFT
0 0 0 0 0 0 /
= MXSFLO MXSTHP NMSVFT
0 0 0 /
= NAQFET NCAMAX
0 0 /
= DAY MONTH YEAR
01 'JAN' 2000 /
= QSOLVE NSTACK
T 4 F F F F /

GRID ================================================== ============
DPGRID
EQUALS
'DX' 75 / PROPERTIES COMMON TO MATRIX AND FRACTURES
'DY' 75 /
'DZ' 30 /
'TOPS' 4000 /
'PERMZ' 0 /
'PORO' 0.19 1 8 1 8 1 1 / MATRIX PROPERTIES
'PERMX' 1 /
'PERMY' 1 /
'PORO' 0.01 1 8 1 8 2 2 / FRACTURE PROPERTIES
'PERMX' 50000 /
'PERMY' 50000 /
/

PSEUDOS

RPTGRID FIELD 18:07 29 NOV 83
1 1 1 1 1 0 1 1 0 1 0 1 1 1 1 1 0
0 0 0 0 0 0 1 /

--SIGMA FOR 10 BY 10 BY 30 FT BLOCKS
SIGMA
0.08 /

PROPS ================================================== ============

--TO OBTAIN CASE WITH NO IMBIBITION, SET CAPILLARY PRESSURES TO ZERO
SWFN
.2500 .0000 4.0000
.3000 .0200 2.9500
.4000 .0550 1.6500
.5000 .1000 .8500
.6000 .1450 .3000
.7000 .2000 .0000
/
.0000 .0000 4.0000
.1000 .0500 1.8500
.2000 .1100 .9000
.2500 .1450 .7250
.3000 .1800 .5500
.4000 .2600 .4000
.5000 .3550 .2900
.6000 .4750 .2000
.7000 .5850 .1600
.8000 .7150 .1100
.9000 .8500 .0500
1.0000 1.0000 .0000
/

SOF2
.3000 .0000
.4000 .1100
.5000 .2400
.6000 .4200
.7000 .7050
.7500 .9200
/
.0000 .0000
.1000 .0210
.2000 .0570
.3000 .1020
.4000 .1730
.5000 .2400
.6000 .3300
.7000 .4500
.7500 .5190
.8000 .5870
.9000 .7700
1.0000 1.0000
/

PVTW 1 TABLES 13 P NODES 5 R NODES FIELD 18:05 29 NOV 83
.0000000 1.00000 3.03E-06 .50000 0.00E-01 /

PVDO
.00000 1.00000 2.00000
8000.00 .92000 2.00000
/

ROCK
4000.00 .30E-05 /

DENSITY
52.0000 64.0000 .04400 /

RPTPROPS FIELD 15:56 29 NOV 83
1 1 1 1 1 1 1 1 /

REGIONS ================================================== ==========

EQUALS FIELD 18:09 29 NOV 83
'SATNUM ' 1 , 1, 8 , 1, 8 , 1, 1 /
'SATNUM ' 2 , 1, 8 , 1, 8 , 2, 2 /
/

RPTREGS FIELD 18:05 29 NOV 83
0 1 0 0 0 0 0 0 0 0 /

SOLUTION ================================================== ===========

EQUIL

4000 3959 6000 0 0 0 0 0 0 /


RPTSOL
1 1 1 0 0 0 0 2 1 0 0 /

SUMMARY ================================================== ========
FPR
FWCT
FOPR

WWCT
/
SCHEDULE ================================================== =========

RPTSCHED FIELD 10:17 13 JUN 85
1 0 1 0 0 0 2 2 2 0 0 2 0 0 2 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 /

WELSPECS
'I' 'G' 1 1 4000 'WAT' /
'P' 'G' 8 8 4000 'OIL' /
/

COMPDAT FIELD 18:10 29 NOV 83
'I ' 1 1 2 2 'OPEN' 0 .000000 .5000 .00000 .0000 0.000E-01/
'P ' 8 8 2 2 'OPEN' 0 .000000 .5000 .00000 .0000 0.000E-01/
/

WCONINJ FIELD 18:10 29 NOV 83
'I ','WAT','OPEN','RATE' 200.000000 1.0E20 .000 'NONE'
1.0E20 1.0E20 0 .00000000 /
/

WCONPROD FIELD 18:11 29 NOV 83
'P ','OPEN','LRAT' 1.0E20 1.0E20 1.0E20
210.000000 1.0E20 .000000 .000000 0 .00000000/
/



TSTEP
6*200
/


END

[vmauricio1]

Thanks. gracias

[Shakespear]

Thank you DAH7542

I checked out the manuals n found a Dual Poro example, all thanks to you..
I tried to model a polymer flood with it, but some how gave me the following error
" NON-LINEAR EQUATION CONVERGENCE FAILURE
ITERATION LIMIT REACHED - TIME STEP CHOPPED FROM 0.53"

I took your data and read it into OFFICE. Then ran the model. There was no ERROR.

Take a look at this work which I think will explain a lot regarding dual-porosity modeling and may give you some ideas.

Simulating water tracer test in naturally fractured reservoirs using discrete fracture and dual porosity models

[link Point to another website Only the registered members can access]
and this

[link Point to another website Only the registered members can access]