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.     

局部大地水准面精化系统软件的设计与实现

给出了局部大地水准面精化系统的总体结构以及各分系统的详细设计与实现方法,开发了旨在对相关重力测量数据实现自动化处理的局部大地水准面精化系统工具软件,并应用于某区域实际大地水准面的精确计算。实践证明,该系统能有效提高相关数据的处理效率,有助于局部大地水准面精化过程的规范化。     

GIS支持下长江口南支河道百年来的演变

利用长江口南支河道1900—2001年的5幅海图资料,在GIS的支持下,建立了百年来长江口南支河道的水下数字高程模型(DEM),在此基础上对其演变特征进行了计算和分析。结果表明:(1)1900—2001年间南支河道表现为微冲,不同时段的冲淤状况不同。(2)南支河道的冲淤演变主要表现为5个水道和4个成形沙体的冲淤变化,各水道的最大水深和成形沙体的总面积和总体积都呈不断增大趋势。(3)南支河道各主要断面以复式河槽为主,主槽位于南岸,且最大水深不断增大。位于南支河道入口处和出口处断面的变化比中段断面的变化剧烈。(4)南支河道的演变是自然演变和人工控制共同作用的结果。     

Plane Symmetric Domain Wall in Lyra Geometry

In this paper general solutions are found for domain walls in Lyra geometry in the plane symmetric spacetime metric given by Taub. Expressions for the energy density and pressure of domain walls are derived in both cases of uniform and time varying displacement field β. It is also shown that the results obtained by Rahaman et al [IJMPD, 10, 735 (2001)] are particular case of our solutions. Finally, the geodesic equations and acceleration of the test particle are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

矢量数据栅格化的一种有效方法——环绕数法

矢量数据栅格化是GIS数据处理的一项重要任务。目前,已存在多种栅格化方法,但是每种方法都有各自的局限性。本文首次提出了环绕数法,它基于计算几何转角理论,通过计算跨越数来确定环绕数,是一种简单、高效、易于实现的栅格化方法。应用环绕数法,不但可以完成简单多边形的栅格化,而且可以完成自交多边形的栅格化。     

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.