fracture.geomp10

Syntax

fp10 = fracture.geomp10(gset<,dfn>)

Get the P10 of edges. The P10 (i.e., the number of fracture intersections with a transect per unit length) of each edge in the geometry set gset is calculated, and the average over all edges is returned. If dfn is not given the calculation is undertaken over all fractures of all DFNs.

Returns:

fp10 - average P10 of geometry set edges

Arguments:

gset - geometry set pointer

dfn - optional DFN pointer

Usage Example

The following example illustrates how fracture.geomp10 can be used to get the P10 along a borehole for various DFNs, or for all fractures (3D).

model new
model random 10000
model domain extent -5 5

;Create the borehole geometry
geometry set 'borehole'
geometry edge create by-position (-5,0,0) (5,0,0)

;Generate fractures using the default template
fracture generate dfn 'p1' p10 0.1 geometry 'borehole'
fracture generate dfn 'p2' p10 0.2 geometry 'borehole'
fracture generate dfn 'p3' p10 0.3 geometry 'borehole'
fracture generate dfn 'p4' p10 0.4 geometry 'borehole'

; fish access to densities
fish define access_p10    
  local geom = geom.set.find('borehole')    
  local dfn1 = dfn.find('p1')    
  local dfn2 = dfn.find('p2')    
  local dfn3 = dfn.find('p3')    
  local dfn4 = dfn.find('p4')    
  io.out('P10 of dfn p1 = '+string(fracture.geomp10(geom,dfn1)))    
  io.out('P10 of dfn p2 = '+string(fracture.geomp10(geom,dfn2)))    
  io.out('P10 of dfn p3 = '+string(fracture.geomp10(geom,dfn3)))    
  io.out('P10 of dfn p4 = '+string(fracture.geomp10(geom,dfn4)))    
  io.out('P10 of all dfns = '+string(fracture.geomp10(geom)))
end
[access_p10]

program return

The same example in 2D:

model new

model domain extent -5 5

; create geometry
geometry set 'borehole'
geometry edge create by-position (-5,0) (5,0)

; generate dfns
fracture template create 'default'
fracture generate dfn 'p1' p10 0.1 geometry 'borehole'
fracture generate dfn 'p2' p10 0.2 geometry 'borehole'
fracture generate dfn 'p3' p10 0.3 geometry 'borehole'
fracture generate dfn 'p4' p10 0.4 geometry 'borehole'

; fish access to densities
fish define access_p10    
    geom = geom.set.find('borehole')    
    dfn1 = dfn.find('p1')    
    dfn2 = dfn.find('p2')    
    dfn3 = dfn.find('p3')    
    dfn4 = dfn.find('p4')    
    oo = io.out('P10 of dfn p1 = '+string(fracture.geomp10(geom,dfn1)))    
    oo = io.out('P10 of dfn p2 = '+string(fracture.geomp10(geom,dfn2)))    
    oo = io.out('P10 of dfn p3 = '+string(fracture.geomp10(geom,dfn3)))    
    oo = io.out('P10 of dfn p4 = '+string(fracture.geomp10(geom,dfn4)))    
    oo = io.out('P10 of all dfns = '+string(fracture.geomp10(geom)))
end
[access_p10]

program return