fp21 = fracture.geomp21(gset<,dfn>) (3D only)

Get the P21 of polygons. The P21 is the cumulative sum of intersection lengths 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.


fp21 - P21 value


gset - geometry set pointer

dfn - optional DFN pointer

Usage Example

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

model new
model domain extent -5 5

;Create the 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_p21
    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('P21 of dfn p1 = '+string(fracture.geomp21(geom,dfn1)))
    io.out('P21 of dfn p2 = '+string(fracture.geomp21(geom,dfn2)))
    io.out('P21 of dfn p3 = '+string(fracture.geomp21(geom,dfn3)))
    io.out('P21 of dfn p4 = '+string(fracture.geomp21(geom,dfn4)))
    io.out('P21 of all dfns = '+string(fracture.geomp21(geom)))


program return