Dear all,
I have borehole data from which I want to interpolate the basis of a 
faulted sandstone formation using kriging.  The structural analysis 
clearly shows the effect of faults with a directional variogram that is 
nonstationary perpendicular to them, but I am a bit at a loss as to how 
these should be integrated in the kriging process. The faults are all more 
or less parallel to one another, so I cannot simply divide up the study 
area in smaller subareas delimited by the faults. I have found one method 
proposed by Maréchal (Marechal, A. (1984), Kriging seismic data in 
presence of fault) to map discontinuities using an external drift function 
with a step function in it, and also a simpler approach consisting of 
"screening" all datapoints that cannot be "seen" from a point across a 
faultline (Chilès and Delfiner, Geostatistics, modeling spatial 
uncertainty), but I am at a loss as to whether either of these are already 
implemented in the R gstat library.
Advice most welcome,
Julien Farlin


Tel: +352 425991 648
Mail: julien.farlin@tudor.lu

Resource Centre for Environmental Technologies
Public Research Centre Henri Tudor
Technoport Schlassgoart
66, rue de Luxembourg
P.O BOX 144
L-4002 Esch-sur -Alzette
Luxembourg
	[[alternative HTML version deleted]]

