土壤吸湿-湿过程中水力传导度的拟合研究

非饱和含水介质在吸湿和脱湿的不同过程中,水力传导度随介质含水率的变化规律是不同的。本文选取大沽河流域两种代表性砂土,对其吸湿-脱湿两个过程分别进行试验观测,并使用VanGenuchten—Mualem模型对其吸湿-脱湿过程进行水力传导度的拟合。研究结果表明：砂土在吸湿过程的非饱和水力传导度与脱湿过程的非饱和水力传导度大小并不相同,在含水率相同的情况下,吸湿过程的水力传导度要比脱湿过程大,在含水率较低时（小于15％）相差可达1个数量级。     

Modelling Positional Uncertainty of Line Features by Accounting for Stochastic Deviations from Straight Line Segments

The assessment of positional uncertainty in line and area features is often based on uncertainty in the coordinates of their elementary vertices which are assumed to be connected by straight lines. Such an approach disregards uncertainty caused by sampling and approximation of a curvilinear feature by a sequence of straight line segments. In this article, a method is proposed that also allows for the latter type of uncertainty by modelling random rectangular deviations from the conventional straight line segments. Using the model on a dense network of sub‐vertices, the contribution of uncertainty due to approximation is emphasised; the sampling effect can be assessed by applying it on a small set of randomly inserted sub‐vertices. A case study demonstrates a feasible way of parameterisation based on assumptions of joint normal distributions for positional errors of the vertices and the rectangular deviations and a uniform distribution of missed sub‐vertices along line segments. Depending on the magnitudes of the different sources of uncertainty, not accounting for potential deviations from straight line segments may drastically underestimate the positional uncertainty of line features.     

Mesozoic-Cenozoic tectonics and geodynamics of Altai, Tien Shan, and Northern Kazakhstan, from apatite fission-track data

Apatite fission-track (AFT) thermochronological modeling as a diagnostic tool for periods of stability (peneplanation) and tectonic activity (orogeny) has been broadly used in tectonic studies of Central Asia in recent years. We discuss more than 100 AFT ages of samples from the Kyrgyz Tien Shan and Altai and compare them with AFT data from northern Kazakhstan. Geological, geomorphological, and AFT data indicate intense activity in the Late Cenozoic Eurasian continental interior. The impact from the India-Eurasia collision on the northern Tien Shan, Altai, and northern Kazakhstan regions showed up at 11, 5, and 3 Ma, respectively, as a result of stress propagation into the continent, with the ensuing reactivation and mountain growth. We hypothesize that a distant effect of the Late Cenozoic India-Eurasia collision was to rejuvenate Paleozoic fault zones and to deform the Mesozoic sedimentary cover north of the collision front as far as the West Siberian Plate. The reactivation facilitated formation of tectonic oil and gas traps. The activity in northern Central Asia under the effect of the Indian indentation into Eurasia appears to continue and may evolve to include uplift of southern West Siberian plate with uplift.     

某基坑工程土钉支护工作性状试验研究

通过对广州天河荟雅苑基坑工程土钉支护进行试验和研究,分析了土钉支护体系的工作性状,并为同类工程的设计提供了一定的依据。     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Numerical simulation of migration–accumulation of oil resources

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Recommended Values of 239Pu, 240Pu and 239+240Pu Concentrations in Reference Material IAEA-315 (Marine Sediment) Estimated by Thermal Ionisation Mass Spectrometry, Inductively Coupled Plasma-Mass Spectrometry and Alpha Spectrometry

The recommended concentrations of ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu in reference material IAEA‐315 (marine sediment) were estimated by three analytical methods: isotope dilution thermal ionisation mass spectrometry (TIMS), isotope dilution inductively coupled plasma‐mass spectrometry (ICP‐MS) and alpha spectrometry. The determination of ²³⁹Pu and ²⁴⁰Pu (^239+240Pu by alpha spectrometry) was carried out with samples from randomly selected bottles using each method. Plutonium‐238 was also measured by alpha spectrometry. A plutonium‐242 reference material was used as a spike for the quantitative analysis. The influence of ²⁴²Pu in the samples was therefore calculated; however, this contribution was less than the range of uncertainty and did not influence the final results. The obtained data were statistically analysed using variance component analysis and paired comparison. The combined standard uncertainties from “method/measurement”, “bottle” and “sub‐sample” were in the order of 3 to 6%. The main contributions to the uncertainty were from the material heterogeneity and from systematic differences between methods. Based on this study with twenty‐seven analyses using 10–14 g sample mass, concentrations of (38 ± 3) Bq kg^?1, (28 ± 3) Bq kg^?1 and (66 ± 4) Bq kg^?1 are proposed as recommended values for ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu, respectively, and (9.5 ± 0.4) Bq kg^?1 for ²³⁸Pu as an information value in reference material IAEA‐315. In mass concentration units, these amount to (16.4 ± 1.2) ng kg^?1, (3.3 ± 0.4) ng kg^?1 and (0.015 ± 0.003) ng kg^?1 for ²³⁹Pu, ²⁴⁰Pu and ²³⁸Pu, respectively. The certified reference materials NIST 4350B and NIST 4354 were also analysed by TIMS for quality assurance of the method used in this study.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。     

顾及DEM误差自相关的坡度计算模型精度分析

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。