Triplecombo+Zonation+Facies

This demonstrates how to plot a three-combo log alongside flui-lithofacies log and user-defined well tops

In [1]:
# Importing the load_table function from the petrolib.file_reader module
from petrolib.file_reader import load_table

# Loading data from the CSV file 'wellA.csv' using the load_table function from the petrolib.file_reader module,
# specifying the delimiter as ','
well_df = load_table('wellA.csv', delimiter=',')

# Displaying the DataFrame containing the loaded data
well_df
Out[1]:
DEPTH GR RT RHOB NPHI DT FLUID-LITHOFACIES
0 5500.0 102.740799 2.5668 2.3228 0.450000 127.483002 1
1 5500.5 100.633003 2.5882 2.3104 0.450000 133.423004 1
2 5501.0 101.032898 2.6036 2.3056 0.450000 133.145996 1
3 5501.5 101.927498 2.6520 2.3144 0.450000 131.022995 1
4 5502.0 101.630699 2.6885 2.3279 0.450000 129.104004 1
... ... ... ... ... ... ... ...
6996 8998.0 91.238297 3.5703 2.4271 0.353446 91.988098 2
6997 8998.5 87.492203 4.6154 2.4087 0.326226 91.277496 2
6998 8999.0 80.359398 4.9716 2.3830 0.315407 90.565697 2
6999 8999.5 75.007797 5.6926 2.3516 0.299560 89.148399 2
7000 9000.0 77.320297 6.0549 2.3249 0.315323 92.421501 2

7001 rows × 7 columns

In [4]:
# Importing the plotFaciesZoneCombo function from the petrolib.plots module
from petrolib.plots import plotFaciesZoneCombo

# Generating a facies zone combo plot using the plotFaciesZoneCombo function, with parameters specifying
# DataFrame (well_df), depth ('DEPTH') and gamma ray ('GR') curves, resistivity ('RT'), neutron porosity ('NPHI'),
# density ('RHOB'), and fluid-lithofacies ('FLUID-LITHOFACIES') columns, depth range (5500 to 9000),
# zone tops (5561.38, 7116.47, 8780.19), zone bottoms (5761.46, 7351.19, 8923.58), zone names ('SAND A', 'SAND B', 'SAND C'),
# facies colors ('#FFC300', '#C80C23', '#120CC8', '#0CC80F'), optional limit, facies labels ('Oil Sand', 'Shale/Shaly-Sand', 'Brine Sand', 'Gas Sand'),
# reservoir threshold (20), fill method ('both'), flag name ('Fluid-Lithofacies'), figure size (13, 30), and title ('WELL A')
plotFaciesZoneCombo(well_df, 'DEPTH', 'GR', 'RT', 'NPHI', 'RHOB', 'FLUID-LITHOFACIES', 5500, 9000,
                     ztops=[5561.38, 7116.47, 8780.19], zbots=[5761.46, 7351.19, 8923.58],
                     zonename=['SAND A', 'SAND B', 'SAND C'],
                     facies_colors=['#FFC300', '#C80C23', '#120CC8', '#0CC80F'],
                     limit=None, facies_labels=['Oil Sand', 'Shale/Shaly-Sand', 'Brine Sand', 'Gas Sand'],
                     res_thres=20, fill='both', flag_name='Fluid-Lithofacies', figsize=(13, 30), title='WELL A')