Giter VIP home page Giter VIP logo

Comments (5)

halbmy avatar halbmy commented on June 23, 2024

Your startmodel is somehow automatically determined with 2.28, which is outside of your given saturation limit [0, 1]. Can you try with a valid starting model invert(limits=[0, 1], startModel=0.5, ...)?

from gimli.

shdlovezxm avatar shdlovezxm commented on June 23, 2024

lake.txt

The range of the Saturation model is [0,1]. No matter how I change the startModel, I cannot obtain the desired result. If I remove the model range, I can obtain the result. However, the range of the Saturation model is [0.81, 5.35], which is the wrong result?

Actual behavior

satERT = ERTPetro.invert(data, lam=20, limits=[0, 1], startModel=0.5, verbose=True)

fop: <pygimli.frameworks.modelling.PetroModelling object at 0x000002DC84675EA0>
Data transformation: <pgcore._pygimli_.RTrans object at 0x000002DC8465FD60>
Model transformation (cumulative):
	 0 <pgcore._pygimli_.RTransLogLU object at 0x000002DC8464D340>
min/max (data): 11.44/87.51
min/max (error): 2%/2.58%
min/max (start model): 0.5/0.5
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
inv.iter 1 ... chi² = 402481.45 (dPhi = 78.89%) lam: 20
--------------------------------------------------------------------------------
inv.iter 2 ... chi² = 132241.06 (dPhi = 67.14%) lam: 20.0
--------------------------------------------------------------------------------
inv.iter 3 ... chi² = 99527.0 (dPhi = 24.7%) lam: 20.0
--------------------------------------------------------------------------------
inv.iter 4 ... chi² = 99372.5 (dPhi = 0.22%) lam: 20.0
################################################################################
#                 Abort criteria reached: dPhi = 0.22 (< 2.0%)                 #
################################################################################

image

satERT = ERTPetro.invert(data, lam=20, verbose=True)

fop: <pygimli.frameworks.modelling.PetroModelling object at 0x000002DC9A609900>
Data transformation: <pgcore._pygimli_.RTrans object at 0x000002DC9A5D6D00>
Model transformation: <pgcore._pygimli_.RTransLog object at 0x000002DC9A609A90>
min/max (data): 11.44/87.51
min/max (error): 2%/2.58%
min/max (start model): 2.28/2.28
--------------------------------------------------------------------------------
--------------------------------------------------------------------------------
inv.iter 1 ... chi² = 94.54 (dPhi = 92.29%) lam: 20
--------------------------------------------------------------------------------
inv.iter 2 ... chi² = 11.3 (dPhi = 87.39%) lam: 20.0
--------------------------------------------------------------------------------
inv.iter 3 ... chi² = 0.6 (dPhi = 87.54%) lam: 20.0


################################################################################
#                  Abort criterion reached: chi² <= 1 (0.60)                   #
################################################################################

image

from gimli.

halbmy avatar halbmy commented on June 23, 2024

You took the lake case (https://www.pygimli.org/_tutorials_auto/3_inversion/plot_8-regionWise.html) which is an underwater case but you don't treat it like that. In the pure resistivity inversion, you obtain subsurface resistivities of about 8 Ohmm, which is 2.5 times below the specified fluid resistivity of 20 Ohmm. If you formally convert 8 Ohmm into saturations you end up with values above 1. So technically all is working well but the data set does not make sense in terms of Archie unless you include the water body. (In this case there is a lot of clay and Archie won't work anyway)

from gimli.

shdlovezxm avatar shdlovezxm commented on June 23, 2024

I understand. Thank you for your patient explanation.

from gimli.

halbmy avatar halbmy commented on June 23, 2024

I guess it can be closed. If there is still a problem, just reopen.

from gimli.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.