See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/328917231

Reservoir Modeling using CMG-2007 for Single Well of radial Flow

Thesis · April 2012

DOI: 10.13140/RG.2.2.19516.28806

CITATION
READS
1
11,399

4 authors, including:

Ali Yahya Jirjees

Ministry of Oil - North Oil Company
34 PUBLICATIONS 10 CITATIONS

SEE PROFILE

All content following this page was uploaded by Ali Yahya Jirjees on 13 November 2018.

The user has requested enhancement of the downloaded file.

Reservoir Modeling using CMG-2007 for
Single Well of radial Flow
TABLE OF CONTENT
acknowledgement

abstract

table of content

list of illustrations

1 Introduction to the Reservoir Modeling

1.1 Reservoir Simulation Objectives and considerations

1.2 Fundamentals of Reservoir Simulation

1.3 Other engineering approaches

1.4 History of Reservoir Modeling

1.5The Basics of Reservoir Simulation

1.6 Simulation Technology

2. Modeling with CMG 2007 (WinBuilder& IMEX)

2.1 Benefit of CMG

2.2 Applications of CMG’s IMEX

2.3 Selection of model grid

2.4 Specification of reservoir rock and fluid description data

2.5 Validity of simulation results

2.6 Model Assumption

3. The project objectives

4. input data of the project

5. Results of the project

6. Comparing CMG results with the companies results of the project

7. Scenarios

8. Discussion and suggestion

Reference
Nomenclature

appendix.A & appendix.B

List of illustrations:

Fig.1.1 Modeling Complexity Vs. physics and Details

Fig.1.2 Flowchart of the CMG’s how to work

Fig.2.1 CMG‘s launcher veiw screen

Fig 5.1. oil rate sc- daily production with time

Fig 5.2. gas rate sc production with time

Fig 5.3. gas oil ratio sc with time

Fig 5.4. water cut sc with time

Fig 5.5. wellbore bottom hole pressure with time

Fig.6.1 comparing oil rate

Fig.6.2 water cut comparing

Fig.6.3 gas oil ratio comparing

Fig.6.4 well bottom hole pressure comparing

Fig.7.1 perforation change affect on oil rate

Fig.7.2 perforation change affect on gas oil ratio

Fig.7.3 perforation change affect on water cut

Fig.7.4 perforation change affect on oil rate

Fig 7.5 BHP change affect on oil rate

Fig 7.6 BHP change affect on gas oil ratio

Fig 7.7 BHP change affect on reservoir bottom hole flowing pressure with time

Fig 7.8 BHP change affect on water cut

Fig 7.9. 3D view of the reservoir ( production and injection well location)

Fig 7.10 The effect of water injection on the daily oil rate

Fig 7.11 The effect of water injection on the water cut

Fig 7.12 The effect of water injection on the gas oil ratio

Fig 7.13 The effect of water injection on the flowing bottom hole pressure

Fig.8.1 oil rate and gas rate with time

Fig.A.1. gas expansion with Pressure plot

Fig.A.2.solution gas oil ratio with Pressure plot

Fig.A.3. oil and gas viscosity with Pressure plot

Fig.A.4. relative permeability with water saturation plot

Fig.A.5. relative permeability with liquid saturation plot

Fig.A.6. gas oil capillary pressure with liquid saturation plot

Fig.A.7. deviation factor with pressure

Fig.A.8. water oil capillary pressure with water saturation plot

Fig.A.9. gas formation volume factor with pressure plot

Fig.A.10. 3D veiw of reservoir grid top

Fig.A.11. 3D veiw of reservoir grid thickness

Fig.A.12. 3D veiw of reservoir permeability

Fig.A.13. 3D veiw of reservoir permeability

Chapter 2

Modeling with CMG 2007 (WinBuilder& IMEX)

1. Benefit of CMG

2. Applications of CMG’s IMEX

3. Selection of model grid

4. Specification of reservoir rock and fluid description data

5. Validity of simulation results

6. Model Assumption
I MEX is one of the world's fastest conventional reservoir simulator. IMEX
models complex, heterogeneous, faulted oil and gas reservoirs, to achieve
accurate predictions and forecasts quickly. Reservoir engineers use IMEX to move
from history-matched, primary production and water floods to enhanced recovery
processes in GEM and STARS quickly and easily.

IMEX is a full-featured three-phase, four-component, black oil reservoir simulator

for accurately modeling complex heterogeneous faulted structures, primary and
secondary (water & polymer flooding) recovery processes using complex
horizontal and multi-lateral wells. IMEX includes features such as local grid
refinement (LGR), comprehensive well management, pseudo-miscible option,
horizontal wells, dual porosity/permeability, flexible grids, gas adsorption and
many more.

IMEX models multiple PVT and equilibrium regions, as well as, multiple rock
types and has flexible relative permeability choices. Regardless of the size or
complexity of your reservoir problem, IMEX is an effective tool for a broad range
of reservoir management issues.

To optimize your reservoir simulation model, use IMEX after history-matching to

screen multiple recovery techniques to obtain the best net present value (NPV).
The speed of IMEX allows engineers to conduct more runs than any other
simulator.
 1. Benefits of CMG
 Achieve simulation results faster than any other conventional simulator
 Ability to quickly screen a variety of recovery mechanisms before moving forward
to a more complex simulation
 Accurate modeling of the matrix-fracture transfer in fractured reservoirs
 Use the speed of IMEX to simulate shale gas modeling
 Seamless integration with CMOST for rapid, accurate history matches and
uncertainty analysis, while leveraging limited engineering time
 Fast and easy transition to enhanced oil recovery (EOR) process modeling in GEM
and STARS

Fig.2.1 CMG‘s launcher veiw screen

2. Applications of CMG’s IMEX

 Unconventional gas &  secondary recovery

Liquids reservoir  Waterflooding
 Coal Bed Methane (CBM)  Surface facilities
 Shale gas modeling
 Shale liquids  Polymer injection
 Tight gas  Dry gas injection
 primary depletion  Pseudo-miscible solvent
 Infill drilling optimization injection
 Horizontal & multi-lateral well  WAG processes
placement  Gas storage fields
 Coning studies  Cycle optimization
 Under-saturated/saturated oil
reservoirs
 Gas condensate reservoirs
 Gas deliverability forecasting
 Reservoir/surface facility
optimization
 Naturally or hydraulically
fractured reservoirs
 Naturally fractured reservoirs
 Gas-oil gravity drainage
3. Selection of model grid
Selection of the x-y-z grid block network involves many factors, including
available budget and the engineer’s judgment and experience . For any type of
model, the arithmetic or computing expense per time step is at least linearly
proportional to the total number of gridblocks employed. The computing expense
of a single model run is proportional to product of the number of gridblocks and
the number of time steps required by the model to cover the total time period of
interest. In the many cases, the time steps size is controlled by the maximum rate
of change (overall gridblocks) in one or more calculated quantities such as pressure
and saturations. This maximum rate of change generally occurs at or near a well or
in the vicinity of a flood front. A doubling of the number of gridblocks can result
approximately in a doubling of this maximum rate of change since each gridblock
is (on the average) one-half as large. The average time steps size, then might
decrease by a factor of two if the number of blocks were double. The final result is
a computing expense per model run which can approach a proportionality to the
square of the total number of gridblocks. This indicates the importance of selecting
the smallest number of gridblocks consistent with reservoir/well description,
recovery process characteristics, and the questions asked regarding reservoir
performance.

The number of gridblocks and the resultant study computing expense are the
lowest in cases where the engineer can justify use of a representative element of
the total field as the basis for the model study. This may be possible in reservoirs
developed with repeated well patterns, for any recovery process---waterflooding,
CO2 injection, steamflooding , etc. In such cases, the representative element
ideally should be a symmetrical element of the reservoir. In strict terms, this
requires a repeated, regular pattern of identically completed and operated wells, a
horizontal, areally homogeneous reservoir formation of uniform thickness, and
areally uniform initial fluid saturation distributions. If these conditions were met,
then questions regarding total field optimization, forecasting and comparative
evaluation of recovery processes could be addressed inexpensively by simulation
of the single pattern (element).

While actual reservoirs never satisfy these conditions exactly, representative-

element simulation studies are frequently performed for repeated pattern processes.
In some cases, a substantial portion of the reservoir may exhibit only moderate
areal heterogeneity and thickness variation. Resultant variation in performance
from one pattern to another may be sufficiently small for engineering purposes to
justify scale-up of single-pattern results to total field performance.

Representative-element simulation is often performed where the study purpose is

comparative evaluation of alternative recovery processes as opposed forecasting of
total field performance for a specific process and operating scheme. The
justification of single-element simulation implied in such cases is that the resultant
ranking of alternative processes is un affected by the variations in pattern (element)
properties over the field. This justification can be and frequently is checked by
repeating the various process simulations for two or more patterns of different
properties representative of different potions of the reservoir.

Finally, the relatively inexpensive single-element simulation applies to design or

optimization studies of a specific recovery process operated in a repeated pattern
mode. For a repeated-pattern steamflood, single-pattern model runs have been
performed to “ optimize” pattern type (e.g., five-, seven- or nine-spot) and size,
injected steam quality and rate, well completion, etc. occasional publications
describe single-pattern simulation studies using a one-quarter five spot or one-
quarter nine-spot as the symmetrical element of the respective pattern. Actually, a
one-eight five-spot or nine –spot(and 1/12 seven-spot) are the smallest symmetrical
elements and should be used to minimize computing expense.

Currently, a major portion of industry-wide effort and computing expense is

simulation studies is associated with total-field forecasting of black-oil reservoir
performance under a sequence of recovery processes. Typically, the engineer must
select a 3D grid for a large reservoir with significant heterogeneity, large areal
variation in dip and thickness, irregular well location and increasing numbers of
wells with successive development stage. The engineer may face a several to
many-year period of historical performance under natural depletion, frequently
with some natural water encroachment. Study objectives may include history
matching, followed by matching and forecasting for awaterflood period, in turn
followed by forecasting for some tertiary scheme such as CO2 injection.
The number of gridblocks is product of the number of areal blocks, N x N y, and
the number of grid layers, N z. Different consideration enter into selection of these
two numbers of spacing. Factors indicating a need for fine areal grid spacing are
high well density and sharp or rapid change (areally) in permeability, porosity,
thickness, and dip. Since these factors frequently vary over the field, the x-and-y-
direction grid spacing are often non uniform. Grid spacing are generally increase
toward the down dip reservoir boundaries and increase greatly with distance into
the aquifer if the latter is present and included in the grid.

In general, of course, the number of areal gridblocks required increases with size
of the reservoir and the number of wells. However, grid spacing ranging from very
fine to very coarse may be appropriate for different reservoirs of comparable size.
The smallest number of areal blocks (coarsest areal spacings) are associated with
reservoir studies limited to natural depletion and crestal or flank gas/or water
injection. In such a case, a coarse grid may result in a number of areal blocks that
include two or more similar type (e.g., production) wells, with little loss in
engineering significance of the simulator results. Large numbers of areal blocks
may be required in case of pattern waterfloods or enhanced recovery processes. A
rough guide in this case is the need for at least two, preferably three or more,
gridblocks separating each injection-production well pair. However, recent studies
describe estimation of pseudo relative-permeability curves, which allow adjacent-
block placement of an injector/producer well pair.

The major factors affecting the number of grid layers (vertical gridblocks) required
are the formation stratification, vertical communication, and total thickness. Many
reservoirs possess a number of formation layers, which correlate from well to well
over much of or all the field. Variations of layer thickness, permeability, and
porosity may be significant areally even greater from one layer to another. The
vertical communication (vertical permeability) between adjacent layer-pairs may
from zero to very high, both areally and from one layer-pair to another. In the
general, at least one grid layer should be used for each correlatable formation layer.
However, common sense and budget constraints argue against definition of a large
number of very thin grid layers. Three dimensional reservoir studies typically
employ 4 to 12 grid layers, and one or more of these grid layers may be a lumped
representation of several thin formation layers. The need for subdivision of one
formation layer into two or more grid layers depends on the layer thickness and
fluid-segregation characteristics of the recovery process and operating rates. Most
recovery processes result in moderate to severe gravity segregation of oil and
injected fluids; injected water or gas tend to under run or override oil, respectively;
many steamflood projects exhibit sever override of oil by the steam. A formation
layer that has significant thickness and zero to poor vertical communication with
layers above and below may exhibit a pronounced phase segregation and require
two or more grid layers. In the idealized example of a field wide, pronounced
gravity override in a vertically homogeneous reservoir, a variable grid spacing
increasing from top to bottom might be specified. That is, four layers of
thicknesses 5,10,20, and 25-ft might give more accurate results than four layers of
equal 15-ft thickness.

A customary approach to determining N z involves use of the simulation model

itself in 2D cross-sectional (x-z slice) mode. For the particular recovery process of
interest, x-z model runs are performed by using different numbers of grid layers.
Pseudorelative-permeability curves reflection phase segregation are calculated
from model runs performed with fine vertical grid spacing. These pseudocurves are
then used in equivalent x-z model runs using fewer grid layers to obtain coarse
(vertical) definition results similar to fine spacing “correct” results. The fewer grid
layers of the coarse definition are then employed in the 3d reservoir study grid.
This concept of generating pseudocurves for coarse vertical grids that reproduce
vertical fine-grid results (using rock or laboratory relative permeabilities) has been
extended to the areal spacing problem, as mentioned earlier.

Obviously, a minimum computing expense follows from use of a single grid layer
representing the entire formation thickness. This results in a 2D x-y areal grid as
opposed to a 3D grid and occasionally is justified in the two extremes of a very
high vertical permeability and a layer formation with zero vertical permeability.
Pseudorealtive permeability and capillary pressure curves are discussed for the
former case in papers describing the vertical equilibrium (VA) concept and for the
latter case by Hearn.
4. Specification of reservoir rock and fluid description data
Geological and petro-physical work based on logs and core analyses yields maps
of structure, net Øh, and kh products for each of the several reservoir layers. The
kh and Øh data are often augment or modified by results of drill-stem, pressure
buildup, and pulse tests. For each layer, the engineer can overlay his areal x-y grid
spacing network on these maps and read off the values of subsea depth, Øh, and kh
at the center of each gridblock. These values along with gross thickness of each
block are then transposed to a data file in a format compatible with that required by
the simulation model. Current research effort is directed toward developing
computer programs that accept digitized core analysis, log and geological data the
selected grid network, and through mapping and interpolation techniques,
automatically prepare the simulation input data file. Laboratory core analysis work
includes measurement of relative-permeability, k r, and capillary-pressure, P c,
curves for a number of field cores. Variations in rock lithology may result in
different sets of k r, an P c, curves for different layers and/or different areal
portions of the reservoir. Most simulation models allow multiple sets of such data
in tubular from with assignment of each set to a user-specified layer/portion of the
reservoir. If the rock water/oil (gas/oil) capillary pressure values are small, the
water/oil (gas/oil) transition zone in the reservoir may be a very small fraction of
total formation thickness. In such cases pseudocapillary-pressure curve(s) should
be used.For black-oil studies, laboratory tests are performed to determine gas
compressibility factor and saturated oil and gas viscosities vs. pressure.
Differential and/or constant-composition expansion tests on oil samples yield the
saturated oil pressure-dependent formation volume factor, Bo (RB/STB), and
solution gas, Rs (scf/STB). The resulting oil and associated gas properties vs.
pressure are entered in the data file in tabular from compatible with simulator input
requirements. For gas condensate depletion studies, constant-volume and constant-
composition expansion tests yield the required pressure-dependent liquid content,
CL (STB/scf), and condensate density values. A wide variety of laboratory tests
are performed for compositional model studies that involve injection of a
nonequilibrium fluid (dry or enriched gas, CO2, N2, etc.). swelling tests yield
relative volumes, saturation pressure, and equilibrium phase composition for each
of a sequence of mixture, say. 1 mole of original reservoir oil and injected fluid.
5. Validity of simulation results
Uncertainties or errors in simulation model results may arise from questionable
assumption or mechanisms not represented in the differential form of the model,
spatial and time truncation error introduced by replacement of the model
differential equations by finite difference approximations, and inadequately known
reservoir rock and/or fluid description data. In addition, the exact solution of the
difference equation is not attained because of round-off error introduced by the
finite word length of the computer. Round-off error is generally negligible
compared with errors from the other three sources. With some exceptions, the
above sources of error are listed in order of increasing importance. However,
successful history matching can reserve the importance of the second and third
sources.

Comparisons of model and laboratory experiment results can indicate model

validity in the absence of the uncertainty 3 above. Several such comparisons show
good model-experiment for gas/oil systems, water/oil coning and fractured-matrix
imbibition.

6. Model Assumption
An assumption common to many black-oil models is complete re-solution of free
gas in accordance with the saturated Rs (p) curve during re-pressurization. This
may be a poor assumption in a case where gridblock thickness is large and gas/oil
gravity (vertical) segregation is pronounced. Prior to re-pressurization in a given
block, the free gas may exist as a high gas saturation in only the upper portion of
the block. This contradicts its representation in the model as a lower saturation
distributed throughout the entire block volume. In the segregated state , the gas will
re-dissolve only in the lower or residual oil saturation in the upper, gas-occupied
portion of the block volume. However, the model will allow re-solution in the
entire rocks oil volume. Pressure hysteresis in the Rs (p) curve has been used to
cope with this problem; an alternative remedy where the computing budget
permits is the use of more grid layers.
An assumption common in early black-oil models was that the reservoir oil obeyed
a single pair of Bo (p) and Rs (p) curves. Some black-oil reservoirs exhibit a
significant variation of oil API gravity and PVT behavior with depth or with depth
and areal location. In some cases, this variation can be represented in a black-oil
model by simply allowing initial solution gas Rsi to vary with depth in the under-
saturated oil column, retaining a single setoff Bo (p) and Rs (p) curves. In the
other cases, multiple sets of these curves and two oil components are necessary and
the single oil-type assumption in a black-oil model can lead to appreciable error.

Mechanisms or phenomena that are significant in some reservoirs and may not be
represented in the model include compaction, hysteresis in wetting and non-
wetting relative permeabilities, and interlayer wellbore cross-flow. The latter is a
particularly difficult modeling problem and the subject of continuing research. A
production well completed in a number of layers may exhibit production from
some layers and, simultaneously, injection (backflow or recirculation) into others.
Factors that promote this possibility are low-pressure drawdown (high PI and/or
low rate) and poor vertical communication between the reservoir layers in the
vicinity of the well. A rigorous treatment of this problem requires modeling of
wellbore multiphase hydraulics and phase segregation combined with calculation
of correct phase mixtures for the layers undergoing injection .
Chapter 3
The project objectives

1. Construct models for a single well radial flow reservoir,

using CMG’s IMEX simulator

2. To assess the differences between the result for recent project with the
companies caring on the same model and data.

3. Apply scenarios to obtain significant differences.

4. To reach the final proper scenario that will be and recommended by using
CMG.
Chapter 4

Input data of the project

1. Basic data

2. Rock and fluid data

3. Initial condotion data

4. Well data

5. Array properties data

6. General pvt data

7. Pvt data

8. Saturation data
 1. Basic Data
Geometry

Radial Extent, ft 2050

Wellbore Radius, ft 0.25

angullar theta division 10

Radial position of first block center,ft 0.84

Number of radial blocks 10

Radial block boundaries, ft

0.365632, 0.90038, 2.21721, 5.45995, 13.4453, 33.1094, 81.5329, 200.777,494.42,
1217.52
dipanngle, degrees 0

Depth to top of formaton,ft 9000

Number of vertical layers 15

 2. Rock And Fluid Data

Pore compressibility, psi-1 0.000004

reference pressure for calculating effect of rock compressibility = 3600 psi

Water Copressibility, psi-1 0.000003

Oil Compressibility for under saturated oil, psi-1 0.00001

Oil viscosity compresssibility for undersaturated oil, psi-1 0

Stock-tank oil density, Ibm/cu ft 45

Stock-tank water density, Ibm/cu ft 63.02

Standard-condition gas density, Ibm/cu ft 0.0702

3. Initial Condition

Presf (referencepressure)= 3600

Depth of gas/oil contact, ft 9035

Oil pressure at gas/oil contact, psi 3600

Capillary pressure at gas/oil coantact, psi 0

Depth of water/oil contact, ft 9209

Capillary pressure at water/oil contact, psi 0

Pb = pi=3600 saturated reservoir

 4. Well Data
Skin 0

Produced well completed in blocks 1,1,7 1,1,8

Minimum bottom hole pressure 3000 psi

Production Schedule

Time Period (Days) Oil Production Rate (STB/D)

1 to 10 1000

10 to 50 100

50 to 720 1000

720 t0 900 100

Layer Thickness Porosity Kx kz

5. Array properties 1 20 0.087 35 3.5
Depth to top of formation 9000 Ft 2 15 0.097 47 4.75
3 26 0.111 148 14.8
4 15 0.16 202 20.2
5 16 0.13 90 9
6 14 0.17 418 41.8
7 8 0.17 775 77.5
8 8 0.08 60 6
9 18 0.14 682 68.2
10 12 0.13 472 47.2
11 19 0.12 125 12.5
12 18 0.105 300 30
13 20 0.12 137 13.75
14 50 0.116 191 19.1
15 100 0.157 350 35
 6. The general pvt data
Reservoir temperature (TRES) 180 F

Stock tank oil density 45 lb/ft3

Gas density 0.0702 lb/ft3

Water density 63.02 lb/ft3

Under-saturated oil compressibility (CO) 0.00001 1/psi

Oil phase viscosity pressure dependence (CVO) 0 cp/psi

Water properties

Water formation volume factor (BWI) 1.01303

Water compressibility (CW) 0.000003 1/psi

Ref. pressure for water (REFPW) 3600 psi

Water viscosity (VWI) 0.96 cp

Pressure dependence of water viscosity (CVW) 0 cp/psi
 7. pvt data:

Oil Oil Solution Water Water Gas Gas

pressure Bo Bw BG
density viscosity GOR density viscosity density viscosity
psia RB/STB IB/FT3 CP SCF/STB RB/STB IB/FT3 CP MCF/STB Ib/ft3 cp
400 1.012 46.497 1.17 165 1.01303 62.212 0.96 5.9 2.119 0.013

800 1.0255 48.1 1.14 335 1.01182 62.286 0.96 2.95 4.238 0.0135

1200 1.038 49.372 1.11 500 1.01061 62.36 0.96 1.96 6.379 0.014

1600 1.051 50.726 1.08 665 1.0094 62.436 0.96 1.47 8.506 0.0145

2000 1.063 52.072 1.06 828 1.0082 62.51 0.96 1.18 10.596 0.015

2400 1.075 53.318 1.03 985 1.007 62.585 0.96 0.98 12.758 0.0155

2800 1.087 54.399 1 1130 1.0058 62.659 0.96 0.84 14.885 0.016

3200 1.0985 55.424 0.98 1270 1.0046 62.734 0.96 0.74 16.896 0.0165

3600 1.11 56.203 0.95 1390 1.00341 62.808 0.96 0.65 19.236 0.017

4000 1.12 56.93 0.94 1500 1.00222 62.883 0.96 0.59 21.192 0.0175

4400 1.13 57.534 0.92 1600 1.00103 62.958 0.96 0.54 23.154 0.018

4800 1.14 57.864 0.91 1676 0.99985 63.032 0.96 0.49 25.517 0.0185

5200 1.148 58.267 0.9 1750 0.99866 63.107 0.96 0.45 27.785 0.019

5600 1.155 58.564 0.89 1810 0.99749 63.181 0.96 0.42 29.769 0.0195
 8. Saturation data

sw krw krow pcow

sg krg krog pcgo
0.22 0 1 7
0 0 1 0
0.3 0.07 0.4 4
0.04 0 0.6 0.2
0.4 0.15 0.125 3
0.1 0.022 0.33 0.5
0.5 0.24 0.0649 2.5
0.2 0.1 0.1 1
0.6 0.33 0.0048 2
0.3 0.24 0.02 1.5
0.8 0.65 0 1
0.4 0.34 0 2
0.9 0.83 0 0.5
0.5 0.42 0 2.5
1 1 0 0
0.6 0.5 0 3
0.7 0.8125 0 3.5
0.78 1 0 3.9

Note : you can see the input data graphs in appendix.A.

Chapter 5
Results of the project

1. CMG calculation results

2. CMG results graphs

 1. CMG Calculation results

After run and validate , CMG calculate the following data :

Total oil in placeSTB 0.28998E+08

Total water in place STB 0.73907E+08

Total gas in place SCF 0.47142E+08

HC. Pore Volume M RBBL 36695

Total Pore Volume. M RBBL 110748

Field Total Fluid

Oil Gas Water Solvent Polymer

(MSTB) (MMSCF) (MSTB) (MMSCF) (MLB)

Cumulative Production 530.94 1156.0 218.32 NA NA
Cumulative Injection NA 0 0 NA NA
Cumulative Gas Lift NA 0 NA NANA
Cumulative Water Influx NA NA 0 NA NA
Current Fluids In Place 28467 45986 73688 NA NA
Production Rates .10000 .13425 .05590 NA NA
Injection Rates NA 0 0 NA NA
 2. CMG Result graphs

Fig 5.1. oil rate sc- daily production with time

Fig 5.2. gas rate sc production with time

Fig 5.3. gas oil ratio sc with time

Fig 5.4. water cut sc with time

Fig 5.5. wellbore bottom hole pressure with time

Note :You can see The changes in oil,water&gas productions with time (in every
days) and also gor ,water cut changes with time in appendix.B.
Chapter .6
Comparing CMG results with the company’s results

1. Comparing of Initial fluid in place

2. Comparing the oil production rate

3. Comparing the water cut

4. Comparing the gas oil ratio

5. Comparing the bottom hole pressure

 1. Comparing of Initial fluid in place

Oil water Gas Declined

Company (106STB) (106STB) (106SCF) (days)

ARCO 28.80 74.03 47.01 257

Chevron 28.88 73.94 47.13 217

CMG 28.998 73.90 47.14

D&S 29.11 74.97 47.11 815

Fran lab 28.89 73.93 47.09 280

Gulf 29.29 73.49 47.63 230

Harwell 28.89 73.96 47.09 232

Intercom 28.92 73.93 47.13 210

McCord-Le 28.68 74.12 46.98 257

Nolen 28.89 73.96 47.08 237

SSC 28.87 74.03 47.04 250

Shell 28.76 74.08 46.94 222

We see that the our CMG’s results are approximately between chevron and ARCO
works.
 1. Comparing the oil production rate

Fig.6.1 comparing oil rate

 3. Comparing the water cut

Fig.6.2 water cut comparing

 4. Comparing the gas oil ratio

Fig.6.3 gas oil ratio comparing

 5. Comparing the bottom hole pressure

Fig.6.4 well bottom hole pressure comparing

Chapter .7. The Scenarios

1. perforation changes

2. bottom hole pressure changes

3. water injection
 1. perforation change
Change perforation layers to (5 & 6) and ( 3&4 ) and (9&10) and then choose the
best perforation layers for the reservoir .

The following figures show the difference between the 4 casesof perforations:

Fig.7.1 perforation change affect on oil rate

Fig.7.2 perforation change affect on gas oil ratio

Fig.7.3 perforation change affect on water cut

Fig.7.4 perforation change affect on oil rate

From the above figures we see that the case of perforation layers 5 & 6 is the best
case .
2. Bottom hole pressure change
Change the bottom hole pressure from 3000 to 2500 , 2800 & 3200 psi

The following figures show the difference between the 4 cases of BHP changing :

Fig 7.5 BHP change affect on oil rate

Fig 7.6 BHP change affect on gas oil ratio

Fig 7.7 BHP change affect on reservoir bottom hole flowing pressure with time
Fig 7.8 BHP change affect on water cut

We see that when we decrease the wellbore bottom hole pressure the oil production
stay constant at 1000 bbl/day more time.
3. water injection
inject water to support pressure of the reservoir
(voidage or replacement)

Use the following information :

Injection well constraints :

BHP bottom hole pressure MAX 6000psi
BHW reservoir water rate MAX 4500bbl/day
Injection well perforation layers :
10,1,7 and 10,1,8

Fig 7.9. 3D view of the reservoir ( production and injection well location)
The following results calculated after run :

Field Total Fluid

Oil Gas Water Solvent Polymer

(MSTB) (MMSCF) (MSTB) (MMSCF) (MLB)

Cumulative Production 701.47 1810.7 390.95 NA NA
Cumulative Injection NA 0 3130.7 NA NA
Cumulative Gas Lift NA 0 NA NANA
Cumulative Water Influx NA NA 0 NA NA
Current Fluids In Place 28211 45507 76742 NA NA
Production Rates .10000 .15305 .09135 NA NA
Injection Rates NA 0 2.8971 NA NA

The following figures show the difference between base case (without
injection) and injection well case

Fig 7.10 The effect of water injection on the daily oil rate
Fig 7.11 The effect of water injection on the water cut

Fig 7.12 The effect of water injection on the gas oil ratio
Fig 7.13 The effect of water injection on the flowing bottom hole pressure

We see that the water injection increase the reservoir pressure and water cut and
oil rate and the gas oil ratio
Chapter .8.
Disscusion and Suggestion

1. Discussion

2. suggestion
1. Discussion :
1. The time of start of decreasing in the oil rate from constraint 1000 bbl/day can
be expressed by breakthrough time because of the gas conning which occur after
this time, such as the below figure , we see that the gas rate increase after oil rate
will decrease .

Fig.8.1 oil rate and gas rate with time

2-When decrease the bottom hole pressure BHP the oil recovery will increase and
then gas breakthrough time delay’s.

3. When change the perforated layer to 9 and 10 the oil recovery increase and the
gas breakthrough time delay but the water conning possibility increase.

4. Increased or reduced density difference delays breakthrough time and

also reduces GOR after breakthrough.

5. A significant effect on coning tendency is observed when there is a decrease in

oil viscosity and an increase in gas viscosity which delays gas breakthrough time
and reduces GOR after breakthrough but increases ultimate recovery.

6. An increased interval thickness results in a long breakthrough time and increase

GOR after breakthrough.
2. suggestion :

1. We suggest to perforate on block 1,1,9 and 1,1,10 because

these layers gives more oil recovery.

2. We suggest to decrease the BHP in order to increase oil

recovery .

3. We suggest to drill an injection well if the above suggestions

is not possible to do.
References :
NOMENCLATURE

A= cross - sectional area of core, ft2

Bo= oil formation volume factor, rb/stb
Bg= gas formation volume factor, rb/stb
GOR = gas-oil ratio
H= oil formation thickness, ft
K= permeability, md
Kh= horizontal permeability, md
Kv= vertical permeability, md
Kro= oil relative permeability, md
Kro= oil relative permeability, md
∆ P= pressure gradients
Pcg= gas capillary pressure
Pcgo= gas-oil capillary pressure
Pco= oil capillary pressure
qo= oil production, STB/cu ft reservoir-day
qg= gas production, STB/cu ft reservoir-day
qi= initial oil production rate
R= recovery
RS= solution gas-oil ratio, stb/stb
re= drainage radius, ft
rw= wellbore radius, ft
Sg= gas saturation, fraction
Swi= initial water saturation
t= time, days
Tbt= breakthrough time, days
= porosity
= oil density, Ib /cu ft
= gas density, Ib /cu ft
= oil viscosity phase, cp
= gas viscosity phase, cp
Appendix.A

1. graphs of the inputpvtdata :

fig.A.1. gas expansion with Pressure plot

Fig.A.2.solution gas oil ratio with Pressure plot

Fig.A.3. oil and gas viscosity with Pressure plot

Fig.A.4. relative permeability with water saturation plot

Fig.A.5. relative permeability with liquid saturation plot

Fig.A.6. gas oil capillary pressure with liquid saturation plot

Fig.A.7. deviation factor with pressure

Fig.A.8. water oil capillary pressure with water saturation plot

Fig.A.9. gas formation volume factor with pressure plot

2. graphs of the geometry and array data :

Fig.A.10. 3D veiw of reservoir grid top

Fig.A.11. 3D veiw of reservoir grid thickness

Fig.A.12. 3D veiw of reservoir permeability

Fig.A.13. 3D veiw of reservoir permeability

Apendix .B. Total Production Table of the project (base case)

zou:o4:ze iooo.ooo ioaa . aaa
zou:o4:zs iooo.ooo i4oi.=z4 i ioaa ie. az
ic.as
o ss s.
e s
is nn
a s a.
iooo.ooo wsa.uw isz.ozss wsa ic.4•s
zou:o4:ze iooo.ooo ins.izzy i4eo ie. s4
zou:o4:ao iooo.ooo isoz.ons isoz
zou:os:oi iooo.ooo isos.v•sa 2OO.Z459 isos is . sa
zou:os:oz iooo.ooo mos.a4a zoi.aoas mos ic.z•s
zou:os:oa iooo.ooo uw.uz uw ia.se
zou:os:o4 iooo.ooo uw.esi uw ic.ei
zou:os:os iooo.ooo iszo.os4 iszi ic.e‹a
zou:os:oc
zou:os:oz iooo.ooo iszs.ozs zou.eisi iszs iz.oz
zou:os:oa iooo.ooo iszs.zza iszs iz.is
zou:os:os iooo.ooo
zou:os:io iooo.ooo isaa
zou:os:ii iooo.ooo u4z.ioi is4z
zou:os:iz iooo.ooo is46. zii.ssas is46 iz.4s
zou:os:ia iooo.ooo as6 i s so . zie.zzea wsi iz. sa
zou:os:i4 iooo.ooo s s z zij.4esa i s s s iz.ee
zou:os:is iooo.ooo isso.szz eis. who
iooo.ooo use.cu szz= iz.so
wo4.rio
zou.os.iz
zou.os.ia iooo.ooo
iooo.ooo zia. aaah ia.o4
zou.os.zo iooo.ooo zzo.wei
zou.os.u iooo.ooo ia.zz
isO4.oaz moi ia. zz
zou.os.za
zou.os.z4 iooo.ooo woo.e4z ia.
iooo.ooo ui4
wei sz
ici4.ioz
woo.sz‹a ia.e4
zou:os:ze
zou.os.zz iooo.ooo
zou:os:za iooo.ooo uss.ome ia . a=
zou:os:zO iooo.ooo i•aso
zou.os.ao iooo.ooo u4z.zw
zou:os:u iooo.ooo use.O=z
io.ii
zou:os:oi iooo.ooo
zou.os.oz iooo.ooo ieao
zoia:os.-oa iooo.ooo
zou:os:o4 iooo.ooo ier=
iooo.ooo wos io.ao
ie‹ae.•ais
wos.u4
zou:oc:os
zoia:oa:oa iooo.ooo
iooo.ooo wie zo.oz
zoia.o€.oz iooo.ooo wie.<ea izzz zo.iz
zoiz.oa.oa iooo.ooo izzz.l‹aa
zou:o€:oO iooo.ooo zs€.oaiG izsz zo.ao
zoia:oc:io iooo.ooo izsz.ais zo.ao
zou.os.ii iooo.ooo wei zo.co
iaoo zo.=o
zoia.oa.iz iooo.ooo wei.ozc
iaoO.uz iazs
zou.o€.ia iooo.ooo iazs.zGo
zou:o€:i< iooo.ooo iasi z o . a a
zoia:oc:is iooo.ooo
zou.os.ie iooo.ooo iaoz
zoia.oa.iz
zou.o€.ia iooo.ooo
iooo.ooo
iacz.
ooi a a zi.ii
zou:o€:iO iooo.ooo iazz.io4 iaeo
zoia:oc:zo iooo.ooo iaai.zie
zoia.o€.u iooo.ooo z zi . aa s iooc
zoiz.oa.zz iooo.ooo zzz.‹aizi
a<s ioi4. s
zoiz:oa:zz iooc.
iooo.ooo zi. zi.
zou:oc:z4 iooo.ooo so
az•s sz
zou:oc:zs iooo.ooo iezz.iG•a zi.az
zou:oc:zo iooo.ooo is<s
issa
zi.=o
zou:oc:zz iooo.ooo isco zi.ze
zou:oc:zs iooo.ooo
zou:oc:zo iooo.ooo issz.szs iscs zi.se
zou:oc:ao iooo.ooo isoo.4sc zsi.4zsz
zou:oz:oi iooo.ooo isos.oom wsa
zou:oz:oz iooo.ooo iszs.coo iosi
zou:oz:oa iooo.ooo issa.eis za4.ai4i uss zz.io
zou:oz:o4 iooo.ooo isoo.som zooü zz.iz
zou:oz:os iooo.ooo isos.aos zai.ia64 zou zz.z4
zz.ai
zou:oz:of iooo.ooo zaa.zz6i
zou:oz:oz iooo.ooo zas. aazi
zou:oz:os iooo.ooo zou.aoi
zou:oz:of iooo.ooo zsi.siis zz. sz
zou:oz:io zz. ss
zou:oz:ii zo3= zz . aa
zou:oz:iz zo4a.isi zo4. saia
zou:oz:ia hoi.os4o zoso.liz
zoso.aio ao=o
zou:oz:i4 zou.oos zoza
zou.4ao zoai
zou.toa
zosi.a4a za.ii
wow:ıo:oa so‹a
wow:ıo:o4 ısas.a4ı
wow:ıo:os sze.ııaz ısao.4zı ıaoa.cao
wow:ıo:oc az4. ıoa4.sos
wow:ıo:oz ıeaz.zz4 az4.aıea
wow:ıo:os sıas mao.s4s
wow:ıo:os ıozs.sza
azı . z aaa
wow:ıo:ıo sis.zsos ıozz.oos azs.cıoı az6.ızaa az6.6zaı azı.ıosz
zow:ıo:ıı
zow:ıo:iz ısza.
zow:ıo:ıa su.csoı asz
zow:ıo:ı4 ıszı.4ss aza.o4aı
zow:ıo:is aıa.6aıs ısıs.czc
zow:ıo:ıc siz.ıısı ısır.zcc
zow:ıo:iz sıo.kızı ısır.su
zora:ıo:ıs sos.ızoo ısıa.ssz
zora:ıo:ıs ısız.ıcs aza. sızo
zora:ıo:zo son.wa4 ıoıo. a44 zz9.ızze
zora:ıo:zı 8O4.?OZ8 ıoos.sza
zora:ıo:zz 8O3.ZB3O woo.toz
soı.saoa ıoo4.ssa
zoız:ıo:zz ıooa.oco
zow:ıo:z4 zzo.izi= aao. huo
zora:ıo:zs ıOoı.z€o
ıs o o .
zora:ıo:zc zzo.saaı zzo.se4ı zzı.oazı
zora:ıo:z= a s = . aaı.zili
zora.ıo.zs A s z ıas
zora:ıo:zo s s s
zoız:ıo:ao
zow:ıo:ai s . aaa
zora:u:oi ıasz. aaı.osoA
zora:u:oz aoz
zora:u:oa
zora.u.o4 =a s . asia ısoo.szo
zora:u:os ıssc.oD‹a zzz. sas
zora:u:o€ ı a as
zoıa:u:o= ısı.w‹aao ıa a z
zora:u:oa =so.zısı . z az zaz.zsıs
zora:u:oo
zora:u:ıo .i s z
zora.u.u za4.oıso
ısso.oıı
zora:u:iz
zora:u:ra
zora:u:ı
zora.u.is
zora.u.ı€ ismi.zzo
zora.u.ız ıa=o.osz
wow.u.is ı aaa .
wow.u.ıG aaa zzs. asıs
wow.u.zo
wow.u.zı zzs.szoı
wow:u:zz ıacc.mız
ıaas.ooo
wow:u:za
wow:u:z
wow:u:zs mai.osza ıaaı.zoa
ıs€o.oso
wow:u:zc asa. zss
zora:u:zz ı a s s . 4 a
zora:u:za ısss.ıoa
zora:u:zo ıasz.4a4
zora.u.zo zza.z4sı
zora.iz.oi
zsz.ıe4z ısso.ı‹aa aaa.oıoz
wow:iz:oz ian-a. siz
zora:ız:oa
zora.iz.oh ıa4s.
zora.iz.os zzı
wow:iz:o‹a
wow:iz:oz ıa4z.ıoo
zoıa.ız.oa ı a 4o
zora.iz.oo ıa z
wow:iz:ıo . so
wow:iz:u a .a s
zora:ız:iz o
zora.iz.ra i
zora:iz:ı4 ıaaz.zOa
wow:iz:is ıaas.zoz
zora:ız:ıc
zora.iz.iz ısa4.ııo
zora.iz.is ıaaz. szı
wow:iz:ı‹a zzs.ıasz ısao.‹aa€
wow:iz:zo
zoıa.ız.zı
ıszc.ıO‹a
=zo.ısıa ıszs.onu

ısır.aos
ısır.z4P
ısız.uc
ısıo.
soa
ısos.osz
ısoz.
za4z zs.ıa a
o
za44 zs.is .
i
s
zası za . za

zwza
zo.ai

zaso zo.ıo
zasz zo.is
z%sı
ao.es
sı.oa
zı.o=
ai.ıo
ai.ı4
sı.ıa
zwoo
ai.zı zı.zs

c
h
o zı.
= zı.
sa
s=

z
e
.
a
i

z
4
ı
L

z
o
.
o
o
3O.OZ
ao.oz
z
o
ı
o
z
o
.
u
z
o
z
ı
ZOL3:JL2:22
zora:iz:za zza.zoıe ıaoı.
ıaoz.a94
aoo ai.e4
ZOL3: 2:24 ız9o.azı
ZOL3: 2:25
ızsa. ai.zo
zora:iz:zo ooo
ZOL3:JL2:2Z
zora:iz:zs ai.sı
ZOL3:JL2:29 w‹as.z4a ai.a4
ZOL3:i2: azı.4z94
3O ız9z.zor ai.sı
zora:iz:ai ızso.ca< ai.o4
ZOL4:OH:OL ma. saz4 ızso.ıoc
zoı4:oi:oz wan.eso
ZOL4:OH:O3 mae.ıze
ZOL4:OH:O4 ıza4. z4sı az.oz
zı4.aziz sam
zoı4:oi:os zıa.szc4 i? z4oz az.ıo
83.IL4 2493 32.L3
zoı4:oi:oz zıa.oosı mao.o9a
ızsı.cos
ZOL4:OH:O8
ZOL4:OH:O9 imza.siz
zoı4:oi:ıo mo.z6za wzz.oaa
was.
wow.oi.u sas
wow.oi.iz zoa.4sıo
wow.oi.ra zsoı az.a€
wow.oi.ı4
wow:oi:is wei.ooo
wow:oi:ı€ m€s.sao
wow:oi:iz zos.Paıo ız€a.oui
wow:oi:ıa w€€.
zoıa:oi:ıO zoz.zısa saz
zoıa:oi:zo w€s.oas
zoıa:oi:zz ız€z.oPs zsıo az.aa
zoıa:oi:zz mao.aoı
zoıa:oi:za msO.ıoO zsız az.zo
zoıa:oi:zs ızsz.aıa zsız
zoı4:oi:ze eee.wo= asia
zoıa:oi:z= s='s.zam ızsa.ıza
ızs4.e4o aao.ıooo zsıs
zoıa:oi:za ızsz.ısı aao.ıs=ı zsıe
wow:oi:zo wsı.ccz aah.zzoı ısı=
wow:oi:ao ı=so.ian asia zz. aa
zoı4.oi.ai
wow.oz.oi •aoı.amos aah. mor zsıe az.sı
wow.oz.oz
wow:oz:oz
zoı<:oz:o oea. sıu
wow:oz:os
zoı4.oz.o€
wow.oz.oz
wow.oz.oa imza. zsa
wow:oz:oO
zoı<:oz:ıo
wow:oz:u aaz.
zoı4.oz.iz zszo
wow.oz.ıa zaıı
wow.oz.ı4
wow:oz:is ızza.•s‹as
zoıa:oz:ıc oaı.ismi zzs.
wow:oz:iz sıs4 zszz aa. zı
zoı4.oz.ıa
wow.oz.is
zoı4:oz:zo wzı.ız•a
zoı4.oz.zı izin.
zoı4.oz.zz eıa azs.
zoı4:oz:za ama.oızs i ma.
zoı4:oz:z4 aso saeı
zoıa:oz:zs
zoı4:oz:zc wu.ın
zoıA:oz:z= zsaı za. s=
zoı4:oz:za
zoı4.os.oi em.eız4
zoı4.oa.oz soo.aazı ı=ıo.‹aas zae.=ı4e
zoı4:oa:oa mor.cıo
zoı4:oa:o4 ı=os.zos
zoı4:oa:os wow.‹acs aao.sıoı
zoı4:oa:oc ızos.eza
zoı4:oa:o= woa.zaa
zoı4:oa:oa eee. sus
zoı4.os.os ses.aıoz ı=oz.ssa
ı=oı. aaa
zoı4.oa.ıo woo.zaı zao.ooh
wow:oa:u ıcos.s=ı
wo=.soo
wsc.u4
ıasz.4oa
ıasz.oaz
zoi4:oa:iż
zoi4:oa:ia ccz.zz3z icas.zsa zssi aa.sa
zoi4:oa:ix •sci.mo4 icaz.szo
zoi4:oa:is aao.zeas iaaa.s4a a4o.osia
zoi4:oa:ie eeo.osei zoo.izas
zoi4:oa:ia
zoi4:oa:io ess.spis ieaz.4zz
ieai.o4s a4.ii
zoi4:oa:zo
zoi4:oa:zi is=a. a4.i6
zoi4:oa:zz zoo a4.ia
zoi4.oa.za •sss.isci iczc.szt
zoi4:oa:zA iszs.
zoi4:oa:zs ssz
zoi4:oa:ze iaz4.iz‹a
zoi4:oa:z= esz. azis iezz.ao4
zoi4:oa:zs esi.e=ss zssi
zoi4:oa:zs iezo.osi
zoi4:oa:ao ie e a.
zoi4:oa:ai i s s=
64s.6izz s o s
zoi4:o4:oi . zs a
zoi4:o4:oz iccs.sis a4o.ssio
zoi4:o4:os iaaa.iaa a4o.sczi
zoi4:o4:oc ich.zas ono. soai
zoi4:o4:oz icco.4oz
zoi4:o4:oa iess.ono
zoi4:o4:oo
zoi4:o4:io a4o.o4is
zoi4:o4:ii s4z.max ics4.sos
zoi4:o4:iż
•s4i.4sss iesz.ise
zoi4:o4:i4 ieso.za4
zoi4:o4:is e4o.iszz ic4s.4i4 a4o.sos
zoi4:o4:ic ic4a.o4s
z4.ai
zoi4:o4:ia •saa.izso icAs.cii
zoi4:o4:io
zoi4:o4:zo ic4a.ono
ic4z.ssc a4o.aoia
zoi4:oA:zi aza.miss ic4i.z4a
zoia:oa:zz icae.eoi zAo.azoa
zoiw.od.zz ea4.eiza i o z a . s s s zao.aasi
zoia:od:zs
zoi4:o4:ze
icas.azo
zoia.oa.zz iez4.
zoiw:od:za em.samo swe
zoi4:o4:zo osi.oraz icaa.no
zoi4:o4:ao
zoi4.os.oi iaao.
zoi4.os.oz nia
eze'.izoz iczo.iaa
iez=.asi
iaza.
szo
iszz.asz
iezz.
sze
i€zi.cii
zoiw:os:oz azz.zizo niż.o=zo i€ie.
zoi4:os:oA eze. a a4o.lou
zoi4:os:os s=9i cii.a‹nz sas isoz.
zoi4.os.os ozc
zoi4:os:oz a icis.Oz uOz.O
zoi4:os:os z oa
zoi4:os:oo i€i4.ai4 isoi.
zoi4.os.io iaia. 4co
isso.
zoi4:os:ii
zoi4:os:iż
zsz
aoa isao. a4o.ez4i a4o.oaii
zoi4:os:ia ieiz.oo= oi=
zoi4:os:i4 azi. icio.=o=
zoi4:os:is wasz i€oO.mo isae.
zoi4:os:ia au.ooio au.oiaz au.oaoi
icoa.u4 sz4
zoi4:os:ie eis.
zoi4.os.iO icoa.az
zoi4:os:zo zzzs o ieos.
zoio.os.zi €io.izoz sza
zoia:os:zz ico4.za au.oasi
z
€iz.z‹azo ieoz.s4z
zoio:os:zz icoi.Oa
€ia.oaOa a i€oo.
zoiA:os:za bis.4sos zaO
isoo.im
zoi4:os:zc ei4.aesz
isoa.
seo uOs.
zoi4:os:za ci4.zzoi
€ia.js4o ai4
iso4.o4•
 zs.iż
zs.i4
zsao
zsao
zsei zs.zi
zsai

zsoi zssi as. aa

zs.=o
ZOU:OK:29
ZOL4:O9: cii.aois isas.z9i z4i.oasz zsez as.ai
3O ust.o4s au.osoi
ZOU:O5:3L
ZOU:OK:O1
mo.izoi isaz.aoa au.os44
isar.sea z4i.oeas
ZOU:OK:O2 isao.zze 34i.LO24
ZOU:OK:O3 isze.o9z 34i.
ZOL4:OK:O4 iszz.asc LOGO
ZOU:OK:O5
ZOU:OK:OK
sos.iazz isze.ezi 34i.
LO94
iszs.
ZOU:OK:OK Ci O €i . O L4 zaa z4i.iiz4
ZOL4:OK:O8 2 isz4.is6 34 . L 32
ZOU:OK:O9 cos.4aiz moż.9z6 34i. Li?6
zoi4:oc:io 34i.
ZOU:OK:1L iszo.4oo LJL98 2B99 3€. LO
ZOL4:O€:i2 isee.z4a a4i.izie 2BOO 3O. L3
ZOL4:O€:JL3 BO3. L22 JL9O8.OL8 z4i.izzz 2 C' O O 3 C' . L 5
zoi4:oc:i4 isee.z94 34 . L243
ZOL4:OD:18 eoi.esae ises.szz 34- . L2.3.3 2€O1 3€.L9
ZOL4:OD:1B eoi.zaa4 ise4.zso 34i. L2'i
ZOU:O‹b:i? issa.izi z4i.izes
coo.zaci i5€L.9L2 34i. L2€
zoi4.oa.ia €
zoi4:oc:io isao.• zcoa aa.4i
zoi4:os:zo as i s ss . 34i.L2'G8
34i.L2'6O
zoi4:os:zi dwa
zoi4:os:zz isse.zeA zm.izsz
zoia:oa:za iss=.oso
zoi4:os:z4 i s s s . a za zm.izza
zoi4.oa.zs
zoi4:os:zs
zoi4:os:zz s s . za= issz.ni=
moż.zos
zoi4:os:za issi.ooh
zoi4:os:zs so4.scio is4o.zoz
zoia:oa:ao isoa. soz zm.ioaz
zoi4:o=:oi
zoi4.oz.oz isoz. zao au.io4o
ssz. au.ioii
isoo.oaa
zzio aai.oszz
zm.oaai
zm.oazz
zoi4:o=:oz
zoi4:oz:o4 soi.ross isoz. sao zm.o=ai
zoi4:o=:os sso . oss is4i.a‹az zeio
zsio ze.si
zoi4:oz:os zm.osso
iszs.ooh zm.oeio ccii
zoia:oz:oz
zoi4.oz.oa isa€.•aiz zai.oaaz zaii
zoi4.oz.oO iszs.aza au.o4i4 zaiz za.=i
zoi4.oz.io zaiz
zoia:oz:ii isaa.on‹a zaiz
zoi4:oz:iż
zoi4:oz:iż sec.stoi iszo.az•s au.oiia zaia
zoia.oz.is sas. o isze.ae:z z4i.oozs
zoia.oz.ia sa4.•aoi‹a isza.
zoe zais
iszz.izz
sa z .s ix s i s zs . st a zais za . aa
zoi4:oz:ia
zoi4:oz:io iszz. zaia
zoi4.oz.zi
zoia.oz.zz sai. saa
zsio
sto.
aiza iszi.taa
iszo.osi
is i a . aaa a4o.oni zaia zz.oi
szO. zzai i s i a . st s zaie
zoi4:oz:zs zaia
zoi4:oz:za mix. zoio
sza.izzs aza
isia.zoz
zoi4:oz:za sca.isei isia.o4 isoo
.zał
zoi4:oz:zs o zao.asis
mu.azs
zoi4:oz:ai
zoi4:oa:oi i s i
zoi4:oa:oz
i s o s ioos.
o . m o a4o.aoie
zoi4:oa:oa . s 4 aas
zoi4:oa:o4 z
zoi4:oa:os isoa. usa.
zoi4:oa:oc zał
zoi4:oa:oz i .s zo zz
asz
zoi4:oa:oa szz.
zoi4:oa:os
zoi4:oa:io ai‹ae
.moocc4 i4si a4o.miss
4
zoi4:oa:ii iso4.sos .io4
zoi4:oa:iż
zoio:oa:iż moż.ss
zoi4:oa:i4 a
zoi4:oa:is isoi.oa
zoi4:oa:ic e
i4sc.aaa
zezo az.ii
zz.iż
zezo zczi zczi
az.is
zezi
zow:oa:iz
zow:oa:is mas.su
zow:oa:is scc.haz ıoaa.azı
wan.mai
zow:oa:zo a4o.
zoı4:oa:zı
zoı4:oa:zz s4zı
zoıo.oa.za
zoıo.oa.z4 ıoaz.izi
zoı4:oa:zs ıoaz.ooh
zoı4:oe:ze mao.arz
zoıo:oa:zz sen. gezi
zoıo:oa:za ı4za.sis gezi
zoı4:oa:zs szı4 gezi
zoı4:oa:ao
zoı4:oa:ai seı. alo.zazı
wow:oe:oi wz4.us
zoı4:os:oz szoı
zoı4:os:oa seı.oısa
zoı4:os:oh seo.
zoıo:os:os size ı4es.ses
zoı4:os:oa
zoı4:os:oz
zoı4:os:oe w€a.zao
zoı4:os:oo w€s.wz s4o.ıPee
zoı4:os:ıo w€4.o€s
zoı4:os:u
zoı4.os.iz ı4€ı.a4a
zoı4.os.ra woo.zao s4o.ıııo
zoı4.os.ı4
wow.oo.is
wow.oo.ıa
wow.oo.iz wsa. s4o.one aa.u
wow.oo.is
wow.oo.ıO zzo
ws4.ın
zoıo.os.zo
wow.oo.zı
zoıo.os.zz aaO.lou aa.zı
zoıo.os.za
wow.oo.z4
wow.oo.zs za. z
zoı4.oo.z€ za a . z a
zoı4.oG.zz zaaı za. zo
zoı4.oO.za
wow.oO.zo s4O.ıaz‹a
wow:oO:ao
wow:ıo:oi s4a.ına
wow:ıo:oz
zoı<:ıo:oz
zoı<:ıo:o< su a . =s ı z
zoıa.ıo.os saa. zam ı<aa.•oaı
zoıa.ıo.o•s
zoıa.ıo.oz ıazs.aoz
zoı4.ıo.oa s4s.
wow.ıo.oo zıza zaas za.sı
wow:ıo:ıo
wow:ıo:u
wow:ıo:iz ıazo.zzz
mzo.we
zoıa:ıo:ra
zoıo:ıo:ı<
zoıa.ıo.is saz. sue imza.
ona
ı<zz.ooo zz‹a. zzız
s4ı .s a4 s ı<zs.‹alo
wow:ıo:ıc
wow:ıo:ı= su.
sızo
wow.ıo.ıa su.o4oı ı4za.aoe
wow.ıo.ıO
zoı4.ıo.zo zawa za.=ı
wow.ıo.zı ıozo.au
zoı4.ıo.zz szo.ıooı ı4ıo.
zoı4.ıo.za sza.
zoı4:ıo:z4 aezı ssa zaso ze.ai
zoı4:ıo:zs ı4ı€.
zoı4:ıo:z€
wow:ıo:zz aza
wow:ıo:za zası za.aa
wow:ıo:zo sz=. wu. aıo zssı za.aa
wow:ıo:ao zezı zcsı
wow:ıo:ai ı4ıa.zoo aa.Oo
wow.u.oi wu.ıo4
wow.u.oz ı4ıo.osa
wow.u.oa mor.oo4 zaa.=sız
wow.u.o4 wo=.ess zaa.=ıOs
wow.ooa
wow.aaı
szz.ızo= wo4.aıa
szz. ama
szz.zız=
szı.=sez
zon:u:os
zow:u:oo
zow:u:oz woz.cao
ıaoı.aaa
265d 30.Od
zow:u:os sza.ızıa
szz.‹szoı woo.‹sas aaa.asaı
265d 30.06
2654 39.O8
zow:u:os szz.zıss wss.los aaa.aooı as.ıo
zow:u:ıo 526. ıass.s-o 3 38 . 4 2 0 as.iz
Ü69T• ıas-.
caz as.ıa
526.32O2
wow:u:u ses.ama ıaso.nsa 338.3603
wow:n:ız za as.is
zow:n:ıa 525.423Ü wss.<sı
ıasa.oza aaa.zsao
ss zs
ss
as.iz
as.is
zow:u:ıa 524. nsa.ass aaa.zssı zs ss as.zı
zon:u:is D.’64 wsz.ama 338.
zow:u:ıo 524. zs ss
zow:u:ız 52DO wsı.aaa Z2d'9
zow:u:is saa.oaaı ıaso.aos aaa.ısoa
ıaas.zzz aaa.ısoz
zow:u:is
sza.ısa‹s ıaas.zso aaa.ızıs 338.O866
522 . 5 OO
S 2 2 . 30 8 0
zow:u:zo 3 zı . ss 3s
wow:u:zı szı.4z4a ıas-.zzs aas.o5ıo zoss as.ai
wow:n:zz 52O. OB3€ mao.zoı aaa.oıoa
zow:u:za 52O. was.ws 33Z.O8OB
wow:u:za 543Cı ıaaa.isa
ıaaa.ıaa
zoıa.ıı.zs szo.ıo<a
wow.u.z€ izah. we
wow:u:zz sı s . a as a mai.ıoz zaz.•sısz z a as us .r ı
zoıa:u:ze ızao.oaa
wow.ıı.zO razı.oui z aa a zs . oı
zoı4.u.ao sıa. zs#a ızza. zaz. soıa zaaı
wow:iz:oi siz.‹aıoz osa zaaı
zoıa:iz:oz ızza.ozs
wow.ız.oa siz.osoa
wow.iz.o4 su.oızo ıazs.ozc
wow.iz.os sı s . z s z
wow:iz:o€ ıaza az.oo
zoıa:iz:oz s ızzz.ooa za.ıa
wow.iz.oa sis. z z z .ı s as
wow.iz.oO zzız ı aaa ao.sı
wow:iz:ıo zz=.ooıo
zoıa:iz:u su.ocıa ao.oo
w€z.ooı
sıa.oaza ıaoc.oos
zoıa.iz.ıa
wow.iz.ı4 siz.ao4z wos.ooo
ıza4.oıa
zoı4:iz:is siz.ısır izah.oza
zoı4.ız.iz su.Gooa ızez.ozz
zoı4:iz:ıa su.azar sa.z4osı
ıaao.ona sı
su.osza . = saso
zoıo:ız:zo sıo.szoı
mo.zoas
ızss.ooo
ızsz.ıoo es.
zoı4:ız:zı ızsa.u= zızee
zoıo:ız:zz
zoıw:iz:zz ızss.soo
zoıo:ız:zz ızss.osz
zoıa:ız:zz
zoıa:iz:zz soa.zzaı ıs ı .
zoıo:ız:zz soa.owzı s s .
zoıa:ız:za ıoo.oooo
zoıa:iz:zs ıoo.oooo
ıoo.oooo z e k i
ıoo.oooo s u s a
ıoo.oooo
ıoo.oooo isa.sıra
zoı4:ız:ze ıoo.oooo ı sz. zass ısza o.oo
zoıa:ız:z= ıoo.oooo ısz. zos=
zoıa.iz.za ıoo.oooo woo 4o.is
ıszz
View publication stats