20世纪末广东省三大城市大气环境质量综合评价——改进密切值法的扩展应用

大气环境质量综合评价是一个城市环境质量考核的重要指标，改进的密切值法是其中一种切实可行的评价方法。通过对1991-2000年10年间广东省三大城市——广州、深圳和珠海的大气环境质量进行综合评价，发现20世纪末，大气环境质量最好的是珠海市，最差的是广州市。与广东省平均水平相比、三大城市中。广州市污染程度高于全省平均水平。从空气环境质量的变化趋势来看，自1997年开始，三大城市大气环境质量趋于好转，为新世纪城市空气环境质量改善提供了良好的基础。     

Tectonic relation between northeastern China and the Korean peninsula revealed by interpretation of GRACE satellite gravity data

The major continental blocks in northeastern Asia are the North China block and the South China block, which have collided starting from the Korean peninsula. Geologic and geophysical interpretations reveal a well defined suture zone in northeastern China from Qinling through Dabie to Jiaodong. The discovery of high-pressure metamorphic rocks in the Hongseong area of the Korean peninsula, prominent evidence for the collision zone, indicates extension of the collision zone in northeastern China into the Korean peninsula. Interpretation of the GRACE satellite gravity dataset shows two prominent structural boundaries in the Yellow Sea. One extends from the Jiaodong Belt in eastern China to the Imjingang Belt in the Korean peninsula. The other extends from near Nanjing, eastern China, to Hongseong. Tectonic movement in or near the suture zone may be responsible for seismic activity in the western Korean peninsula and the development of the Yellow Sea sedimentary basin.     

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).     

Measurement of the total nitrogen in lake sediments： Comparison and its environmental significance

Nitrogen cycle is an important bio-geochemical process in the environment. Measurement of the total nitrogen （TN） is a routine experiment in agriculture, biology and environmental sciences. The Kjeldahl method （KM） and elemental analyzer method （EA） are both commonly used to determine TN. Total nitrogen by EA is the sum of nitrate （NO3）, nitrite （NO2）, organic nitrogen and ammonia. Total nitrogen by KM （TKN） is made up of both organic nitrogen and ammonia. A comparative study focused on the two methods is conducted by analysis of TN in 97 samples from the sediment sequence of Gouchi, a salt lake in North China. KM presents a higher degree of accuracy than EA with a standard deviation of 0.007 vs. 0.024. With the presence of nitrate and/or nitrite nitrogen, however, measurement by KM is considerably lower than that by EA. Therefore, for samples from lake sediment sequences or soils in North China, KM is inapplicable to determining TN because of usually high contents of nitrous salt. Despite the inconsistency, use of both methods to the same samples makes sense in reconstructions of climatic and environmental changes from lake sediments. In Lake Gouchi, the contents of nitrate and nitrite nitrogen vary from 1.40% in the lower part of the sequence to 14.77% in the uppermost part, suggesting a gradual evolution process from a fresh water lake to the present-day salt lake.     

Engineering geological studies and tunnel support design at Sulakyurt dam site, Turkey

This paper presents the results of engineering geological investigations and tunnel support design studies, carried out at the Sulakyurt dam site, northeast of Ankara, Turkey. The Sulakyurt dam will be used for flow control and water storage for irrigation projects. Studies were carried out both in the field and the laboratory. Field studies include engineering geological mapping, intensive discontinuity surveying, core drilling and sampling for laboratory testing. The diversion tunnel will be driven in rock mass, consisting of granite and diorite. Empirical, analytical and numerical methods were combined for safe tunnel design. Rock mass rating (RMR), Rock mass quality (Q) and Geological strength index (GSI) systems were used for empirical rock mass quality determination, site characterization and support design. The convergence–confinement method was used as analytical method and software called Phase², a 2D finite element program, was utilized as numerical method. According to the results acquired from the empirical, analytical and numerical methods, tunnel stability problems were expected in both granite and diorite rock masses. The support system, suggested by empirical methods, was applied and the performance of suggested support system was evaluated by means of numerical modelling. It was concluded that the suggested support systems were adequate, since after applying the suggested support system to granite and diorite, tunnel deformation and the yielded elements around the tunnel decreased significantly. Thus, it is suggested that for more reliable support design empirical, numerical and analytical methods should be combined.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Analysis of a rigid footing lying on three-layered soil using the finite difference method

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.     

Performance of Mortar and Chemical Grout Injection into Surrounding Soil When Slurry Pipe-jacking Method is Used

Small-diameter shallow tunnels are often being built by using the slurry pipe-jacking method. This system involves the pushing or thrusting of a drivage machine and concrete pipes into the ground. Chemical grout injection into the surrounding soil around the tunnel is carried out after the drivage and pushing processes are finished. The purpose of the chemical grout injection is to maintain permanent stability of the surrounding soil. However, the behavior of the chemical grouting material in the surrounding soil around the tunnel and the amount of optimum injection is not clearly understood. From these points of view, this paper discusses the performance of the chemical grouting material, when it is injected into the surrounding soil around the tunnel, by means of 2-D Eulerian–Lagrangian seepage analysis. Moreover, the effectiveness of the chemical grout injection was evaluated by using the non-linear finite element method. This investigation show when the range of the grouted zone is designed; it is necessary that the relationship between Young’s modulus of the soil/grouted zone and the confining stress be taken into consideration in order to establish effective, economical and safe chemical grout injection system. Understanding the performance of the seepage/dispersion behavior of the chemical grout and the characteristics of soil/ grouted zone is also important.     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.