Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

Preliminary 3–D geological model of the Kalgoorlie region,Yilgarn Craton,Western Australia,based on deep seismic‐reflection and potential‐field data

The granite‐greenstone terranes of the Eastern Goldfields Province, Yilgarn Craton, Western Australia, are a major Australian and world gold and nickel source. The Kalgoorlie region, in particular, hosts several world‐class gold deposits. To attempt to understand why these deposits occur where they do, it is important to understand the crustal architecture in the region and how the major mineral systems operate in this architecture. One way to understand these relationships is to develop a detailed 3–D geological model for the region. The best method to map the 3–D geometry of major geological structures is by acquisition and interpretation of seismic‐reflection profiles. To contribute to this aim, a grid of deep seismic‐reflection traverses was acquired in 1999 to examine the 3–D geometry of the region in an area including the Kalgoorlie mineral region and mineral fields to the north and west. This grid was tied to the 1991 regional deep seismic traverse and 1997 high‐resolution seismic profiles in the same region. The grid covers an area measuring approximately 50 km wide by 50 km long and extended to a depth of approximately 50 km (below the base of the crust in this region). The resulting 3–D geological model was further constrained by both surface geological data and geophysical interpretations, with the seismic interpretations themselves also constrained by gravity and magnetic modelling. The 3–D model was used to investigate the geometric relationships between the major faults and shear zones in the area, the relationship between the granite‐greenstone succession and the basement, and the spatial relationships between the greenstones and the granites. Interpretation of the grid of seismic lines and construction of the 3–D geological model confirmed the existence of the detachment surface and led to the recognition that the granite‐greenstone contact usually occurs at a much shallower level than the detachment. Also, west‐dipping faults in the vicinity of the Golden Mile, including the Abattoir Shear through to Boulder‐Lefroy Fault, appear to be more important than previously thought in controlling the structure of that area. An antiformal thrust stack occurs beneath a triangle zone centred on the Golden Mile. The Black Flag Group was deposited in a probable extensional setting, and late extension was also probably more important than previously thought. The granite‐gneiss domes were uplifted by the formation of antiformal thrust stacks at depth beneath them.     

Gas generation and expulsion characteristics of Middle–Upper Triassic source rocks,eastern Kuqa Depression,Tarim Basin,China: implications for shale gas resource potential

Frontier exploration in the Kuqa Depression, western China, has identified the continuous tight-sand gas accumulation in the Lower Cretaceous and Lower Jurassic as a major unconventional gas pool. However, assessment of the shale gas resource in the Kuqa Depression is new. The shale succession in the Middle–Upper Triassic comprises the Taliqike Formation (T₃t), the Huangshanjie Formation (T₃h) and the middle–upper Karamay Formation (T_2–3k), with an average accumulated thickness of 260 m. The high-quality shale is dominated by type III kerogen with high maturity and an average original total organic carbon (TOC) of about 2.68 wt%. An improved hydrocarbon generation and expulsion model was applied to this self-contained source–reservoir system to reveal the gas generation and expulsion (intensity, efficiency and volume) characteristics of Middle–Upper Triassic source rocks. The maximum volume of shale gas in the source rocks was obtained by determining the difference between generation and expulsion volumes. The results indicate that source rocks reached the hydrocarbon expulsion threshold of 1.1% VR and the hydrocarbon generation and expulsion reached their peak at 1.0% VR and 1.28% VR, with the maximum rate of 56 mg HC/0.1% TOC and 62.8 mg HC/0.1% TOC, respectively. The volumes of gas generation and expulsion from Middle–Upper Triassic source rocks were 12.02 × 10¹² m³ and 5.98 × 10¹² m³, respectively, with the residual volume of 6.04 × 10¹² m³, giving an average gas expulsion efficiency of 44.38% and retention efficiency of 55.62%. Based on the gas generation and expulsion characteristics, the predicted shale gas potential volume is 6.04 × 10¹² m³, indicating a significant shale gas resource in the Middle–Upper Triassic in the eastern Kuqa Depression.     

Generation of correlated properties in heterogeneous porous media

The spatial distribution of rock properties in porous media, such as permeability and porosity, often is strongly variable. Therefore, these properties usefully may be considered as a random field. However, this variability is correlated frequently on length scales comparable to geological lengths (for example, scales of sand bodies or facies). To solve various engineering problems (for example, in the oil recovery process) numerical models of a porous medium often are used. A need exists then to understand correlated random fields and to generate them over discretized numerical grids. The paper describes the general mathematical methods required to do this, with one particular method (the nearest neighbor model) described in detail. How parameters of the mathematical model may be related to rock property statistics for the nearest neighbor model is shown. The method is described in detail in one, two, and three dimensions. Examples are given of how model parameters may be determined from real data.     

上海地区深部构造特征和地震活动趋势探讨

本文介绍由重磁资料计算上海地区居里面、莫霍面、均衡重力异常及地幔流应力场的方法和成果,并在对成果初步分析的基础上,探讨了深部地质构造特征及其与地震活动趋势的关系,其中着重讨论了居里面及其与地震活动的关系。文章最后提出了对上海地区地震活动趋势的几点认识。     

A numerical scheme for calculating stresses and strain rates in glaciers

The stress and strain-rate fields in a glacier or other type of rock with known rheological properties can be calculated by using a numerical solution technique. Calculations are based on force-balance equations and the constitutive relation for polar ice, and do not involve any mathematical approximations. The geometry is prescribed. Basal velocities also are prescribed, either by specifying their values or by relying on a sliding relation. All other quantities are determined numerically.     

In situ yield behaviour of a coal pillar

Summary An indication of the complete stress-strain behaviour of coal was obtained through instrumenting a headgate yield pillar and monitoring the changes in vertical and horizontal pillar stress, pillar strain, and roof-floor convergence as a function of time and distance from the moving longwall face. A semiautomated monitoring program was conducted over a 5 month period as the longwall face advanced toward and eventually passed the instrumented pillar. The development of a confined core and yield zone within the pillar was confirmed by the vertical and horizontal stress changes. Thein situ pre-failure (208 to 661 MPa) and post-failure moduli (55 to 208 MPa) of the coal were determined from the change in pillar stress versus pillar strain. Thein situ moduli are an order of magnitude lower than the 3450 MPa, commonly assumed for coal.     

渤海陆缘地震的力学特征和成因

渤海陆缘是华北强震活动区,如1966年的邢台震群,1967年的河间地震,1969年的渤海地震,1975年的海城地震,1976年的唐山震群和1983年的菏泽地震,都是在该区发生的。显然,深入研究该地区的构造运动和地震成因是有重要科学意义的。本文拟从应力场角度进行研究。