Patent EA-005195-B1 (EA)
Title: Method for predicting the presence of
IPC: G01V9/00
Filed: 2003-04-01
Abstract: 1. A method for predicting mineral resources and determining spatial values of new mineral deposits of i different types, mainly areal, characterized in that an investigated Sp [belongs to] r region r of the Earth is divided into sites DeltaSp [belongs to] Sp and for each site multi-parameter evaluations of the Earth interior Z<R>(DeltaSp,T) are calculated for different i using a vector of mapped interior parameters F<R> = (Fj), j = 1, 2, ..., Mi, as a functional Q from the multi-parameter evaluations of probability the presence WH<R>(F<R>,T) and absence Wo<R>(F<R>,T)i of minerals in sites of type Z<R>(DeltaSp,T)i = Q[WH<R>(F<R>,T)i, Wo<R>(F<R>,T)], which, using an analog criterion omegak<R>[(DeltaSa,T)i], taking into account the interior in the DeltaSa sites, located on Si area compare for correspondence the interior DeltaSp with that in sites of standard analogs DeltaSa and if the correspondence provisions are met Z<R>(DeltaSp,T)i [belongs to] omegak<R>[Z<R>(DeltaSa,T)i] assuming hypothesis about deposit presence with K analog criterion in the DeltaSp and spatial coordinates of local zones with predictable deposits on the Sp. area are calculated using coordinates of the DeltaSp sites, in which deposits were predicted. 2. The method of claim 1, characterized in that calculation of evaluations of the interior state is carried out only for informative parameters Fj, j = 1, ..., m number of which differs for i different types of mineral deposits. 3. The method of claim 1, characterized in that calculation of evaluations of the interior state in the DeltaSp sites is carried out using models of the multi-parameter evaluations of WH<R>(F,T)i, Wo<R>(F,T)i, precalculated for the DeltaSa sites, belonging to the Si areas discovered on the step of regional adaptation in the region using r standard analogs of deposits. 4. The method of claim 3, characterized in that the F = (Fj), j = 1, 2, ..., Mi mapped parameters of the region interior are determined for each coordination mesh point, discrete interval Delta of which is determined by the DeltaSa sites area as for adapting and by the DeltaSp sites area as for predicting. 5. The method of claim 4, characterized in that values of the omegak<R>[Z<R>(DeltaSa,T)i] creterion analog are determined by intervals when calculating the distribution of of a number of standard analog sites of discovered deposits Nz(DeltaSa,T)i = f(K,T)i, being in different intervals of the K zone of the interior state. 6. The method of claim 5, characterized in that the intervals of the interior state, combining the state of different sites of the DeltaSa standard analogs are determined as a zone of crossing intervals of variation states within each DeltaSa site remembering its spatial position and observed therewithin a variation of Fj parameters. 7. The method of claim 6, characterized in that the variation zone within the DeltaSa standard analogs sites are approximated by a circle with DeltaR radius, the center of which is determined by the interior state in the central point of the DeltaSa site. 8. The method of claim 7, characterized in that the optimum DeltaR* radius is determined providing minimum deviation of contour of predicted deposit from actual contour, which is ensured at the adaptation step in the r region performing test analysis of areas with discovered deposits with known contours. 9. The method of claim 8, characterized in that the radius of interior state variation within the standard analog sites is assigned as DeltaRpes<DeltaR* for pessimistic prediction of new deposits with confidence probability more than 0.9 reflecting contours of predicted deposits and as DeltaRopt>DeltaR* for optimistic prediction with confidence probability at least 0.8. 10. The method of claim 9, characterized in that local bases of the Fj mapped parameters are formed for all or part j parameters having maps of 1:500000 scale and less, that allows to assign DeltaSp<DeltaSa, wherein DeltaSa = 1 km<2>. 11. A method for predicting mineral resources and determining coordinates thereof, comprising the steps of: a) forming a regional database of Fj<R> interior parameters, where F is parameter value, r is region index, j is parameter index, j = 1, 2, ..., M, representing in existing maps parameters of corresponding regions interior, the interior parameters are formed as F<R> = F1, F2, ..., FM = {Fj} vectors of parameters values in the DeltaSp sites of the coordinate scale for Sp area of corresponding r region; b) forming database about discovered different types of mineral deposits on the territory of corresponding regions areas Sai, where i = 1, 2, ..., N, indicating geologic age of T formations and about areas of proved deposit absence, wherein data about discovered deposits are formed in the form of a plurality of DeltaSai<H> sites and DeltaSai standard analogs correspondingly of deposit presence or absence; c) calculating one-parameter evaluations of probability of presence WH<R>(Fj,T)i and absence Wo<R>(Fj,T)i of deposits of i type, i = 1, 2, ..., N in formations of T geologic age in r region for DeltaSa, sites based on the regional database of the F<R> = F1, F2, ..., FM = {Fj} parameters and the DeltaSai<H> sites and DeltaSai standard analog of deposit presence/absence; d) forming an informative set of parameters F1...Fm for assigned geologic age and i type deposit and calculating multi-parameter evaluations of probabilities of WH<R>(F1...Fm,T)i presence and Wo<R>(F1...Fm,T)i absence of i type deposits in formations of T geologic age in r region by the the DeltaSai<H> sites and DeltaSai standard analog of deposit presence/absence; e) calculating multi-parameter evaluations of the interior state Z<R>(DeltaSa,T))i = Q<R>[(WH<R>(F1...Fm,T)i, Wo<R>(F1...Fm,T)i] by the DeltaSai<H> sites and DeltaSai standard analog of deposit presence/absence in formations of T geologic age in r region; f) determining local zones omegaH<R>, omegao<R> of variations of values in the DeltaSai<H> sites and DeltaSai standard analog of deposit presence/absence and carrying out optimization and classification of said zones by the K = 1, ..., k of analog criterion according to the number of Nz<R>(K,T) being within said local zones; g) calculating multi-parameter evaluations of the interior state Z<R>(DeltaSa,T))i = Q<R>[(WH<R>(F1...Fm,T)i, Wo<R>(F1...Fm,T)i] by the DeltaSai<H> sites based on calculated multi-parameter evaluations of probabilities of WH<R>(F1...Fm,T)i presence and Wo<R>(F1...Fm,T)i absence of i type deposits in formations of T geologic age in r region; h) taking hypothesis HH about the presence mineral deposit in the DeltaSp site if evaluation of the interior state Z<R>(DeltaSp,T)i belongs to a local zone omegaH<R>of value variations Z<R>(DeltaSp)i [belongs to] omegaH<R>(DeltaSa) setting the K level of analog criterion. 12. The method of claim 11, characterized in that the Earth terrain is used as regions on the area of which observe the change of the Fj<R> parameter values in all field of their values. 13. The method of claim 11, characterized in that geological, geophysical, structural, tectonic, morphologic, geochemical, thermal, seismic, magnetic, stratigraphic, gravimetric maps of the Fj<R> interior parameters are used as existing maps of the interior parameters. 14. The method of claim 11, characterized in that a number Nz<R>(K,T)i of standard-analog sites is taken as analog criterion when classifying local zones of interior state variations by levels of analog criterion according to the table: Levels of K analog criterion Number of analog criterion sites Nz<R>(K,T), related to the interior local zones in corresponding sites of deposit presence 1 Nz<R>(K,T)>=2000 2 1000<=NZ<R>(K,T)<2000 3 500<=NZ<R>(K,T)<1000 4 200<=NZ<R>(K,T)<500 5 50<=NZ<R>(K,T)<200 6 NZ<R>(K,T)<50