乌鲁木齐河山区流域360年径流量的重建

本文利用在乌鲁木齐河山区流域所建立的树轮年表，重建了乌鲁木齐河山区流域３６０年上年７月～当年３月的平均径流量。校准方程的相关系数为０．６７０，交叉检验的误差缩减值达０．３６６。在重建的３６０年径流量的变化中，有４个偏丰期和３个偏枯期，第２，３两个偏丰期与乌鲁木齐河河源１号冰川的两次冰进期相对应。对重建径流量的丰枯频率分析发现，平水年份出现最多，偏枯水年份多于偏丰水年份约５．８％，特枯水年没有出现，特丰水年出现６次，约占１．７％。这表明３６０年来径流量变化基本上稳定。     

流体在同心环状空间中流动时的压力损失

同心环状空间是工程上较常见到的流体流动空间，在此空间中流体流动情况的分析较复杂。工程上对其流动情况的分析常常是引用流体力学中的分析结果，但往往由于对环状空间结构的条件重视不够，引用时很容易出毛病。经过仔细推导分析，指出其应用条件，并有助于对同心环状空间中流体流动时压力损失的分析与计算     

树轮──流域径流变化的记录

根据树轮稳定碳同位素组成（δ１３Ｃｔ）序列可较灵敏地记录降水变化的事实，利用长白山红松树轮δ１３Ｃｔ序列与松花江年径流量之间的平行变化关系，重建了近２００年来松花江年径流量的变化：年径流量变化有明显的准２２年周期；１９００年之前，径流量在周期波动中呈上升趋势，１９００年之后，波动中呈下降趋势     

Are open fractures necessarily aligned with maximum horizontal stress?

Fluid flow in fractured rock is an increasingly central issue in recovering water and hydrocarbon supplies and geothermal energy, in predicting flow of pollutants underground, in engineering structures, and in understanding large-scale crustal behaviour. Conventional wisdom assumes that fluids prefer to flow along fractures oriented parallel or nearly parallel to modern-day maximum horizontal compressive stress, or S_Hmax. The reasoning is that these fractures have the lowest normal stresses across them and therefore provide the least resistance to flow. For example, this view governs how geophysicists design and interpret seismic experiments to probe fracture fluid pathways in the deep subsurface. Contrary to these widely held views, here we use core, stress measurement, and fluid flow data to show that S_Hmax does not necessarily coincide with the direction of open natural fractures in the subsurface (>3 km depth). Consequently, in situ stress direction cannot be considered to predict or control the direction of maximum permeability in rock. Where effective stress is compressive and fractures are expected to be closed, chemical alteration dictates location of open conduits, either preserving or destroying fracture flow pathways no matter their orientation.     

湖北省土壤侵蚀景观空间异质性及驱动因子分析

地域分异是地球表层大小不等、内部具有一定相似性地段之间的相互分化以及由此产生的差异。为了研究不同区位土壤侵蚀问题，从土壤生态景观及系统论出发，运用地质学、地理学、景观生态学、环境学的理论和研究方法，研究湖北省土壤侵蚀景观空间格局及其驱动因子，使土壤侵蚀问题研究提高到一个新的水平。湖北省土壤侵蚀景观具有南北分带、东西分区，为一不对称的断块一环组合，土壤流呈现向长江、江汉盆地中心轴带辐聚、单流向特点。景观空间异质性形成的首要驱动因子是大地构造背景，以房县一襄樊一广济断裂带为界，南北两侧地壳物质组成和构造发展史存在较明显的差异，现代气候带、降雨量、温热程度及土地利用等差异，造成了湖北省区域土壤地理、土壤生态的分异，形成湖北省土壤生态带、区具有南北分带，东西分区的宏观格局；其次大兴安岭一武陵山深部构造陡变带两侧新构造运动强度差异、大别造山带构造强烈隆升，导致土壤侵蚀强度的西强东弱、南北强中间弱的态势；成土母岩差异性决定了土壤可蚀性的多变；空间上“土壤侵蚀内城区”分布在湖北省的周边地区，经济贫困、管理落后，这一地区的经济水平与水土流失间形成“自反馈作用”，这一现象在我国水土保持、生态建设工作中应该引起重视。     

基于二维小波变换的FMI图象分割

为了从FMI资料中定量提取参数，一个重要的步骤是从实际FMI资料中分离出反映溶孔、溶洞、裂缝的子图像。本文给出的方法，考虑图像像元邻域的特征，应用二维小波变换求出目标与背景边缘的点集，按这个边缘点集的坐标点所对应的原图像像素灰度值的平均值作为分割阈值进行图像分割。实际资料处理表明，应用这种方法可以从实际的FMI资料中准确地分割出孔洞、裂缝的子图像并且可以按深度段连续自动处理，为后续定量计算参数奠定了良好基础。     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

Penny-shaped fractures revisited

All methods of seismic characterization of fractured reservoirs are based on effective media theories that relate geometrical and material properties of fractures and surrounding rock to the effective stiffnesses. In exploration seismology, the first-order theory of Hudson is the most popular. It describes the effective model caused by the presence of a single set of thin, aligned vertical fractures in otherwise isotropic rock. This model is known to be transversely isotropic with a horizontal symmetry axis (HTI). Following the theory, one can invert the effective anisotropy for the crack density and type of fluid infill of fractures, the quantities of great importance for reservoir appraisal and management.Here I compute effective media numerically using the finite element method. I deliberately construct models that contain a single set of vertical, ellipsoidal, non-intersecting and non-interconnected fractures to check validity of the first-order Hudson’s theory and establish the limits of its applicability. Contrary to conventional wisdom that Hudson’s results are accurate up to crack density e ≈ 0.1, I show that they consistently overestimate the magnitudes of all effective anisotropic coefficients ε^(V), δ^(V), and γ^(V). Accuracy of theoretically derived anisotropy depends on the type of fluid infill and typically deteriorates as e grows. While the theory gives | ε^(V)|, |δ^(V)|, |γ^(V)| and close to the upper bound of the corresponding numerically obtained values for randomly distributed liquid-filled fractures, theoretical predictions of ε^(V), δ^(V) are not supported by numerical computations when the cracks are dry. This happens primarily because the first-order Hudson’s theory makes no attempt to account for fracture interaction which contributes to the final result much stronger for gas- than for liquid-filled cracks. I find that Mori-Tanaka’s theory is superior to Hudson’s for all examined crack densities and both types of fluid infill.The paper was presented at the 11th International Workshop on Seismic Anisotropy (11IWSA) held in St. John’s, Canada in 2004.     

Anisotropic damage of rock joints during cyclic loading: constitutive framework and numerical integration

This work describes a constitutive framework for modeling the behavior of rough joints under cyclic loading. Particular attention is paid to the intrinsic links between dilatancy, surface degradation, and mobilized shear strength. The framework also accounts for the important effect of shear‐induced anisotropy. The resulting approach is fully three‐dimensional and is not restricted to plane‐displacement kinematics. Both the governing formulation and an algorithm for implicit numerical integration are presented. While the proposed methods are general, we also postulate a specific model that is compared with experimental data. It employs relatively few free parameters but shows good agreement with laboratory tests. Copyright © 2013 John Wiley & Sons, Ltd.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.