allow manual vsh_gr in hc correction and rename notebooks

Author: imranfadhil

.gitignore

@@ -111,4 +111,4 @@ notebooks/data/02_intermediate/*
 notebooks/data/03_processed/*
 notebooks/data/04_project/*
 notebooks/local/*
-notebooks/03_qpp_ssc_single_*
+notebooks/03_qpp__*

notebooks/03b_qpp_carb_batch.ipynb notebooks/03_qpp_carb_batch.ipynb

@@ -105,27 +105,27 @@
     "        SWT=1\n",
     "    ),\n",
     "}\n",
+    "vsh_gr = estimate_vsh_gr(well_data['GR'], min_gr=0, max_gr=50)\n",
     "carb_model = Carbonate(**args)\n",
     "vclw, vcalc, vdolo = carb_model.estimate_lithology(\n",
-    "    well_data['NPHI'], well_data['RHOB'], gr=well_data['GR'], pef=well_data['PEF'],\n",
+    "    well_data['NPHI'], well_data['RHOB'], gr=well_data['GR'], vsh_gr=vsh_gr,  # pef=well_data['PEF'],\n",
     "    model=model, method=method, normalize= True, carbonate_type=carbonate_type\n",
     ")\n",
     "args.update(carb_model.__dict__)\n",
     "well.update_config(args)  # Save lithology model to well\n",
     "\n",
     "# Implement hydrocarbon correction\n",
     "nphihc, rhobhc, hc_flag = neu_den_xplot_hc_correction(\n",
-    "    well_data['NPHI'], well_data['RHOB'], gr=well_data['GR'],\n",
+    "    well_data['NPHI'], well_data['RHOB'], gr=well_data['GR'], vsh_gr=vsh_gr,\n",
     "    dry_sand_point=args['dry_sand_point'],\n",
     "    dry_clay_point=args['dry_clay_point'],\n",
     "    corr_angle=args['hc_corr_angle'], buffer=args['hc_buffer']\n",
     ")\n",
     "\n",
     "# Estimate lithology\n",
-    "vsh_gr = estimate_vsh_gr(well_data['GR'], min_gr=0, max_gr=110)\n",
     "carb_model = Carbonate(**args)\n",
     "vclw, vcalc, vdolo = carb_model.estimate_lithology(\n",
-    "    nphihc, rhobhc,  gr=well_data['GR'], pef=well_data['PEF'], vsh_gr=vsh_gr,\n",
+    "    nphihc, rhobhc,  gr=well_data['GR'], vsh_gr=vsh_gr,  # pef=well_data['PEF'], \n",
     "    model=model, method=method, normalize= True, carbonate_type=carbonate_type\n",
     ")\n",
     "\n",
@@ -168,7 +168,7 @@
     "from quick_pp.rock_type import estimate_vsh_gr\n",
     "\n",
     "vsh_gr = estimate_vsh_gr(well_data['GR'], alpha=0.1)\n",
-    "vsh_gr2 = estimate_vsh_gr(well_data['GR'], min_gr=0, max_gr=110)\n",
+    "vsh_gr2 = estimate_vsh_gr(well_data['GR'], min_gr=0, max_gr=50)\n",
     "\n",
     "fig, axs = plt.subplots(2, 1, figsize=(20, 3), sharex=True)\n",
     "\n",

notebooks/03b_qpp_carb_single.ipynb notebooks/03_qpp_carb_single.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -22,15 +22,15 @@
     "from quick_pp.objects import Project\n",
     "\n",
     "# Load well from saved file\n",
-    "project = \"MOCK\"\n",
+    "project = \"MOCK_clastic\"\n",
     "project_path = rf\"data\\04_project\\{project}.qppp\"\n",
     "project = Project().load(project_path)\n",
     "project.get_well_names()"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 71,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -50,7 +50,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 72,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -75,7 +75,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 73,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -179,9 +179,12 @@
    "source": [
     "from quick_pp.rock_type import estimate_vsh_gr\n",
     "vsh_gr = estimate_vsh_gr(well_data['GR'], alpha=0.1)\n",
+    "vsh_gr_2 = estimate_vsh_gr(well_data['GR'], min_gr=10, max_gr=135)\n",
+    "\n",
     "fig, axs = plt.subplots(2, 1, figsize=(20, 3), sharex=True)\n",
     "\n",
     "axs[0].plot(well_data.DEPTH, vsh_gr, label='vsh_gr')\n",
+    "axs[0].plot(well_data.DEPTH, vsh_gr_2, label='vsh_gr_2')\n",
     "axs[0].legend()\n",
     "\n",
     "axs[1].plot(well_data.DEPTH, well_data['GR'], label='GR')\n",
@@ -190,7 +193,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 78,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -235,7 +238,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 80,
+   "execution_count": 13,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -247,6 +250,9 @@
     "    dry_clay_point=args['dry_clay_point'],\n",
     ")\n",
     "\n",
+    "rho_ma = rho_matrix(vsand, vsilt, vcld)\n",
+    "phid = density_porosity(rhob_corr, rho_ma)\n",
+    "\n",
     "# Normalize lithology\n",
     "vsand, vsilt, vcld = normalize_volumetric(vsand, vsilt, vcld, phit)\n",
     "\n",
@@ -255,8 +261,6 @@
     "vclb = vcld * clay_phit\n",
     "vclw = vcld + vclb\n",
     "\n",
-    "rho_ma = rho_matrix(vsand, vsilt, vcld)\n",
-    "phid = density_porosity(rhob_corr, rho_ma)\n",
     "phie = phit - vclb\n",
     "\n",
     "# Estimate rock types\n",
@@ -405,7 +409,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 85,
+   "execution_count": 18,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -438,7 +442,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 86,
+   "execution_count": 19,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -482,7 +486,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 89,
+   "execution_count": 22,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -494,7 +498,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 90,
+   "execution_count": 23,
    "metadata": {},
    "outputs": [],
    "source": [

notebooks/03_qpp_ssc_single.ipynb

@@ -90,7 +90,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "core_df = pd.read_csv(r\"C:\\Users\\imran\\OneDrive\\Desktop\\Petrophysics\\Data\\Core\\M1_M3_RAW_SCAL.csv\")\n",
+    "core_df = pd.read_csv(r\"xxx\")\n",
     "core_df = core_df[core_df.WELL_NAME == well_name].sort_values('DEPTH')"
    ]
   },

notebooks/shf-examples.ipynb

@@ -160,7 +160,7 @@ def porosity_correction_averaging(nphi, dphi, method='weighted'):
 
 
 def neu_den_xplot_hc_correction(
-        nphi, rhob, gr=None,
+        nphi, rhob, gr=None, vsh_gr=None,
         dry_sand_point: tuple = None,
         dry_clay_point: tuple = None,
         fluid_point: tuple = (1.0, 1.0),
@@ -185,7 +185,7 @@ def neu_den_xplot_hc_correction(
     D = fluid_point
     rocklithofrac = length_a_b(A, C)
 
-    frac_vsh_gr = estimate_vsh_gr(gr)
+    frac_vsh_gr = vsh_gr if vsh_gr is not None else estimate_vsh_gr(gr)
     nphi_corrected = np.empty(0)
     rhob_corrected = np.empty(0)
     hc_flag = np.empty(0)