COMPOSITIONAL MODELING OF MULTIPHASE FLOW AND ENHANCED OIL RECOVERY PROSPECTS IN LIQUID-RICH
UNCONVENTIONAL RESERVOIRS
by Najeeb S. Alharthy © Copyright by Najeeb S. Alharthy, 2015 All Rights Reserved A thesis submitted to the Faculty and the Board of Trustees of the Colorado School of Mines in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Petroleum Engineering).
Golden, Colorado Date
Signed: Najeeb S. Alharthy
Signed: Dr. Hossein Kazemi Thesis Advisor
Signed: Dr. Ramona M. Graves Thesis Advisor
Golden, Colorado Date
Signed: Dr. Erdal Ozkan Professor and Interim Department Head Department of Petroleum Engineering
ii ABSTRACT
Production of tight oil from shale reservoirs in North America reduces oil imports and provides better economics than natural gas. Thus, many companies have directed their e↵orts to liquids production from Bakken, Eagle Ford, Niobrara, etc. Bakken recoverable reserves is estimated to be 7.4 billion barrels. Despite advances in technology, the oil recovery factor remains low (4% to 6%) (Energy Information Administration, 2013). To produce these liquid-rich shale reservoirs e ciently, a thorough understanding of flow mechanisms, reservoir properties, and rock and fluid interactions is necessary. This work will present two areas for investigation. First, this research work presents compositional modeling of liquid-rich unconventional reservoirs using volume balance method. A 2D three-phase single and dual-porosity models using volume balance method are developed and presented. Due to the explicit nature of the phase saturation calculations, a discrepancy in the number of moles in the system was observed and a “mole correction term” was introduced to rectify the material balance error for the system. Since the volume balance approach uses partial molar volume as weighting factors, a robust partial molar volume algorithm is presented and validated against published experimental data by Wu and Ehrlich (1973). The volume balance dual-porosity model aforementioned is used to model depletion of liquid-rich unconventional reservoir going below saturation pressure and the model results are validated with CMG GEM compositional simulator. Finally, the model is used to study multiphase flow regimes observed in liquid-rich reservoirs in the field. The analysis helps decipher multiphase bilinear and multiphase linear flow regimes using compositional flow rate–normalized pressure analysis from the volume balance method. From the analysis, the e↵ective permeability and hydraulic fracture permeability is calculated and presented.
iii Second,theenhancedoilrecoverypotentialofliquid-richshalereservoirswasevaluated using laboratory data from experiments conducted at Energy and Environmental Research Center (EERC) on several Bakken core samples of di↵erent size. We present both laboratory and numerical modeling of carbon dioxide (CO2)oilrecoveryfromtheseBakkencores.We also evaluate the EOR potential of using produced associated gas for injection. In laboratory experiments CO2 injection recovers higher than 90% of oil from several Middle Bakken cores and up to 40% from Lower Bakken cores. To decipher the oil recovery mechanisms in these experiments, a numerical compositional model was used to match laboratory results. We concluded that CO2 injection mobilizes matrix oil by miscibility and solvent extraction – leading to counter-current flow of oil from the matrix instead of oil displacement in the matrix (which is the conventional EOR wisdom). Specifically, the controlling factors include re-pressurization, di↵usion-advection mass transfer, oil swelling, and viscosity reduction. Laboratory results were scaled to field application in a North Dakota Middle Bakken well. The primary oil depletion period was history matched and oil production was forecasted for 10 years, recovering 6.2% oil. Then, we devised an EOR protocol using hu↵-and-pu↵ supercritical CO2 injection and natural gas liquids (NGL). Approximately 5% additional oil was produced by CO2 solvent and 6.25% by NGL solvent for fracture spacing of 500 feet. We believe oil recovery will increase further with closer fracture spacing.
iv TABLE OF CONTENTS
ABSTRACT ...... iii
LISTOFFIGURES ...... x
LISTOFTABLES...... xiii
LISTOFSYMBOLS...... xiv
LISTOFABBREVIATIONS ...... xxi
ACKNOWLEDGMENTS ...... xxii
DEDICATION ...... xxiv
CHAPTER1 INTRODUCTION ...... 1
1.1 BACKGROUND AND PROBLEM STATEMENT ...... 1
1.2 OBJECTIVES ...... 4
1.2.1 Compositional Rate Transient Analysis in Liquid-Rich Shale Reservoirs . 4
1.2.2 Appraisal of EOR potential in Liquid-Rich Shale Reservoirs ...... 4
1.3 CONTRIBUTIONOFTHESTUDY ...... 5
1.4 THE ORGANIZATION OF THE THESIS ...... 5
CHAPTER2 LITERATUREREVIEW ...... 7
2.1 COMPOSITIONALMODELING ...... 7
2.2 COMPOSITIONAL RATE TRANSIENT ANALYSIS ...... 8
2.3 ENHANCED OIL RECOVERY IN UNCONVENTIONAL RESERVOIRS . . . 9
CHAPTER 3 COMPOSITIONAL MODELING ...... 13
3.1 VOLUME BALANCE METHOD (VBM) ...... 13
v 3.2 FORMULATION AND IMPLEMENTATION OF VBM IN DUAL-POROSITYSYSTEMS ...... 15
3.3 THERMODYNAMIC MODEL ...... 17
3.3.1 PhaseEquilibriaandFlash...... 19
3.3.2 Equation of State (EOS) (PR 1976) ...... 19
3.4 VALIDATION OF THE THERMODYNAMIC MODEL ...... 20
3.5 VALIDATION OF THE VBM COMPOSITIONAL MODEL ...... 21
3.5.1 Model Parameters and Setup ...... 25
3.5.2 Results and Comparison ...... 25
CHAPTER 4 COMPOSITIONAL RATE TRANSIENT ANALYSIS ...... 31
4.1 MODIFICATION OF RATE TRANSIENT ANALYSIS ...... 31
4.2 RATE-NORMALIZED PRESSURE EQUATION ...... 31
4.2.1 Single-PhaseFlow...... 32
4.2.2 Multi-phaseFlowBlackoilCase ...... 34
4.2.3 Multi-phaseFlowCompositionalCase ...... 36
4.2.4 SummaryoftheAnalyticalSolutions ...... 38
4.3 LIQUID-RICH UNCONVENTIONAL RESERVOIR CASE STUDY . . . . . 39
4.3.1 Model Parameters and Setup ...... 40
4.3.2 FluidParameters ...... 42
4.3.3 Rock-FluidParameters ...... 42
4.4 CASE STUDY RESULTS AND ANALYSIS ...... 42
CHAPTER 5 ENHANCED OIL RECOVERY - LABORATORY AND FIELD STUDY...... 49
5.1 SUPERCRITICAL FLUID EXTRACTION ...... 49
vi 5.2 BAKKEN CO2 SOAKING EXPERIMENTS ...... 50
5.2.1 Laboratory Experiments and Experimental Procedures ...... 50
5.2.2 FluidSystemProperties ...... 51
5.2.3 Bakken Core Description ...... 52
5.2.4 Laboratory Results ...... 57
5.3 MODELINGEXPERIMENTS...... 58
5.3.1 Laboratory Model: Grid System ...... 58
5.3.2 LaboratoryModel: FluidSystem ...... 60
5.3.3 LaboratoryModel: Rock-FluidSystem ...... 60
5.3.4 LaboratoryModel: HistoryMatching ...... 60
5.3.5 DiscussionofLaboratoryResults ...... 65
5.4 MODELINGFIELD ...... 65
5.4.1 FieldModel:GridSystem ...... 67
5.4.2 FieldModel:FluidSystem...... 68
5.4.3 FieldModel: Rock-FluidSystem ...... 70
5.4.4 FieldModel:HistoryMatching ...... 73
5.4.5 Field CO2 Enhanced Oil Recovery Scheme ...... 73
5.4.6 Field NGL Enhanced Oil Recovery Scheme ...... 77
5.4.7 DiscussionofFieldResults...... 77
CHAPTER 6 MASS TRANSFER MECHANISMS ...... 81
6.1 TRANSPORT MEANS ...... 81
6.2 ADVECTIVE FLOW ...... 81
6.3 MOLECULAR DIFFUSION FLUX ...... 82
vii 6.3.1 Maxwell-Stephan Model ...... 83
6.3.2 GeneralizedFick’sLaw ...... 85
6.3.3 ClassicalFick’sLaw ...... 86
6.3.4 Di↵usion Coe cients Correlations ...... 86
6.3.4.1 Wilke(1950) ...... 87
6.3.4.2 WilkeandChang(1955) ...... 87
6.3.4.3 Sigmund (1976a, 1976b) ...... 88
6.3.4.4 HaydukandMinhas(1982) ...... 89
6.3.4.5 Renner (1988) ...... 90
6.3.4.6 RiaziandWhitson(1993) ...... 90
6.3.4.7 Maxwell-Stefan (MS) Multicomponent Molecular Di↵usion Coe cients ...... 90
6.3.5 Di↵usion Coe cients Calculations (Bakken Oil) ...... 92
6.4 GRAVITYDRAINAGE ...... 92
6.5 UNDERLYING EFFECTS OF TRANSPORT PRINCIPLES ...... 94
6.5.1 OilSwellingandViscosityReduction ...... 95
6.5.2 Reduction of Interfacial Tension (IFT) at the matrix-fracture interface ...... 95
6.5.3 Better CO2 Miscibility with Lower Temperature at Matrix-Fracture Interface ...... 96
6.5.4 SummaryofUnderlyingTransportPrinciples ...... 96
CHAPTER 7 CONCLUSIONS, RECOMMENDATIONS AND FUTURE WORK . . . 97
7.1 MULTIPHASE TRANSIENT ANALYSIS IN LIQUID-RICH SHALES . . . . 97
7.2 ENHANCED OIL RECOVERY IN LIQUID-RICH SHALES ...... 97
viii 7.3 RECOMMENDATIONS AND FUTURE WORK ...... 98
REFERENCESCITED ...... 100
APPENDIX A - COMPOSITIONAL MODELING USING VOLUME BALANCE APPROACH...... 105
A.1 Volume Balance Formulation for Single-Porosity and Dual-Porosity . . . . . 105
A.2 Derivation of Compositional Equation and Pressure Equation ...... 105
APPENDIX B - THERMODYNAMICS ...... 109
B.1 Peng-Robinson Equation of State ...... 109
B.2 Fugacity ...... 110
B.3 Derivative of Fugacity with respect to Pressure and Composition ...... 110
B.4 Derivative of Compressibility Factor with respect to Pressure and Composition ...... 111
B.5 Partial Molar Volume ...... 111
B.6 FluidCompressibility ...... 112
APPENDIX C - SATURATION EQUATIONS ...... 113
C.1 LiquidandVaporEquations ...... 113
C.2 Derivation of Saturation Equations ...... 113
ix LIST OF FIGURES
Figure 1.1 Significance of tight oil production (Energy InformationAdministration,2013)...... 2
Figure3.1 Generalvolumebalanceimplementation...... 18
Figure 3.2 Phase envelope for C1 =0.70 ,C4 =0.20 ,and C10 =0.10...... 21
Figure 3.3 Thermodynamic validation between developed routine for density and z factor calculations with CMG PVT Package (WinProp)...... 22
Figure 3.4 Thermodynamic validation between developed routine for fugacity and viscosity calculations with CMG PVT Package (WinProp)...... 23
Figure 3.5 Thermodynamic validation of partial molar volume and fluid compressibility calculations...... 24
Figure3.6 Relativepermeabilitycurves ...... 26
Figure 3.7 Validation of pressure and saturation profile ( VBM vs CMG GEM simulator)...... 28
Figure 3.8 Validation of cummulative oil and cummulative gas (VBM vs CMG GEM simulator)...... 29
Figure 3.9 Comparison of material balance error (VBM vs CMG GEM simulator). . 30
Figure 4.1 Refined gridding and well dimensions for multiphase depletion model. . . 41
Figure 4.2 Phase envelope and component specification for tight oil system...... 43
Figure 4.3 Case study validation for multiphase flow depletion run...... 45
Figure4.4 Di↵erent flow regimes in stimulated horizontal well...... 46
Figure 4.5 Deciphered flow regimes and dual porosity feature...... 47
Figure 5.1 Enhanced oil recovery experiments on Bakken Cores (performed at EERC)...... 51
x Figure 5.2 Compositions of separator samples and produced streams for Middle Bakken...... 53
Figure 5.3 Produced composition stream for Lower Bakken core...... 54
Figure 5.4 Thin sections for Middle Bakken core at di↵erent resolutions, mineralogy composed of abundant monocrystalline quartz grains (white color) with non-skeletal calcerous grains, minor calcite and Fe-Dol (tan and brown color), and some K-spar, Plagioclase, and Pyrite (black color)...... 55
Figure 5.5 Thin sections for Lower Bakken core at di↵erent resolution, mineralogy composed of quartz and calcite dominated (white and tan color), minor amount of clays such as illite (dark brown color), and kerogen patches (black color)...... 56
Figure 5.6 Oil recovery factor for Middle Bakken core soaking experiment...... 57
Figure 5.7 Oil recovery for Lower Bakken core soaking experiment...... 58
Figure 5.8 Single-porosity radial grid system used in Bakken core CO2 soaking experiments...... 59
Figure 5.9 Middle Bakken synthetic lumped fluid composition and phase envelope. . 61
Figure 5.10 Lower Bakken synthetic lumped fluid composition and phase envelope. . . 62
Figure5.11 Relativepermeabilitycurves...... 63
Figure5.12 Fracturerelativepermeabilitycurves...... 64
Figure 5.13 History match results for Middle Bakken and Lower Bakken CO2 core flooding experiments...... 66
Figure 5.14 Reservoir dimensions (single-stage HF) for a North Dakota Bakken well model...... 67
Figure5.15 MiddleBakkenreservoirfluiddescription...... 69
Figure 5.16 Equation of state (EOS) model tuning of Gas-Oil Ratio (GOR) and oil densitywithPVTlaboratorydata...... 71
Figure 5.17 Equation of State (EOS) model tuning of oil viscosity and swelling factorwithPVTlaboratorydata...... 72
xi Figure5.18 Historymatchingprocesswithoilratescontrol ...... 73
Figure 5.19 History match of bottom hole pressure and gas rates ...... 74
Figure 5.20 Bottom hole pressure and oil rates during EOR scheme1...... 76
Figure5.21 GasratesandcomparisonofallEORschemes...... 78
Figure 5.22 Comparison of two solvent types and e↵ect of molecular di↵usion . . . . 79
Figure 6.1 Flowchart for calculating molecular di↵usion coe cients using Leahy-DiosandFiroozabadi (2007)approach...... 93
xii LIST OF TABLES
Table 3.1 Three-component fluid system used for thermodynamic validation . . . . . 21
Table3.2 Testcasereservoirparameters ...... 25
Table 3.3 Three-component fluid system used for simulation run ...... 27
Table 4.1 Summary of Bilinear solutions for single-phase, multi-phase black oil, and multi-phasecompositionalmodels ...... 39
Table 4.2 Summary of Linear solutions for single-phase, multi-phase black oil, and multi-phasecompositionalmodels ...... 39
Table4.3 Multiphasecasestudyreservoirparameters...... 40
Table 4.4 Three-component fluid system used for multiphase depletion run...... 42
Table4.5 Bilinearmultiphaseflowanalysisfordepletionrun ...... 48
Table4.6 Linearmultiphaseflowanalysisfordepletionrun ...... 48
Table5.1 XRDanalysisofMiddleBakkenCore ...... 55
Table5.2 XRDanalysisofLowerBakkenCore ...... 56
Table5.3 RadialcaseforMiddleBakkencore ...... 59
Table 5.4 NorthDakotaBakkenwellreservoirparameters...... 68
Table 5.5 Lumped-component Middle Bakken fluid system used for field case . . . . . 70
Table 5.6 CO2 Enhanced oil recovery schemes ...... 75
Table 5.7 NGL Enhanced oil recovery schemes ...... 77
Table 5.8 Summary of results for Enhanced oil recovery schemes ...... 80
Table 6.1 Molecular Di↵usion Calculations for Middle Bakken fluid system ...... 94
Table 6.2 Summary of swelling tests for Middle Bakken fluid system ...... 95
xiii LIST OF SYMBOLS
A ...... PR-EOS coe cient a ...... attraction parameter for EOS
B ...... PR-EOS coe cient b ...... repulsive parameter for EOS
psi bblSINGLE [ ]...... y-intercept of bilinear single-phase flow equation (RB/d)cp
psi blSINGLE [ ]...... y-intercept of linear single-phase flow equation (RB/d)cp
psi bblMULTI [ ]...... y-intercept of bilinear multi-phase flow equation (RB/d)cp
psi blMULTI [ ]...... y-intercept of linear multi-phase flow equation (RB/d)cp
RB B [ STB ]...... formation volume factor
RB Bo [ STB ]...... oil formation volume factor
RB Bg [ STB ]...... gas formation volume factor
RB Bw [ STB ]...... water formation volume factor
C1 ...... methane
C4 ...... butane
C7 ...... heptane
C10 ...... decane
1 cv [psi ]...... fluid compressibility
1 c [psi ]...... rock compressibility
1 ct [psi ]...... total compressibility
1 ct,f [psi ]...... total fracture compressibility
xiv 1 ct,m [psi ]...... total matrix compressibility f [psi]...... fugacity
FcD [mdft]...... fracture conductivity h [ft]...... gravity head between fracture and matrix kc [fraction]...... equilibrium ratio of component c kB ...... Boltzmann constant khf [md]...... hydraulic fracture permeability kf,eff [md]...... e↵ective fracture permeability km [md]...... matrix permeability kf [md]...... fracture permeability krg [fraction]...... gas relative permeability kro [fraction]...... oil relative permeability krw [fraction]...... water relative permeability
Lx,Ly,Lz [ft]...... fracture spacing in x,y, and z
psi mblSINGLE [ ]...... slope of bilinear single-phase flow equation (RB/d)cppd
psi mlSINGLE [ ]...... slope of linear single-phase flow equation (RB/d)cppd
psi mblMULTI [ ]...... slope of bilinear multi-phase flow equation (RB/d)cppd
psi mlMULTI [ ]...... slope of linear multi-phase flow equation (RB/d)cppd
1 MoleCorr [ day ]...... mole correction term nf ...... total number of fractures nc ...... total number of components
N [lb mole]...... total number of moles in hydrocarbon vapor phase g
xv N [lb mole]..total number of moles of component c in hydrocarbon vapor phase gc N [lb mole]...... total number of moles in hydrocarbon oil phase o N [lb mole]....total number of moles of component c in hydrocarbon oil phase oc N [lb mole]...... total number of moles in aqueous phase w N [lb mole]...... total number of moles of component c in aqueous phase wc N [lb mole]...... total number of moles t N [lb mole]...... total number of moles of component c tc po [psia]...... oil pressure pcog [psia]...... oil and gas capillary pressure pcow [psia]...... oil and water capillary pressure pc [psia]...... critical pressure pwf [psia]...... flowing bottom hole pressure pwell [psia]...... well flowing bottom pressure (constraint)
ft3 q [ day ]...... flow rate
1 lb mole qˆ [ day ][ fy 3day ]...... reservoir flow rate per rock volume
ft3 qo [ day ]...... oil flow rate
ft3 qg [ day ]...... gas flow rate
ft3 qw [ day ]...... water flow rate
3 ft ⇧psia R [ oRlb mole ]...... gas constant (10.731) shf ...... apparent skin sg ...... gas saturation so ...... oil saturation
xvi sr ...... residual saturation sw ...... water saturation t [day]...... time
T [oR]...... temperature (oF +459.67)
ft2md Tx [ ft ]...... Single phase transmissibility in x-direction
o o Tc [ R]...... critical temperature ( F +459.67)
lb mole Uc [ day ]...... net molar flux of component c
3 VR [ft ]...... rock volume
3 Vt [ft ]...... total system volume
ft3 v [ lb mole ]...... specific volume
ft3 vt [ lb mole ]...... total specific volume
ft3 vtc [ lb mole ]...... total partial molar volume with respect to component c wc [fraction]...... mole fraction of component c in aqueous phase whf [ft]...... hydraulic fracture width (pseudoized)
whforiginal [ft]...... hydraulic fracture width (original) xc [fraction]...... mole fraction of component c in hydrocarbon liquid phase yc [fraction]...... mole fraction of component c in hydrocarbon vapor phase yf [ft]...... fracture half length xf [ft]...... fracture half length z [ratio]...... deviation factor zc [fraction] ...... overall mole fraction of component c in hydrocarbon and aqueous phases
↵ ...... EOS coe cient of attractive term
xvii 1 ↵o [ day ]...... oil product term used in multiphase di↵usivity
1 ↵g [ day ]...... gas product term used in multiphase di↵usivity
1 ↵w [ day ]...... water product term used in multiphase di↵usivity