seismicreservoirmodeling / serempy Goto Github PK

GeostatisticsDiscreteDriver.py : yields an error at line 45, because the parameter l (correlation range) is int, while function IndicatorKriging checks it is a float (line 127). The authors should check for int or float

GeostatisticsDiscreteDriver.py: figure 4, subplot 1: the title should be "Indicator Kriging Probability for facies 0". Without it the user might think it is the probability of facies 1 because the higher values are yellow, just like yellow is the color for facies 1.

GeostatisticsDiscreteDriver.py: figure 4, subplot 2: the title should be "Indicator Kriging most likely facies"

RockPhysicsInversionDriver.py: there should be titles in each figure to help the user understand what method is being used

SeismicInversionDriver.py: lines 82 and 83: substitute the function np.matlib.repmat with np.tile. At least my numpy version did not have this function (maybe it is deprecated), while np.tile worked.

ESSeisInversionDriver.py: there should be titles on the figures 2 and 3 to help the user understand what he is looking at (prior realizations and posterior, respectively)

ESPetroInversionDriver.py: include legend in the figures to help understand what is being plotted.

Line 39 in the script "Inversion.py" there is the following mathematical expression:

np.sqrt(alpha * sigmaerr)

This is not usual mathematical operation for matrices (numpy.sqrt takes the element-wise square root of the matrix), and maybe what you want to use here is the function numpy.linalg.cholesky:

np.linalg.cholesky(alpha * sigmaerr)

Cheers

Great code btw