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.
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-positions (-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-positions (-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
Was this helpful? ... | Itasca Software © 2024, Itasca | Updated: Oct 31, 2024 |