def test_get_faciesfraction_some_wells(testpath, tmp_path):
"""Import some wells and get the facies fractions per zone, for
wells < 70 degrees inclination.
"""
wlist = [
Well(wpath, zonelogname="Zonelog")
for wpath in glob.glob(str(testpath / WFILES2))
]
mypoints = Points()
facname = "Facies"
fcode = [1]
mypoints.dfrac_from_wells(wlist,
facname,
fcode,
zonelist=None,
incl_limit=70)
# rename column
mypoints.zname = "FACFRAC"
myquery = 'WELLNAME == "OP_1" and ZONE == 1'
usedf = mypoints.dataframe.query(myquery)
assert abs(usedf[mypoints.zname].values[0] - 0.86957) < 0.001
mypoints.to_file(
tmp_path / "ffrac_per_zone.rmsasc",
fformat="rms_attr",
attributes=["WELLNAME", "ZONE"],
pfilter={"ZONE": [1]},
)
def test_get_faciesfraction_some_wells():
"""Import some wells and get the facies fractions per zone, for
wells < 70 degrees inclination.
"""
wlist = []
for w in sorted(glob.glob(str(WFILES2))):
wlist.append(Well(w, zonelogname="Zonelog"))
logger.info("Imported well {}".format(w))
mypoints = Points()
facname = "Facies"
fcode = [1]
nwell = mypoints.dfrac_from_wells(
wlist, facname, fcode, zonelist=None, incl_limit=70
)
# rename column
mypoints.zname = "FACFRAC"
logger.info("Number of wells is %s, DATAFRAME:\n, %s", nwell, mypoints.dataframe)
myquery = 'WELLNAME == "OP_1" and ZONE == 1'
usedf = mypoints.dataframe.query(myquery)
assert abs(usedf[mypoints.zname].values[0] - 0.86957) < 0.001
mypoints.to_file(
"TMP/ffrac_per_zone.rmsasc",
fformat="rms_attr",
attributes=["WELLNAME", "ZONE"],
pfilter={"ZONE": [1]},
)
def test_get_faciesfraction_some_wells():
"""Import some wells and get the facies fractions per zone, for
wells < 70 degrees inclination.
"""
wlist = []
for w in sorted(glob.glob(wfiles2)):
wlist.append(Well(w, zonelogname='Zonelog'))
logger.info('Imported well {}'.format(w))
mypoints = Points()
facname = 'Facies'
fcode = [1]
nwell = mypoints.dfrac_from_wells(wlist,
facname,
fcode,
zonelist=None,
incl_limit=70)
# rename column
mypoints.zname = 'FACFRAC'
logger.info('Number of wells is %s, DATAFRAME:\n, %s', nwell,
mypoints.dataframe)
myquery = 'WELLNAME == "OP_1" and ZONE == 1'
usedf = mypoints.dataframe.query(myquery)
assert abs(usedf[mypoints.zname].values[0] - 0.86957) < 0.001
mypoints.to_file('TMP/ffrac_per_zone.rmsasc',
fformat='rms_attr',
attributes=['WELLNAME', 'ZONE'],
pfilter={'ZONE': [1]})
