Kinematic model of the Gippsland Basin

We analyse here uncertainties in a kinematic model of the Gippsland Basin, Victoria. The model is built on the basis of a detailed analysis of the geological history and adjusted to observed structures in the area.

Here, we first evaluate the model and analyse the setting a little bit and then proceed to investigate how errors in the parameterised geological history (i.e. the parameters of the kinematic events) propagate into model uncertainties.



In [12]:

    
!pwd









    



/Users/flow/Documents/02_work/10_Geomodels/06_Noddy



In [11]:

    
from IPython.core.display import HTML
css_file = 'pynoddy.css'
HTML(open(css_file, "r").read())









    



---------------------------------------------------------------------------
IOError                                   Traceback (most recent call last)
<ipython-input-11-64a794565859> in <module>()
      1 from IPython.core.display import HTML
      2 css_file = 'pynoddy.css'
----> 3 HTML(open(css_file, "r").read())

IOError: [Errno 2] No such file or directory: 'pynoddy.css'



In [1]:

    
%matplotlib inline



In [2]:

    
import sys, os
import matplotlib.pyplot as plt
import pynoddy.history
import pynoddy.output
import copy
import pickle



In [3]:

    
# os.chdir(r"/Users/Florian/Documents/10_Geomodels/Noddy/GBasin/GBasin")
os.chdir(r"/Users/flow/Documents/02_work/10_Geomodels/06_Noddy/")

Model analysis

First step: simple evaluation of the model structure for the new version of the model.



In [4]:

    
# the following reloads are only required during development phase
reload(pynoddy.history)
reload(pynoddy.events)
PH = pynoddy.history.NoddyHistory("GBasin_Ve1_V3.his")









    



 STRATIGRAPHY
 FOLD
 UNCONFORMITY
 FAULT
 FAULT
 UNCONFORMITY
 FAULT
 FAULT
 UNCONFORMITY
 FAULT
 FOLD
 UNCONFORMITY



In [5]:

    
# get some basic model information
print PH.get_extent()
print PH.get_origin()









    



(26630.0, 19291.0, 1500.0)
(0.0, 0.0, 1500.0)



In [6]:

    
# Let's have a look at the defined events
PH.events









    Out[6]:





{1: <pynoddy.events.Stratigraphy instance at 0x110258878>,
 2: <pynoddy.events.Fold instance at 0x110258dd0>,
 3: <pynoddy.events.Unconformity instance at 0x110258518>,
 4: <pynoddy.events.Fault instance at 0x1102580e0>,
 5: <pynoddy.events.Fault instance at 0x110258fc8>,
 6: <pynoddy.events.Unconformity instance at 0x110258758>,
 7: <pynoddy.events.Fault instance at 0x110258050>,
 8: <pynoddy.events.Fault instance at 0x110258f38>,
 9: <pynoddy.events.Unconformity instance at 0x110258830>,
 10: <pynoddy.events.Fault instance at 0x1102584d0>,
 11: <pynoddy.events.Fold instance at 0x110258b48>,
 12: <pynoddy.events.Unconformity instance at 0x1102581b8>}



In [7]:

    
# Compute the model
reload(pynoddy)
his = 'GBasin_V3.his'
PH.write_history(his)
out = 'GBasin_V3_out'
pynoddy.compute_model(his, out)
print os.getcwd()









    



/Users/flow/Documents/02_work/10_Geomodels/06_Noddy

And now we visualise the output:



In [8]:

    
reload(pynoddy.output)
PO = pynoddy.output.NoddyOutput(out)



In [9]:

    
# create section plots in axes diretions:
PO.plot_section('x', ve = 5.)
PO.plot_section('y', position = 100, ve = 5.)



In [10]:

    
# Export to VTK for 3-D visualisation and analysis
PO.export_to_vtk(vtk_filename = "GBasin_V3")

$$ \nabla^2 T = 0$$



In [10]: