equinor / flownet

FlowNet - Data-Driven Reservoir Predictions
GNU General Public License v3.0
63 stars 29 forks source link

Support I,J,K-based definition of fault segment orientation #426

Open olwijn opened 3 years ago

olwijn commented 3 years ago

Item 8 of each record under the Eclipse keyword FAULTS indicates the orientation of the face of a fault segment. This should be one of X, Y, Z or I, J, K. FlowNet currently only supports X, Y, Z. The I, J, K orientation is used in the current field case. The script _from_flow.py should be modifed as follows:

if row["FACE"] == "X" or row["FACE"] == "X+" or row["FACE"] == "I" or row["FACE"] == "I+": points.append((row["NAME"], i + 1, j, k)) elif row["FACE"] == "Y" or row["FACE"] == "Y+" or row["FACE"] == "J" or row["FACE"] == "J+": points.append((row["NAME"], i, j + 1, k)) elif row["FACE"] == "Z" or row["FACE"] == "Z+" or row["FACE"] == "K" or row["FACE"] == "K+": points.append((row["NAME"], i, j, k + 1)) elif row["FACE"] == "X-" or row["FACE"] == "I-": points.append((row["NAME"], i - 1, j, k)) elif row["FACE"] == "Y-" or row["FACE"] == "J-": points.append((row["NAME"], i, j - 1, k)) elif row["FACE"] == "Z-" or row["FACE"] == "K-": points.append((row["NAME"], i, j, k - 1)) else: raise ValueError( f"Could not interpret '{row['FACE']}' while reading the FAULTS keyword." )