search

andymcdgeo / Petrophysics-Python-Series

A series of Jupyter notebooks showing how to load well log and petrophysical data in python.
MIT License
283 stars 132 forks source link

Create code for sort of a MultiMin or Elan Optimization approach to log analysis #1

Open Philliec459 opened 3 years ago

Philliec459 commented 3 years ago

Hi Andy, do you know how to create an Optimization example by minimize the difference between the actual TNPH and RHOB vs. the Theoretical TNPH and RHOB calculated from the log response equations for each one solving for Mineralogy?

I am thinking a simple Calcite/Dolomite/PHIT bulk volume.

RHOB_THEORETICAL = 2.65 VOL_QUARTZ + 2.71 VOL_CALCITE + 2.847 VOL_DOLO + 1.1 PHIT TNPH_THEORETICAL = -0.05 VOL_QUARTZ + 0 VOL_CALCITE + 0.005 VOL_DOLO + 1 PHIT

and VOL_QUARTZ = 0 and PHIT could be estimated from a chart book solution based on TNPH and RHOB in this clean rock.

RHOB - RHOB_THEORETICAL = DIFF_RHOB THPH - TNPH_THEORETICAL = DIFF_TNPH

DIFF_RHOB + DIFF_TNPH = 0

by varying VOL_CALCITE and minimize on DIFF_RHOB + DIFF_TNPH to approach zero.

VOL_DOLO = 1 - VOL_CALCITE - PHIT