面向稀疏散布数据集的时空Kriging优化

时空Kriging法通过将变异函数向时空域进行扩展得到时空变异函数,有效地利用时空邻近的采样点综合进行插值,由于时空稀疏散布数据集具有单一时刻下样本点数量少以及时空分布不规律的特点,难以满足使用时空Kriging插值法的基本条件,导致插值精度不高,据此本文提出了优化方法：通过多时段叠置拟合空间变异函数的方法,综合利用时空邻域内的采样点以解决单一时刻下空间邻域内数量不足情况;控制时间变异对空间变异函数拟合的误差影响;采用积合式模型构建时空变异函数进行插值。最后使用Argo海温数据进行插值实验,在相同条件下与时空Kriging法以及时空权重法的交叉验证结果对比得出,该方法在保证拟合所需采样点数量要求的同时,有效削减了一般时空Kriging法中时间变异对空间变异函数拟合结果的干扰,插值结果的绝对误差均值从0.5降低至0.2以内,稳定性进一步增强,改善了时空Kriging法在稀疏散布数据条件下精度上的不足。     

Geostatistical estimation of chemical contamination in stream sediments: The case study of Vale das Gatas mine (northern Portugal)

Based on an environmental geochemistry case study carried out in the neighbourhood of a W–Sn abandoned mine, the pollution in stream sediments was modelled through a Global Contamination Index. Such an index permits one to summarize the combination of deleterious elements in a single variable, obtained by the projection of samples onto the first axis of a PCASD (Principal Components Analysis of Standardized Data) applied to the entire n × p matrix containing the available concentrations of p = 16 elements in the set of n = 220 collected samples.In order to provide a sound basis for a coherent planning of the remediation process which will be put in operation in the affected area, it is necessary to balance the costs of reclaiming with the probabilities of exceeding the upper limits accepted for concentrations of environmentally harmful elements in sediments. Given these limits, they are back-transformed in the index values, providing a practical threshold between ‘clean’ and ‘contaminated’ samples. On the other hand, the minimum dimension of the cell to be reclaimed is restrained by the selected remediation process to be applied in the affected area. Hence, to meet the constraints of such a remediation process, it is required to estimate the probabilities of exceeding the index threshold in technologically meaningful sub-areas. For this end, the Indicator Block Kriging technique was applied, producing a series of maps where sub-areas to be reclaimed can be spotted for different probability levels. These maps, on which the decision making remediation agency can rely for its cost-benefit analysis, take into account both the spatial structure of ‘clean’ vs. ‘contaminated’ samples and the constraints of the reclaiming process.     

Spatial Variability of Rock Depth in Bangalore Using Geostatistical, Neural Network and Support Vector Machine Models

Geospatial technology is increasing in demand for many applications in geosciences. Spatial variability of the bed/hard rock is vital for many applications in geotechnical and earthquake engineering problems such as design of deep foundations, site amplification, ground response studies, liquefaction, microzonation etc. In this paper, reduced level of rock at Bangalore, India is arrived from the 652 boreholes data in the area covering 220 km². In the context of prediction of reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth, Geostatistical model based on Ordinary Kriging technique, Artificial Neural Network (ANN) and Support Vector Machine (SVM) models have been developed. In Ordinary Kriging, the knowledge of the semi-variogram of the reduced level of rock from 652 points in Bangalore is used to predict the reduced level of rock at any point in the subsurface of the Bangalore, where field measurements are not available. A new type of cross-validation analysis developed proves the robustness of the Ordinary Kriging model. ANN model based on multi layer perceptrons (MLPs) that are trained with Levenberg–Marquardt backpropagation algorithm has been adopted to train the model with 90% of the data available. The SVM is a novel type of learning machine based on statistical learning theory, uses regression technique by introducing loss function has been used to predict the reduced level of rock from a large set of data. In this study, a comparative study of three numerical models to predict reduced level of rock has been presented and discussed.     

Triple diagram method for the prediction of wave height and period

Many formulations have been developed so far to predict the wave height and period from fetch length and wind blowing duration for a constant wind speed. This study aimed to predict wave parameters from fetch length and meteorological factors by using triple diagram methodology based on Kriging principles. Proposed model results were compared with Joint North Sea Wave Project (JONSWAP) model which is used so commonly in the ocean and coastal engineering studies. For the implementation of the methodology hourly wave and wind data were obtained from a buoy located in Lake Ontario. Numerical and graphical comparisons demonstrated that the proposed method outperforms the classical formulation.     

实验变差函数计算方法的研究与运用

本文系统、详尽地分析讨论了不同样品数据条件下实验变差函数的计算方法。在实验变差函数计算流程的设计和有效样品搜索域的确定方面提出了自己独特的见解。最后用一个矿山的控矿工程教据对文中介绍的方法进行了测试。结果表明,利用本文的方法,容易得到稳健、准确、信息量高的实验变差函数。     

Structure sismique du socle paléozoı̈que du bassin des Doukkala,Môle côtier,Maroc occidental. Indication en faveur de l'existence d'une phase éo-varisqueSeismic structure of the Doukkala basin,Palaeozoic basement,western Morocco: a hint for an Eovariscan fold-and-thrust belt

Seismic profiles and well data from the Doukkala basin unravel the structure of the Palaeozoic basement and suggest that this coastal zone of western Morocco was affected by a compressive phase during the Frasnian. This resulted in the formation of upright, plurikilometric folds associated with reverse faults (North Doukkala), and of asymmetrical folds associated with mostly west verging ramps (South Doukkala). Folding involved all pre-Upper Frasnian formations and caused partial or total hiatus of Upper Frasnian–Strunian strata. This event can be correlated with the orogenic phase reported from more internal domains of the Morocco Hercynian belt, where it is referred to as the ‘Bretonne’ or ‘Eovariscan’ phase. To cite this article: H. Echarfaoui et al., C. R. Geoscience 334 (2002) 13–20     

Caractérisation des paragenèses et des paléocirculations fluides dans l'indice d'or de Bleı̈da (Anti-Atlas,Maroc)Characterization of the paragenesis and the fluid palaeo-circulations in the gold showings of Bleida (Anti-Atlas,Morocco)

The gold showings at Bleida are hosted in Late Pan-African N50–80 °E quartz–hematite–chlorite 1 tension lenses that are related to the activity of major sinistral sub-east–west thrusts. Ores result from three superimposed stages of fluid migration. Gold occurs in microcracks offsetting the earlier minerals. Fluids evolved from COHN compositions with a saline component to boiling aqueous fluids. Pressure and temperature decreased from 50 MPa and 300 °C to less than 4 MPa and 150 °C. Thus, the gold showings at Bleida were formed in a typical geothermal (epithermal) setting, likely controlled by the Late Pan-African magmatism. To cite this article: A. Barakat et al., C. R. Geoscience 334 (2002) 35–41     

Tectonique transpressive en terre Adélie au Paléoprotérozoı̈que (Est Antarctique)Palaeoproterozoic transpression in Terre Adélie (East Antarctica)

The Palaeoproterozoic units of Terre Adélie show two types of structural domains associated with HT–LP metamorphic conditions: domes and NS–N340° striking vertical shear zones. Shear zones reflect dextral transpressive motions. Domes reflect sub-vertical shortening and principal stretching subparallel to shear zones. They could partly result from longitudinal flow coeval with transpression. Deformations are comparable to those described along the eastern and western boundaries of the Archean Gawler Craton (South-East Australia), which underlines the continuity between these two areas before opening of the Austral Ocean. To cite this article: A. Pelletier et al., C. R. Geoscience 334 (2002) 505–511.     

Les dykes basiques du massif ancien de l'Ourika (Atlas de Marrakech,Maroc) : géochimie et significationThe basic dykes of the Ourika old massif (High Atlas of Marrakech): Geochemistry and significance

The Precambrian massif of Ourika is crosscut by two systems of basic dykes, striking N40°E and N90–120°E. Using incompatible trace elements, the two systems form two distinct chemical groups, displaying a continental tholeiitic affinity. The composition variations between the two defined groups can be due to heterogeneities of mantle sources and to contamination, during the magma ascent, by the continental crust. The emplacement of these basic dykes, before the late-PIII formations, can be related to the Neoproterozoic distension generalised to the Anti-Atlas chain. To cite this article: A. Barakat et al., C. R. Geoscience 334 (2002) 827–833.