A level-set corrector to an adaptive multiscale permeability prediction

We consider the inverse problem of permeability estimation for two-phase porous-media flow. The novel approach is based on regularization by zonation, where the geometry and size of the regions are chosen adaptively during the optimization procedure. To achieve this, we have utilized level-set functions to represent the permeability. The available data are sparsely distributed in space; hence, it is reasonable to confine the estimation to coarse-scale structures. The level-set approach is able to alter the boundaries between regions of different permeability without strict restrictions on their shape; however, when the data are sparse, a reasonable initial guess for the permeability is required. For this task, we use adaptive multiscale permeability estimation, which has the potential of identifying main permeability variations. These are described by a piecewise constant function, where the constant values are attained on rectangular zones. In the current work, we develop a level-set corrector strategy, assuming adaptive multiscale permeability estimation as a predictor.     

Roughness control on hydraulic conductivity in fractured rocks

The influence of joint roughness on the typologies of fluid flow inside fractures is well known and, thanks to experiences in the field of hydraulics, it has been studied from both a physical and mathematical point of view. Nevertheless, the formulations adopted by traditional hydraulic models are hardly applicable in the geological field, because of the difficulty encountered in the roughness parameter estimation. Normally this parameter can be estimated using the joint roughness coefficient (JRC), which considers both the asperity height and its regularity and directional trend. The main advantage in using the JRC arises from the fact that it can easily be obtained from geological-technical surveys and from comparison with the standard Barton profiles. Some relationships have been built up that allow for the estimation of the hydraulic conductivity tensor (an essential parameter for understanding water flow in fractured rock masses), not only as a function of traditional parameters like aperture, spacing, dip and dip direction, etc., but also of joint roughness, precisely expressed in terms of the roughness coefficient. These relationships have been studied initially from a theoretical point of view and then practically, through laboratory investigations.

---

Resumen  Se conoce muy bien la influencia de la rugosidad de las grietas en las tipologías del flujo de fluidos a lo interior de las fracturas y gracias a las experiencias en el campo de hidráulica ha sido posible estudiarla desde puntos de vista matemáticos y físicos. Sin embargo, las formulaciones adoptadas por los modelos hidráulicos tradicionales tienen poca aplicabilidad en el campo geológico debido a la dificultad relacionada con la estimación del parámetro de rugosidad. Normalmente este parámetro puede estimarse usando el coeficiente de rugosidad de grieta (JRC) el cual considera tanto la altura de la aspereza como su regularidad y tendencia direccional. La principal ventaja de utilizar el JRC se deriva del hecho que puede obtenerse fácilmente de levantamientos técnico-geológicos y de la comparación con los perfiles Standard Barton. Se han construido algunas relaciones que permiten la estimación del tensor de conductividad hidráulica (un parámetro esencial para el entendimiento del flujo de agua en masas de roca fracturadas), no solo en función de parámetros tradicionales como apertura, espaciado, buzamiento y dirección de buzamiento, etc., sino también en función de la rugosidad de la grieta estimada con precisión en términos del coeficiente de rugosidad. Estas relaciones se han estudiado inicialmente desde un punto de vista teórico y luego de modo práctico a través de investigaciones de laboratorio.

---

Résumé  L’influence de la rugosité des joints sur les types d’écoulement de fluide dans les fractures est bien connue et a été étudiée aussi bien du point de vue physique que mathématique grace à des expériences menées dans le domaine de l’hydraulique. Cependant les formulations adoptées dans les modèles hydrauliques traditionnels sont difficilement applicables dans le domaine de la géologie à cause de la difficulté rencontrée pour estimer la rugosité. Ce paramètre peut normalement être apprécié grace au coefficient de rugosité du joint (JRC), lequel prend en compte à la fois la hauteur de l’aspérité ainsi que sa régularité et sa direction. Le principal avantage dans l’utilisation du JRC réside dans le fait qu’il peut facilement être obtenu à partir d’études techniques-géologiques et par comparaison avec la classification de Barton. Des relations qui permettent une estimation du tenseur de conductivité hydraulique (un paramètre essentiel pour comprendre l’écoulement de l’eau dans les masses rocheuses fracturées) ont été élaborées, pas seulement en fonction de paramètres traditionnels tels que l’ouverture, l’espacement, l’inclinaison et la direction d’inclinaison, etc , mais aussi en prenant en compte la rugosité des joints à travers le coefficient de rugosité. Ces relations ont initialement été étudiées d’un point de vue théorique puis expérimentalement à travers des recherches en laboratoire.

     

磁铁矿采空区处理施工工艺与技术方法探讨

文章系统地介绍了采用地质钻机成孔,高压注浆机进行充填注浆的施工方法。为了将水泥砂浆及水泥粉煤灰浆填入采空区,采用能施加高压的输送泵及高压注浆设备、铸铁注浆管、水泥封孔的注浆工艺处理采空区铁路路基。在信息化施工过程中,通过试验不断地调整设计方法、改进施工工艺,在处理巨型磁铁矿采空区工程实践中达到了预期的目标。施工中,在控制充填边界、封孔方法及施工工艺等方面有一定的创新。提出的采空区施工技术控制、质量检测与评价方法,对同类工程勘察、设计及施工具有一定的参考价值。     

利用水源热泵开采浅层地热能若干问题的探讨

地下水源热泵采能技术作为一种浅层低温地热能的可持续开发方式,在我国推广应用已有近十年的时间。尽管推广时间不长,但发展速度非常之快,尤其是近3～5年。由于多方面的原因,我国部分地下水源热泵采能工程逐渐暴露出一些问题,需要给予足够的重视。本文介绍了地下水源热泵采能系统的运行模式及特点,简要地对国内外热泵采能技术的应用现状进行了综述,并对该项技术在我国推广应用过程中存在的一些问题,进行了一定的分析和讨论。     

多普勒雷达探测“艾利”台风风场不对称结构

2004年8月25日, “艾利”台风在靠近我国台湾东北部海域后不同于正常路径, 先后发生两次左折, 先是西北转偏西, 再偏西转西南, 这种路径在历史上极为罕见。通过“艾利”台风的雷达观测事实, 从单多普勒雷达和双多普勒雷达均发现在台风前进方向的右前方, 即西北象限风速极值大于其他象限, 进一步证明了“艾利”台风的风场分布遵循波数为1的非对称性, 从位于台风不同象限的福州长乐和厦门两站雷达资料, 以及常规的地面观测、高空探测资料和NCEP数值预报再分析资料都得到证实。研究表明, “艾利”台风处在大陆高压东南侧, 台风外围存在一支7 m/s左右的环境风场引导气流。因此, 西北象限风速极值大于其他象限的不对称分布可能是由环境风场和台风本身两方面共同造成的, 这种不对称分布有利于维持“艾利”台风向西南方向移动。在我国东南沿海地区, 由于单部多普勒雷达不可能同时探测到台风的西北、东北和西南3个象限的风速极值, 而双多普勒雷达或多多普勒雷达联合探测台风, 则可以同时探测到3个象限的风速极值, 根据风速极值可能存在的不对称分布情况来预测台风的路径转折趋势, 因此重视双多普勒雷达或多多普勒雷达联合探测对于提高监测预报水平具有重要的意义。     

Remote sensing research issues of the National Land Use Change Program of China

In 1999, the Ministry of Land and Resources (MLR) of China launched the National Land Use Change Program especially to monitor the scale and distribution of urban expansion and the decrease in cultivated land through remote sensing technology. This Program has been carried out annually and continuously for seven years since then and played an important role in the policy-making of MLR about land management and planning. This paper gives an overview about this Program and discusses several research issues. First, the remote sensing data sources and other ancillary data used in this Program are presented. The approaches for image preprocessing, i.e. radiometric normalization, image geometric rectification and image fusion are then introduced with an emphasis on the algorithm development for image registration. Second, land use change detection technique is the most critical and complex aspect of the Program. The methodologies for change detection using either bi-temporal image pair or one existing land use map and one remotely sensed image are detailed. Third, since the data of land use changes derived from remote sensing will be operationally used for local and central government, field validation and accuracy assessment are crucial to ensure the reliability of change detection results. The strategy of field work and the resulting accuracy evaluations is presented. The land use and change information derived from remotely sensed data has wide applications for land management, including land use database updating, verification of land use planning and monitoring of national high-tech parks. Last, suggestions on how to make full use of the images and change detection result, to improve the consistency of land use classification and to develop change detection algorithms for diverse and complex remote sensing data are given.     

一种扩展的基于结构检测的最优多视极化SAR相干斑滤波算法

全极化数据的可得,使得利用极化散射矩阵以减少多视全极化SAR图像中的相干斑成为了可能。其中,被称为最优的多视极化白化滤波器（MPWF）就是一个典型例子。然而,为提高滤波器参数估计的精度,同时能够自适应的检测图像中的结构特征,在滤波时采用适当的窗算法是必要的。本文提出了两个有效的窗算法,并用1994年NASA SIR-C/X-SAR,L band,经过四视处理的天山森林的极化数据进行了仿真。实验结果表明,这两种算法不仅很好的达到了降噪效果,而且对图像的纹理信息以及结构特征（即边缘特征、线性特征、强散射体特征）的保持效果是最佳的。     

1480型多普勒天气雷达中尺度气旋识别模块设计

借鉴国内外中尺度气旋识别算法的研究成果特别是NSSL(National Severe Storm Laboratory)的中尺度气旋识别算法,结合多普勒天气雷达二次产品软件的特点,提出了改进的中尺度气旋识别算法,以动态链接库的形式,设计完成了中尺度气旋识别模块.对2002年5月27日安徽省北部地区的一次强对流天气过程的雷达资料进行了分析,并与北京敏视达(METSTAR)雷达有限公司Build 10.8软件系统生成的中尺度气旋产品进行了对比分析,结果表明该模块有比较好的识别能力,预警效果较好.     

On the flow structure and turbulence during sweep and ejection events in a wind-tunnel model canopy

Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.     

Analysis of the impact of conservation measures on stream flow regime in catchments of the Loess Plateau,China

Catchments in the Loess Plateau have been under the influence of human activities for centuries. In the last four decades, soil conservation measures have accelerated and intensified. These measures were designed to reduce soil erosion, improve agricultural productivity, and enhance environmental quality. It is important to evaluate the effects of these measures on hydrology in order to develop sustainable catchment management plans in the region. This study evaluated changes in stream flow data for four selected catchments in the Loess Plateau following large‐scale soil conservation measures. The non‐parametric Mann–Kendall test was used to identify trends in annual stream flow and the results showed significant downward trends in three of the four catchments. The Pettitt test indicated that a change point occurred in 1978 in these three catchments. Annual precipitation in all the catchments showed no significant trend during the period of record. Comparison of daily flow duration curves for two 20‐year periods (1957–1978) and (1979–2003) showed significant changes in stream flow regime. Reduction in most percentile flows varied between 20 and 45%, and the reduction in low flows was greatest. Overall, the reductions in daily flow were increasing with time, with significant changes occurring in the 1990s. However, it is not clear whether these catchments have seen the full effects of the soil conservation measures, so the results of this study might underestimate the final impact of soil conservation on stream flow regime. Copyright © 2007 John Wiley & Sons, Ltd.