Dynamic assessment of pollution risk of groundwater source area in Northern China

Based on the dynamic analysis and research of pollution risk of groundwater sources, this paper creates the dynamic assessment method of pollution risk of groundwater source area under the theory of “source-pathway-receptor”, and applies this method to one typical fissure karst groundwater source area in northern China. Following the 30-year petroleum pollutant migration simulation and pollution risk assessment of groundwater source area, this study finds that the very high risk zone is mainly located in Q Petrochemical Company and the surrounding area and the area adjacent to River Z. Within this period of thirty years, the pollution risk of groundwater source area has showed a dynamic trend that features an inverted “V” shape. The ratio of very high risk zone to the total area will be 18.1%, 17.47% and 16.62% during the tenth year, the twentieth year and the thirtieth year separately, and will reach the highest level of 19.45% during the fifteenth year. Meanwhile, the vertical migration distance of pollutant centre concentration changed from the surface soil at the outset to the deepest point of about 250 meters underground during the tenth year. The results of this risk assessment indicate the dynamic feature of pollution risk. The dilution, degradation and migration of petroleum pollutants in groundwater system contribute to an ultimate decline in pollution risk.     

Nano-scale pore structure and fractal dimension of organic-rich Wufeng-Longmaxi shale from Jiaoshiba area,Sichuan Basin: Investigations using FE-SEM,gas adsorption and helium pycnometry

Much attention have been recently paid to the upper Ordovician Wufeng shale (O₃w) and lower Silurian Longmaxi shale (S₁l) in the Jiaoshiba area of Sichuan Basin, which is now the largest producing shale gas field in China. Field emission scanning electron microscopy (FE-SEM), low pressure gas (N₂ and CO₂) adsorption, helium pycnometry, X-ray diffraction and geochemical analyses were performed to investigate the pore structure and fractal dimension of the pores in O₃w-S₁l shale formation in the Jiaoshiba area. FE-SEM images show that organic matter (OM) pores are dominant in the organic-rich samples and these pores are often irregular, bubble-like, elliptical and faveolate in shape, while in organic-poor samples, limited and isolated interparticle (interP), intraparticle (intraP) and OM pores are observed. Reversed S-shaped isotherms obtained from nitrogen adsorption are type Ⅱ, and hysteresis loops indicate that the shape of micropore in the samples is slit-or plate-like. BET surface areas and total pore volume vary from 12.2 to 27.1 m²/g and from 1.8 × 10⁻² to 2.9 × 10⁻² cm³/g, with an average of 19.5 m²/g and 2.3 × 10⁻² cm³/g, respectively. Adsorption volume from both N₂ and CO₂ adsorption increases with respect to TOC contents. Porosities obtained from helium porosimetry are comparable with these from gas (CO₂ and N₂) adsorption in O₃w-S₁l shale. However, porosity determined by quantitative FE-SEM analysis is much smaller, which is mainly related to limited resolution and the small areas of investigation.Based on the Frenkel-Halsey-Hill (FHH) model of low-pressure N₂ adsorption, fractal dimensions of the pores varied from 2.737 to 2.823. Relationships between pore structure parameters and TOC content, mineral composition and fractal dimension reveal that the fractal dimension is mainly associated with micropores. Samples with higher TOC content, higher quartz content and lower clay content tend to contain more heterogeneous micropores, resulting in higher fractal dimensions and more complicated pore structure in shales. Therefore, fractal dimension is an effective parameter to reflect the complexity of pore structure and the degree of micropore development in O₃w-S₁l shale.     

Early Mesozoic Southward Subduction of the Eastern Mongol–Okhotsk Oceanic Plate: Evidence from Zircon U–Pb–Hf Isotopes and Whole‐rock Geochemistry of Triassic Granitic Rocks in the Mohe Area,NE China

In this paper, we present new U–Pb zircon ages, Hf isotope data and major and trace elements for Early Mesozoic granitic rocks in Mohe area in the Erguna Massif of northeast China to elucidate the southward subduction of the eastern Mongol–Okhotsk Oceanic plate in Early Mesozoic. Zircons from two representative intrusions, syenogranites and monzogranites, in the Mohe area are euhedral–subhedral in shape, display oscillatory growth zoning in cathodoluminescence (CL) images, and have Th/U ratios of 0.10–0.72, and in combination these features indicating that the zircons are of igneous origin. U–Pb zircon dating results demonstrate that the syenogranites formed at 245.1 ± 1.4 Ma and monzogranites formed at 212.2 ± 1.7 Ma. These granitic rocks are characterized by high SiO₂, Al₂O₃ and (Na₂O + K₂O), low TFeO, MgO, TiO₂ and P₂O₅ concentrations, belonging to the high‐K calc‐alkaline series. They are enriched in LREE and large ion lithophile elements (e.g., Rb, K, and Sr), depleted in HREE and high field strength elements (e.g., Nb, Ta, Th, and Ti), as well as very weak negative Eu anomalies (Eu/Eu* = 0.48 ~ 1.01). Their zircon εHf(t) values range from −7.9 to −2.0 and range from 0.20 to 0.49, in response to their two‐stage Hf model ages (T_DM2) range from 1.40 Ga to 1.77 Ga range from 0.94 Ga to 1.24 Ga, respectively, indicating that primary magmas of syenogranites were derived from partial melting of newly accreted juvenile crustal material that formed from the enriched mantle during the Mesoproterozoic, monzogranites are generated by partial melting of newly accreted juvenile crustal material that formed from the depleted mantle during the Meso‐ to Neoproterozoic. We conclude, therefore, that the early Mesozoic granitic rocks of the Mohe area are associated with the continuous southward subduction of the Mongol–Okhotsk oceanic plate rather than the Paleo‐Asian and circum‐Pacific tectonic regimes.     

The simulation of LUCC based on Logistic-CA-Markov model in Qilian Mountain area,China

The Qilian mountain area was examined for using the Logistic-CA-Markov coupling model combined with GIS spatial analyst technology to research the transformation of LUCC, driving force system and simulate future tendency of variation. Results show that:(1) Woodland area decreased by 12.55%, while grassland, cultivated land, and settlement areas increased by 0.22%, 7.92%, and 0.03%, respectively, from 1986 to 2014. During the period of 1986 to 2000, forest degradation in the middle section of the mountain area decreased by 1,501.69 km~2. Vegetation cover area improved, with a net increase of grassland area of 38.12 km~2 from 2000 to 2014.(2) For constructing the system driving force, the best simulation scale was 210m×210m. Based on logistic regression analysis, the contribution(weight) of composite driving forces to land use and cover change was obtained, and the weight value was more objectively compared with AHP and MCE method.(3) In the natural scenarios, it is predicted that land use and cover distribution maps of Qilian mountain area in 2028 and 2042, and the Lee-Sallee index test was adopted. Over the next 27 years(2015–2042), farmland, woodland, grassland, settlement areas show an increasing trend, especially settlements with an obvious change of 0.56%. The area of bare land will decrease by 0.89%. Without environmental degradation, tremendous structural change of LUCC will not occur, and typical characteristic of the vertical zone of the mountain would remain. Farmland and settlement areas will increase, but only in the vicinity of Qilian and Sunan counties.     

Discussion on applying an analytical method to optimize the anti-freeze design parameters for underground water pipelines in seasonally frozen areas

Adopting the quasi-three-dimensional (Quasi-3D) numerical method to optimize the anti-freeze design parameters of an underground pipeline usually involves heavy numerical calculations. Here, the fitting formulae between the safe conveyance distance (SCD) of a water pipeline and six influencing factors are established based on the lowest water temperature (LWT) along the pipeline axis direction. With reference to the current widely used anti-freeze design approaches for underground pipelines in seasonally frozen areas, this paper first analyzes the feasibility of applying the maximum frozen penetration (MFP) instead of the mean annual ground surface temperature (MAGST) and soil water content (SWC) to calculate the SCD. The results show that the SCD depends on the buried depth if the MFP is fixed and the variation of the MAGST and SWC combination does not significantly change the SCD. A comprehensive formula for the SCD is established based on the relationships between the SCD and several primary influencing factors and the interaction among them. This formula involves five easy-to-access parameters: the MFP, buried depth, pipeline diameter, flow velocity, and inlet water temperature. A comparison between the analytical method and the numerical results based on the Quasi-3D method indicates that the two methods are in good agreement overall. The analytic method can be used to optimize the anti-freeze design parameters of underground water pipelines in seasonally frozen areas under the condition of a 1.5 safety coefficient.     

Uncertainty analysis of space–time prisms based on the moment-design method

Space–time prism (STP) is an important concept for the modeling of object movements in space and time. An STP can be conceptualized as the result of the potential path of a moving object revolving around in the three-dimensional space. Though the concept has found applications in time geography, research on the analysis and propagation of uncertainty in STPs, particularly under high degree of nonlinearity, is scanty. Based on the efficiency and effectiveness of the moment-design (M-D) method, this paper proposes an approach to deal with nonlinear error propagation problems in the potential path areas (PPAs) of STPs and their intersections. Propagation of errors to the PPA and its boundary, and to the intersection of two PPAs is investigated. Performance of the proposed method is evaluated via a series of experimental studies. In comparison with the Monte Carlo method and the implicit function method, simulation results show the advantages of the M-D method in the analysis of error propagation in STPs.     

Bias estimation and correction for triangle-based surface area calculations

The calculation of surface area is meaningful for a variety of space-filling phenomena, e.g., the packing of plants or animals within an area of land. With Digital Elevation Model (DEM) data we can calculate the surface area by using a continuous surface model, such as by the Triangulated Irregular Network (TIN). However, just as the triangle-based surface area discussed in this paper, the surface area is generally biased because it is a nonlinear mapping about the DEM data which contain measurement errors. To reduce the bias in the surface area, we propose a second-order bias correction by applying nonlinear error propagation to the triangle-based surface area. This process reveals that the random errors in the DEM data result in a bias in the triangle-based surface area while the systematic errors in the DEM data can be reduced by using the height differences. The bias is theoretically given by a probability integral which can be approximated by numerical approaches including the numerical integral and the Monte Carlo method; but these approaches need a theoretical distribution assumption about the DEM measurement errors, and have a very high computational cost. In most cases, we only have variance information on the measurement errors; thus, a bias estimation based on nonlinear error propagation is proposed. Based on the second-order bias estimation proposed, the variance of the surface area can be improved immediately by removing the bias from the original variance estimation. The main results are verified by the Monte Carlo method and by the numerical integral. They show that an unbiased surface area can be obtained by removing the proposed bias estimation from the triangle-based surface area originally calculated from the DEM data.     

基于层次分析法的德州地面沉降易发区的可拓学划分应用

基于物元模型的可拓学理论分析,提出华北平原德州地面沉降易发区的可拓学划分方法。利用地面沉降易发性分级、分类标准和影响因子分析,构造出经典域物元和节域物元,应用物元和可拓集合中的关联函数,建立了易发性等级综合评判的可拓评价模型,通过基于层次分析法的可拓学评价分析,采用正方形等间距(2km×2km)剖分了135个网格,再将剖分单元转化成面元,经过等值差分,得到了地面沉降易发区划结果。计算结果显示,可拓学方法能够实现定量化、多因子评价地面沉降易发性等级问题,从而科学合理地指导地面沉降综合分区防治。     

青海省祁连县大拉洞柏杨沟矿区煤岩层对比研究

通过研究太原组的沉积特征及聚煤变化规律,采用了标志层、煤岩层组合特征、层间距、测井曲线形态等方法,对矿区内煤层、岩层进行了综合对比。确定了该区主要煤层的赋存层位、形态、分布及变化特征,对周边地区找煤工作具有一定的指导作用。     

20世纪40年代初以来中国大陆沿海主要海湾形态变化

受气候变化和人类活动的影响,中国海岸带正在经历各种剧烈的变化,海湾形态变化是其中的一个重要方面.基于地形图资料与遥感影像提取20世纪40年代初至2014年7个时相的中国大陆沿海主要海湾的空间信息,进而分析海湾岸线类型结构,开发利用程度和空间位置的变化特征,以及海湾面积,形状和重心分布的变化特征.结果表明,20世纪40年代初以来,中国大陆沿海主要海湾的形态变化显著且复杂,主要体现为:海湾岸线的结构变化显著,自然岸线的长度和比例急剧减少;海湾岸线的开发利用程度持续增强;海湾岸线普遍背陆向海运动;海湾面积普遍大幅萎缩;海湾形状总体趋向复杂但个体间的差异在缩小;海湾重心分布位置普遍背陆向海迁移;海湾变化具有显著的时空差异性等.总的来说,以围填海为主的人类活动是导致海湾发生上述变化特征的主要原因.通过研究,指出应该对中国的海湾进行功能定位划分和分类管理,通过强化以海湾为基本单元的管理,促进中国海岸带综合管理能力的提高和管理目标的实现.