Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

基于蚁群算法的城市可持续发展综合评价

在可持续发展体系分类的基础上，应用蚁群算法对可持续发展各分类指数公式和综合指数公式进行优化，建立可持续发展评价模型。将该模型应用于乌鲁木齐市可持续发展水平评价，评价结果与现状分析一敛。结果表明该模型简单适用、有良好的通用性。     

分组密码的安全性与强化设计

分组密码是数据通讯中最常用的数据加密方式，以DES为例分析现有分组加密算法的安全隐患，并提出了可变密钥加密和变长密文输出两个新思路，可应用于所有现有分组加密算法以提高安全性，并就该方法的安全性、效率、具体应用做出了分析。     

用遗传算法对自动控制系统数学模型寻优

探讨用遗传算法对数学模型进行优化。考虑到控制系统稳健性的要求 ,用遗传算法寻找出控制系统最佳稳定域 ,实现控制系统数学模型的寻优。阐述了用遗传算法求解问题的步骤和参数的取值 ,并用仿真实验对优化结果进行了检验。结果表明 :在控制系统数学模型的优化中 ,遗传算法具有其他算法无可比拟的优越性。     

天文(重力)水准的新手段及其计算方法

介绍当今确定天文大地垂偏差的新仪器－CCD自动天体测量仪和确定两点间高程异常差的新方法，该方法用单极坐标代替繁琐的双极坐标进行计算，讨论了天文重力水准的误差及垂线偏差非线性影响等问题，估计在不久的将来，用这一手段施测山区似大地水准面的精度可望达到厘米级。     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

Migration of the Ganga river and its implication on hydro-geological potential of Varanasi area, U.P., India

Borehole data reveals that during Late Quaternary, the Ganga river was non-existent in its present location near Varanasi. Instead, it was flowing further south towards peripheral craton. Himalayan derived grey micaceous sands were being carried by southward flowing rivers beyond the present day water divide of Ganga and mixed with pink arkosic sand brought by northward flowing peninsular rivers. Subsequently, the Ganga shifted to its present position and got incised. Near Varanasi, the Ganga river is flowing along a NW-SE tectonic lineament. The migration of Ganga river is believed to have been in response to basin expansion caused due to Himalayan tectonics during Middle Pleistocene times. Multi-storied sand bodies generated as a result of channel migration provide excellent aquifers confined by a thick zone of muddy sediments near the surface. Good quality potable water is available at various levels below about 70 m depth in sandy aquifers. Craton derived gravelly coarse-to-medium grained sand forms the main aquifer zones of tens of meter thickness with enormous yield. In contrast, the shallow aquifers made up of recycled interfluve silt and sandy silt occur under unconfined conditions and show water-level fluctuation of a few meters during pre-and post-monsoon periods.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.