当前国际水文学科技术领域的研究方向

根据联合国有关机构发布的公报和世界各国近年来的学术研究情况，整理出了近期水文学科在地表水、地下水、水质、水资源系统、大气～土壤～植被、大陆侵蚀、冰雪、遥感技术和示踪技术等领域的热门研究课题，在一定程度上反映了当前国际上水文学科的研究和发展方向，对国内同行有一定的参考价值。     

An invariance in the carbazoles of petroleum

A remarkable invariance in the ratio of 1,3-dimethylcarbazole （DMC） to 1,6-DMC was discovered in crude oils from the Pearl River Mouth Basin, South China Sea. The remarkably invariant ratio is kept at a constant of about unity regardless of their concentrations, sources or maturities for the sampies. In combination with the molecular structures of 1,3- and 1,6-DMCs, the invariance might indicate that the nitrogen compounds share a common precursor with a skeleton of 1-methylcarbazole and are formed through methylation at C3 and C6 with an essentially identical rate.     

Gas accumulation from oil cracking in the eastern Tarim Basin: A case study of the YN2 gas field

Oil and gas exploration in eastern Tarim Basin, NW China has been successful in recent years, with several commercial gas accumulations being discovered in a thermally mature to over-mature region. The Yingnan2 (YN2) gas field, situated in the Yingnan structure of the Yingjisu Depression, produces gases that are relatively enriched in nitrogen and C₂₊ alkanes. The δ¹³C₁ (−38.6‰ to −36.2‰) and δ¹³C₂ values (−30.9‰ to −34.7‰) of these gases are characteristic of marine sourced gases with relatively high maturity levels. The distributions of biomarkers in the associated condensates suggest close affinities with the Cambrian–Lower Ordovician source rocks which, in the Yingjisu Sag, are currently over-mature (with 3–4%Ro). Burial and thermal maturity modeling results indicate that paleo-temperatures of the Cambrian–Lower Ordovician source rocks had increased from 90 to 210 °C during the late Caledonian orogeny (458–438 Ma), due to rapid subsidence and sediment loading. By the end of Ordovician, hydrocarbon potential in these source rocks had been largely exhausted. The homogenization temperatures of hydrocarbon fluid inclusions identified from the Jurassic reservoirs of the YN2 gas field suggest a hydrocarbon emplacement time as recent as about 10 Ma, when the maturity levels of Middle–Lower Jurassic source rocks in the study area were too low (<0.7%Ro) to form a large quantity of oil and gas. The presence of abundant diamondoid hydrocarbons in the associated condensates and the relatively heavy isotopic values of the oils indicate that the gases were derived from thermal cracking of early-formed oils. Estimation from the stable carbon isotope ratios of gaseous alkanes suggests that the gases may have been formed at temperatures well above 190 °C. Thus, the oil and gas accumulation history in the study area can be reconstructed as follows: (1) during the late Caledonian orogeny, the Cambrian–Lower Ordovician marine source rocks had gone through the peak oil, wet gas and dry gas generation stages, with the generated oil and gas migrating upwards along faults and fractures to form early oil and gas accumulations in the Middle–Upper Ordovician and Silurian sandstone reservoirs; (2) since the late Yanshanian orogeny, the early oil accumulations have been buried deeper and oil has undergone thermal cracking to form gas; (3) during the late Himalayan orogeny, the seals for the deep reservoirs were breached; and the gas and condensates migrated upward and eventually accumulating in the relatively shallow Jurassic reservoirs.     

“探宝号”船可变阵列震源的实现及应用

2001年，广州海洋地质调查局对“探宝号”船原有震源更新升级为阵列可变的BOLT震源系统。随着更新改造工程的成功，该船已成为具备承担高分辨率、常规、低频深水油气地震勘探的多功能地震调查船。本文介绍了地震勘探震源的基本工作原理以及三种阵列的分布、波形、频谱等特性。     

An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.     

Mixing of event and pre‐event water in a shallow Entisol in sloping farmland based on isotopic and hydrometric measurements,SW China

Water percolation and flow processes in subsurface geologic media play an important role in determining the water source for plants and the transport of contaminants or nutrients, which is essential for water resource management and the development of measures for pollution mitigation. During June 2013, the dynamics of the rainwater, soil water, subsurface flows and groundwater in a shallow Entisol on sloping farmland were monitored using a hydrometric and isotopic approach. The results showed that effective mixing of rainwater and soil water occurred in hours. The rebound phenomenon of δD profiles in soils showed that most isotope‐depleted rainwater largely bypassed the soil matrix when the water saturation in the soil was high. Preferential‐flow, which was the dominant water movement pattern in the vadose zone, occurred through the whole soil profile, and infrequent piston‐flow was mainly found at 20–40 cm in depth. The interflow in the soil layer, composed of 75.2% rainwater, was only generated when the soil profile had been saturated. Underflow in the fractured mudrock was the dominant flow type in this hillslope, and outflow was dominated by base flow (groundwater flow) with a mean contribution of 76.7%. The generation mechanism of underflow was groundwater ridging, which was superimposed upon preferential‐flow composed mainly of rainwater. The quick mixing process of rainwater and soil water and the rapid movement of the mixture through preferential channels in the study soil, which shows a typical bimodal pore size distribution, can explain the prompt release of pre‐event water in subsurface flow. Water sources of subsurface flows at peak discharge could be affected by the antecedent soil water content, rain characteristics and antecedent groundwater levels. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.