琼东南盆地深水区新近系海底扇沉积特征与资源潜力

综合利用钻井、岩心、薄片及分析化验资料研究了琼东南盆地深水区新近系海底扇沉积特征,并利用最新的三维地震资料,通过井震精细标定、多属性融合技术、方差体切片、三维地貌砂体镂空等综合技术手段,精细刻画了海底扇砂体的空间分布特征。研究结果表明,深水区新近系海底扇是由陆架区的砂体滑塌并二次搬运形成,形成过程具有多期次性。受不同物源的影响,海底扇岩性和物性存在较大的差异。海底扇岩性及沉积构造具有砂质滑塌、碎屑流、浊流和深水底流改造的特征。海底扇的沉积微相、厚度、砂泥比和砂泥岩空间配置关系直接控制了地震振幅反射强度和频率的变化。砂体纵向叠置,横向连片,并被后期泥质水道切割分块形成多个岩性圈闭。综合分析认为,深水区海底扇砂体发育区烃源条件优越,储盖配置关系和圈闭条件良好,具备形成大中型岩性油气藏的有利条件,具有较大的油气勘探潜力。     

Numerical simulation of groundwater and early warnings from the simulated dynamic evolution trend in the plain area of Shenyang,Liaoning Province (P.R. China)

Groundwater level is the most direct factor reflecting whether groundwater is in a virtuous cycle. It is the most important benchmark for deciding whether a balance can be struck between groundwater discharge and recharge and whether groundwater exploitation will trigger problems pertinent to environment, ecology and environmental geology. According to the borehole and long-term monitoring wells data in the plain area of Shenyang, a numerical groundwater model is established and used to identify and verify the hydrogeological parameters and balanced items of groundwater. Then the concept of red line levels, the control levels of groundwater is proposed, the dynamic evolution trend of groundwater under different scenarios is analyzed and predicted and groundwater alerts are given when groundwater tables are not between the lower limit and the upper limit. Results indicated: (1) The results of identification and verification period fitted well, and the calculation accuracy of balanced items was high; (2) with the implementation of shutting wells, groundwater levels in urban areas of Shenyang would exceed the upper limit water level after 2020 and incur some secondary disasters; (3) under the recommended scenario of water resources allocation, early-warnings for groundwater tables outside the range would occur in the year of 2020, 2023, 2025 respectively for successive wet, normal and dry years. It was imperative to reopen some groundwater sources and enhance real-time supervision and early-warning to prevent the occurrence of potential problems.     

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个网格,再将剖分单元转化成面元,经过等值差分,得到了地面沉降易发区划结果。计算结果显示,可拓学方法能够实现定量化、多因子评价地面沉降易发性等级问题,从而科学合理地指导地面沉降综合分区防治。