黄土沟壑区小流域土壤pH值的空间分布及条件模拟

土壤pH值是影响土壤养分有效性和化学物质在土壤中行为的主要因素,研究土壤pH值的空间分布特征对于土壤养分管理和土壤污染预测具有重要意义。该文用地统计学方法研究了环境因素复杂的黄土高原小流域土壤pH值空间分布特征。结果表明,黄土沟壑区小流域土壤pH值具有球形—指数套合模型的空间结构特征,其空间异质性主要来源于流域内土地利用和土壤侵蚀等随机因素。与有机质协同的Kriging法能较好地对土壤pH值进行估值,其估值范围小于实测数据,估值误差来源于复杂的环境因素。序贯高斯条件模拟的土壤pH值范围与实测数据接近,模拟的平均值低于实测数据,模拟误差来源于模拟过程中独特的Kriging算法及高斯特性。

油田高含水期剩余油分布规律研究

After systemic investigation of the techniques，route lines and mechanisms about the remaining oil，the dynamic migration and congregation behavior of the remaining oil are discussed on base of interaction between flowing and enriching of water and oil.After the micro-scope modeling of the fluid flow in porous media and the changes in petrol-physical properties of the flowing system, the characters of fluid fields and the dynamic distribution of oil are discussed, among which the preference-flowing is focused on. Based on the preference-flowing in porous media, the concept of the preference-flowing channels is developed. According to above, heterogeneous distribution of water and oil in the field and dynamic mechanism of remaining oil are all obvious. media can be divided into three kinds, directional, stochastic, arbitrary porous media. The main research results are as following: 1. Treating the characteristic parameters such as permeability, porosity and wettability as regional parameter, the fluid field with high water-cut has been established by geostatistical method, among which the difference of flowing pores and the changes of its petrol-physical properties during flooding are studied. 2. The flow process of water and oil are recurrent in physical simulation experiments, in which the mechanisms and phenomena are caught and analyzed. Fluid flow mechanics in porous media with preference-flowing channels have been studied. 3. The mutual coupling between water and oil is induced and the mathematical evolution equations including this interaction were built. . 4. Through coupling effect between flowing water and oil, the dynamic migration and congregation behavior of remaining oil depend upon this coupling. 5. Coupling between water and oil act as driving force and trapping force for the remaining oil. The coupling model of thesis has been verified by simplified the numerical model and compared results with Ng35 oil reservoir in Gudao oil field, it has important theoretical and application values for improving precision of remaining oil and production performance prediction, and is a new method for studying the mechanics of remaining oil in channeled porous media has been established. Key words：flow field，high water-cut，coupling，dominant flow in porous media，remaining oil

Groundwater quality, vulnerability and potential assessment in Kobo Valley development project, Ethiopia

This study deals with investigating the groundwater quality for irrigation purpose, the vulnerability of the aquifer system to pollution and also the aquifer potential for sustainable water resources development in Kobo Valley development project. The groundwater quality is evaluated up on predicting the best possible distribution of hydrogeochemicals using geostatistical method and comparing them with the water quality guidelines given for the purpose of irrigation. The hydro geochemical parameters considered are SAR, EC, TDS, Cl-, Na+, Ca++, SO4 2- and HCO3 -. The spatial variability map reveals that these parameters falls under safe, moderate and severe or increasing problems. In order to present it clearly, the aggregated Water Quality Index (WQI) map is constructed using Weighted Arithmetic Mean method. It is found that Kobo-Gerbi sub basin is suffered from bad water quality for the irrigation purpose. Waja Golesha sub-basin has moderate and Hormat Golena is the better sub basin in terms of water quality. The groundwater vulnerability assessment of the study area is made using the GOD rating system. It is found that the whole area is experiencing moderate to high risk of vulnerability and it is a good warning for proper management of the resource. The high risks of vulnerability are noticed in Hormat Golena and Waja Golesha sub basins. The aquifer potential of the study area is obtained using weighted overlay analysis and 73.3% of the total area is a good site for future water well development. The rest 26.7% of the area is not considered as a good site for spotting groundwater wells. Most of this area fall under Kobo-Gerbi sub basin.

Assessing the uncertainty associated with intermittent rainfall fields

[1] In many practical situations where spatial rainfall estimates are needed, rainfall occurs as a spatially intermittent phenomenon. An efficient geostatistical method for rainfall estimation in the case of intermittency has previously been published and comprises the estimation of two independent components: a binary random function for modeling the intermittency and a continuous random function that models the rainfall inside the rainy areas. The final rainfall estimates are obtained as the product of the estimates of these two random functions. However the published approach does not contain a method for estimation of uncertainties. The contribution of this paper is the presentation of the indicator maximum likelihood estimator from which the local conditional distribution of the rainfall value at any location may be derived using an ensemble approach. From the conditional distribution, representations of uncertainty such as the estimation variance and confidence intervals can be obtained. An approximation to the variance can be calculated more simply by assuming rainfall intensity is independent of location within the rainy area. The methodology has been validated using simulated and real rainfall data sets. The results of these case studies show good agreement between predicted uncertainties and measured errors obtained from the validation data.

Geochemical characterization of heavy metal contaminated area using multivariate factorial kriging

This paper describes a geostatistical method, known as factorial kriging analysis, which is well suited for analyzing multivariate spatial information. The method involves multivariate variogram modeling, principal component analysis, and cokriging. It uses several separate correlation structures, each corresponding to a specific spatial scale, and yields a set of regionalized factors summarizing the main features of the data for each spatial scale. This method is applied to an area of high manganese-ore mining activity in Amapa State, North Brazil. Two scales of spatial variation (0.33 and 2.0 km) are identified and interpreted. The results indicate that, for the short-range structure, manganese, arsenic, iron, and cadmium are associated with human activities due to the mining work, while for the long-range structure, the high aluminum, selenium, copper, and lead concentrations, seem to be related to the natural environment. At each scale, the correlation structure is analyzed, and regionalized factors are estimated by cokriging and then mapped.

Geoelectrical mapping of contamination in the cemeteries: the case study in Piracicaba, São Paulo/Brazil

This paper presents the results of electrical resistivity methods in the area delineation that was potentially contaminated by liquefaction products, which are also called putrefactive liquids in Vila Rezende municipal cemetery, Piracicaba, So Paulo, Brazil. The results indicate a depth of water table between 3.1 and 5.1 m, with two groundwater direction flows, one to the SW and another to the SE. Due to the contamination plumes, which have the same groundwater direction flow, as well the conductive anomalies observed in the geoelectric sections, the contamination suspicions in the area were confirmed. The probable plume to the SE extends beyond the limits of the cemetery. The location of the conductive anomalies and the probable contamination plumes showed that the contamination is linked with the depth of the water table and the burial time. Mapping using the geostatistical method of ordinary kriging applied to the work drew structural characteristics of the regional phenomenon and spatial behavior of the electrical resistivity data, resulting in continued surfaces. Thus, this method has proved to be an important tool for mapping contamination plumes in cemeteries.

Estudo do nivel piezometrico por analise geoestatistica

160 groundwater bores in the Ribeirao Preto region were used to study the behaviour of the piezometric surface of the Botucatu Formation aquifer by geostatistical methodology. The semi-variogram showed the following results: nugget = 50 m2; sill = 450 m2 and range = 4.5 km. The piezometric surface, obtained by punctual kriging in a 100 x 100 knot grid, revealed two lowered zones: one in a southwest to north-northeastern trend, indicating the groundwater flow north toward the Rio Pardo; and the other in the central part due to the high concentration of holes in the urban region of Ribeirao Preto city. -from English summary

Suscetibilidade magnética e sua associação com atributos do solo e do relevo em uma topossequência de Latossolos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Interpolação geoespacial da evapotranspiração de referência (ETo) em regiões com escassez de dados: estudo de caso no Sul de Minas Gerais, Brasil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Geoelectrical mapping of contamination in the cemeteries: the case study in Piracicaba, Sao Paulo/Brazil

This paper presents the results of electrical resistivity methods in the area delineation that was potentially contaminated by liquefaction products, which are also called putrefactive liquids in Vila Rezende municipal cemetery, Piracicaba, So Paulo, Brazil. The results indicate a depth of water table between 3.1 and 5.1 m, with two groundwater direction flows, one to the SW and another to the SE. Due to the contamination plumes, which have the same groundwater direction flow, as well the conductive anomalies observed in the geoelectric sections, the contamination suspicions in the area were confirmed. The probable plume to the SE extends beyond the limits of the cemetery. The location of the conductive anomalies and the probable contamination plumes showed that the contamination is linked with the depth of the water table and the burial time. Mapping using the geostatistical method of ordinary kriging applied to the work drew structural characteristics of the regional phenomenon and spatial behavior of the electrical resistivity data, resulting in continued surfaces. Thus, this method has proved to be an important tool for mapping contamination plumes in cemeteries.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Global compilation of benthic data sets I

Approaches to quantify the organic carbon accumulation on a global scale generally do not consider the small-scale variability of sedimentary and oceanographic boundary conditions along continental margins. In this study, we present a new approach to regionalize the total organic carbon (TOC) content in surface sediments (<5 cm sediment depth). It is based on a compilation of more than 5500 single measurements from various sources. Global TOC distribution was determined by the application of a combined qualitative and quantitative-geostatistical method. Overall, 33 benthic TOC-based provinces were defined and used to process the global distribution pattern of the TOC content in surface sediments in a 1°x1° grid resolution. Regional dependencies of data points within each single province are expressed by modeled semi-variograms. Measured and estimated TOC values show good correlation, emphasizing the reasonable applicability of the method. The accumulation of organic carbon in marine surface sediments is a key parameter in the control of mineralization processes and the material exchange between the sediment and the ocean water. Our approach will help to improve global budgets of nutrient and carbon cycles.

Spatially stratified sampling using auxiliary information for geostatistical mapping

This paper presents a method of spatial sampling based on stratification by Local Moran’s I i calculated using auxiliary information. The sampling technique is compared to other design-based approaches including simple random sampling, systematic sampling on a regular grid, conditional Latin Hypercube sampling and stratified sampling based on auxiliary information, and is illustrated using two different spatial data sets. Each of the samples for the two data sets is interpolated using regression kriging to form a geostatistical map for their respective areas. The proposed technique is shown to be competitive in reproducing specific areas of interest with high accuracy.

A Geostatistical Study pf Alumina-Content in Ores of An Indian Iron-Ore Deposit

Mining and blending operations in the high grade iron ore deposit under study are performed to optimize recovery with minimal alumina content while maintaining required levels of other chemical component and a proper mix of ore types. In the present work the regionalisation of alumina in the ores has been studied independently and its effects on global and local recoverable tonnage as well as on alternatives of mining operations have been evaluated. The global tonnage recovery curves for blocks (20m x 20m x 12m) obtained by simulation closely approximated the curves obtained theoretically using a change of support under the discretised gaussian model. Variations in block size up to 80m x 20m x 12m did not affect the recovery as the horizontal dimensions of the blocks are small in relation to the range of the variogram. A comparison of the local tonnage recovery curves obtained through multiple conditional simulations made with that obtained by the method of uniform conditioning of block grades on an estimate of panel 100m x 100m x 12m panel grade reveals comparable results only in panels which have been well conditioned and possesing an ensemble simulation mean close to the ordinary kriged value for the panel. Study of simple alternative sequence of mining on the conditionally simulated deposit shows that concentration of mining operations simultaneously on a single bench enhances the fluctuation in alumina values of ore mined.

Nassau grouper (Epinephelus striatus) spawning aggregations: hydroacoustic surveys and geostatistical analysis

With the near extinction of many spawning aggregations of large grouper and snapper throughout the Caribbean, Gulf of Mexico, and tropical Atlantic, we need to provide baselines for their conservation. Thus, there is a critical need to develop techniques for rapidly assessing the remaining known (and unknown) aggregations. To this end we used mobile hydroacoustic surveys to estimate the density, spatial extent, and total abundance of a Nassau grouper spawning aggregation at Little Cayman Island, Cayman Islands, BWI. Hydroacoustic estimates of abundance, density, and spatial extent were similar on two sampling occasions. The location and approximate spatial extent of the Nassau grouper spawning aggregation near the shelf-break was corroborated by diver visual observations. Hydroacoustic density estimates were, overall, three-times higher than the average density observed by divers; however, we note that in some instances diver-estimated densities in localized areas were similar to hydroacoustic density estimates. The resolution of the hydroacoustic transects and geostatistical interpolation may have resulted in over-estimates in fish abundance, but still provided reasonable estimates of total spatial extent of the aggregation. Limitations in bottom time for scuba and visibility resulted in poor coverage of the entire Nassau grouper aggregation and low estimates of abundance when compared to hydroacoustic estimates. Although the majority of fish in the aggregation were well off bottom, fish that were sometimes in close proximity to the seafloor were not detected by the hydroacoustic survey. We conclude that diver observations of fish spawning aggregations are critical to interpretations of hydroacoustic surveys, and that hydroacoustic surveys provide a more accurate estimate of overall fish abundance and spatial extent than diver observations. Thus, hydroacoustics is an emerging technology that, when coupled with diver observations, provides a comprehensive survey method for monitoring spawning aggregations of fish.