Re: AI-GEOSTATS: variogram and fractals

2009-03-18 Thread Syed Shibli
lags. I think that for small > lags (compared to geological layers or laminae thickness) you still > have a positive correlation.or not? > Sebas > > > > > At 02.14 18/03/2009, Syed Shibli wrote: >>Borehole logs of porosity or density can sometimes exhibit >>a

Re: AI-GEOSTATS: variogram and fractals

2009-03-17 Thread Syed Shibli
Borehole logs of porosity or density can sometimes exhibit anti-persistent behaviour, with more zero crossings around the normalised trace. On 3/17/09, seba wrote: > Dear List members > I received a lot of references from you and I was also lucky because > last week I had the possibility to be at

Re: AI-GEOSTATS: variogram and fractals

2009-03-11 Thread Syed Shibli
Try Jens Feder's book "Fractals". You can calculate the (omni) 2D variogram as per usual and plot it on a log-log scale. The fractal dimension should be D=2-(S/2) where S is the slope, and assuming the variable can be modelled using fractional Brownian motion. Syed On Wed, Mar 11, 2009 at 8:34 A

Re: AI-GEOSTATS: Re: Numerical method to solve kriging equations

2008-03-31 Thread Syed Shibli
Sorry to be nitpicky, but is that real CPU time or elapsed time (which is what you've indicated)? The latter will be quite misleading if the CPU is supporting other loads. Try this instead clock_t t1=clock(); do_whatever_it_is_youre_doing(); clock_t t2=clock(

Re: AI-GEOSTATS: kriging or IDW in case study of hydrology?

2008-02-20 Thread Syed Shibli
A strong nugget may be the exception rather than the norm for thickness data. You can try cross-validating kriged thickness results based on some a priori variogram model to see whether your estimates of thicknesses can be improved using a spatial correlation model. Syed On 2/20/08, Peter Bossew

Re: AI-GEOSTATS: DACE implementation in C

2008-02-11 Thread Syed Shibli
Try http://gstl.sourceforge.net/ which is the geostatistics template library written in C++ and developed at Stanford. If you're not familiar with objects and methods then it may be a struggle initially since C is procedural. But the algorithms lend themselves well to the OOP paradigm. Syed On F

Re: AI-GEOSTATS: upscaling methods

2007-06-01 Thread Syed Shibli
Generally arithmetic for flow parallel to the bed, harmonic for flow in series; and geometric for homogeneous formations. These are all specific cases of the power law averaging method for single phase flow. For multiphase flow upscaling, God help you. + + To post a message to the list, send it

Re: AI-GEOSTATS: nugget effect plus linear model Gslib

2007-05-29 Thread Syed Shibli
A fractal model with exponent 1 reverts to just a linear interpolation. No need for fancy kriging unless you want to derive the kriging variance. Syed + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and "unsubscr

Re: AI-GEOSTATS: FW: cosim files

2007-05-22 Thread Syed Shibli
I think COSIM was compiled using an IBM supplied TIMER routine (Fortran IV). Here is a link with code that calls F90/F95 timer routines. http://orion.math.iastate.edu/burkardt/f_src/timer/timer.html He'll have to modify it and rename the routine to TIMER or whatever is convenient for his compil

Re: AI-GEOSTATS: kriging variance and accuracy

2006-10-02 Thread Syed Shibli
application. Every technique is a double edged sword, so I guess it is a question of "know what you're doing and when"! Regards, Syed Shibli On 2 Oct 2006, at 23:42, Mario Rossi wrote: Dear all, The purpose of the paper (and also of this discussion) is to make people understand about a

Re: AI-GEOSTATS: Log versus nscore transform -- hydraulic conductivity fork

2006-08-12 Thread Syed Shibli
Michael, Carmen-Kozeny is the only attempt I know to relate particle size to permeability, backed by numerous experimental data. If you like glass beads then their formula is for you. There is also numerous empirical data looking at (log normal) distributions of pore sizes, which effectively relat

Re: AI-GEOSTATS: Log versus nscore transform

2006-08-11 Thread Syed Shibli
Michael, Looking at the plots you mentioned, I doubt that you'd ever get a straight line at high saturations purely because of the presence of the capillary transition. I am thiking that straight line (relative) conductivities will only be discernable at 100% saturation of any one phase, and w

Re: [ai-geostats] Pareto vs Lognormal distribution

2005-09-01 Thread Syed Shibli
May I suggest that you look at other analogous datasets with n > 25, e.g. the North Sea basin or Gulf of Mexico, before making some firm conclusions about whether Pareto or Lognormal works best. A lot of this information is in the public domain, one can browse the Websites of the UK DTI, Norweg

Re: [ai-geostats] Simulating an autocorrelated field

2005-06-13 Thread Syed Shibli
Have a look at Alexis Brandeker's pages (specifically his MSc thesis), among many on the Internet with sources for generating multidimensional fBm (power law) distributions, including multifractals. The link is http://www.astro.su.se/~alexis/. You may want to write him for the C code. The 1D

Re: [ai-geostats] 10th anniversary of AI-GEOSTATS

2005-04-26 Thread Syed Shibli
have learnt a tremendous amount from the various participants around the world (an have co-authored papers with colleagues I have yet to meet in real life, thanks to the power of the Net). Here's to another decade (or two, or three) of geostatistical discussion! Cheers, Syed Shibli Abe

Re: AI-GEOSTATS: Anisotropy with varying Range and Sill

2001-09-10 Thread Syed Shibli
On Monday, September 10, 2001, at 12:19 PM, Isobel Clark wrote: > Any basic Geostatistics book will show you models with > varying ranges and sills with direction. > > Except mine. We believe that sill varying with > direction is a symptom of deeper problems such as > non-stationarity, trend, d