Petkov K. et al.

3D GEOSTATISTICAL MODEL OF THE ORE BODY …

50 years University of Mining and Geology “St. Ivan Rilski”
Annual, vol. 46, Part I, Geology and Geophysics, Sofia, 2003, pp.113-118

3D GEOSTATISTICAL MODEL OF THE ORE BODY

IN ELATSITE PORPHYRY COPPER DEPOSIT, PANAGYURISHTE ORE REGION

Kamen Petkov Popov, Kalin Ivanov Ruskov, Georgi Iliev Georgiev

University of Mining and Geology "St. Ivan Rilski", Sofia 1700, Bulgaria
[email protected], [email protected], [email protected]

ABSTRACT
The results of three-dimensional geostatistical modeling of copper ore body in Elatsite porphyry copper deposit are present in this work. The deposit is located in
NNW part of Elatsite-Chelopech Ore Field from Panagyurishte Ore Region, Bulgaria. Paleozoic granodioritic rocks are developed in northern part of the deposit and
Pre-Cambrian metamorphites are located in south. Upper Cretaceous granodioritic to quartzdioritic or monzodioritic porphyrites are intruded near and along the
contact between granodiorites and metamorphites.
Drillhole data from the open pit area, collected during preliminary and detail exploration of deposit and located between levels 1390 and 1045 are used for the model
development. Copper contents in composite samples from each drill for each level represent primary data. Three dimensional variogram analysis, based on spherical
model is used for determination of properties of spatial variability in data. Automatic routines of Least Squares method are applied for approximation of experimental
variogram values, which aims precise determination of variogram models. Cross-validation for the variogram model effectivity is done and it shows that chosen model
represent adequately natural variability of source data.
Smooth variation of copper content differences is typical for the Elatsite deposit, which is premise for precise determination and interpretation of variogram model.
Derived ore body anisotropy is oriented along fault structures, dykes and host rock contacts, confirming their main ore-controlling role in deposit.
Three-dimensional model of ore body is designed on basis of variogram analysis results, as ordinary block kriging interpolation is used. The deposit is modeled with
discrete blocks with dimension 15х15х15 m. The ore body is stockwork type, prolonged in NE direction. Digital terrain model is constructed for the modeling purposes
and it is used as boundary to exclude eroded part. The constructed 3D model is proximal and generalized as the ore body development is not studied in depth and
spatial characteristics of rock types are not separately analyzed.

GEOLOGICAL NOTES part of the deposit. The contact turns to the southeast (140о)
and then rapidly turns to the northeast (60о) as shown in fig.1.
The Elatsite porphyry copper deposit is located at the The strike of the contact in easternmost part is about 80о. The
northern slope of Chelopeshka Baba pick, about 55-60 km east general dip is 20-25о to the south near the surface, up to 45-
from Sofia City and about 6 km from Etropole town. It is formed 55о in depth. Tectonic brecciation is observed along the
within the frame of Elatsite-Chelopech Ore Field, as it is contact as a result of younger tectonic deformations.
located about 6 km NW from the center of Chelopech volcano
(Popov et al, 2001). 114.78 M.t. ore with 0.37 % copper
content were extracted during the period from 1981 to 1995.
Calculated in 01.01.2001 reserves are 139.57 M.t. ore with
contents of 0.32 % Cu, 0.108 g/t Au, 0.004 % Mo, as well as
the resources are 425.98 M.t. with 0.28 % copper content.

The Lower Paleozoic metamorphic rocks, Lower

Carboniferous rocks of the Vejen pluton and Late Cretaceous
rocks of Elatsite dyke-like intrusion with the associated dykes
are spread across the region of deposit (fig. 1). The Lower
Paleozoic metamorphic rocks are found along the southern
deposit’s boundary and phyllite, chlorit-sericitic schist and
quartz-sericitic schist represent them. Various interpretations
about their lithostratigraphic affiliation exist (Angelov et al.,
1995; Cheshitev et al., 1995; Antonov and Jelev, 2002). They
are alterated into hornfelse and knotted schist along the Figure 1.Geological map of the Elatsite deposit deposit (by
contact zone with the rocks of the Vejen pluton. Their main G. Georgiev, unpubl.)
orientation is from 90 to 110o with 15-45o dip to south. 1 - microdiorite, micro-gabbro-diorite, microgabbro, 2 –
granodiorite porphyrite, granite porphyrite, 3 - quartz-
The Vejen pluton is developed in northern part of the monzodiorite porphyrite, quartz-diorite porphyrite, 4 – Vejen
deposit, mainly represent by granodiorite and granite. Von pluton - granodiorite, granite, quartz-diorite, aplite, 5 -
Quadt et al. (2002) determine its absolute age as 314+4 Ma by hornfelse, knotted schist, schist, 6 – Kashana reverse slip -
206Pb/238U zircon method. The strike in the contact between the
overthrust, 7 - faults: а) Elatsite, b) Murgana, 8 – open pit
granodiorite and the horfelse is about 100о in the westernmost horizons, 9 – boundary of designed 3D model

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
113
Petkov K. et al. 3D GEOSTATISTICAL MODEL OF THE ORE BODY …
The Elatsite intrusion (The big dyke) consists of quartz- probably with Laramian age are established as well. It should
monzonitie intruded in the contact zone between be mentioned that Georgiev (unpubl.) determines radial-
metamorphites and Vejen pluton. The wide fracturing zone of concentric development of the jointing and small faults within
the Kashana reverse slip - overthrust obviously defines its the deposit.
position. This intrusion is represented by a long plate-like body
with over 3.3 km length, as its width varying from several The Elatsite deposit is spatially and genetically related to
meters up to 200-400 m. In the western part of the deposit the the Upper Cretaceous Elatsite intrusion and associated dykes
strike of that body is about 90о, in the central part - 110-120о, (Hadjyiski et al., 1970; Kalaydjiev et al., 1984; Dimitrov, 1988;
than it turns sharply to the northeast (50о), and finally turns to Popov et al., 2001). This is determined according to
the southeast. The intrusion’s dip is 20-45о to the south in circumstance that it is intensively affected together with the
upper parts, while in depth it’s probably increasing up to 55о- Paleozoic rocks by the post-ore hydrothermal alterations and
70о. by the ore mineralisation. The hydrothermal alterations are
A great number of dykes represent by quartz-monzodiorite represented mainly by K-silicate (Ksp+bi+qtz+il+pl+calc or
porphyrite, granodiorite porphyrite, granite-porphyries and Ksp+chl+qtz+il+ab), propyllitic (Ep+act+ab+bi+qtz+calc or
aplite associate with Elatsite intrusive. Microdioritic, Ep+chl+ab+bi+qtz) and sericitic (Il+qtz+py) alterations
micromonzonitic, diorite porphyritic and quartz-diorite (Strashimirov et al., 2002). The copper is primary economic
porphyritic dykes are also observed (Trashliev and Trashlieva, element in the Elatsite deposit and the Au is secondary one.
1964; Von Quadt et al., 2002). The molybdenite concentrate was extracted during some
exploitation periods as well. Several consecutively formed
The U/Pb zircon method dating determines the age of mineral associations represent the ore mineralization:
igneous rock bodies between 92.1+0.3 Ma – 91.84+0.31 Ma magnetite-bornite-chalcopyrite, quartz-pyrite-chalcopyrite,
(Von Quadt et al., 2002), which is very close to the quartz-molybdenite, quartz-pyrite, quartz-galena-sphalerite and
Cenomanian – Turonian boundary. Velichkova et al. (2001) quartz-calcite-zeolitic associations (Strashimirov et al., 2002).
determines 88-90 Ма age of the subvolcanic bodies near The development of the magnetite-bornite-chalcopyrite
Elatsite and Chelopech villages by 40Ar/39Ar analysis on biotite. mineralisation is quite irregular, as it includes native gold and
Besides, these intrusions are most probably formed by the end some minerals, containing elements from the PGE (Dragov
of Turonian or in the beginning of Coniasian ages, according to and Petrunov, 1996). The major economic interest takes the
geological data (Popov et al., 2001). quartz-pyrite-chalcopyrite association, which comprises the
main content of copper, gold (electrum) and molybdenum. The
Examining the structure of Elatsite deposit, Kalaydjiev et al. quartz-pyrite and quartz-galena-sphalerite associations are
(1984) suggest the development of a dense net of rectilinear developed mainly in outlying parts of the deposit.
faults with E-SE, N-NW and NE direction. The recent study
(Georgiev, unpubl.) does not confirm this assumption. The ore mineralisation forms a big stockwork, which is
Practically, there are three main fault structures, traced across developed within intensively jointed Paleozoic granodiorite and
the deposit: Elatsite fault, Murgana fault and Kashana reverse hornfelse, as well as within the rocks from Elatsite intrusive
slip - overthrust (fig. 1). The Elatsite fault posses a strike, and associated dykes (fig. 1). This ore body posses an ellipse-
varying from 85о to 125о. It lengthwisely crosses Elatsite like shape in plan, elongated to the northwest. Its length is
intrusion in the western part of the deposit, while to the east it about 1200 m and the width varies from 200 to 750 m. The
passes along the intrusion’s southern contact. Its dip is 70-85о area of its horizontal projection is about 0.616 km2 (Hadjyiski et
to the south. Movements, typical for right normal-slip fault are al., 1970). The stockwork is outcropped by the erosion on level
accomplished. This fault basicly coincides “The first Elatsite 1400 m (before the beginning of exploitation), as the ore
fault”, defined by Kalaydjiev et al. (1984). The Murgana fault mineralisation reaches level 550 m in depth. Its main axis
possesses a strike about 125-150о. It diagonally intersects the posses a general strike of 40-50о to south, as in general it
deposit as it’s dip is subvertical 80-90 to SW mainly. It is
o follows the contact between Lower Paleozoic rocks and
nominated by Kalaydjiev et al. (1984) as “Central Elatsite fault”. granodiorite of the Vejen pluton. The ore stockwork is not
Right-slip movements with amplitude about 200-250 m and outlined in depth, as most of drills do not reach its boundaries.
collaps of 50-200 m are accomplished. It is turning in It has continuous strike and dip. The ore body boundary within
southeastern part, where it is noted by Popov et al. (2001), as host rocks is gradual and it is defined by copper content of
the strike is 150-1550. The Kashana fault is registered for the 0.18 %. The transition is more sharply-outlined within
first time by Trashliev (1961), and later it is traced out and metamorphites, along their contact with granodiorite. The ore
described by Kouykin and Milanov (1970) and Kouykin et al. mineralisation is represented by multiple veinlets and
(1971). These authors determine its age as Austrian in first impregnations within intrusive rocks, while veinlets
mentioned paper, while in second paper they assign Laramian predominate in metamorphites. Small quartz-pyrite veins are
(or Post-Turonian) age to it. A thrust fault zone with width of rarely observed.
80-150 m, composed by anastomosing subordinate fault
surfaces is observed in southwestern part of the Elatsite
deposit, within the Late Paleozoic metamorphic rocks. It DATA USED
passes along the contact between the Early Paleozoic and
Triassic rocks to the west. The zone probably represents the Drillhole data from preliminary and detail exploration of
appearance of Kashana thrust in this part of the deposit. The Elatsite deposit (Hadjyiski et. al., 1970), between horizons
partial intrusion within this zone of Elatsite intrusive, as well as 1390 and 1045, are used in this study. The preliminary
the intrusion at Kashana and some other smaller bodies, point exploration is realised on 100x100 m grid, down to level 805.
to its Austrian age. Younger post-ore movements most The drill grid is chessmate condensed up to 71x71 m in the

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
114
 Petkov K. et al. 3D GEOSTATISTICAL MODEL OF THE ORE BODY …
central parts and up to 50x100 m in the periphery during detail defined as the relation of average differences between
exploration. Outlining boreholes placed at 200 to 300 m from measured values towards distances between samples:
the base exploration network are drilled around ore body. This r ∑ h (xi − xi +h )2
n

exploration network geometry is preserved down to level 1000, ()

γ h = i =1
2nh
,
as in deeper levels it is destored due to borehole inclination.
The deposit is studied in deep below level 802 with 12 drills, where xi-xi+h are nh of number differences between values
located in central part at intervals of 200 m. The drill core is (e.g. element’s contents), measured in samples at average
sampled by section samples with average length of 2.38 m and distance h between samples and situated along given
the copper content is determined by chemical analyses. direction. The main advantage of variogram analysis is
Composite core samples are used for the developing of average estimation of different natural features, as anisotropy,
geostatistical model, where every composite sample represent character and degree of correlation between neighbor samples

a) b)

c) d)
Figure 3. Experimental variograms and respective spherical models, determined for vertical sections in directions 0о (а), 45о (b),
90о (c) and 135о (d). Four variograms, calculated at dips 0о, 45о, 90о and 135о, are shown for every direction.

the copper content, which is averaged from section samples of

every drill at current level. Total number of 2819 composite
samples from 24 horizons is used. Copper contents in these
samples posses asymmetric statistical distribution shown on
figure 2, as the statistical parameters describing their
distribution low are represented below in table 2.

VARIOGRAM ANALYSIS

The variogram analysis is an essential part from the

geostatistical modeling of deposits. Its aim is a determination
of the spatial variability nature for studied geological feature
(Matheron, 1967; Rendu, 1981). The variogram function is Figure 2.Histogram presenting copper distribution in used
composite core samples.

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
115
Petkov K. et al. 3D GEOSTATISTICAL MODEL OF THE ORE BODY …
at different distances, value of general variability, continuity, are presented in table 1 and figure 3, which illustrate existing
etc., which are typical for the object studied. After calculation of anisotropy of copper ore body. Geometrical anisotropy
representing elongated ore morphology is examined, as well
as zonal anisotropy, which is determined by ore body location
around the contact between granodiorites and metamorphic
rocks. 18 variograms with angular lag 20о (half-angle solution)
are used for anisotropy analysis in horizontal direction with aim
to describe the structure of copper spatial variability (fig. 4).
The variogram surface shows that smoothest changes of
copper contents are examined in NE-SW direction, along the
prolongation of ore body, while highest values of copper
contents dispersion are examined in NW-SE direction.
Interesting fact is that major anisotropy axis is oriented either in
average differences in separate directions, represented by NE-SW or NW-SE direction due to chosen automatic
experimental variograms, it is very important to select an approximation method, which can be explained with the strike

a) b)
Figure 4. Variogram structure in horizontal direction, determined by 18 variograms (half-angle solution), shown as experimental
variograms (a) and variogram surface (b).

appropriate theoretical model, to describe adequately natural variation of host rocks contact as well as orientation of major
structure in variability of data used. faults and dykes.

Cross-validation routine is performed for estimation of

averaged variogram model significance, which is applied later
The experimental variograms in plane and in four vertical in kriging interpolation modeling. This validation compares true
directions with azimuth 0о, 45о, 90о and 135о are calculated in copper contents measured in composite samples with their
this variogram modeling of Elatsite deposit. Sixteen kriging estimations by neighbor samples. Popular geostatistical
experimental variograms are calculated in general – four modeling software GSLIB (Deutsch and Journel, 1998) is used.
variograms at 0о, 45о, 90о and 135о directions in plane, as well The results are shown on table 2 and figure 5, where a high
as at four variograms with 0о, 45о, 90о and 135о slopes in every similarity between true values and estimates, with correlation
vertical direction. Variograms are calculated by three- of 0.86, is illustrated. Similarity between real and estimated
dimensional searching of neighbor samples, in following values confirms the good representation of natural copper
conditions: from 20 to 500 m distance between samples with variability by variogram model, which is obvious condition for
lag 20 m, horizontal and vertical bandwidth 200 m, horizontal precise reserve calculation trough kriging method. Smoothing
and vertical angle tolerance 22.5о. effect of kriging interpolation, espatialy for samples with higher
Experimental variograms are approximated with spherical contents, is illustrated on figure 5.
model, represented in general form as:
 3h h 3 
γ (h) = C  − 3  for h<a, Table 1. Spherical variogram model parameters, calculated for
 2a 2a  individual directions and vertical angles.
γ ( h) = C for h>a, Azimuth Dip Nugget effect Sill Range
where C and a are variogram’s sill and range respectively, 0 0 0.00592432 0.0233568 285
while h is distance between samples. 0 45 0.000786367 0.0226083 180
The variogram analysis is completely done by UNCERT 0 90 0.00535023 0.0195682 305
0 135 0.00332431 0.0234292 190
software (Wingle, et al., 1997), as the approximation of
45 0 0.00415641 0.0215626 270
theoretical spherical model is automatically calculated by least 45 45 0.00273686 0.0203063 220
squares method. Resulting variogram models and parameters 45 90 0.00765853 0.0161969 300

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
116
 Petkov K. et al. 3D GEOSTATISTICAL MODEL OF THE ORE BODY …
45 135 0.0028656 0.0198173 170 represents central co-ordinates, estimated copper contents
90 0 0.00543918 0.0249688 315 and the kriging estimation variance for each block. Detailed
90 45 0.00573138 0.0249761 295 analysis of variability and reserves calculation conditions for
90 90 0.00277493 0.0210528 275 separate host rocks and ore types is presented by Todorov et
90 135 0.00713603 0.0196509 225
al. (2002).
135 0 0.00709468 0.0233214 295
135 45 0.00673393 0.0247538 250 Specialized module to visualize three-dimensional model of
135 90 0.00408765 0.026182 495 copper ore body is developed in OpenDX software, which is
135 135 0.00609529 0.0213959 285 based on IBM Visualization Data Explorer software (IBM,
Average: 0.00487 0.0221 1997). The ore body morphology is illustrated by individual
isosurfaces on figure 6, as the relief surface is also shown for
Table 2. Statistical parameters from cross-validation, higher reality. Constructed model could be used for reserves
describing distribution lows of copper contents in composite calculation for separate horizons as well, if the geological-
core samples (Cu), their kriging estimations (Cu*), estimation economic and mining-technical conditions in the deposit would
variance (SCu*) and the differences between estimations and be taken into account.
true contents (Cu* - Cu).
Cu Cu* SCu* Cu* - Cu
Mean. 0.282811 0.283932 0.002157 0.001121 CONCLUSIONS
Median 0.26 0.27 0.002 0.009
Minimum 0.007 0.01 0.001 -0.818 Designed three-dimensional model describes the spatial
Maximum 1.37 0.788 0.024 0.324
features of copper contents distribution and ore body
Variance 0.028384 0.019908 4.16E-07 0.007325
Std. Deviation 0.168476 0.141096 0.000645 0.085585
morphology in Elatsite deposit. The variogram analysis shows
Skewness 1.163632 0.568922 14.986 -1.5891 smooth copper contents variance, which is precondition for the
Kurtosis 2.70775 0.283673 469.9832 9.997335 precise model design. Invoked anisotropy posses a
complicated morphology due to complex geological setting and
different fault directions, as well as the variation of host rock
contacts. The major anisotropy axis is oriented either in
northeast or northwest direction, parallel or orthogonal to ore
body length respectively, in depends of applied method for
automated approximation of variogram model. Both geometric
and zonal components are observed in invoked anisotropy, as
the second is most obviously manifested in W-E direction. The
geometrical anisotropy component reflects the differences in
ore-forming process intensity at central and outern parts of
stockwork and it is main controlling factor for ore body
morphology. The zonal component is caused by rapid changes
in lithological features of host rocks. It could be concluded that
copper content variances are smoothest in NE-SW direction
generally, while highest variance values are observed in NW-
SE direction.
Figure 5. Correlation between true copper contents in core
samples (Cu) and their kriging estimations (Estimate).

KRIGING BLOCK MODEL

The geostatistical block model of copper ore body is

constructed for the volume with dimensions 1560x1080 m in
plan and 345 m in depth, between horizons 1390 and 1045.
Outline areas with lack of data are cutted, as the contour in
which the copper ore body is designed is shown on figure 1.
Digital elevation model is constructed by topographic map and
it is used to remove eroded parts. The studied volume is
disintegrated to small blocks with dimension 15x15x15 m,
which are estimated by block kriging method with 4x4x4
discretization points for every block. Kriging modeling is
performed with kt3d program from GSLIB software (Deutsch
and Journel, 1998), where 4 to 10 neighbor samples are used. Figure 6. 3D model of copper ore body represented by
Octant type searching within maximum radius is performed for isosurfaces on 0.25, 0.4 and 0.55 % Cu contents. Relief model
declustering of drill data. Applied average variogram model, is overlaid as transparent light-gray surface. View from
determined by variogram analysis (tabl. 1) posses nugget Northeast.
effect 0.00487 and sill 0.01723, as the anisotropy is described
with maximum, minimum and vertical range respectively 315
(azimuth 110о), 270 and 300 m. The resulting digital model

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
117
 Petkov K. et al. 3D GEOSTATISTICAL MODEL OF THE ORE BODY …
The contact between Vejen pluton and Lower Paleozoic metamorphites in this area. This contact controls the
metamorphites as well as the sub-parallel to contact Kashana development of higher metal contents as well. It should be
reverse slip – overthrust posses main magma-controlling and mentioned that the ore mineralisation within Paleozoic
ore-controlling role within Elatsite deposit area. The metamorphites is developed mainly within the zone, where
intersection area of this contact, Elatsite intrusive and Elatsite they are alterated into hornfelse. The knotted schist and
and Murgana faults defines particular ore stockwork location. phyllite located out of this zone posses screening role during
The last two faults obviously posses an ore-conducting role the ore forming, due to their plastic features.
within studied depth interval. The area of intersection between Designed model describes in general the upper levels of
separate structures is characterized by the highest deposit, as the highest contents are observed between
permeability, which is marked by the highest copper horizons 1400 and 1200. The cross-validation performed on
concentration in ore body. The ore stockwork stretch in variogram model shows a high correlation between estimated
northeastern direction is obviously determined by orientation of and true contents, which represents adequate describing of
the contact between Paleozoic granodiorites and real geological situation by the created model.

REFERENCES: Kouykin, S., Gercheva, Y., Milanov, L., Hristov, St. 1971.
Angelov, B., Iliev, K., Haydutov, I., Sapounov, I., Geological Setting of Stara Planina Mnt., between Zlatitsa
Choumachenko, P., Chounev, D.,. Tsankov, Ts., Marinova, and Troyan passage. Jub. Ann. Com. Geol. 179-200. (in
R., Rousanov, Iv., Yanev, Sl. 1995. - In: Explanatory Note Bulgarian).
to the Geological Map of Bulgaria in Scale 1:100 000, Popov, P., Raditchev, R., Dimovski, S. 2001. Geology and
Botevgrad map sheet. Geol. I-t BAS, “Geology and evolution of the Elatsite-Chelopech Porphyry Copper -
Geophysics”, Sofia, 117 p. (in Bulgarian). Massive Sulphide Ore Field. Ann. Univ. Min. Geol., Sofia,
Antonov, M., Jelev, V. 2002. Ductile Shear Zone and Brittle 43-44, I, 31-43.
Faults in the Southwestern Slope of Zlartitsa-Teteven Rendu, J.-M. 1981. An Introduction to Geostatistical Methods
Mountain (Central Bulgaria). Ann. Univ. Min. Geol., Sofia, of Mineral Evaluation. South African Institute of Mining and
45, 1, 13-20. Metallurgy, Printpack Ltd.
Cheshitev, G., Milanova, V., Sapounov, I., Choumachenko, P. Strashimirov, S., Petrunov, R., Kanazirski, M. 2002. Porphyry-
1995. - In: Explanatory Note to the Geological Map of copper mineralisation in the central Srednogorie zone,
Bulgaria in Scale 1:100 000, Teteven map sheet. Geol. I-t Bulgaria. Mineralium Deposita, 37: 587-598.
BAS, “Geology and Geophysics”, Sofia, 94 p. (in Todorov, J., Popov, K., Shanov, S., Boykova, A. 2002.
Bulgarian). Geological Conditions for a Correct Geostatistical
Deutsch, C., Journel, A. 1998. GSLIB: Geostatistical Software Evaluation: Example from the Elatsite Copper Deposit in
Library and User's Guide. Oxford University Press. Bulgaria. In: Geostatistics Rio 2000, Kluwer Academic
Dimitrov, S. 1988. Mineral Composition of Elatsite Plutogen- Publishers, 177-189.
Impregnated Copper-Molybdenum Deposit. Ann. Kom. Trashliev, S. 1961. On the Genesis and Age of Kashana Barite
Geol., vol. 28, 67-84. (in Bulgarian). Deposit, Pirdop area. Bull. Bulg. Geol. Soc., v. 22, 3, 245-
Dragov, P., Petrunov, R. 1996. Elatsite porphyry copper - 252. (in Bulgarian).
precious metal (Au and PGE) deposit. In: Plate tectonic Trashliev, S., Trashlieva, J. 1961. On the Young Intrusions in
aspects of the Alpine metallogeny in the Carpatho-Balkan Zlatitsa part of Stara Planina Mnt. In: Coll. Hon. Acad. Y. S.
region. Proceeding of the annual meeting - Sofia,. Unesco Yovchev. Publ. GUGOZN and NIGI, 545-558. (in
- IGCP Project № 356. Volum 1, 171-174. Bulgarian).
Hadjyiski, G., Angelkov, K., Nedkova, Ts., Tsvetkova, H. 1970. Velichkova, S., Handler, R., Neubauer, F. Ivanov, Z. 2001.
Report for the results from the geological exploration of Preliminary 40Ar/39Ar mineral ages from the Central
Elatsite copper-ore deposit – Etropole and realized during Srednogorie Zone, Bulgaria: Implications for Cretaceous
1959-1968 with copper reserves calculation for veinlet geodynamics. Rom. J. Mineral Deposits, Abstracts volume
impregnated ore, in state at 01.07.1968. Ministry of of ABCD – GEODE 2001 workshop, Vata Bai, Romania,
Environment and Water, National Geofond, I-744. (in 112-113.
Bulgarian). Von Quadt, A., Peytcheva, I., Kamenov, B., Fanger, L.,
IBM. 1997. IBM Visualization Data Explorer. Seventh Edition. Heinrich, C., Frank, M. 2002. The Elatsite porphyry copper
IBM Corporation, NY 10598-0704. deposit in the Panagyurishte ore district, Srednogorie zone,
Matheron, G. 1967. Traite de Geostatistique Applique. Bulgaria: U-Pb zircon geochronology and isotope-
Technic, Paris. geochemical investigation of magmatism and ore genesis.
Kalaydjiev, S., Hadjyiski, G., Angelkov, K. 1984. Structural (in print).
Conditions for Localization of Elatsite Porphyry Copper Wingle, W., Poeter, E., McKenna, S. 1997. UNCERT User’s
Deposit. Bull. Bulg. Geol. Soc., v. 45, 2, 189-196. (in Guide. Colorado School of Mines, Golden, Colorado
Bulgarian). 80401.
Kouykin, S., Milanov, L. 1970. Notes on the Geology of a part
of Zlatitsa Stara Planina. - Bull. Bulg. Geol. Soc. No 1, 120-
126 (in Bulgarian).

ANNUAL University of Mining and Geology “St. Ivan Rilski”, vol. 46 (2003), part I GEOLOGY AND GEOPHYSICS
118
View publication stats
 Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …
50 години Минно-геоложки университет “Св. Иван Рилски”
Годишник, том 46, свитък І, Геология и геофизика, София, 2003, стр.147-153

ТРИМЕРЕН ГЕОСТАТИСТИЧЕСКИ МОДЕЛ НА РУДНОТО ТЯЛО В МЕДНОПОРФИРНО

НАХОДИЩЕ ЕЛАЦИТЕ, ПАНАГЮРСКИ РУДЕН РАЙОН

Камен Попов, Калин Русков, Георги Георгиев

Минно-Геоложки Университет “Св. Иван Рилски”, София, 1700, България;

[email protected], [email protected], [email protected]

РЕЗЮМЕ
В настоящата работа са представени резултатите от тримерното геостатистическо моделиране на медното рудно тяло в меднопорфирно находище
Елаците. То е разположено в ССЗ част на Елашко-Члопешкото рудно поле от Панагюрския руден район. В района на находището са развити палеозойски
гранодиоритови скали от север и докамбрийски филити от юг, като в близост и по контакта им са внедрени горнокредни гранодиоритови до кварцдиоритови
порфирити или монцодиоритови порфирити.
За съставянето на модела са използвани данните от сондажите от детайлното проучване на находището и разположени между хоризонти 1390 и 1045
включително. Първичните данни представляват съдържанията на мед в групови проби от всеки сондаж за всеки експлоатационен хоризонт. За определяне
характера на пространствената променливост на медта е използван тримерен вариограмен анализ на базата на сферичен модел. За апроксимация на
експерименталните вариограмни стойности са използвани автоматизирани процедури по метода на най-малките квадрати, целящи прецизното извеждане
на вариограмните модели. Направена е проверка на адекватността на вариограмния модел, която показва че избрания модел добре описва природната
променливост на първичните данни.
Характерно за находище Елаците е плавното изменение на разликите в съдържанията на мед, което е предпоставка за прецизното съставяне и
интерпретация на вариограмния модел. Изведената анизотропия на рудното тяло е ориентирана по посоката на разпространение на разломните
структури, дайките и контактите между вместващите скали, което потвърждава тяхната главна рудоконтролираща роля в находището.
Въз основа на резултатите от вариограмния анализ е съставен тримерен модел на медното рудно тяло, чрез интерполация по метода на блоков крайгинг.
За целта находището е моделирано чрез дискретни блокове с размери 15х15х15 m. Рудното тяло е от щокверков тип, удължено в СИ посока. Също така,
за нуждите на модела е съставен и цифров модел на релефа, който е използван като ограничителна повърхност за изключване на отнесените от ерозията
части. Съставеният тримерен модел има приблизителен и обобщаващ характер, тъй като не е проследено развитието на рудното тяло в дълбочина и не са
анализирани пространствените характеристики по отделните типове вместващи скали.

ГЕОЛОЖКИ БЕЛЕЖКИ

Меднопорфирното находище Елаците е разположено по

северния склон на вр. Челопешка баба, на около 55-60 km
източно от гр. София и на около 6 km от гр. Етрополе. То е
формирано в рамките на Елашко–Челопешкото рудно
поле, като се намира на около 6 км северозападно от
центъра на Челопешкия вулкан (Popov et al, 2001). За
периода от 1981 до 1995 г. са добити 114 783.9 х.т. руда
със съдържание на мед 0.37 %. Към 01.01.2001 г. в
находището има запаси от 139 571 х.т. руда със
съдържание на мед 0.32 %, Au 0.108 g/t, Mo 0.004 % и
резерви от 425 980 х.т. руда със съдържание на мед
0.28 %.

В ареала на находище Елаците са развити старопалео-

зойски нискометаморфни скали, долнокарбонският Ве- Фигура 1. Геоложка карта на находище Елаците
женски плутон, къснокредният Елашки дайковиден (по Г. Георгиев, непубл.)
интрузив и многобройни дайки, асоцииращи с него (фиг. 1). 1 - микродиорити, микро-габродиорити, микрогабро,
Старопалеозойските нискометаморфни скали са развити в 2 – гранодиоритови порфирити, гранит порфири,
южната окрайнина на находището, като тяхната литостра- 3 - кварц-монцодиоритови порфирити, кварц-диоритови
тиграфска принадлежност се интерпретира различно порфирити, 4 - Веженски плутон - гранодиорити,
(Ангелов и др., 1995; Чешитев и др., 1995; Antonov and гранити, кварц-диорити, аплити, 5 - хорнфелзи, пъпчиви
Jelev, 2002). Те са представени от филити, хлорит- шисти, шисти, 6 -Кашански възсед-навлак, 7 - разломи:
серицитови шисти, кварц-серицитови шисти и др. В а) Елашки, b) Мургански, 8 - хоризонти на кариерата,
екзоконтакта на Веженския плутон те са контактно 9 - граници на съставения тримерен модел
променени в хорнфелзи и пъпчиви шисти. Тези скали са с
посока 90-1100, като затъват на юг с наклон 15-450.

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА

147
Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …
Веженският плутон заема северната част на находището, Отбелязан е от Калайджиев и др. (1984) като
като е представен предимно от гранодиорити и по-малко – Централноелашки разлом. По този разлом е осъществено
гранити. Според определението на Von Quadt et al. [2002] дясно отсядане с амплитуда 200-250 m и пропадане 50-
по 206Pb/238U цирконов метод, неговата абсолютна възраст 200 m. На югоизток, където е отбелязан от Popov et al.
е 314+4 Ma. Посоката на контакта между гранодиоритите и (2001), той завива, като посоката му е 150-1550.
хорнфелзите в най-западните части е около 100о, след Кашанският навлак е маркиран за първи път от Трашлиев
това завива на югоизток (140о), после рязко на североизток (1961), а след това е проследен и описан от Куйкин и
(60о) и в най-източните части е около 80о (фиг. 1). Като Миланов (1970) и Куйкин и др. (1971). В първата статия на
цяло той затъва на юг с наклон от 20-25о близо до последните автори той е приет за австрийски, а във
повърхността, а в дълбочина - до 45-55о. По контакта се втората – за ларамийски (или следтуронски). В
наблюдава тектонско брекчиране, което е резултат от по- югозападната част на находище Елаците, в
младите тектонски деформации. старопалеозойските метаморфозирани скали е развита
една възсед-навлачна разломна зона с ширина 80-150 m,
Елашкият интрузив (Голямата дайка) е с кварц-монцо- съставена от серия субпаралелни второстепенни
диоритов състав. Той е внедрен в зоната на контакта разломни повърхности. На запад тя преминава по контакта
между старопалеозойските метаморфити и Веженския между старопалеозойските и триаските скали. Най-
плутон. Очевидно неговата позиция в голяма степен се вероятно тя представлява проява на Кашанския навлак в
определя от широката зона на разломяване на Кашанския този участък. Обстоятелството, че в тази първоразрядна
възсед–навлак. Този интрузив представлява силно структура е внедрен частично Елашкия интрузив, както и
удължено пластинообразно тяло, чиято дължина интрузива при Кашана и други по-малки тела, указва, за
надхвърля 3.3 km, а ширината му варира от няколко метра неговата австрийска възраст. Същевременно по него са
(западно от находището) до 200-400 m. Посоката му осъществени и по-млади, следрудни възседни движения,
варира значително, като взападната част на находището е които най-вероятно са ларамийски. Следва да се
около 90о, в централната 110-120о, след това рязко завива отбележи, че Георгиев (непубл.) установява радиално-
на североизток (50о) а накрая също рязко завива на концентрично развитие на напукаността и малките
югоизток. В горните нива интрузивът затъва на юг с наклон разломи в участъка на находището.
20-45о, а в дълбочина наклонът му вероятно достига 55о-
70о. Находище Елаците пространствено и генетично е свър-
зано с горнокредния Елашки интрузив и асоцииращите с
С Елашкия интрузив асоциират редица дайки от кварц- него дайки (Хаджийски и др., 1970; Калайджиев и др., 1984;
монцодиоритови порфирити, гранодиоритови порфирити, Димитров, 1988; Popov et al., 2001). Това се определя от
гранитови порфири и аплити. Освен това се срещат дайки обстоятeлството, че заедно с палеозойските скали той е
от микродиорити, микромонцодиорити, диоритови засегнат интензивно от пострудните хидротермални
порфирити и кварцдиоритови порфирити (Трашлиев и метасоматити и рудната минерализация. Хидротермал-
Трашлиева, 1964; Von Quadt et al., 2002). ните метасоматити са представени предимно от К-
силикатни (Ksp+bi+qtz+il+pl+calc или Ksp+chl+qtz+il+ab),
Датирането по U/Pb цирконов метод определя възрастта пропилитови (Ep+act+ab+bi+qtz+calc или Ep+chl+ab+bi+qtz)
на магмените тела в диапазона 92.1+0.3 Ma – и серицитови (Il+qtz+py) изменения (Strashimirov et al.,
91.84+0.31 Ma (Von Quadt et al., 2002), което отговаря 2002). Главен промишлен елемент в находище Елаците е
приблизително на границата между ценомана и турона. медта, а второстепенен – златото. В отделни периоди от
Velichkova et al. (2001), по 40Ar/39Ar анализ на биотита, експлоатацията е добиван и молибденитов концентрат.
определят възрастта на субвулканските тела при Елаците Рудната минерализация е представена от няколко
и Челопеч на 88-90 Ма. Същевременно, по геоложки последователно формирани минерални асоциации:
данни, тези интрузиви най-вероятно са формирани към магнетит-борнит-халкопиритова, кварц-пирит-халкопирито-
края на турона и началото на кониаса (Popov et al., 2001). ва, кварц-молибденитова, кварц-пиритова, кварц-галенит-
сфалеритова и кварц-калцит-зеолитова (Strashimirov et al.,
Изследвайки структурата на находище Елаците, 2002). Магнетит-борнит-халкопиритовата асоциация е
Калайджиев и др. (1984) предполагат развитието на гъста неравномерно развита, като в нея присъстват самородно
мрежа от праволинейни разломи с изток-югоизточна, злато и минерали, носители на елементи от платиноид-
север-северозападна и североизточна посока. Извършени- ната група (Dragov and Petrunov, 1996). Главно промишле-
те по-нови проучвания (Георгиев, непубл.) не потвържда- но значение има кварц-пирит-халкопиритовата асоциация,
ват тази картина. Практически, в участъка на находището която съдържа основното количество от медта, златото
се проследяват три главни разломни структури: Елашки (електрум) и молибдена. Кварц-пиритовата и кварц-
разлом, Мургански разлом и Кашански възсед-навлак (фиг. галенит-сфалеритовата асоциации са развити основно в
1). Елашкият разлом е с посока варираща от 85о до 125о. В окрайните части на находището.
западната част на находището той пресича надлъжно
Елашкия интрузив, а на изток преминава по неговия южен Рудната минерализация изгражда един голям щокверк,
контакт. Той затъва на юг с наклон 70-85о, като по него са който е развит в интензивно напуканите палеозойски
реализирани десни разсед-отседни движения. В общи гранодиорити и хорнфелзи, както и в скалите на Елашкия
линии той се покрива с Първия Елашки разлом по интрузив и асоцииращите с него дайки (фиг. 1). В хоризон-
Калайджиев и др. (1984). Мурганският разлом е с посока тален план този щокверк е с елипсовидна форма, като е
125-150о, като разсича диагонално находището. Наклонът удължен в североизточно направление. Дължината му е
му е субвертикален 80-90о и затъва предимно на югозапад. около 1200 m, а ширината му варира от 200 до 750 m като

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА
148
 Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …
площта на хоризонталната му проекция е около 0.616 km2 ВАРИОГРАМЕН АНАЛИЗ
(Хаджииски и др., 1970). На повърхността рудния щокверк
е разкрит от ерозията на кота 1400 m (преди началото на Вариограмният анализ е съществена част от геостатис-
експлоатацията), като в дълбочина рудната минерали- тическото моделиране на находищата. Неговата цел е
зация достига до кота 550 m. Оста му е с генерален наклон изявата на закономерностите в пространствената промен-
40-50о на юг, като в основни линии следва контакта между ливост на изучавания геоложки показател (Matheron, 1967;
скалите на стария палеозой и гранодиоритите на Rendu, 1981). Вариограмата се дефинира от зависимостта
Веженския плутон. В дълбочина рудният щокверк не е на средните разлики между измерваните стойности спрямо
окунторен, тъй като повечето сондажи не достигат разстоянията между пробите:
неговите граници. Той е издържан по посока и наклон.
Границата на рудното тяло с вместващите скали е (xi − xi+h )2 ,
() ∑
nh
r
постепенна, като се определя по съдържание на мед γ h = i =1 (1)
0.18 %. Преходът е значително по-бърз в метаморфитите, 2nh
с отдалечаване от техния контакт с гранодиоритите. В
интрузивните скали рудната минерализация е представена където xi-xi+h са nh на брой разлики между стойностите
от многочислени жилки и впръслеци, докато в метамор- (напр. съдържания на полезен компонент), измерени в
фитите – предимно от жилки. Рядко се срещат малки проби, отстоящи на разсояние h помежду си в дадено
кварц-пиритови жили. направление.

Основно предимство на вариограмният анализ е, че чрез

ДАННИ, ИЗПОЛЗВАНИ ПРИ МОДЕЛИРАНЕТО него се отчитат като средно различни природни свойства,
като анизотропия, характер и степен на зависимост между
В настоящото изследване са използвани сондажните съседните проби при различни разстояния между тях, ниво
данни от предварителното и детайлното проучване на на общата променливост, прекъснатост и др., които са
находище Елаците (Хаджийски и др., 1970), от хоризонт характерни за изследвания обект. След пресмятане на
1390 до 1045 включително. Предварителното проучване е средните разлики в различни направления, представени от
извършено по мрежа 100x100 m до кота 805. По време на експерименталните вариограми, особено важно е да се
детайлното проучване сондажната мрежа е шахматно подбере подходящ теоретичен модел, който адекватно да
сгъстена в централните части до 71x71 m, а в описва природната структура на променливостта в
периферията – до 50x100 m. Встрани от рудното тяло са използваните данни.
прокарани оконтурващи сондажи, разположени от 200 до
300 m от основната проучвателна мрежа. Така описаната При настоящото вариограмно моделиране на находище
геометрия на сондажната мрежа се запазва до кота 1000, Елаците са пресметнати експерименталните вариограми в
докато в по-ниските нива се размества, поради хоризонтално и в четири вертикални направления с
изкривяването на сондажите. Находището е проучено в азимути 0о, 45о, 90о и 135о. Общо са пресметнати 16
дълбочина, под кота 805, с 12 сондажа, разположени в експериментални вариограми - четири вариограми с
централните части през 200 m. Сондажите са опробвани азимути 0о, 45о, 90о и 135о в хоизонтално направление,
със секционни проби със средна дължина 2.38 m, като както и по четири вариограми под наклони 0о, 45о, 90о и
съдържанието на мед е определяно с химичен анализ. За 135о за всяко от вертикалните направления. Вариограмите
съставянето на геостатистическия модел са използвани са пресмятани чрез тримерно търсене на съседни проби,
групови проби от сондажите, като всяка групова проба при следните условия: разстояния между пробите от 20 до
представя съдържанието на мед, усреднено по 500 m през стъпка 20 m, хоризонтална и вертикална ивица
секционните проби от всеки сондаж за даден хоризонт. на търсене 200 m, хоризонтален и вертикален ъглов
Използвани са общо 2819 групови проби от 24-те хоризонта. толеранс на търсене 22.5о.
Съдържанията на мед в тези проби показват асиметрично За апроксимация на експерименталните вариограми е
статистическо разпределение, илюстрирано на фигура 2, а използван сферичен модел, задаван в общ вид:
статистическите параметри, описващи закона на
разпределение са представени по-долу в таблица 2.  3h h 3 
γ (h) = C  − 3  за h<a, (2)
 2a 2a 
γ ( h) = C за h>a,

където C и a са съответно прага и ранга на вариограмата,

а h е разстоянието между пробите.

Вариограмният анализ е извършен изцяло с програмата

UNCERT (Wingle, et al., 1997), като апроксимацията на
теоретичния сферичен модел е извършена
автоматизирано по метода на най-малките квадрати.
Получените вариограмни модели и параметри са
представени в таблица 1 и фигура 3, които илюстрират и
Фигура 2. Хистограма, представяща разпределението съществуващата анизотропия на медното рудно тяло.
на медта в използваните сондажни групови проби Установява се както геометрична анизотропия, отчитаща

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА

149
Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …

a) b)

c) d)

Фигура 3. Експериментални вариограми и съответните им сферични модели, изведени за вертикални направления с

ориентация 0о (а), 45о (b), 90о (c) и 135о (d). За всяко направление са показани четирите вариограми, пресмятани за
наклони 0о, 45о, 90о и 135о
удължената морфология на орудяването, така и зонална което е необходимо условие за прецизното пресмятане на
анизотропия, обусловена от разположението на запасите по метода на крайгинг. Фиг. 5 илюстрира и
орудяването около контакта между гранитоидните и изглаждащият ефект на крайгинг интерполацията, особено
метаморфните скали. За анализ на анизотропията в план за пробите с високи съдържания.
са построени 18 вариограми през стъпка 20о (half-angle
solution), характеризиращи структурата на пространстве- Таблица 1. Параметри на сферичните вариограмни
ната променливост на медта (фиг. 4). От вариограмната модели, изведени за отделните направления и наклони
повърхност се вижда, че най-плавни изменения на съдър- Азимут Наклон Ефект на Праг Ранг
жанията се наблюдават в ЮЗ-СИ направление, по посока включенията
на удължението на рудното тяло, докато в СЗ-ЮИ направ- 0 0 0.00592432 0.0233568 285
ление се наблюдават най-високи стойности на диспер- 0 45 0.000786367 0.0226083 180
сията на съдържанията. Интересен факт е, че в зависи- 0 90 0.00535023 0.0195682 305
0 135 0.00332431 0.0234292 190
мост от избрания метод за автоматична апроксимация на
45 0 0.00415641 0.0215626 270
вариограмния модел, главната ос на анизотропията се 45 45 0.00273686 0.0203063 220
ориентира или в ЮЗ-СИ или в СЗ-ЮИ направление, което 45 90 0.00765853 0.0161969 300
се обяснява с промяната на посоката на контакта между 45 135 0.0028656 0.0198173 170
вместващите скали и с ориентацията на по-главните 90 0 0.00543918 0.0249688 315
разломи и дайкови скали. 90 45 0.00573138 0.0249761 295
90 90 0.00277493 0.0210528 275
Проверката за адекватността на усреднения вари- 90 135 0.00713603 0.0196509 225
ограмен модел (cross-validation), използван по-късно при 135 0 0.00709468 0.0233214 295
135 45 0.00673393 0.0247538 250
крайгинг интерполацията, е извършена чрез сравняване на 135 90 0.00408765 0.026182 495
истинските измерени съдържания на мед в груповите 135 135 0.00609529 0.0213959 285
проби от сондажите и техните оценки по метода на Средно: 0.00487 0.0221
крайгинг. Използван е популярния пакет за геостатис-
тическо моделиране GSLIB (Deutsch and Journel, 1998).
Резултатите са представени в таблица 2 и фигура 5,
където се илюстрира доброто съответствие между
истинските и оценените стойности с коефициент на
корелация 0.86. Сходството между реалните и оценените
съдържания показва, че вариограмния модел добре
описва съществуващата природна променливост на медта,

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА
150
 Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …

a) b)

Фигура 4. Вариограмна структура в хоризонтално направление, изведена по 18 вариограми (half-angle solution),

представена във вид на експериментални вариограми (a) и вариограмна повърхност (b)

Таблица 2. Статистически параметри, характеризиращи дискретизация 4x4x4 точки за всеки блок. Крайгинг
законите на разпределение на съдържанията на мед в моделирането е извършено с програмата kt3d от GSLIB
сондажните проби (Cu), техните крайгинг оценки (Cu*), (Deutsch and Journel, 1998), като за всеки блок са
дисперсията на оценките (SCu*) и разликите между използвани от 4 до 10 съседни проби, търсени по октант
оценките и истинските съдържания (Cu* - Cu) при максимален радиус, за да се избегне клъстерирането
Cu Cu* SCu* Cu* - Cu на данните по сондажите. Зададеният усреднен сферичен
Средно-аритм. 0.282811 0.283932 0.002157 0.001121 модел, определен чрез вариограмния анализ (табл. 1) е с
Медиана 0.26 0.27 0.002 0.009 ефект на включенията 0.00487 и праг 0.01723, като
Минимум 0.007 0.01 0.001 -0.818 анизотропията е характеризирана с максимален, мини-
Максимум 1.37 0.788 0.024 0.324
мален и вертикален ранг съответно 315 (с азимут 110о),
Дисперсия 0.028384 0.019908 4.16E-07 0.007325
Станд. отклонение 0.168476 0.141096 0.000645 0.085585 270 и 300 m. Полученият цифров модел представя за
Асиметрия 1.163632 0.568922 14.986 -1.5891 всеки блок неговите координати, оцененото съдържание на
Ексцес 2.70775 0.283673 469.9832 9.997335 мед и дисперсията на крайгинг оценката. По-детайлен
анализ на променливостта и условията за пресмятане на
запасите по отделни блокове вместващи скали и типове
руда е представен от Todorov et al. (2002).

За визуализирането на получения тримерен модел на

медното рудно тяло е създаден специализиран модул,
разработен чрез програмата OpenDX, която е базирана на
IBM Visualization Data Explorer (IBM, 1997). На фигура 6 е
илюстрирана морфологията на рудното тяло, визуализи-
рана чрез отделни изоповърхнини, като за по-голяма
реалност е показана и повърхността на релефа. Ако бъдат
отчетени геолого-икономическите и минно-техническите
условия в находището, съставеният модел може да бъде
използван и за пресмятане на запасите по отделните
хоризонти.
Фигура 5. Взаимовръзка между реалните съдържания на
мед в сондажните проби (Cu) и техните крайгинг оценки
(Estimate)

БЛОКОВ МОДЕЛ ПО МЕТОДА НА КАРЙГИНГ

Блоковият геостатистически модел на медното рудно

тяло е конструиран за обем с размери 1560x1080 m в план
и 345 m в дълбочина, между хоризонти 1390 и 1045.
Окрайните участъци, в които липсват достатъчно данни са
изрязани, като контурът, в който е моделирано рудното
тяло е показан на фигура 1. По топографска карта е Фигура. 6. Тримерен модел на медното рудно тяло,
създаден цифров модел на релефа, който е изполван за представено чрез изоповърхности по съдържания 0.25,
изключване на ерозиралите части от модела. Изслед- 0.4 и 0.55 % Cu. Моделът на релефа е наложен като
ваният обем е разбит на блокове с размери 15x15x15 m, прозрачна светлосива повърхност. Поглед от
които са оценявани по метода на блоков крайгинг при североизток.

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА

151
Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …
ЗАКЛЮЧЕНИЕ ЛИТЕРАТУРА

Направеният тримерен модел изявява пространствените Ангелов, В., Илиев, К., Хайдутов, И., Сапунов, И.,
закономерности в разпределението на съдържанията на Чумаченко, П., Чунев, Д., Цанков, Ц., Маринова, Р.,
мед и морфологията на рудното тяло в находище Елаците. Русанов, Ив., Янев, Сл. 1995. Обяснителна записка към
Вариограмният анализ показва плавно изменение на геоложката карта на България М 1:100000, картен лист
съдържанията на мед, което е предпоставка за прецизното Ботевград. София, Геол. и-т БАН, Геология и
съставяне на модела. Изведената анизотропия има геофизика АД, 117 с.
сложна морфология, поради комплицираният геоложки Димитров, С. 1988. Минерален състав на
строеж и различните посоки на тектонските нарушения, плутогенновпръснатото медно-молибденово находище
както и изменението в посоката на контактите между Елаците. – Год. КГ, 28, 67-84.
вместващите скали. В зависимост от избрания метод на Калайджиев, С., Хаджийски, Г., Ангелков, К. 1984.
автоматична апроксимация на вариограмния модел, Структурни условия за локализацията на медно-
дългата ос на анизотропията се ориентира в североиз- порфирното находище Елаците. – Сп. БГД, 45, 2, 189-
точно или северозападно направление, съответно пара- 196.
лелно или напречно на удължението на рудното тяло. В Куйкин, С., Миланов, Л. 1970. Бележки за геоложкия
изведената анизотропия се наблюдава както геометрична, строеж на част от Златишка Стара планина. – Сп.
така и зонална анизотропия, като последната се изявява БГД, 31, 1, 120-126.
най-ясно в И-З направление. Геометричната компонента в Куйкин, С., Герчева, Я., Миланов, Л., Христов, Ст. 1971.
анизотропията отразява разликите в интензитета в Геологичен строеж на Стара планина между Златишкия
процесите на рудоотлагане в централните и окрайните и Троянския проход. – Юбил. Год. КГ, 179-200.
части на щокверка и е основният фактор, влияещ върху Трашлиев, С. 1961. Върху генезиса и възрастта на
морфологията на рудното тяло. Зоналната компонента се баритовото находище Кашана, Пирдопско. – Сп. БГД,
дължи на рязката смяна на литоложките особености на 22, 3, 245-252.
вместващите скали. Може да се каже, че като цяло Трашлиев, С., Трашлиева, Ж. 1964. Върху младите
измененията в съдържанията на медта са най-плавни в интрузии в Златишка Стара планина. – В: Сб. в чест
СИ-ЮЗ направление, докато в СЗ-ЮИ направление се на акад. Й. С. Йовчев. Изд. ГУГОЗН и НИГИ, 545-558.
наблюдават най-големи стойности на променливостта. Хаджийски, Г., Ангелков, К., Недкова, Ц., Цветкова, Х. 1970.
Доклад за резултатите от геоложките проучвания на
Основна магмоконтролираща и рудоконтролираща роля медно-рудното месторождение “Елаците” – Етрополе,
в ареала на находище Елаците имат котакта между извършени през 1959-1968 г. с изчисление на запасите
Веженския плутон със старопалеозойските метаморфити и от мед на прожилково впръсната руда, по състояние
субпаралелния на контакта Кашански възсед-навлак. към 01.07.1968 г. Министерство на Околната Среда и
Конкретното местоположение на рудния щокверк се Водите, Национален Геофонд, I-744.
определя от участъка на пресичане между този контакт, Чешитев, Г., Миланова, В., Сапунов, И., Чумаченко, П.
Елашкия интрузив, Елашкия и Мурганския разломи. 1995. Обяснителна записка към геоложката карта на
Последните два разлома очевидно имат рудопроводяща България М. 1:100 000, картен лист Тетевен. София,
роля в изследвания дълбочинен интервал. Именно Геол. и-т БАН, Геология и геофизика АД, 94 с.
участъка на отбелязаното пресичане на различни струк-
турни елементи се отличава с най-висока проницаемост, Antonov, M., Jelev, V. 2002. Ductile Shear Zone and Brittle
което се маркира от най-голямата концентрация на мед в Faults in the Southwestern Slope of Zlartitsa-Teteven
рудното тяло. Удължението на рудния щокверк в северо- Mountain (Central Bulgaria). – Ann. Univ. Min. Geol.,
източна посока явно се определя от ориентацията на Sofia, 45, 1, 13-20.
контакта между палеозойските гранодиорити и метамор- Deutsch, C., Journel, A. 1998. GSLIB: Geostatistical Software
фити в този участък. Този контакт също контролира Library and User's Guide. Oxford University Press.
развитието на по-високи съдържания на метал. Следва да Dragov, P., Petrunov, R. 1996. Elatsite porfiry copper –
се отбележи, че в палеозойските метаморфити рудната
precious metal (Au and PGE) deposit. In: Plate tectonic
минерализация е развита основно в зоната където те са
aspects of the Alpine metalogeny in the Carpatho-Balkan
променени в хорнфелзи. Пъпчивите шисти и филитите
region. Proceeding of the Annual Meeting, Sofia, Unesco –
извън тази зона, поради тяхната пластичност, имат
IGCP Project № 356. Volum 1, pp. 171-174.
екранираща роля при рудоотлагането.
IBM. 1997. IBM Visualization Data Explorer. Seventh Ed., IBM
Corporation, New York.
Като цяло построеният модел характеризира горните
Matheron, G. 1967. Traite de Geostatistique Applique.
нива на находището, като най-високи съдържания се
Technic, Paris.
установяват между нива 1400 и 1200. Направената
Popov, P., Raditchev, R., Dimovski, S. 2001. Geology and
проверка на вариограмния модел отчита висока корелация
между оценените и реалните съдържания, което показва, evolution of the Elatsite-Chelopech Porphyry Copper –
че построеният модел отразява адекватно реалната Massive Sulphide Ore Field. – Ann. Univ. Min. Geol.,
геоложка ситуация. Sofia, 43-44, I, 31-43.
Rendu, J.-M. 1981. An Introduction to Geostatistical Methods
of Mineral Evaluation. South African Institute of Mining and
Metallurgy, Printpack Ltd.

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА
152
 Попов К. и др. ТРИМЕРЕН ГЕОСТАТИЧЕСКИ МОДЕЛ …
Strashimirov, S., Petrunov, R., Kanazirski, M. 2002. Porphyry- Volume of ABCD – GEODE 2001 Workshop, Vata Bai,
copper mineralization in the central Srednogorie zone, Romania, 112-113.
Bulgaria. –Mineralium Deposita, 37, 587-598. Von Quadt, A., Peytcheva, I., Kamenov, B., Fanger, L.,
Todorov, J., Popov, K., Shanov, S., Boykova, A. 2002. Heinrich, C., Frank, M. 2002. The Elatsite porphyry copper
Geological Conditions for a Correct Geostatistical deposit in the Panagyurishte ore district, Srednogorie zone,
Evaluation: Example from the Elatsite Copper Deposit in Bulgaria: U-Pb zircon geochronology and isotope-
Bulgaria. – In: Geostatistics Rio 2000. Kluwer Academic geochemical investigation of magmatism and ore genesis
Publishers, 177-189. (in print).
Velichkova, S., Handler, R., Neubauer, F. Ivanov, Z. 2001. Wingle, W., Poeter, E., McKenna, S. 1997. UNCERT User’s
Preliminary 40Ar/39Ar mineral ages from the Central Guide. Colorado School of Mines, Golden, Colorado.
Srednogorie Zone, Bulgaria: Implications for Cretaceous
geodynamics. – Rom. J. Mineral Deposits, Abstracts

Препоръчана за публикуване от
катедра “Геология и проучване на полезните изкопаеми”, ГПФ

ГОДИШНИК на Минно-геоложкия университет “Св. Иван Рилски”, том 46 (2003), свитък I, ГЕОЛОГИЯ И ГЕОФИЗИКА

153