Numerical Model of Two-Phase Flow in Dissolvable Porous Media and Simulation of Reservoir Acidizing

Mathematical model of porous media dissolution coupled with two-phase flow is proposed. The model is based on the conception of dissolvable porous medium with deformable mass-variable porous skeleton. Model can be used for simulation of coupled chemo- and hydrogeomechanical processes which are difficult to examine experimentally. Acidizing of calcite oil reservoir is used as an example of the process. Water solution of hydrochloric acid and oil are two fluid phases of the model with several components. Dissolvable porous media is treated as deformable mass-variable solid phase. Change in mass of the solid phase is caused by hydrochloric acid dissolving the calcite part of the solid phase. Dissolution is supposed to be congruent; kinetics is governed by the Nernst law. Software for numerical solution of the model is developed. It uses AmgCL parallel library for high-performance computing in order to deal with large algebraic systems on the each time step of calculations. The library uses algebraic multigrid methods for preconditioning and parallel iterative solvers. NVidia CUDA framework is used as a backend to perform GPGPU calculations, because it proved to be faster than OpenCL framework on this problem. Numerical experiments on the basis of data set from real reservoirs are conducted with the developed software. Good correlation between field and calculated data is achieved. Numerical experiments for different configurations of heterogeneous layer are performed. Acidizing of layers with highly permeable conduit and with random distribution of permeability is modeled.     

Potential promoted productivity and spatial patterns of medium- and low-yield cropland land in China

With a continuously increasing population and better food consumption levels, improving the efficiency of arable land use and increasing its productivity have become fundamental strategies to meet the growing food security needs in China. A spatial distribution map of medium- and low-yield cropland is necessary to implement plans for cropland improvement. In this study, we developed a new method to identify high-, medium-, and low-yield cropland from Moderate Resolution Imaging Spectroradiometer (MODIS) data at a spatial resolution of 500 m. The method could be used to reflect the regional heterogeneity of cropland productivity because the classification standard was based on the regionalization of cropping systems in China. The results showed that the proportion of high-, medium-, and low-yield cropland in China was 21%, 39%, and 40%, respectively. About 75% of the low-yield cropland was located in hilly and mountainous areas, and about 53% of the high-yield cropland was located in plain areas. The five provinces with the largest area of high-yield cropland were all located in the Huang-Huai-Hai region, and the area amounted to 42% of the national high-yield cropland area. Meanwhile, the proportion of high-yield cropland was lower than 15% in Heilongjiang, Sichuan, and Inner Mongolia, which had the largest area allocated to cropland in China. If all the medium-yield cropland could be improved to the productive level of high-yield cropland and the low-yield cropland could be improved to the level of medium-yield cropland, the total productivity of the land would increase 19% and 24%, respectively.     

Data-Driven Predictive Modeling of Mineral Prospectivity Using Random Forests: A Case Study in Catanduanes Island (Philippines)

The Random Forests (RF) algorithm is a machine learning method that has recently been demonstrated as a viable technique for data-driven predictive modeling of mineral prospectivity, and thus, it is instructive to further examine its usefulness in this particular field. A case study was carried out using data from Catanduanes Island (Philippines) to investigate further (a) if RF modeling can be used for data-driven modeling of mineral prospectivity in areas with few (i.e., <20) mineral occurrences and (b) if RF modeling can handle predictor variables with missing values. We found that RF modeling outperforms evidential belief (EB) modeling of prospectivity for hydrothermal Au–Cu deposits in Catanduanes Island, where 17 hydrothermal Au–Cu prospects are known to exist. Moreover, just like EB modeling, RF modeling allows analysis of the spatial relationships between known prospects and individual layers of predictor data. Furthermore, RF modeling can handle missing values in predictor data through an RF-based imputation technique whereas in EB modeling, missing values are simply represented by maximum uncertainty. Therefore, the RF algorithm is a potentially useful method for data-driven predictive modeling of mineral prospectivity in regions with few (i.e., <20) occurrences of mineral deposits of the type sought. However, further testing of the method in other regions with few mineral occurrences is warranted to fully determine its usefulness in data-driven predictive modeling of mineral prospectivity.     

Capacity of soil loss control in the Loess Plateau based on soil erosion control degree

The capacity of soil and water conservation measures, defined as the maximum quantity of suitable soil and water conservation measures contained in a region, were determined for the Loess Plateau based on zones suitable for establishing terraced fields, forestland and grassland with the support of geographic information system (GIS) software. The minimum possible soil erosion modulus and actual soil erosion modulus in 2010 were calculated using the revised universal soil loss equation (RUSLE), and the ratio of the minimum possible soil erosion modulus under the capacity of soil and water conservation measures to the actual soil erosion modulus was defined as the soil erosion control degree. The control potential of soil erosion and water loss in the Loess Plateau was studied using this concept. Results showed that the actual soil erosion modulus was 3355 t?km^-2?a^-1, the minimum possible soil erosion modulus was 1921 t?km^-2?a^-1, and the soil erosion control degree was 0.57 (medium level) in the Loess Plateau in 2010. In terms of zoning, the control degree was relatively high in the river valley-plain area, soil-rocky mountainous area, and windy-sandy area, but relatively low in the soil-rocky hilly-forested area, hilly-gully area and plateau-gully area. The rate of erosion areas with a soil erosion modulus of less than 1000 t?km^-2?a^-1 increased from 50.48% to 57.71%, forest and grass coverage rose from 56.74% to 69.15%, rate of terraced fields increased from 4.36% to 19.03%, and per capita grain available rose from 418 kg?a^-1 to 459 kg?a^-1 under the capacity of soil and water conservation measures compared with actual conditions. These research results are of some guiding significance for soil and water loss control in the Loess Plateau.     

Degrading river network due to urbanization in Yangtze River Delta

Evolution of river systems under the background of human activities has been a heated topic among geographers and hydrologists. Spatial and temporal variations of river systems during the 1960s–2010s in the Yangtze River Delta (YRD) were investigated based on streams derived from the topographic maps in the 1960s, 1980s and 2010s. A list of indices, drainage density (Dd), water surface ratio (WSR), ratio of area to length of main streams (R), evolution coefficient of tributaries (K) and box dimension (D), were classified into three types (quantitative, structural, and complex indices) and used to quantify the variations of stream structure. Results showed that: (1) quantitative indices (Dd, WSR) presented decreasing trend in the past 50 years, and Dd in Wuchengxiyu, Hangjiahu and Yindongnan have decreased most, about 20%. Structurally, the Qinhuai River basin was characterized by significant upward R, and K value in Hangjiahu went down dramatically by 46.8% during the 1960s–2010s. Decreasing tendency in D was found dominating across the YRD, and decreasing magnitude in Wuchengxiyu and Hangjiahu peaks for 7.8% and 6.5%, respectively in the YRD. (2) Urbanization affected the spatial pattern of river system, and areas with high level of urbanization exhibited least Dd (2.18 km/km²), WSR (6.52%), K (2.64) and D (1.42), compared to moderate and low levels of urbanization. (3) Urbanization also affected the evolution of stream system. In the past 50 years, areas with high level of urbanization showed compelling decreasing tendency in quantitative (27.2% and 19.3%) and complex indices (4.9%) and trend of enlarging of main rivers (4.5% and 7.9% in periods of the 1960s–1980s and the 1980s–2010s). In the recent 30 years, areas with low level of urbanization were detected with significant downward trend in Dd and K. (4) Expanding of urban land, construction of hydraulic engineering and irrigation and water conservancy activities were the main means which degraded the river system in the YRD.     

Spatiotemporal distribution and historical evolution of polders in the Dongting Lake area,China

Dongting Lake is the largest lake in the middle reaches of the Yangtze River in China. For centuries, people inhabiting the Dongting Lake area have been reclaiming land and constructing dams for flood resistance, agricultural production, and rural settlement, forming geographical entities known as polders. In this study, the regional spatial distributions of polders in the Dongting Lake area in 1949, 1998, and 2013 were obtained using historical maps and modern remotely sensed data, revealing changes since the establishment of the People’s Republic of China. Nanxian County was then selected to demonstrate polder changes at the county level, because it has undergone the most dramatic changes in the area. Different polder change models for the Datonghu, Yule, and Renhe polders were analyzed for eight periods: 1644 (the early Qing Dynasty), 1911 (the late Qing Dynasty), 1930 (the Republic of China), 1949 (the People’s Republic of China), 1963, 1970, 1998, and 2013. Three resulting polder evolution models are: 1) reclaiming polders from lakes, 2) integrating polders by stream merging, and 3) abandoning polders for flood release. The polder evolution models demonstrate the wisdom of local people in using land resources according to the specific regional conditions. Throughout their long-term historical evolution, the spatial distribution of polders in the Dongting Lake area tended to be homogeneous, and the degree of human disturbance tended to be stable. However, a shift occurred, from pure polder area growth or removal to more comprehensive management and protection of the regional environment.     

山东沂沭河流域2000BC前后古文化兴衰的环境考古

Neolithic culture series in the Yishu River Basin developed in the order of Beixin culture-Dawenkou culture-Longshan culture-Yueshi culture. During the early and middle stage of Longshan culture （4600-4300 cal. yr BP）, the climate in the Yishu River Basin was warm and wet. Paddy-oriented agriculture planted paddy was very developed. The society was flourishing with great amount of archaeological sites. The cooling starting in 4200 cal. yr BP made the paddies shortfall in output or even no seeds were gathered. This situation intensified the discrepancy between population and resource. The scarcity in natural resource led to substantial decrease in population and subsequent drop in archaeological sites. About 4000 cal. yr BP Longshan culture was displaced by Yueshi culture which was relatively un- derdeveloped, simple and unsophisticated.     

第四纪环境研究新理论

On the basis of extensive survey to the Quatemary paleosols, soils and weathering layers in the vast northern China and the Loess Plateau, we found some unusually special phenomena of chemical components such as unsuccessive illuvium, twin illuviums, unusually thick illuvium and multi-illuviums etc. According to the analysis on the content of CaCO3 and the data of penetrating experiment, a new theory of illuvial depth or removal depth of CaCO3 in weathering and leaching belt was constructed, which indicates that the time of CaCO3 removing to the illuvial depth is very short,and the effect of time factor on illuvial depth can be ignored. The theory can be taken as a credible foundation for studying many geological and geographical problems in weathering and leaching belt.When the illuvial depth of CaCO3 iS bigger than the thickness of developing belt of soil or paleosol, it can be determined that the paleosol has turned into weathering crust. When the illuvial depth of CaCO3 is bigger than the thickness of paleosol, paleosol is leaching moderately acid soil. When two, three layers or unusually thick CaCO3 illuviums exist in the same weathering section or at the bottom of the same paleosol, there were two or more periods forming paleosol and corresponding climatic stages at that time. On the basis of the equation of relationship between mean annual rainfall (y) and illuvial depth of CaCO3 (x) (y = 305.5x 168.5) determined in the paper, mean annual precipitation during the development of paleosol can be calculated.     

An approach to adaptation to climate changes in Poland

Adopted by COP 10 (Dec 1/CP.10) and approved by the MOP1, the Buenos Aires programme of adaptation and response measures opens doors to intensify preparations for expected climate change. By this decision the COP, requested the SBSTA to develop a structured 5-year programme of work of the SBSTA on the scientific, technical and socio-economic aspects of impacts of, and vulnerability and adaptation to, climate change. Consequently, the COP, by its decision 2/CP.11, adopted the “Five-year programme of work of the Subsidiary Body for Scientific and Technological Advice on impacts, vulnerability and adaptation to climate change” Finally during COP12 this programme was approved as “Nairobi Work Programme on impacts, vulnerability and adaptation to climate change”. This programme has fundamental significance not only for developing countries, but also for industrialized nations in which some sectors of the or social life are particularly vulnerable to climate change, specifically, inter alia EIT countries and new EU Member States. Further development of this adaptation programme economy should contain steps that provide optimum economic and social effectiveness, risk management, identification of vulnerable sectors and gaps in knowledge, preparation of a list of policy options, including an analysis of cost effectiveness, selection of the most effective policies, and a preparedness implementation plan. In Poland the preliminary adaptation programme covered agriculture, water management, and coastal zone management. For the time being, gaps in knowledge and preparedness measures have been identified. An estimation of possible impact on these areas was based on chosen GCMs, and sea level rise IPCC scenarios. In conclusion, it was stated that the results achieved should be seen as a first step forward and a more comprehensive study is necessary to update the results and cover other sectors of the economy, such as health protection, spatial planning, ecosystems and forestry, and to develop specific guidelines and recommendations for policy-makers.     

Observation of hydrogen peroxide concentrations in a Japanese red pine forest

In order to study the concentrations of hydrogen peroxide (H₂O₂) and the factors controlling its concentrations, we monitored concentrations of H₂O₂ and other gases such as sulfur dioxide, ozone, and NO _x as well as meteorological factors such as air temperature, relative humidity, and wind direction/speed during eight measurement periods from 2000 to 2002 in a Japanese red pine forest in Japan. The H₂O₂ concentrations ranged from below 0.01 to 1.64 ppb, and analysis of the diurnal variation in H₂O₂ concentration showed high concentrations around noon, and low concentrations in the morning and late afternoon. The H₂O₂ concentrations were high in early summer, when O₃ concentration, temperature, and solar radiation were high, and were low in fall, when O₃ concentration, temperature, and solar radiation were low. We propose that O₃ concentration affects the production of H₂O₂ in the monitored region during the period under study, but that high H₂O₂ concentrations were sometimes caused by the transport of polluted air from urban regions. H₂O₂ concentrations decreased remarkably when SO₂ concentrations increased by transported volcanic emission on Miyake Island. In the absence of the effects of SO₂, H₂O₂ concentrations increased with increasing O₃ concentration and temperature.