fracture.geomp20

Syntax

fp20 = fracture.geomp20(gset<,dfn>) (3D only)

Get the P20 of polygons. The P20 is the number of fracture intersections with the geometry polygons per unit area of the geometry polygons. If dfn is not given the calculation is undertaken over all fractures of all DFNs.

Returns:

fp20 - average P20 of geometry set edges

Arguments:

gset - geometry set pointer

dfn - optional DFN pointer

Usage Example

The following example illustrates how fracture.geomp20 can be used to get the P20 on an outcrop for various DFNs, or for all fractures.

model new
model random 10000
model domain extent -5 5

;Create outcrop geometry
geometry set 'outcrop'
geometry polygon create by-position (-5,-5,0) (-5,5,0) (5,5,0) (5,-5,0)

;Generate fractures using the default template
fracture generate dfn 'p1' fracture-count 100 
fracture generate dfn 'p2' fracture-count 200
fracture generate dfn 'p3' fracture-count 300 
fracture generate dfn 'p4' fracture-count 400

; fish access to densities
fish define access_p20
    local geom = geom.set.find('outcrop')
    local dfn1 = dfn.find('p1')
    local dfn2 = dfn.find('p2')
    local dfn3 = dfn.find('p3')
    local dfn4 = dfn.find('p4')
    io.out('P20 of dfn p1 = '+string(fracture.geomp20(geom,dfn1)))
    io.out('P20 of dfn p2 = '+string(fracture.geomp20(geom,dfn2)))
    io.out('P20 of dfn p3 = '+string(fracture.geomp20(geom,dfn3)))
    io.out('P20 of dfn p4 = '+string(fracture.geomp20(geom,dfn4)))
    io.out('P20 of all dfns = '+string(fracture.geomp20(geom)))
end

[access_p20]

program return