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.     

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.     

第四纪环境研究新理论

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.     

Investigation on semi-direct and indirect climate effects of fossil fuel black carbon aerosol over China

A Regional Climate Chemistry Modeling System that employed empirical parameterizations of aerosol-cloud microphysics was applied to investigate the spatial distribution, radiative forcing (RF), and climate effects of black carbon (BC) over China. Results showed high levels of BC in Southwest, Central, and East China, with maximum surface concentrations, column burden, and optical depth (AOD) up to 14 μg?m^?3, 8 mg?m^?2, and 0.11, respectively. Black carbon was found to result in a positive RF at the top of the atmosphere (TOA) due to its direct effect while a negative RF due to its indirect effect. The regional-averaged direct and indirect RF of BC in China was about +0.81 and ?0.95 W?m^?2, respectively, leading to a net RF of ?0.15 W?m^?2 at the TOA. The BC indirect RF was larger than its direct RF in South China. Due to BC absorption of solar radiation, cloudiness was decreased by 1.33 %, further resulting in an increase of solar radiation and subsequently a surface warming over most parts of China, which was opposite to BC’s indirect effect. Further, the net effect of BC might cause a decrease of precipitation of ?7.39 % over China. Investigations also suggested large uncertainties and non-linearity in BC’s indirect effect on regional climate. Results suggested that: (a) changes in cloud cover might be more affected by BC’s direct effect, while changes in surface air temperature and precipitation might be influenced by BC’s indirect effect; and (b) BC second indirect effect might have more influence on cloud cover and water content compared to first indirect effect. This study highlighted a substantial role of BC on regional climate changes.     

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.     

Unique isoprene oxidation products demonstrate chlorine atom chemistry occurs in the Houston,Texas urban area

As part of the 2000 Texas Air Quality Study (TexAQS), we studied the isoprene oxidation process under ambient conditions to discern the presence of chlorine atom (Cl) chemistry in the Houston, Texas urban area. By measuring chloromethylbutenone (CMBO) and an isomer of chloromethylbutenal (CMBA), we clearly observed sixteen episodes of active Cl chemistry during the 24-day experiment. Estimated median Cl concentration during each of these episodes was between the detection limit of ~10² atoms cm⁻³ and 50 _{- 30} ^{+ 70} ×10⁴ 50_{ - 30}^{ + 70} \times {10^4} atoms cm⁻³. Cl concentration during all the episodes averaged 7.6 _{- 2.0} ^{+ 4.7} ×10⁴ 7.6_{ - 2.0}^{ + 4.7} \times {10^4} atoms cm⁻³ and thus amounted to less than 3% of the OH concentration during the same periods. During the episodes, the fraction of oxidation chemistry initiated by Cl ranged from 3–43% and was strongly dependent on the quantity and type of hydrocarbons present in the atmosphere. Because of its intermittent presence and low concentration, Cl is not a broadly influential oxidant in the Houston, Texas urban area.     

南海地区赤道辐合带的动力学与能量学特征分析

本文利用三角形区域法,对南海地区ITCZ的21个个例分析计算了有关动力学特征量随高度的分布,从而得出了平均的结果,这与蒋全荣和余志豪(1984)的分析结论基本一致。分析还发现在平均和分类情况下,南海地区ITCZ的高层(约300hpa)可能有一个相对的干燥层或“干盖”存在。能量学特征的计算表明,南海地区的ITCZ低层具有大量水汽、潜热和全能量的水平辐合,然后垂直向上输送并在高层向副热带地区水平辐散开去。因此,它在热带大气环流能量、水汽等方面的传输过程中起着重要的作用。      

深圳海岸复杂地形气流与湍流特征的数值模拟

利用三维非静力能量闭合（Ｅ－ε）的边界层模式，以深圳海岸复杂地形进行了实际模拟。结果给出了该地区海陆风情形下气流和湍流孤变化特征。在海风发展盛期，气流方向由大面积的水域和内陆的位置决定，不规则海岸线对局地气流影响不大，夜间陆风时，不规则海岸线对局地气流影响较大，湍流能量高值颁在陆地上空的不稳定层内，水面上湍能很小。模拟结果与实测结果上比较吻合。     

集合预报及其在中期天气预报中的应用

作者概述了集合预报的基本概念，集合预报中期天气预报中的意义与作用，现有的主要集合预报产品以及集合预报在美国美国国家气象中心和欧洲中期天气预报的业务运行。简要介绍了我国国家气象中心在动力延伸预报和集合预报方面的试验研究情况。