经历剪切变形历史的岩石节理表面粗糙特性分析

基于两类岩石节理剪切试验，采用非接触式激光测距仪，实测了剪切前后节理表面粗糙数据，统计分析了岩石节理表面粗糙特性的三维粗糙参数的物理含义，并将这些参数与相应的二维参数进行了充分比较。结果表明，三维参数与二维参数间存在的近似线性关系，而多个三维参数相对比单一的二维参数更全面有效地描述了节理表面粗糙特性。     

昌马水库右岸山体稳定性研究

昌马水库导流排砂洞塌方导致右岸山体发生松动 ,本文在大量地质资料的基础上 ,对右岸岩体裂隙网络进行模拟 ,同时 ,采用极限分析能量法和有限元数值模拟方法等对右岸山体的稳定性进行分析研究 ,并提出相应的加固方案。     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Texture and petrology of modern river, beach and shelf sands in a volcanic back-arc setting, northeastern Japan

Abstract The focus in the present study is on characterizing spatial patterns of textural and petrological variabilities, and on evaluating mechanisms influencing the textural and petrological components of modern river, beach and shelf sands in a volcanically active back‐arc tectonic setting. Abashiri Bay and the surrounding area in eastern Hokkaido, Japan, has volcanic source land within a back‐arc region associated with subduction of the Pacific Plate beneath the Okhotsk (North American) Plate. A total of 41 river, beach and shelf sands were obtained for grain‐size and modal composition analyses. Multivariate analytical techniques of hierarchical cluster and principal component analyses were performed on the textural and petrological data for investigating relations among quantitative variables. On the basis of grain‐size data, four sedimentary zones were identified: zone I, palimpsest zone; zone II, relict zone; zone III, modern (proteric) zone; zone IV, coastal sedimentary zone. All sands are feldspatholithic and quartz‐deficient. The framework (quartz, feldspar and rock fragment) modal compositions were also classified into four clusters, A–D. The characteristic components of each cluster are as follows: cluster A, felsic volcanic rock fragments; cluster B, andesitic–basaltic volcanic rock fragments; cluster C, mixed or plagioclase; cluster D, sedimentary rock fragments. Almost all sands in western and central Abashiri Bay belong to cluster A, where the original compositions are influenced by Kutcharo pyroclastic flow deposits. Andesitic–basaltic lava and Neogene volcaniclastic and sedimentary rocks have a major influence on the compositions of shelf sands in eastern Abashiri Bay. The modal compositions are basically controlled by the source lithology. Compositional maturity (percentage of quartz to feldspar and rock fragments; Q/FR%) slightly increased, in order, from river (1.2), zone IV (coastal, 1.7), zone II (relict, 2.2), zone I (palimpsest, 3.6), to zone III (modern proteric, 7.0). Greater maturity in the recycled sediments is indicative of weathering under the sea or abrasion by transportation induced by sea‐level fluctuations, waves, or sea currents. Several controlling factors – (i) source lithological; (ii) mineralogical; (iii) climatic; and (iv) geomorphological controls – might still cause low maturity through all sedimentary zones other than the continental margin sands previously reported.     

Saprolite Remnants as Indicators of Pre-Glacial Landform Genesis in Southeast Sweden

Twenty‐six sites with remnants of gravelly saprolites (grus) have been located in southeast Sweden. Joint block hills (castle kopjes) and steep rock walls with weathered joints as well as rounded boulders are documented to have an origin in deep weathering and subsequent stripping of saprolites. The saprolite remnants and landforms result from the fragmentation of the re‐exposed sub‐Cambrian peneplain along fracture systems. Only shallow saprolites occur on the elevated intact parts of the sub‐Cambrian peneplain, while saprolites up to 20 m thick are encountered in areas where the sub‐Cambrian peneplain is fractured and dissected. Neogene uplift with reactivation of the weathering system is thought to be the main cause of saprolite formation. Deep weathering is thus judged to have been the major agent of landform formation in the study area, while glacial and glaciofluvial erosion has contributed mainly by stripping saprolites, detaching corestones, and plucking joint blocks along weathered joints.     

Conditional Simulation of Random Fields with Bivariate Gamma Isofactorial Distributions

This work focuses on a random function model with gamma marginal and bivariate isofactorial distributions, which has been applied in mining geostatistics for estimating recoverable reserves by disjunctive kriging. The objective is to widen its use to conditional simulation and further its application to the modeling of continuous attributes in geosciences. First, the main properties of the bivariate gamma isofactorial distributions are analyzed, with emphasis in the destructuring of the extreme values, the presence of a proportional effect (higher variability in high-valued areas), and the asymmetry in the spatial correlation of the indicator variables with respect to the median threshold. Then, we provide examples of stationary random functions with such bivariate distributions, for which the shape parameter of the marginal distribution is half an integer. These are defined as the sum of squared independent Gaussian random fields. An iterative algorithm based on the Gibbs sampler is proposed to perform the simulation conditional to a set of existing data. Such ‘multivariate chi-square’ model generalizes the well-known multigaussian model and is more flexible, since it allows defining a shape parameter which controls the asymmetry of the marginal and bivariate distributions.     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

水库群水沙调控的单-多目标调度模型及其应用

为缓解内蒙古河段"二级悬河"形势,以黄河上游沙漠宽谷河段为研究对象,以龙羊峡水库、刘家峡水库为调控主体,开展黄河上游水沙调控研究。建立了输沙量、发电量最大的单目标模型以及多目标模型;分别采用自迭代模拟算法、逐次逼近动态规划算法(DPSA)和改进的非支配排序遗传优化算法(NSGA-Ⅱ)求解模型;设置了初始、常规、优化和联合优化4种方案。通过实例计算,联合优化调度方案的区间总冲刷量达到了0.38亿t,梯级发电量148.22亿kW·h。该方案以较小的电量损失换来了输沙量的大幅度增加,水沙调控效果显著,推荐为最优方案。研究成果量化了水沙调控效果和各目标间的转化规律,为开展黄河上游水沙调控提供了决策依据,具有重要的应用价值和实际指导意义。