A new approach for the prediction of unexposed fractured reservoirs: a case study in Millas Granite (French Pyrenees)
|
|
- Osborne Pope
- 6 years ago
- Views:
Transcription
1 Hydrological Sciences -Journal- des Sciences Hydrologiques,40,3, June A new approach for the prediction of unexposed fractured reservoirs: a case study in Millas Granite (French Pyrenees) S. PISTRE & L. M. BANGOY Laboratoire d'hydrogéologie, Place Eugène Bataillon, F Montpellier Cedex 05, France T. RIVES ELF Aquitaine, CSTJF, Avenue Larribeau, F Pau Cedex, France Abstract A fracture map is an essential aid in the study of fissured aquifers even though it often includes covered areas where the fractures are not visible. Is it possible to predict fracturing or permeable zones under the cover? This problem is dealt with using a stochastic approach based on fracture density. Classical methods of fracture density assessment involve several biases and are not adapted to hydrogeology. A proposed "fracture density index" (FDI) eliminates these biases and gives weight to two important parameters influencing the hydraulic conductivity of a zone: length and proximity of long fractures. The FDI is a regionalized variable and allows one to make predictions by kriging over covered zones. The proposed method has been applied to the Campoussy site (Pyrénées Orientales, France). In contrast to previous methods, the FDI predictions were in good agreement with the apparent hydraulic transmissivity and electrical anisotropy. The results were promising but the method needs to be tested on other sites before being considered a useful and efficient hydrogeological tool. Une nouvelle approche pour la prévision des réservoirs fissurés sous couverture: exemple du granite de Millas (Pyrénées françaises) Résumé Les cartes de fracturation sont des documents indispensables pour l'étude des aquifères fissurés, mais elles présentent souvent des zones couvertes où les fractures ne sont pas visibles. Est-il possible de prévoir la fracturation ou les zones perméables sous cette couverture? Ce problème est discuté grâce à une approche stochastique, à partir de la densité de fracturation. Les méthodes classiques de calcul de la densité de fracturation induisent plusieurs biais et ne sont pas adaptées à Fhydrogéologie. L'"index de densité de fracturation" (FDI) proposé élimine ces biais et donne beaucoup de poids à deux paramètres importants pour la perméabilité d'une zone: longueur et proximité des grandes fractures. Le FDI correspond à une variable régionalisée et permet des prévisions par krigeage sur des zones couvertes. La méthode est appliquée au site de Campoussy (Pyrénées Orientales, France). Contrairement aux méthodes existantes, les prédictions montrent une bonne corrélation avec la transmissivité hydraulique apparente et 1'anisotropic électrique. Ces résultats sont prometteurs mais la méthode doit être testée sur d'autres sites pour être considérée utile et efficace en hydrogéologie. Open for discussion until I December 1995
2 352 S. Pistre et al. INTRODUCTION Crystalline rocks have low matrix porosity so their hydraulic behaviour is mainly controlled by their fracture network geometry. A hydrogeologieal study in this type of environment provides a fracture network analysis. Fracture maps derived from direct observations (field studies and aerial photographs) are often the only reliable data for the characterization of the fracture pattern in a groundwater reservoir. However, fractures are not always visible at the ground surface, in particular where they are buried under natural cover (farmland, wooded areas) or artificial cover (inhabited areas). The characterization of fracture networks in such unexposed areas may be of economic importance and so extrapolations are then necessary. Previous studies have shown that geostatistics (Matheron, 1965, 1973), which analyses regionalized phenomena using a probabilistic approach, is a mathematical tool well adapted to rock fracture characterization at the scale of outcrops (Rossier, 1986), aerial views (Razack, 1986) or SPOT photos (Castaing et al, 1989). A fracture map is transformed into a map of fracture density values, and kriging then allows the estimation of this density parameter in unexposed areas from values in the surrounding zones. The fracture density value is assumed to be the total fracture length in each square unit (Brière, 1982) or circle unit (Kiraly, 1969) on the fracture map. The classical statistical methods involve several geometric biases and do not account for water circulation data, so a new method is proposed herein for the estimation of fracture density values specifically adapted for hydrogeology. THE FRACTURE DENSITY INDEX METHOD Fracture density maps are spatially gridded fracture maps; the fracture density value assigned to the centre of each grid unit corresponds to the total fracture length in that unit. In order to produce a continuous field of fracture density values, this gridding technique involves several biases: in particular, where a long fracture crosses a unit but not an adjacent unit, then two neighbouring fracture density values may be very different. This could involve a discontinuity of data corresponding to an accentuation of the local variability of the fracturing phenomenon, difficult to consider from a geostatistic point of view, and which could be conditioned by the accuracy of the map. Furthermore, for hydrogeologieal applications, such a discontinuity is inconsistent with the spatial variability of the hydrodynamic parameters which are rather continuous at this scale, owing to the small fracture network. This bias may be partially avoided by shifting the grid a half unit length successively in two orthogonal directions and then recalculating the fracture density values. Nevertheless, the parameters of most influence on productivity in drilling are the proximity and length of important fractures which are not taken into account in these methods.
3 Prediction of unexposed fractured reservoirs 353 The fracture density index method (Pistre, 1992a), specifically developed for hydrogeology, aims to predict the fracturing of unexposed zones without any geometric bias. It involves three successive stages: Fracture map In order to consider the entire aquifer, the fracture map must cover a much wider area than the zone of interest. The scale is determined from the minimum fracture length influencing global circulations in the reservoir (generally ten to several tens of metres). Aerial views are generally the means most adopted. Fractures are drawn as linear segments and covered zones, where fractures are unexposed, are shaded. Fracture density index map The fracture map is then transformed into a fracture density index map. A regular grid is placed on the map and the fracture density index is calculated for the centre of each unit. The size of unit is sufficient to allow accurate kriging approximately ten times the size of the unexposed zone, but this has to be limited by time considerations. The fracture density index FDI(x,y) at a point P(x,y), for a single fracture i, is defined as: FDIjfXy) = [fir) dr (1) where r u and r 2i are the distances between P and the two extremities of the fracture. The index FDI depends on the fracture length and on the relative position of P(x,y) and of the fracture (Fig. 1(a)). Fig. 1 Fracture density index estimation for one fracture, at one point.
4 354 S. Pistre et al. In order to adapt the FDI to hydrogeology, one needs to give a large weight to long or proximal fractures and a small weight to short or distal fractures whose hydrodynamic influence on P(x,y) is negligible. This weighting can be made by assuming/(r) to be a function derived from the Gaussian law fir) = e" r2/2t2 where T has the dimension of length. The graph of this function (Fig. 1(b)) shows the influence of fracture proximity and length on FDI ; (x,y), which can be modified by varying T. This latter parameter is derived from the estimated average distance for which a fracture will have only a minor influence on the hydrogeological characteristics at P(x,y), i.e. for which exp( r 2 /!-! 2 ) becomes insignificant (Fig. 2(a)). The effect of a fracture is considered to be negligible at r = AT (Fig. 1(b)). For a set of N f fractures FDI(x,y) is defined as: N f r n FDI(x,y) = f e^^dr ( 2 ) where FDI(J,)>) is in metres of fractures. The method obviously implies an edge effect at map boundaries or close to shaded zones where the "visible" fracture field is reduced. This bias remains restricted (Fig. 2(b)) but is likely to alter the validity of interpolation calculations close to hidden areas. A correcting factor A(x,y) is therefore introduced which takes this bias into account: A(x,y) = ^fir)ds/ N, J>)d 5 -ÇJ/(r)d5 (3)
5 Prediction of unexposed fractured reservoirs 355 where N z is the number of hidden areas, r represents the distance between P(x,y) and the surface differential ds, 5, represents the area of each hidden zone (Fig. 1) and S is the whole area of the field considered. The corrected fracture density index FDI corr (x,y) becomes: N f hi FDI corr (x,y) = J2 f e~ r2,2t2 dr.a(x,y) (4) Calculation of FDI corr (x,y) at each point of the grid leads to the fracture density index map. Prediction of fracture networks by kriging The fracture density index map is an incomplete field of numerical values corresponding to the fracture field of the aquifer. The value FDI corr (x,y) may be considered as a regionalized variable z(x) = FDI corr (x,y) (X represents simultaneously the two coordinates of a point studied) which corresponds to the realization of a random function at each point of the zone studied (Delhomme, 1978). The spatial variability of z(x) = FDI corr (x,y) on the map can be described by an experimental variogram and characterized in numerical terms by adjustment to a theoretical variogram (De Marsily, 1986). The knowledge of the values FDI corr (x,y) at points close to hidden areas and the spatial variability of this parameter allows estimation of the fracture density index FDI(x,y) in unexposed zones by interpolation (point kriging). The more the variable is structured and regular, the more this interpolation is relevant or accurate (deviation estimations are very small). The kriging results in a complete fracture density index map. APPLICATION TO THE MILLAS GRANITE The Campoussy experimental site is located in the late Variscan calc-alkaline granitoids of the Quérigut-Millas massif (Pyrénées Orientales, France, Fig. 3). This granitic massif, located along the North Pyrenean Fault, has a complex geological history (Choukroune, 1976). The development of the observed fracture network can be related to two important tectonic events: the brittle pyrenean north-south compression and the NW-SE Oligo-Miocene extension (Soliva et al., 1992; Arthaud & Pistre, 1993). Intensely fractured, the massif is, in a general way, an aquifer. Moreover, the relief and large faults create a pattern of reservoirs with an extent of several square kilometres which are often in elevated positions with free, semi-captive or captive characters according to the presence or not of an alluvial cover. The Campoussy site is an area where the granitic massif is covered by alluvial deposits. The problem is to characterize its fracturing.
6 356 S. Pistre et al. Fig. 3 Location of the Querigut-Millas massif in the Pyrenees Axial Zone, and of the area studied (black square): 1- gneissic massif; 2- granite massifs; and 3- sedimentary formations. Data collection and corresponding fracture density index map The field selected corresponds to a 5 km 2 rectangle (the entire aquifer) roughly centred on the experimental site. Fracture traces longer than 50 m were drawn from high quality aerial photographs at 1: scale, complemented with field observations (Pistre, 1989). A map of the fractures was obtained with shaded areas representing the unexposed zones (Fig. 4). The fracture network comprised two dominant sets: N170 to N010, N030 to N050, and two subordinate sets: N080 to N100, N120 to N135 (Fig. 5). These four fracture sets occurred over the entire massif (Pistre, 1992b; Arthaud & Pistre, 1993). The values of FDI corr (jc,y) were calculated for each point on a regular grid (50 m unit length) and contours were drawn for different r values (25, 50, 100 and 200 m) to cover the range of the possible average distances ( m) for which the influence of a fracture is not negligible (see above) (Bangoy, 1992). Figure 6 shows the map for T = 50 m. Prediction of fracture network by kriging Maps of the selected field were calculated for the four T values. With increasing T, the map contours became smoother and showed less detail (e.g. alignments, fracture intersections). The selection of r, which reflects the influence of fracture closeness and length, therefore determines the behaviour
7 Prediction of unexposed fractured reservoirs 357 Fig. 4 Fracture map drawn from aerial views: location of the hydrogeological site at Campoussy (star), Fig. 5 Fractures directional histogram, from the fractures cumulated lengths, with azimuthal classes of 10. of the variogram which indicates the spatial variability of FDI and which guides kriging. Thus, this selection is based not only on hydrodynamic considerations (see above), but also, for technical reasons (kriging), on the variograms. They show oscillations (with more instabilities when T decreases) at approximately identical distances for different T values (Fig. 7), which correspond to "hole effects" ("rich" areas surrounded by "poor" ones) except at 1400 m (not visible for T = 25 m) where the regionalization of the phenomenon is proved. The T chosen needs to lead to a reliable map in total
8 358 S. Pistre et al. Ill Recent deposits mm Campoussy site Fig. 6 Map of the fracture density index deduced from Fig. 2, with r = 50 m. T=25m 1=50 m Spherical + Linear Spherical + Linear 0 -,MII.. I. 1 I 1 I.. I MI 1 I.... IM! 0- I...1.1". L h, m h, m 60i m 2 : 40- T= =100m /? so, 2 : m T=200m 20-2 x Spherical i «Se ibbo"" iebb h, m h, m Fig. 7 Experimental and adjusted variograms for different values of T. accordance with the preferential alignments of the fracture map. The value T = 25 m provides the most contrasted map, but is not representative because the distinctness is only illusory according to the resolution of the aerial views. All the variograms were adjusted (Fig. 7) to spherical and linear models as follows:
9 Prediction of unexposed fractured reservoirs 359 u> 1 h + u) 2 [3f2(hfa)-l/2(.h/a) 3 ] + C Q h < a a)j/z + co 2 + C() h > a where a is range and co 1; w 2 are positive coefficients. The adjustments with linear and spherical models were satisfactory for all the variograms (the average difference between the experimental and analytical variograms was less than 3% of the variance values) with the prospect of further calculations during interpolation. The value T = 50 m was selected considering its variogram: - FDI was structured and regular for distances below 300 m. The presence of this 300 m large regionalization was not visible for higher r, whereas it was visible for T = 25 m. For distances beyond 300 m the variogram was quite stable but quite flat (more than for r = 25 m): FDI was no more regular. This structure was important for interpolation under the cover because its dimension was about that of the covered zones: the 300 m distance values, which were barely correlated, had a weak influence on the kriging. - The behaviour at the origin did not show any nugget-effect, T = 50 m seemed to be well adapted to this fracture map, and the nearest values were correlated (the "erratic" component of the FDI was small). This would allow a more stable or more "accurate" kriged map to be obtained, reducing the standard deviation of the estimation. The kriging allowed a prediction of the fracture density index below the Campoussy alluvial parcel. Maps of predicted FDI values and associated standard deviations (Fig. 8) were prepared. The FDI values obtained were larger in the eastern and central parts of the site than in the western part. Fig. 8 (a) kriged map of FDI at the Campoussy site; and (b) map of the standard deviation of estimation (m). PREDICTION VALIDITY The fracture density index method is one of the methods providing a regionalized variable which shows fracture density on a fracture map, and is then able
10 360 S. Pistre et al. to predict fracture densities of unexposed zones. In order to test such predictions, the coherence of the FDI method with other approaches is investigated below. Comparison with previous methods The classical grid methods differ from the fracture index method in the construction of fracture density maps. In the simple grid method (Kulatilake & Wu, 1984) a regular grid is superimposed on the fracture map and the total j 1 Recent deposits ^M Campoussy site Fig. 9 Maps of fracture density, with a 100 m x 100 m unit (isodensities in m per unit), deduced from Fig. 2: estimated from (a) the simple grid method; and (b) the multiple grid method.
11 Prediction of unexposed fractured reservoirs 361 a ;' t ;' ; ' ' / / / / / / // ; / : : i ' / /' / ; /, - 4 y ','' / ^ 60 /,, 0 ' t,--' -r 20m Fig. 10 Kriged map of fracture density (m), with a 50 m x 50 m unit, on the Campoussy site for (a) the simple grid method; and (b) the multiple grid method. fracture length of each unit square is measured. The calculated value of the fracture density d(x,y) is expressed in metres of fracture for the chosen unit and is assigned to the centre of each unit. The unit size is determined by the aerial photograph resolution (100 m x 100 m was chosen). The resulting map (Fig. 9(a)) was quite similar to the fracture density index map (Fig. 6), the only differences being the amplitudes of spatial variations and many details close to the unexposed zones. Nevertheless, the kriged map of the Campoussy site (Fig. 10(a)) was very different from that produced by using the fracture density index method (Fig. 8). In the multiple grid method, the grid is translated half a unit length in N-S and E-W directions in order to multiply the points of determination of the fracture density. This overlapping of squares makes the edge effect disappear. The map obtained (Fig. 9(b)) presented the same characteristics as the fracture density maps drawn with the previous methods. However, it was more contrasted and appeared to point out more details than Fig. 9(a), but it needed a longer time of calculation and bulkier files (x4 for a similar fractures field). The interpolated map of the Campoussy site (Fig. 10(b)) was also very different from Fig. 8. The three methods led to predictions of the fracture density of the site. In spite of biases inherent to the grid methods - edge effect (linked to the position of each square unit) and shape effect (influence of fracture orientations) one could not determine which predicted map was correct. Onlycomparison with geophysical and hydrological data from the Campoussy site could allow the best method to be determined. Validation by other approaches A series of short pumping tests in the seven wells of the Campoussy site (F1-F7 on Fig. 11) was made in order to investigate the hydraulic properties of the captive aquifer located in the deep fissured granite. The hydraulic
12 362 S. Pistre et al. behaviour of the fissured granite appeared to be equivalent to a porous medium (Pistre, 1992b): among the classical models for pumping-test interpretation, the Theis model gave the best results (except for the beginning of pumping) for each borehole. The presence of many small-sized fractures, as observed at the surface (Arthaud & Pistre, 1993), can explain this behaviour. Assuming the reservoir as a porous medium without any leakage, a map of "apparent transmissivities" may be drawn (Fig. 11). The apparent transmissivity is the experimental value which represents the conductive properties of a volume of granite functioning at the time of pumping. Because the circulations in fissured media occur in fractures, the apparent transmissivity is an increasing function of fracture density. An east-west profile across the Campoussy site showed an apparent transmissivity curve (Fig. 12), which suggested a greater fracture density in the centre of the site. A series of electrical surveys were made at the Campoussy site in order to determine the electric properties of the medium and to characterize its structure (Aslanian, 1989). Schlumberger soundings across the Campoussy site showed three successive strata interpreted as alterites, altered granite and unaltered \ \ '. / F2 / F3 # / F4 \ 1 * C1 V \ D C4 0 20m \ / Fig. 11 Map of the apparent transmissivity ( x l(r 5 m 2 s" 1 ) at the Campoussy site. Drill locations (Fl to F7); Schlumberger locations (SI to S4) and square electrical soundings locations (CI to C4). \ ' /CM*?'' / F V\V>\ \ F6 S3 C2y ; *^.\ \ \ t 3-i 1 i i i i i i i i i i i i i i i i i i i " " m Fig. 12 Comparisons at the site between electrical anisotropy coefficient (\: squares), apparent transmissivity (Ta: stars = data; dashed curve = interpolation) and interpolated fracture densities computed from simple grid (iil), multiple grids (dz) and fracture density index (FDI) methods.
13 Prediction of unexposed fractured reservoirs 363 granite with an apparent resistivity (p am ) of about 700 ohm m. Schlumberger drag soundings revealed increasing thickness of the altered zone to the west (from 15 to 25 m). Three multidirectional square soundings (C1-C3 on Fig. 11) showed an anisotropy of the fresh granite along the NE-SW direction. This anisotropy could be interpreted as indicating the presence of N050 fractures at depth. The values of p am and X, the elongation coefficient (i. e. the anisotropy coefficient), are indicative of fracture density (Haberjam, 1972; Darboux-Afouda, 1985). The value of p am was constant across the site and X was greater in the centre (Fig. 12), suggesting a higher fracture density in the centre of the Campoussy site. The hydraulic and electrical methods showed a variation in the east-west structure of the Campoussy site; in particular, a probably more densely fractured zone in its central part. The east-west profiles of the estimated fracture densities obtained by the three methods (Fig. 12) showed a constant fracture density in the simple grid method, a regularly increasing fracture density from west to east in the multiple grid method and a maximum fracture density in the centre of the site in the fracture density index method. The fracture density index profile was the most coherent with the direct hydraulic measurements and electrical surveys, which indicated that the fracture density index method was the most valid fracture density approach, at least for hydrogeological considerations. In the western part of the profile, the apparent transmissivity decreased more quickly than the FDI probably because the connectivity decreased more quickly than the density of fractures (Robinson, 1984). DISCUSSION AND CONCLUSIONS The fracture density index (FDI) is a new variable with a hydrodynamic connotation. It can be considered as a realization of a regionalized variable of a studied area and so allows predictions by kriging in unexposed zones. The fracture density index method avoids geometric biases (effects of the position and the orientation of each unit) and hydrodynamic bias (cutting of long fractures) which are encountered in the simple and multiple grid methods: thus it corrects hidden area effects. Like the other methods, it implies partial interpretations to select the parameter values, in particular for the selection of the fractures to be taken into account. The minimum fracture length of the fracture map corresponds to the minimum fracture length which could significantly influence the hydrodynamic behaviour of groundwater reservoirs. At the Campoussy site, this was assumed to be fractures of tens of metres. As shown by hydraulic tests, the FDI variable seemed to be linked with the hydraulic transmissivity of the reservoir, so the kriged maps indicated the poorly and highly fractured zones. The reliability of predictions depends on the parameter r which reflects the distance (AT) from which a fracture will have only a minor hydraulic influence on a possible well. To provide useful results, kriging preferably necessitates a continuous variable. So, for this influence distance, several r values need to be tested. The behaviours at the origin of the corresponding
14 364 S. Pistre et al. variograms permit selection of the best r value (in the example given herein T = 50 m). The fracture density index method necessitates more complex and longer calculations than those of the other methods. Moreover, in order to validate the FDI method, in particular the link between FDI and hydraulic transmissivity, other applications will have to be carried out at other sites. FDI allows characterization of three dimensional hydraulic behaviour of a groundwater reservoir from two dimensional surface (fracture) data. A remaining problem is the choice of the fracture sets to be taken into account since in a fissured network not all fractures contribute equally to fluid circulation. For example, at the Campoussy site the N170-N010 and N030-N050 fractures, which were the most represented, could only be used following a tectonic analysis of the domain (Arthaud & Pistre, 1993) which attributed the N030-N050 set formation and opening of the N170-N010 fractures to the Oligo-Miocene extension. Moreover, once may also consider that all the fractures are potential drains. The N080-N100 and N120-N135 fracture sets insure local drainage through the main fractures and could not be neglected. The fracture density index method constitutes a promising new predictive technique for hydrogeology. It can be improved by other applications on other sites and by taking into account fracture organization and hierarchy. It is necessary, therefore, to understand fracture development mechanisms via field analysis in well exposed areas (Rawnsley et al., 1992), analogue modelling (Rives & Petit, 1990a, 1990b) and numerical simulations (Rives et al., 1992) to guide further extrapolations. REFERENCES Arthaud, F. & Pistre, S. (1993) Les fractures et les paléocontraintes du granite hercynien de Millas (Zone Axiale des pyrénées). Un cas d'étude de la tectonique cassante d'un aquifère de socle. Geodin. Acta, 6(3), Aslanian, D. (1989) Approche électrique de la géométrie d'un magasin aquifère en milieu cristallin. DEA, Montpellier II, France. Bangoy, L. M. (1992) Hydrodynamique d'un site expérimental en aquifère de socle fissuré; nouvelle méthode d'interprétation des essais hydrauliques. Thèse Doctorat es Sciences, Université de Montpellier II, France. Brière, G. (1982) Introduction à l'analyse géostatistique sur clichés aériens de la fracturation des magasins aquifères en roches fissurées. Thèse 3 Cycle, Montpellier II, France. Castaing, C, Dutartre, P., Gouyet, J.F., Loiseau, P., Martin, P. & Pointet, T. (1989) Etude pluridisciplinaire d'un réseau de discontinuités-images SPOT en milieu granitique couvert. Implication en hydrogéologie des milieux fissurés. Bull. BRGM, no. 1, Choukroune, P. (1976) Structure et évolution tectonique de la zone nord-pyrénéenne. Mém. Soc. Géol. France, no. 127, Darboux-Afouda, 8. (1985) Recherches hydrogéologiques par les méthodes géophysiques dans les régions de socle en Afrique Occidentale. Thèse 3 Cycle, Montpellier II, France. Delhomme, J. P. (1978) Application de la théorie des variables régionalisées dans les sciences de l'eau. Bull. BRGM, série 2, section III, De Marsily, G. (1986) Quantitative Hydrogeology. Groundwater Hydrogeology for Engineers. Academic Press, Orlando, USA. Haberjam, G. M. (1972) The effects of anisotropy on square array resistivity measurements. Geophys. Prosp. 20, Kiraly, L. (1969) Statistical analysis of fractures (orientation and densities). Geol. Rundschau, 59(2),
15 Prediction of unexposed fractured reservoirs 365 Kulatilake, P. H. S. W. & WU, T. H. (1984) The density of discontinuity traces in sampling windows. J. Rock Mech. & Geomech. Abstr. 21(6), Matheron, G. (1965) Les Variables Régionalisées et Leur Estimation. Masson, Paris, France. Matheron, G. (1973) The intrinsic random functions and their applications. Adv. in Appl. Prob. 5, Pistre, S. (1989) Approche tectonique et microtectonique de la géométrie d'un aquifère fissuré en zone de socle: exemple du site expérimental de Campoussy dans le massif de Millas. DEA, Montpellier II, France. Pistre, S. (1992a) Méthode d'évaluation d'un indice de densité de fracturation en hydrogéologie et prévision des zones perméables par une approche probabiliste. C. R. Acad. Se. Paris, 314(11), Pistre, S. (1992b) Structure fissurale et modélisation hydrodynamique. Thèse Doctorat es Sciences, Université de Montpellier II, France. Rawnsley, K. D., Rives, T., Petit, J. P., Hencher, S. R. & Lumsden, A. C. (1992) Joint development in perturbed stress fields near faults. /. Struct. Geol. 14, Razack, M. (1986) Approche probabiliste de l'étude en subsurface de la géométrie des réservoirs fissurés. Effet de l'échelle d'investigation. Bull. BRGM Hydrogéologie no. 2, Rives, T. & Petit, J. P. (1990a) Diaclases et plissement: une approche expérimentale. C. R. Acad. Se. Paris, 310(11), Rives, T. & Petit, J. P. (1990b) Experimental study of jointing during cylindrical and non-cylindrical folding. In: Mechanics of Jointed and Faulted Rock, ed. H. P. Rossmanith, Balkema, Rotterdam, The Netherlands. Rives, T., Razack, M., Petit, J. P. & Rawnsley, K. D. (1992) Joint spacing: analogue and numerical simulations. /. Struct. Geol. 14, Robinson, P. C. (1984) Connectivity, flow and transport in network models of fracture media. PhD dissertation, St Catherine's College, Oxford, UK. Rossier, Y. (1986) Exemple d'applications d'une méthode d'analyse de la fissuration. Illustration en pays granitique et calcaire. Bull. BRGM Hydrogéologie no. 2, Soliva, J., Pistre, S., Brunei, M. & Arthaud, F. (1992) Les zones de cisaillement mylonitiques dans le granite hercynien de Millas (Pyrénées Orientales): déformation ductile d'âge hercynien et/ou pyrénéen. C. R. Acad. Se. Paris 314(11), Received 24 January 1994; accepted 9 December 1994
16
Estimation of the transmissmty of the Santiago aquifer, Chile, using different geostatistical methods
Groundwater Management: Quantity and Quality (Proceedings of the Benidorm Symposium, October 1989). IAHS Publ. no. 188,1989. Estimation of the transmissmty of the Santiago aquifer, Chile, using different
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationThe presentation is entitled ANISOTROPY OF THE ROCK QUALITY DESIGNATION & ITS GEOSTATISTICAL EVALUATION».
It is my pleasure to present the work done jointly by my collaborator Doctor Serge Séguret from France, and myself, on data provided by my team and by my colleague Claudio Rojas The presentation is entitled
More informationCase Study: University of Connecticut (UConn) Landfill
Case Study: University of Connecticut (UConn) Landfill Problem Statement:» Locate disposal trenches» Identify geologic features and distinguish them from leachate and locate preferential pathways in fractured
More informationSome elements for improving interpretation of concrete electrical resistivity
Some elements for improving interpretation of concrete electrical resistivity Jean-François LATASTE 1, Stéphane LAURENS 2, Andrzej MOCZKO 3 1 Université Bordeaux 1, Talence, France, jf.lataste@ghymac.u-bordeaux1.fr
More informationAzimuthal Square Array Configuration and Groundwater Prospecting in Quartzite Terrian at Edaikkal, Ambasamudram, Tirunelveli
Research Journal of Earth Sciences 4 (2): 49-55, 2012 ISSN 1995-9044 IDOSI Publications, 2012 DOI: 10.5829/idosi.rjes.2012.4.2.2937 Azimuthal Square Array Configuration and Groundwater Prospecting in Quartzite
More informationStructural and Petrophysical Characterization of Mixed Drain/Barrier Fault Zones in Carbonates: Example from the Castellas Fault (SE France)*
Structural and Petrophysical Characterization of Mixed Drain/Barrier Fault Zones in Carbonates: Example from the Castellas Fault (SE France)* Christophe Matonti 1, Juliette Lamarche 1, Yves Guglielmi 1,
More informationDETERMINATION OF A SALT-WATER INTERFACE BY ELECTRIC RESISTIVITY DEPTH SOUNDINGS
Hydrological SciencesBulletin des Sciences Hydrologiques, XXI, 4 2/976 DETERMINATION OF A SALTWATER INTERFACE BY ELECTRIC RESISTIVITY DEPTH SOUNDINGS A. GINZBURG Department of Environmental Sciences, Tel
More informationNumerical analysis for an interpretation of the pressuremeter test in cohesive soil
Numerical analysis for an interpretation of the pressuremeter test in cohesive soil J. Monnet Joseph Fourier University L3S-R Domaine Universitaire, BP n 53, 3841, Grenoble, Cedex 9, France jmonnet@ujf-grenoble.fr
More informationNOTICE CONCERNING COPYRIGHT RESTRICTIONS
NOTICE CONCERNING COPYRIGHT RESTRICTIONS This document may contain copyrighted materials. These materials have been made available for use in research, teaching, and private study, but may not be used
More informationAgeostrophic instabilities of a front in a stratified rotating fluid
8 ème Congrès Français de Mécanique Grenoble, 27-3 août 27 Ageostrophic instabilities of a front in a stratified rotating fluid J. Gula, R. Plougonven & V. Zeitlin Laboratoire de Météorologie Dynamique
More informationFracture characterization from scattered energy: A case study
Fracture characterization from scattered energy: A case study Samantha Grandi K., Sung Yuh, Mark E. Willis, and M. Nafi Toksöz Earth Resources Laboratory, MIT. Cambridge, MA. Total Exploration & Production.
More informationEAS 233 Geologic Structures and Maps Winter Miscellaneous practice map exercises. 1. Fault and separation:
Miscellaneous practice map exercises 1. Fault and separation: With respect to Map 1, what are (a) the orientation of the fault, and (b) the orientation of bedding in the units cut by the fault. (c) Mark
More informationTh E Natural Fracture Prediction for Discrete Fracture Modelling
Th E103 01 Natural Fracture Prediction for Discrete Fracture Modelling H. Phillips* (Schlumberger), J.P. Joonnekindt (Schlumberger) & L. Maerten (Schlumberger) SUMMARY There are many uncertainties in the
More informationTu D Understanding the Interplay of Fractures, Stresses & Facies in Unconventional Reservoirs - Case Study from Chad Granites
Tu D201 04 Understanding the Interplay of Fractures, Stresses & Facies in Unconventional Reservoirs - Case Study from Chad Granites D. Lirong (Chinese National Petroleum Company Ltd. (Chad)), C. Shrivastava*
More informationINTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY
INTEGRATION OF BOREHOLE INFORMATION AND RESISTIVITY DATA FOR AQUIFER VULNERABILITY ANDERS V. CHRISTIANSEN 1, ESBEN AUKEN 1 AND KURT SØRENSEN 1 1 HydroGeophysics Group, Aarhus University, Finlandsgade 8,
More informationAzimuthal Resistivity to Characterize Fractures in a Glacial Till. Mark Boris, University of Saskatchewan Jim Merriam, University of Saskatchewan
Azimuthal Resistivity to Characterize Fractures in a Glacial Till Mark Boris, University of Saskatchewan Jim Merriam, University of Saskatchewan Abstract Azimuthal resistivity was used to characterize
More informationHybrid Modelling using Neuro Fractal for Fractured Reservoirs
Hybrid Modelling using Neuro Fractal for Fractured Reservoirs NAM H. TRAN *, KAREN VALENCIA, KIEN TRAN AND SHEIK S. RAHMAN (*) School of Petroleum Engineering The University of New South Wales Sydney NSW
More informationL'importance de la présence de dykes et filons de quartz fracturés pour la recherche d'eau dans les roches métamorphiques argileuses
L'importance de la présence de dykes et filons de quartz fracturés pour la recherche d'eau dans les roches métamorphiques argileuses The importance of the presence of fractured quartz dykes and veins in
More informationEstimation of direction of increase of gold mineralisation using pair-copulas
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Estimation of direction of increase of gold mineralisation using pair-copulas
More informationDETERMINING HIGH VOLTAGE CABLE CONDUCTOR TEMPERATURES. Guy Van der Veken. Euromold, Belgium. INVESTIGATIONS. INTRODUCTION.
DETERMINING HIGH VOLTAGE CABLE CONDUCTOR TEMPERATURES. Guy Van der Veken. Euromold, Belgium. INTRODUCTION. INVESTIGATIONS. Type tests on MV cable accessories are described in CENELEC HD68 and HD69 documents.
More informationA. Baiocchi 1, W. Dragoni 2, F. Lotti 1, S.M. Piacentini 2, V. Piscopo 1
A. Baiocchi 1, W. Dragoni 2, F. Lotti 1, S.M. Piacentini 2, V. Piscopo 1 (1) Department of Ecological and Biological Sciences, University of Tuscia, Viterbo, Italy (2) Department of Physics and Geology,
More informationKey words groundwater; hard-rock aquifer; weathered profile; major ions; residence time; numerical model; Basque Country, France
386 Evolving Water Resources Systems: Understanding, Predicting and Managing Water Society Interactions Proceedings of ICWRS2014, Bologna, Italy, June 2014 (IAHS Publ. 364, 2014). Coupling hydrology, geochemistry
More informationGeologic Considerations of Shallow SAGD Caprock; Seal Capacity, Seal Geometry and Seal Integrity, Athabasca Oilsands, Alberta Canada
Geologic Considerations of Shallow SAGD Caprock; Seal Capacity, Seal Geometry and Seal Integrity, Athabasca Oilsands, Alberta Canada Gordon T. Stabb, Michael Webb Durando Resources Corp, Suncor Energy
More informationB01 Surface-wave Dispersion Stacking on a Granitemicaschists Contact at Ploemeur Hydrological Observatory, France
B01 Surface-wave Dispersion Stacking on a Granitemicaschists Contact at Ploemeur Hydrological Observatory, France S. Pasquet* (UMR CNRS 7619 Sisyphe), L. Bodet (UMR CNRS 7619 Sisyphe), L. Longuevergne
More informationAdebayo O. Ojo, M.Sc. 1* and Martins O. Olorunfemi, Ph.D *
A Graphical and Semi-Quantitative Technique for Investigating Vertical Electrical Sounding (VES) Curves for Indices of Confined Fractured Basement Column. Adebayo O. Ojo, M.Sc. 1* and Martins O. Olorunfemi,
More informationCombining geological surface data and geostatistical model for Enhanced Subsurface geological model
Combining geological surface data and geostatistical model for Enhanced Subsurface geological model M. Kurniawan Alfadli, Nanda Natasia, Iyan Haryanto Faculty of Geological Engineering Jalan Raya Bandung
More informationSimultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations
Simultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations C. Rentier,, A. Dassargues Hydrogeology Group, Departement
More informationSTOCHASTIC CONTINUUM ANALYSIS OF GROUNDWATER FLOW PATHS FOR SAFETY ASSESSMENT OF A RADIOACTIVE WASTE DISPOSAL FACILITY
STOCHASTIC CONTINUUM ANALYSIS OF GROUNDWATER FLOW PATHS FOR SAFETY ASSESSMENT OF A RADIOACTIVE WASTE DISPOSAL FACILITY K. Chang*, C.L. Kim, E.Y. Lee, J.W.Park, H.Y.Park, C.G. Rhee, M.J. Song Nuclear Environment
More informationHYDROGEOLOGICAL PROPERTIES OF THE UG2 PYROXENITE AQUIFERS OF THE BUSHVELD COMPLEX
R. Gebrekristos, P.Cheshire HYDROGEOLOGICAL PROPERTIES OF THE UG2 PYROXENITE AQUIFERS OF THE BUSHVELD COMPLEX R. Gebrekristos Digby Wells Environmental P. Cheshire Groundwater Monitoring Services Abstract
More informationThe failure of the sounding assumption in electroseismic investigations
The failure of the sounding assumption in electroseismic investigations F.D. Fourie, J.F. Botha Institute for Groundwater Studies, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
More informationFRACTURE-DISTRIBUTION MODELING IN ROCK MASS USING BOREHOLE DATA AND GEOSTATISTICAL SIMULATION
FRACTURE-DISTRIBUTIO MODELIG I ROCK MASS USIG BOREHOLE DATA AD GEOSTATISTICAL SIMULATIO Katsuaki Koike 1, Kazuya Komorida 2, and Yuichi Ichikawa 3 1: Department of Civil Engineering, Faculty of Engineering,
More informationPRINCIPAUX ENJEUX LIÉS AUX AQUIFÈRES DE SOCLE
PRINCIPAUX ENJEUX LIÉS AUX AQUIFÈRES DE SOCLE John M. Sharp, Jr. Department of Geological Sciences Jackson School of Geosciences The University of Texas Austin, Texas, USA Uwe Troeger* Lehrstuhl fuer Hydrogeologie
More informationDeveloping rules of thumb for groundwater modelling in large open pit mine design
Developing rules of thumb for groundwater modelling in large open pit mine design Jim Hazzard, Branko Damjanac, Christine Detournay & Loren Lorig Itasca Consulting Group, Minneapolis, MN, USA ABSTRACT
More informationMesurer des déplacements ET des contraintes par correlation mécanique d images
Mesurer des déplacements ET des contraintes par correlation mécanique d images J. RÉTHORÉ a, A. LEYGUE a, M. CORET a, L. STAINIER a, E. VERRON a a. Institut de Recherche en Génie Civil et Mécanique (GeM)
More informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationHYDRAULIC PROPERTIES OF UNCONSOLIDATED SEDIMENTARY AQUIFERS IN KRIBI, CAMEROON: ENVIRONMENTAL IMPLICATIONS
HYDRAULIC PROPERTIES OF UNCONSOLIDATED SEDIMENTARY AQUIFERS IN KRIBI, CAMEROON: ENVIRONMENTAL IMPLICATIONS Vol. 13, n 1, 12 19 Anatole Djieto Lordon 1*, Nenita Ntumba Bukalo 1 and Christopher M. Agyingi
More informationOutils de Recherche Opérationnelle en Génie MTH Astuce de modélisation en Programmation Linéaire
Outils de Recherche Opérationnelle en Génie MTH 8414 Astuce de modélisation en Programmation Linéaire Résumé Les problèmes ne se présentent pas toujours sous une forme qui soit naturellement linéaire.
More informationA CONSTITUTIVE MODEL TO PREDICT THE HYDROMECHANICAL BEHAVIOUR OF ROCK JOINTS
OttawaGeo27/OttawaGéo27 A CONSTITUTIVE MODEL TO PREDICT THE HYDROMECHANICAL BEHAVIOUR OF ROCK JOINTS Dominic Tremblay, Richard Simon and Michel Aubertin Department of civil, geological & mining engineering
More informationModelling effective permeability of fracture networks in permeable rock formation by singular integral equations method
Author manuscript, published in "6th International Conference on Computational and Experimental Methods in Multiphase and Complex Flow., Kos : Greece (2011)" DOI : 10.2495/MPF110241 6th International Conference
More informationPETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR
PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR APPLIED GRADUATE STUDIES Geology Geophysics GEO1 Introduction to the petroleum geosciences GEO2 Seismic methods GEO3 Multi-scale geological analysis GEO4
More informationIs there still room for new developments in geostatistics?
Is there still room for new developments in geostatistics? Jean-Paul Chilès MINES ParisTech, Fontainebleau, France, IAMG 34th IGC, Brisbane, 8 August 2012 Matheron: books and monographs 1962-1963: Treatise
More informationConcept of an Integrated Workflow for Geothermal Exploration in Hot Sedimentary Aquifers
Concept of an Integrated Workflo for Geothermal Exploration in Hot Sedimentary Aquifers J. Florian Wellmann, Franklin G. Horoitz, Klaus Regenauer-Lieb Western Australian Geothermal Centre of Excellence,
More informationGEOELECTRICAL INVESTIGATIONS OF MARINE INTRUSIONS ON THE ROMANIAN BLACK SEA SHORE
GEOELECTRICAL INVESTIGATIONS OF MARINE INTRUSIONS ON THE ROMANIAN BLACK SEA SHORE P. GEORGESCU, D. IOANE, B.M. NICULESCU, F. CHITEA University of Bucharest, Department of Geophysics, 6 Traian Vuia Street,
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationOn the magnetic structure of DyCro3
On the magnetic structure of DyCro3 B. Van Laar, Jacqueline B. A. A. Elemans To cite this version: B. Van Laar, Jacqueline B. A. A. Elemans. On the magnetic structure of DyCro3. Journal de Physique, 1971,
More informationMaximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale
Maximize the potential of seismic data in shale exploration and production Examples from the Barnett shale and the Eagle Ford shale Joanne Wang, Paradigm Duane Dopkin, Paradigm Summary To improve the success
More informationP11 A Methodology to Rescale the Vertical Coordinate into Geological Time Using a Geostatistical Analysis of Log Data
P A Methodology to Rescale the Vertical Coordinate into Geological Time Using a Geostatistical Analysis of Log Data M. Lefranc* (Andra / EMP), B. Beaudoin (Ecole des Mines de Paris), J.P. Chilès (Ecole
More informationAPPLICATION OF NUCLEAR MAGNETIC RESONANCE SOUNDINGS TO GROUNDWATERRESERVES MAPPING IN WEATHERED HARD-ROCK AQUIFERS (BRITTANY, FRANCE)
APPLICATION OF NUCLEAR MAGNETIC RESONANCE SOUNDINGS TO GROUNDWATERRESERVES MAPPING IN WEATHERED HARD-ROCK AQUIFERS (BRITTANY, FRANCE) Francis Mathleu', Robert Wvns 1, Jean-Micbel Baltassar', Anatoly Legchenko't',
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationActivity Submitted by Tim Schroeder, Bennington College,
Structural Analysis of a Hot Dry Rock Geothermal Energy System Activity Submitted by Tim Schroeder, Bennington College, tschroeder@bennington.edu Description: This project applies basic geologic skills
More informationQuantifying uncertainty of geological 3D layer models, constructed with a-priori
Quantifying uncertainty of geological 3D layer models, constructed with a-priori geological expertise Jan Gunnink, Denise Maljers 2 and Jan Hummelman 2, TNO Built Environment and Geosciences Geological
More informationHighly Robust Variogram Estimation 1. Marc G. Genton 2
Mathematical Geology, Vol. 30, No. 2, 1998 Highly Robust Variogram Estimation 1 Marc G. Genton 2 The classical variogram estimator proposed by Matheron is not robust against outliers in the data, nor is
More informationNEW GEOLOGIC GRIDS FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS
FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS EMMANUEL GRINGARTEN, BURC ARPAT, STANISLAS JAYR and JEAN- LAURENT MALLET Paradigm Houston, USA. ABSTRACT Geostatistical modeling of reservoir
More informationA STATIC 3D MODELING OF HYDROCARBONIC RESERVOIR WITH THE HELP OF RMS CASE study: THE SOUTH EAST ANTICLINE OF KHUZESTAN IRAN
:43-48 www.amiemt.megig.ir A STATIC 3D MODELING OF HYDROCARBONIC RESERVOIR WITH THE HELP OF RMS CASE study: THE SOUTH EAST ANTICLINE OF KHUZESTAN IRAN Hamid reza samadi 1,mohammad hadi Salehi 2 1 PH.D
More informationtechnical article Satinder Chopra 1*, Kurt J. Marfurt 2 and Ha T. Mai 2
first break volume 27, October 2009 technical article Using automatically generated 3D rose diagrams for correlation of seismic fracture lineaments with similar lineaments from attributes and well log
More informationPOTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE
POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE OCTOBER 2012 CONTENT Page I INTRODUCTION 1 II FIELD
More informationEnhanced Subsurface Interpolation by Geological Cross-Sections by SangGi Hwang, PaiChai University, Korea
Enhanced Subsurface Interpolation by Geological Cross-Sections by SangGi Hwang, PaiChai University, Korea Abstract Subsurface geological structures, such as bedding, fault planes and ore body, are disturbed
More informationHDR/EGS POTENTIAL OF THE VILKYCIAI AREA, WEST LITHUANIA
HDR/EGS POTENTIAL OF THE VILKYCIAI AREA, WEST LITHUANIA Prepared by Saulius Sliaupa, Dainius Michelevicius, Gediminas Motuza, Larisa Korabliova, Kirstina Ciuraite West Lithuania is situated in the central
More informationInverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study Hermans T. 1, Scheidt C. 2, Caers J. 2, Nguyen F. 1 1 University of Liege, Applied Geophysics
More informationSummary. Introduction
Xian Qiang*, Liu Yonglei, Lv Dong, An Haiting, He Xiaosong, Li Haiyin, Xiao Yong, Zhou Chenguang, Xu Jianyang, Dong Lei,and Mao Xianyu,BGP,CNPC Summary Although high density, FAZ(Full azimuth) seismic
More informationCOAL EXPLORATION AND GEOSTATISTICS. Moving from the 2003 Guidelines to the 2014 Guidelines for the Estimation and Classification of Coal Resources
COAL EXPLORATION AND GEOSTATISTICS Moving from the 2003 Guidelines to the 2014 Guidelines for the Estimation and Classification of Coal Resources Britmindo Australia Pty Ltd 35 Bottlebrush Street, Heathwood
More informationAnomaly effects of arrays for 3d geoelectrical resistivity imaging using orthogonal or parallel 2d profiles
African Journal of Environmental Science and Technology Vol. 4(7), pp. 446-454, July 2010 Available online at http://www.academicjournals.org/ajest ISSN 1991-637X 2010 Academic Journals Full Length Research
More informationGEOMETRIC MODELING OF A BRECCIA PIPE COMPARING FIVE APPROACHES (36 th APCOM, November 4-8, 2013, Brazil)
GEOMETRIC MODELING OF A BRECCIA PIPE COMPARING FIVE APPROACHES (36 th APCOM, November 4-8, 2013, Brazil) Serge Antoine Séguret; Research Engineer in Geostatistics; Ecole des Mines; Center of Geosciences
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationPhilippe LAGACHERIE, Jean Claude FAVROT et Marc VOLTZ. UFR science du sol INRA, 2 place Viala Montpellier cedex 1
Scientific registration n : 461 Symposium n : 17 Presentation : poster Using information from reference areas for mapping soils and their properties over physiographic regions Utilisations de secteurs
More informationPressure Grouting of Fractured Bedrock to Control Acid Mine Drainage
WATER RESOURCES AT RISK May 14-18, 1995 Denver American Institute of Hydrology Pressure Grouting of Fractured Bedrock to Control Acid Mine Drainage S. A. Effner, G. D. Vandersluis, and V. Straskraba Hydro-Geo
More informationP026 Outcrop-based reservoir modeling of a naturally fractured siliciclastic CO 2 sequestration site, Svalbard, Arctic Norway
P026 Outcrop-based reservoir modeling of a naturally fractured siliciclastic CO 2 sequestration site, Svalbard, Arctic Norway K. Senger (University of Bergen / CIPR), K. Ogata* (University Centre in Svalbard),
More information18 Single vertical fractures
18 Single vertical fractures 18.1 Introduction If a well intersects a single vertical fracture, the aquifer s unsteady drawdown response to pumping differs significantly from that predicted by the Theis
More informationTRUNCATED GAUSSIAN AND PLURIGAUSSIAN SIMULATIONS OF LITHOLOGICAL UNITS IN MANSA MINA DEPOSIT
TRUNCATED GAUSSIAN AND PLURIGAUSSIAN SIMULATIONS OF LITHOLOGICAL UNITS IN MANSA MINA DEPOSIT RODRIGO RIQUELME T, GAËLLE LE LOC H 2 and PEDRO CARRASCO C. CODELCO, Santiago, Chile 2 Geostatistics team, Geosciences
More informationGeostatistics: Kriging
Geostatistics: Kriging 8.10.2015 Konetekniikka 1, Otakaari 4, 150 10-12 Rangsima Sunila, D.Sc. Background What is Geostatitics Concepts Variogram: experimental, theoretical Anisotropy, Isotropy Lag, Sill,
More informationDFNModeler: AN EFFICIENT DISCRETE FRACTURE NETWORK MODELER
1 DFNModeler: AN EFFICIENT DISCRETE FRACTURE NETWORK MODELER 0709 Guohai Jin and Jack C. Pashin Geological Survey of Alabama, P. O. Box 869999, Tuscaloosa, AL 35486 ABSTRACT DFNModeler is a Windows application
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationD.NECHEM 1, M.HANNOUCHE 1, N. BOUGHERIRA 2 ABSTRACT. KEYWORDS: thermal stand - characterization of water - sources - chemical nature.
Seventeenth International Water Technology Conference, IWTC17 Istanbul, 5-7 November 213 CARACTEREZATION OF THERMAL SRPING WATERS DJEBEL SAFIA NORTH EASTERN ALGERIA - D.NECHEM 1, M.HANNOUCHE 1, N. BOUGHERIRA
More informationSupplement A. Effective Coulomb Stress Analysis We assessed changes in Coulomb stress beneath reservoirs by the general theory of decoupled linear poroelasticity [Biot, 1956; Roeloffs, 1988; Wang 2000;
More informationJohn E. Gale 1 and Eunjeong Seok 2
Field and Laboratory Coupled Fracture Deformation-Pore Pressure-Permeability Experiments That Provide Insight for Depressurization of Fractured Rock Slopes John E. Gale 1 and Eunjeong Seok 2 1 Fracflow
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More information11301 Reservoir Analogues Characterization by Means of GPR
11301 Reservoir Analogues Characterization by Means of GPR E. Forte* (University of Trieste) & M. Pipan (University of Trieste) SUMMARY The study of hydrocarbon reservoir analogues is increasing important
More information3D geostatistical porosity modelling: A case study at the Saint-Flavien CO 2 storage project
3D geostatistical porosity modelling: A case study at the Saint-Flavien CO 2 storage project Maxime Claprood Institut national de la recherche scientifique, Québec, Canada Earth Modelling 2013 October
More informationSINGULARITY ANALYSIS FOR IMAGE PROCESSING AND ANOMALY ENHANCEMENT. Qiuming CHENG
SINGULARITY ANALYSIS FOR IMAGE PROCESSING AND ANOMALY ENHANCEMENT Qiuming CHENG Department of Earth and Atmospheric Science, Department of Geography, York University, Toronto, 4700 Keele Street, Ont. M3J
More informationElectron beams magnetic field is not a result of electron motion but of their intrinsic magnetic moment.
Electron beams magnetic field is not a result of electron motion but of their intrinsic magnetic moment. (May 2014) Jean de Climont jeandeclimont@yahoo.ca Abstract : This paper proposes an experiment intended
More informationRESISTIVITY IMAGING IN EASTERN NEVADA USING THE AUDIOMAGNETOTELLURIC METHOD FOR HYDROGEOLOGIC FRAMEWORK STUDIES. Abstract.
RESISTIVITY IMAGING IN EASTERN NEVADA USING THE AUDIOMAGNETOTELLURIC METHOD FOR HYDROGEOLOGIC FRAMEWORK STUDIES Darcy K. McPhee, U.S. Geological Survey, Menlo Park, CA Louise Pellerin, Green Engineering,
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationMICRO LEVEL GEO-RESISTIVITY SURVEY THROUGH V.E.S. TEST FOR GROUNDWATER FEASIBILITY STUDY AND SELECTION OF BORE WELL SITES IN PIPILI BLOCK OF PURI
Review of Environment and Earth Sciences 2014 Vol.1, No.2, pp. 37-45 ISSN(e): 2313-8440 ISSN(p): 2409-2150 DOI: 10.18488/journal.80/2014.1.2/80.2.37.45 MICRO LEVEL GEO-RESISTIVITY SURVEY THROUGH V.E.S.
More informationThe Snap lake diamond deposit - mineable resource.
Title Page: Author: Designation: Affiliation: Address: Tel: Email: The Snap lake diamond deposit - mineable resource. Fanie Nel Senior Mineral Resource Analyst De Beers Consolidated Mines Mineral Resource
More informationAutomatic Determination of Uncertainty versus Data Density
Automatic Determination of Uncertainty versus Data Density Brandon Wilde and Clayton V. Deutsch It is useful to know how various measures of uncertainty respond to changes in data density. Calculating
More informationModélisation des écoulements dans un réseau discret de fractures par une approche continue
Réserves prolongées Raffinage propre Véicules économes Carburants diversifiés CO maîtrisé Modélisation des écoulements dans un réseau discret de fractures par une approce continue A. Fourno, C. Grenier,
More informationGeothermEx, Inc. GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION HOLE PROGRAM, KILAUEA EAST RIFT ZONE, HAWAII TASK 1 REPORT
(415) 527 9876 CABLE ADDRESS- GEOTHERMEX TELEX 709152 STEAM UD FAX (415) 527-8164 Geotherm Ex, Inc. RICHMOND. CALIFORNIA 94804-5829 GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION
More informationZONAL KRIGING. 5.1: Introduction and Previous Work CHAPTER 5
CHAPTER 5 ZONAL KRIGING Kriging and conditional indicator simulation are valuable tools for evaluating uncertainty in the subsurface, but are limited by the assumption of stationarity (i.e. the mean and
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationMEMORANDUM. Interpretation of Magnetic and Volterra-3DIP survey MC Claims
11966 95A Avenue, Delta, BC V4C 3W2 Canada Tel +1 (604) 582-1100 www.sjgeophysics.com MEMORANDUM Date: September 13, 2017 From: To: SUBJECT: E. Trent Pezzot Bonanza Mining Corporation Interpretation of
More informationGeoelectricity. ieso 2010
Geoelectricity ieso 2010 1 RESISTIVITY SURVEY AT VENETO VILLA GRITTI AT THE TOWN OF TREVISO (VENETO REGION) The survey was carried out to verify the underground presence of the fondations of a rustic building.
More informationChapter 8 Fetter, Applied Hydrology 4 th Edition, Geology of Groundwater Occurrence
Chapter 8 Fetter, Applied Hydrology 4 th Edition, 2001 Geology of Groundwater Occurrence Figure 8.42. Alluvial Valleys ground-water region. Fetter, Applied Hydrology 4 th Edition, 2001 Fetter, Applied
More informationRole of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation
Role of lithological layering on spatial variation of natural and induced fractures in hydraulic fracture stimulation Vincent Roche *, Department of Physics, University of Alberta, Edmonton roche@ualberta.ca
More informationMultiscale Seismic Signature of a Small Fault Zone in a Carbonate Reservoir: Relationships Between V P Imaging, Fault Zone Architecture and Cohesion*
Multiscale Seismic Signature of a Small Fault Zone in a Carbonate Reservoir: Relationships Between V P Imaging, Fault Zone Architecture and Cohesion* Pierre Jeanne 1,2, Yves Guglielmi 2, and Frédéric Cappa
More informationSimulation of Mineral Precipitation and Dissolution in the 5-km Deep Enhanced Geothermal Reservoir at Soultz-sous-Forêts, France
Proceedings World Geothermal Congress 2005 Antalya, Turkey, 24-29 April 2005 Simulation of Mineral Precipitation and Dissolution in the 5-km Deep Enhanced Geothermal Reservoir at Soultz-sous-Forêts, France
More informationPollution des sols et interaction avec la biosphère
Pollution des sols et interaction avec la biosphère Mazen SAAD Ecole Centrale de Nantes Laboratoire de Mathématiques Jean Leray FRANCE Mathématiques Pour La Planète La Rochelle 17 octobre 2013 Plan 1 Transport
More informationUnderstanding the Mechanical Behavior of Drilling-induced Tensile Fractures through Photoelasticity Lab Tests Conducted on Glass Cubes
Understanding the Mechanical Behavior of Drilling-induced Tensile Fractures through Photoelasticity Lab Tests Conducted on Glass Cubes Qing Jia, Douglas R. Schmitt, Randy Kofman and Xiwei Chen University
More informationSTRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION TESTS
Studia Geotechnica et Mechanica, Vol. 37, No. 1, 2015 DOI: 10.1515/sgem-2015-0003 STRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION
More information