Difference between revisions of "Category: OptiFrac"

From wiki.pengtools.com
Jump to: navigation, search
(→‎Fracturing Optimization Software)
 
(71 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
__TOC__
 
__TOC__
−
== Brief ==
+
== Fracturing Optimization Software ==
  
 
[[File:OptiFrac_i.png|thumb|right|200px|link=https://www.pengtools.com/optiFrac|pengtools optiFrac]]
 
[[File:OptiFrac_i.png|thumb|right|200px|link=https://www.pengtools.com/optiFrac|pengtools optiFrac]]
  
−
[[:Category:optiFrac | optiFrac]] is a single well fracture design optimization tool.  
+
[[:Category:optiFrac | optiFrac]] is a single well fracture design optimization software.  
  
 
For the given set of reservoir and proppant properties [[:Category:optiFrac | optiFrac]] calculates maximum achievable well productivity index ([[JD]]) and required fracture geometry.
 
For the given set of reservoir and proppant properties [[:Category:optiFrac | optiFrac]] calculates maximum achievable well productivity index ([[JD]]) and required fracture geometry.
Line 10: Line 10:
 
In production optimization of hydraulic fracturing, it is important to recognize that for a fixed proppant volume, fracture length and fracture width compete to each other. Consequently, there is an optimum fracture length and width that would maximize the dimensionless productivity index of the fractured well <ref name = pushing/>. '''Unified Fracture Design''' provides a methodology whereby for any given proppant volume the best combination of width and length could be determined <ref name=UFD2002/>.  
 
In production optimization of hydraulic fracturing, it is important to recognize that for a fixed proppant volume, fracture length and fracture width compete to each other. Consequently, there is an optimum fracture length and width that would maximize the dimensionless productivity index of the fractured well <ref name = pushing/>. '''Unified Fracture Design''' provides a methodology whereby for any given proppant volume the best combination of width and length could be determined <ref name=UFD2002/>.  
  
−
[[:Category:optiFrac | optiFrac]] is available online at [https://www.pengtools.com www.pengtools.com] and [https://itunes.apple.com/us/app/optifrac/id1064023255 in AppStore] for iPad.
+
[[:Category:optiFrac | optiFrac]] is available online at [https://www.pengtools.com www.pengtools.com].
  
 
== Typical applications ==
 
== Typical applications ==
  
−
* Optimization of hydraulic fracturing with '''Unified Fracture Design'''<ref name=UFD2002/>
+
* Hydraulic fracture design in oil and gas reservoirs
−
* Understanding post-fracturing production performance
+
* Optimization of hydraulic fracturing with '''Unified Fracture Design'''<ref name=UFD2002/> method
−
* Design Sensitivity Studies and Benchmarking
+
* Calculating '''optimum''' fracture design parameters:
−
* Optimum fracture design parameters determination:
 
 
** Dimensionless productivity index, '''J<sub>D</sub>''' .
 
** Dimensionless productivity index, '''J<sub>D</sub>''' .
−
** Dimensionless Fracture conductivity, '''C<sub>fD</sub>''' .
+
** Dimensionless fracture conductivity, '''C<sub>fD</sub>''' .
−
** Fracture half length, '''X<sub>f</sub>''' .
+
** Fracture half length, '''x<sub>f</sub>''' .
−
** Fracture width, '''''w''''' .
+
** Fracture width, '''''w<sub>f</sub>''''' .
 
** Fracture penetration, '''I<sub>x</sub>''' .
 
** Fracture penetration, '''I<sub>x</sub>''' .
−
 
+
* Understanding post-fracturing production performance
−
== Main features ==
+
* Sensitivity studies
−
* Plot of '''J<sub>D</sub>''' as a function of '''C<sub>fD</sub>''' and '''I<sub>x</sub>''' as parameter.
 
−
* Plot of '''J<sub>D</sub>''' as a function of  '''C<sub>fD</sub>''' and  '''N<sub>p</sub>''' as parameter.
 
−
* Design optimization curve which corresponds to the maximum '''J<sub>D</sub>''' values for different '''N<sub>p</sub>'''.
 
−
* Design Optimum Point at which the dimensionless productivity index, '''J<sub>D</sub>''', is maximized for the given proppant, fracture and reservoir parameters.
 
−
* Physical constraints envelope.
 
−
* Proppant library with predefined proppant properties.
 
−
* Users Data Worksheet for benchmarking vs actual.
 
−
 
 
−
== Interface features ==
 
−
 
 
−
* Save and share models with colleagues
 
−
* Last saved model on current computer and browser is automatically opened
 
−
* Metric and US oilfield units
 
−
* Save as image and print plots by means of chart context menu (button at the upper-right corner of chart)
 
−
* Download '''pdf''' report with input parameters, calculated values and plots
 
−
* Select and copy results to Excel or other applications
 
  
 
==Math & Physics==
 
==Math & Physics==
Line 50: Line 33:
 
:<math>C_{f_D}=\frac{k_f w_f}{k x_f} </math> - dimensionless fracture conductivity <ref name=UFD2002/>,  
 
:<math>C_{f_D}=\frac{k_f w_f}{k x_f} </math> - dimensionless fracture conductivity <ref name=UFD2002/>,  
  
−
:<math>N_p={I_x}^2 C_{f_D} = \frac{4 k_f x_f w_f}{k {x_e}^2} = \frac{4 k_f x_f w_f h}{k {x_e}^2 h} =  \frac{2 k_f V_f}{k V_r} </math> - proppant numer <ref name=UFD2002/>,
+
:<math>N_p={I_x}^2 C_{f_D} = \frac{4 k_f x_f w_f}{k {x_e}^2} = \frac{4 k_f x_f w_f h}{k {x_e}^2 h} =  \frac{2 k_f V_f}{k V_r} </math> - proppant number <ref name=UFD2002/>,
  
−
:<math>{\bar{P}_D}^{pss} = \frac{1}{2} ln{\left (\frac{16}{\gamma C_A {I_x}^2} \right )} + f </math> - pseudo-steady state equation for finite-conductivity fractured wells <ref name = pushing/>,
+
:<math>{\bar{P}_D}^{pss} = \frac{1}{2} ln{\left (\frac{16}{1.78 C_A {I_x}^2} \right )} + f </math> - pseudo-steady state equation for finite-conductivity fractured wells <ref name = pushing/>,
  
−
:<math>{\bar{P}_D}^{ss} = \frac{1}{2} ln{\left (\frac{26.4}{\gamma C_A {I_x}^2} \right )} + f </math> - steady state equation for finite-conductivity fractured wells <ref name = pushing/>,
+
:<math>{\bar{P}_D}^{ss} = \frac{1}{2} ln{\left (\frac{26.4}{1.78 C_A {I_x}^2} \right )} + f </math> - steady state equation for finite-conductivity fractured wells <ref name = pushing/>,
  
 
:<math>C_{A}=7.7327 {I_x}^3 - 25.204 {I_x}^2 - 0.7211 I_x + 30.555 </math> - shape factor <ref name = pushing/>.
 
:<math>C_{A}=7.7327 {I_x}^3 - 25.204 {I_x}^2 - 0.7211 I_x + 30.555 </math> - shape factor <ref name = pushing/>.
  
−
==Type Curves==
+
==Fracturing Design Optimization Type Curves==
  
−
The Type Curves shows the dimensionless productivity index, '''J<sub>D</sub>''',  at steady and pseudo-steady state as function of  '''C<sub>fD</sub>''', using '''I<sub>x</sub>''' as parameter (red curves) and overlapping with the type curve with '''N<sub>p</sub>''' as parameter (black curves).
+
The Type Curves show the dimensionless productivity index, '''J<sub>D</sub>''',  at steady and pseudo-steady state as a function of  '''C<sub>fD</sub>''', using '''I<sub>x</sub>''' as parameter (red curves) and overlapping with the type curve with '''N<sub>p</sub>''' as parameter (black curves).
  
 
The green curve along the maximum points for different '''N<sub>p</sub>''' values is “Design Optimization Curve”<ref name = pushing/>. This curve represents the target of the designs of the fracture treatments in a dimensionless form.
 
The green curve along the maximum points for different '''N<sub>p</sub>''' values is “Design Optimization Curve”<ref name = pushing/>. This curve represents the target of the designs of the fracture treatments in a dimensionless form.
Line 70: Line 53:
 
Type Curves were obtained, through seven hundred runs with a numerical simulator, modeling a fractured well in a closed square reservoir <ref name = pushing/>. For infinite and finite fracture conductivities, the shape factors, '''C<sub>A</sub>''', can be calculated if the '''P<sub>D</sub>''' is known for a specific value of '''I<sub>x</sub>'''. The '''P<sub>D</sub>''' value obtained by numerical simulations. After knowing '''C<sub>A</sub>''' (which would be a function of  '''I<sub>x</sub>'''), f-function values can be calculated if  '''P<sub>D</sub>''' is known for a specific dimensionless fracture conductivity '''C<sub>fD</sub>''' and fracture penetration  '''I<sub>x</sub>'''.
 
Type Curves were obtained, through seven hundred runs with a numerical simulator, modeling a fractured well in a closed square reservoir <ref name = pushing/>. For infinite and finite fracture conductivities, the shape factors, '''C<sub>A</sub>''', can be calculated if the '''P<sub>D</sub>''' is known for a specific value of '''I<sub>x</sub>'''. The '''P<sub>D</sub>''' value obtained by numerical simulations. After knowing '''C<sub>A</sub>''' (which would be a function of  '''I<sub>x</sub>'''), f-function values can be calculated if  '''P<sub>D</sub>''' is known for a specific dimensionless fracture conductivity '''C<sub>fD</sub>''' and fracture penetration  '''I<sub>x</sub>'''.
  
−
==Flow Diagram==
+
==Fracturing Optimization Flow Diagram==
  
 
[[File:optiFrac flow diagram.png|400px]]
 
[[File:optiFrac flow diagram.png|400px]]
  
−
==Workflow==
+
==Fracturing Optimization Workflow==
 +
 
 +
[[File:Fracturedimensions.png|thumb|right|500px| Fracture Plane Dimensions]]
  
 
1. Calculate the '''N<sub>p</sub>''':
 
1. Calculate the '''N<sub>p</sub>''':
Line 82: Line 67:
 
:<math>k_f=k_{prop} * Gel Damage</math> the fracture permeability
 
:<math>k_f=k_{prop} * Gel Damage</math> the fracture permeability
  
−
:<math>M_f=\frac{M_{prop}}{Out Zone Growth} \frac{h_{net}}{h_{gross}}</math> the proppant mass in the pay zone
+
:<math>M_f=M_{prop} \frac{h_{net}}{h_f}</math> the proppant mass in the pay zone
  
 
:<math>V_f=\frac{M_f}{SG_{prop} (1 - \phi_{prop})}</math> the fracture volume in the pay zone
 
:<math>V_f=\frac{M_f}{SG_{prop} (1 - \phi_{prop})}</math> the fracture volume in the pay zone
Line 96: Line 81:
 
:<math>{w_f}^{opt}=\frac{{C_{fD}}^{opt} {x_f}^{opt} k}{k_f}</math>
 
:<math>{w_f}^{opt}=\frac{{C_{fD}}^{opt} {x_f}^{opt} k}{k_f}</math>
  
−
:<math>h_f=h_{gross} * Out Zone Growth</math> the fracture height
+
==Fracturing Physical Constraints==
 +
It is important to mention that the Design Optimization Curve could give unrealistic fracture geometry depending on the reservoir permeability, reservoir mechanical properties and target '''N<sub>p</sub>'''. The two most common scenarios are <ref name = pushing/>:
 +
 
 +
# The required net pressure for the fracture geometry is too high - “maximum net pressure curve”,
 +
# Fracture width is too small (fracture too narrow) - “minimum width curve”.
  
−
==Physical Constraints==
 
−
It is important to mention that the Design Optimization Curve could give unrealistic fracture geometry depending on the reservoir permeability, reservoir mechanical properties and target '''N<sub>p</sub>'''. The two most common scenarios are <ref name = pushing/>:
 
  
−
# The required net pressure for the fracture geometry is too high,
+
The area between the “minimum width curve” and the “maximum net pressure curve” is the “working area” (highlighted in yellow) of the whole type curve for the specific rock mechanical properties, reservoir and proppant properties used. Any fracture design for this specific case should be located on the “optimum design curve” anywhere in this working area depending on the desired '''N<sub>p</sub>'''<ref name = pushing/>.
−
# Fracture width is too small (fracture too narrow).
+
 
 +
[[File:Physical Constraints on a PSS type curve.png  | link=https://www.pengtools.com/optiFrac | Open in optiFrac]]
  
 
====Maximum net pressure====
 
====Maximum net pressure====
Line 114: Line 102:
 
Fracture propped width should be greater than N times mean proppant diameter (to provide at least N proppant layers in the fracture after closure)<ref name = pushing/>. N=3 in the [[:Category:optiFrac | optiFrac]].
 
Fracture propped width should be greater than N times mean proppant diameter (to provide at least N proppant layers in the fracture after closure)<ref name = pushing/>. N=3 in the [[:Category:optiFrac | optiFrac]].
  
−
==Discussion==
+
== Main features ==  
 +
* Fracture design Type Curves (Plot of '''J<sub>D</sub>''' as a function of '''C<sub>fD</sub>''' using '''I<sub>x</sub>''' and '''N<sub>p</sub>''' as parameter).
 +
* Design Optimization Curve which corresponds to the maximum '''J<sub>D</sub>''' values for different '''N<sub>p</sub>'''.
 +
* Design Optimum Point at which '''J<sub>D</sub>''' is maximized for the given proppant, fracture and reservoir parameters.
 +
* Physical constraints envelope.
 +
* [[Hydraulic fracturing proppant catalog]] with predefined proppant properties.
 +
* Users Data Worksheet for benchmarking vs actual.
 +
* "Default values" button resets input values to the default values
 +
* Switch between Metric and Field units
 +
* Save/load models to the files and to the user’s cloud
 +
* Share models to the public cloud or by using model’s link
 +
* Export pdf report containing input parameters, calculated values and plots
 +
* Continue your work from where you stopped: last saved model will be automatically opened
 +
* Download the chart as an image or data and print (upper-right corner chart’s button)
 +
* Export results table to Excel or other application
  
 
== Nomenclature ==
 
== Nomenclature ==
−
:<math>J_D</math> =  
+
:<math>C_{fD}</math> = dimensionless fracture conductivity, dimensionless
−
:<math>C_{fD}</math> =  
+
:<math>C_A</math> = shape factor, dimensionless
−
:<math>C_{fD}</math> =  
+
:<math>E</math> = Young's Modulus, psia
−
:<math>x_f</math> =  
+
:<math>f</math> = f-function, dimensionless
 +
:<math>Gel Damage</math> = proppant permeability reduction due to gel damage, %
 +
:<math>h</math> = reservoir thickness, ft
 +
:<math>I_x</math> = penetration ratio, dimensionless
 +
:<math>J_D</math> = dimensionless productivity index, dimensionless
 +
:<math>k</math> = permeability, md
 +
:<math>M</math> = mass, lbm
 +
:<math>N_p</math> = dimensionless proppant number, dimensionless
 +
:<math>\bar{P}_D</math> = dimensionless pressure (based on average pressure), dimensionless
 +
:<math>P_{net}</math> = net pressure, psia
 +
:<math>SG</math> = specific gravity, dimensionless
 +
:<math>V</math> = volume, ft<sup>3</sup>
 +
:<math>w</math> = width, ft
 +
:<math>x_e</math> = drainage width, ft
 +
:<math>x_f</math> = fracture half-length, ft
 +
 
 +
===Greek symbols===
 +
:<math>\delta</math> = dry to wet width ratio at the end of pumping, usually 0.5-0.7
 +
:<math>\gamma</math> = geometric factor in vertical direction, 0.75 for PKN model, 1 for KGD model
 +
:<math>\nu</math> = Poisson's ratio, dimensionless
 +
:<math>\phi</math> = porosity, fraction
 +
:<math>\pi</math> = 3.1415
 +
 
 +
===Superscripts===
 +
:opt = optimal
 +
:pss = pseudo-steady state
 +
:ss = steady state
 +
 
 +
===Subscripts===
 +
:e = external
 +
:f = fracture
 +
:gross = gross
 +
:max = maximum
 +
:net = net
 +
:prop = proppant
 +
:r = reservoir
  
 
== References ==
 
== References ==
Line 149: Line 186:
 
</references>
 
</references>
 
[[Category:pengtools]]
 
[[Category:pengtools]]
 +
 +
{{#seo:
 +
|title=Hydraulic Fracture Design Optimization Software
 +
|titlemode= replace
 +
|keywords=hydraulic fracturing, hydraulic fracturing formulas, hydraulic fracturing proppant, optimization, petroleum engineering
 +
|description=Hydraulic fracturing formulas for fracture design optimization in optiFrac software at pengtools.com
 +
}}

Latest revision as of 08:20, 16 March 2020

Fracturing Optimization Software

pengtools optiFrac

optiFrac is a single well fracture design optimization software.

For the given set of reservoir and proppant properties optiFrac calculates maximum achievable well productivity index (JD) and required fracture geometry.

In production optimization of hydraulic fracturing, it is important to recognize that for a fixed proppant volume, fracture length and fracture width compete to each other. Consequently, there is an optimum fracture length and width that would maximize the dimensionless productivity index of the fractured well [1]. Unified Fracture Design provides a methodology whereby for any given proppant volume the best combination of width and length could be determined [2].

optiFrac is available online at www.pengtools.com.

Typical applications

  • Hydraulic fracture design in oil and gas reservoirs
  • Optimization of hydraulic fracturing with Unified Fracture Design[2] method
  • Calculating optimum fracture design parameters:
    • Dimensionless productivity index, JD .
    • Dimensionless fracture conductivity, CfD .
    • Fracture half length, xf .
    • Fracture width, wf .
    • Fracture penetration, Ix .
  • Understanding post-fracturing production performance
  • Sensitivity studies

Math & Physics

Fracture Notation
I_x=\frac{2x_f}{x_e} - penetration ratio [2],
C_{f_D}=\frac{k_f w_f}{k x_f} - dimensionless fracture conductivity [2],
N_p={I_x}^2 C_{f_D} = \frac{4 k_f x_f w_f}{k {x_e}^2} = \frac{4 k_f x_f w_f h}{k {x_e}^2 h} = \frac{2 k_f V_f}{k V_r} - proppant number [2],
{\bar{P}_D}^{pss} = \frac{1}{2} ln{\left (\frac{16}{1.78 C_A {I_x}^2} \right )} + f - pseudo-steady state equation for finite-conductivity fractured wells [1],
{\bar{P}_D}^{ss} = \frac{1}{2} ln{\left (\frac{26.4}{1.78 C_A {I_x}^2} \right )} + f - steady state equation for finite-conductivity fractured wells [1],
C_{A}=7.7327 {I_x}^3 - 25.204 {I_x}^2 - 0.7211 I_x + 30.555 - shape factor [1].

Fracturing Design Optimization Type Curves

The Type Curves show the dimensionless productivity index, JD, at steady and pseudo-steady state as a function of CfD, using Ix as parameter (red curves) and overlapping with the type curve with Np as parameter (black curves).

The green curve along the maximum points for different Np values is “Design Optimization Curve”[1]. This curve represents the target of the designs of the fracture treatments in a dimensionless form.

Open in optiFrac

Open in optiFrac

Type Curves were obtained, through seven hundred runs with a numerical simulator, modeling a fractured well in a closed square reservoir [1]. For infinite and finite fracture conductivities, the shape factors, CA, can be calculated if the PD is known for a specific value of Ix. The PD value obtained by numerical simulations. After knowing CA (which would be a function of Ix), f-function values can be calculated if PD is known for a specific dimensionless fracture conductivity CfD and fracture penetration Ix.

Fracturing Optimization Flow Diagram

OptiFrac flow diagram.png

Fracturing Optimization Workflow

Fracture Plane Dimensions

1. Calculate the Np:

V_r=h_{net} {x_e}^2 the volume of the reservoir
k_f=k_{prop} * Gel Damage the fracture permeability
M_f=M_{prop} \frac{h_{net}}{h_f} the proppant mass in the pay zone
V_f=\frac{M_f}{SG_{prop} (1 - \phi_{prop})} the fracture volume in the pay zone
N_p=\frac{2 k_f V_f}{k V_r} the proppant number

2. Read CfDopt, Ixopt, JDopt from the Design Optimization Curve of the Type Curve

3. Calculate optimum fracture half-length and width:

{x_f}^{opt}=0.5 x_e {I_x}^{opt}
{w_f}^{opt}=\frac{{C_{fD}}^{opt} {x_f}^{opt} k}{k_f}

Fracturing Physical Constraints

It is important to mention that the Design Optimization Curve could give unrealistic fracture geometry depending on the reservoir permeability, reservoir mechanical properties and target Np. The two most common scenarios are [1]:

  1. The required net pressure for the fracture geometry is too high - “maximum net pressure curve”,
  2. Fracture width is too small (fracture too narrow) - “minimum width curve”.


The area between the “minimum width curve” and the “maximum net pressure curve” is the “working area” (highlighted in yellow) of the whole type curve for the specific rock mechanical properties, reservoir and proppant properties used. Any fracture design for this specific case should be located on the “optimum design curve” anywhere in this working area depending on the desired Np[1].

Open in optiFrac

Maximum net pressure

The maximum net pressure during the fracturing treatment should provide a surface pressure less than a certain value (which is surface pressure operational limit) [1].

w_{max}=\frac{2 P_{net} h_f (1 - \nu^2)}{E}
w=w_{max} \frac{\pi}{4} \gamma \delta

Minimum fracture width

Fracture propped width should be greater than N times mean proppant diameter (to provide at least N proppant layers in the fracture after closure)[1]. N=3 in the optiFrac.

Main features

  • Fracture design Type Curves (Plot of JD as a function of CfD using Ix and Np as parameter).
  • Design Optimization Curve which corresponds to the maximum JD values for different Np.
  • Design Optimum Point at which JD is maximized for the given proppant, fracture and reservoir parameters.
  • Physical constraints envelope.
  • Hydraulic fracturing proppant catalog with predefined proppant properties.
  • Users Data Worksheet for benchmarking vs actual.
  • "Default values" button resets input values to the default values
  • Switch between Metric and Field units
  • Save/load models to the files and to the user’s cloud
  • Share models to the public cloud or by using model’s link
  • Export pdf report containing input parameters, calculated values and plots
  • Continue your work from where you stopped: last saved model will be automatically opened
  • Download the chart as an image or data and print (upper-right corner chart’s button)
  • Export results table to Excel or other application

Nomenclature

C_{fD} = dimensionless fracture conductivity, dimensionless
C_A = shape factor, dimensionless
E = Young's Modulus, psia
f = f-function, dimensionless
Gel Damage = proppant permeability reduction due to gel damage, %
h = reservoir thickness, ft
I_x = penetration ratio, dimensionless
J_D = dimensionless productivity index, dimensionless
k = permeability, md
M = mass, lbm
N_p = dimensionless proppant number, dimensionless
\bar{P}_D = dimensionless pressure (based on average pressure), dimensionless
P_{net} = net pressure, psia
SG = specific gravity, dimensionless
V = volume, ft3
w = width, ft
x_e = drainage width, ft
x_f = fracture half-length, ft

Greek symbols

\delta = dry to wet width ratio at the end of pumping, usually 0.5-0.7
\gamma = geometric factor in vertical direction, 0.75 for PKN model, 1 for KGD model
\nu = Poisson's ratio, dimensionless
\phi = porosity, fraction
\pi = 3.1415

Superscripts

opt = optimal
pss = pseudo-steady state
ss = steady state

Subscripts

e = external
f = fracture
gross = gross
max = maximum
net = net
prop = proppant
r = reservoir

References

  1. ↑ 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 Rueda, J.I.; Mach, J.; Wolcott, D. (2004). "Pushing Fracturing Limits to Maximize Producibility in Turbidite Formations in Russia"Free registration required (SPE-91760-MS). Society of Petroleum Engineers. 
  2. ↑ 2.0 2.1 2.2 2.3 2.4 Economides, Michael J.; Oligney, Ronald; Valko, Peter (2002). Unified Fracture Design: Bridging the Gap Between Theory and Practice. Alvin, Texas: Orsa Press. 

Pages in category "OptiFrac"

The following 6 pages are in this category, out of 6 total.

4

6

H

J

O

Retrieved from "https://wiki.pengtools.com/index.php?title=Category:OptiFrac&oldid=6626"