非均匀陆面条件下区域蒸散量计算的遥感模型

非均匀陆面条件下的区域蒸散计算是一个复杂的问题。文中首先在利用遥感资料求取地表特征参数 (如植被覆盖度、地表反照率等 )的基础上 ,建立了裸露地表条件下的裸土蒸发和全植被覆盖条件下植被蒸腾计算模型 ,然后结合植被覆盖度 (植被的垂直投影面积与单位面积之比 )给出非均匀陆面条件下的区域蒸散计算方法。实测资料验算表明该模型具有较高的计算精度。文章最后利用该模型对中国北方地区的蒸散量进行了计算 ,并对该研究区蒸散的特点进行了分析     

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.     

非均匀介质热蠕变流动的数值求解

针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.     

200hPa非地转风和非均匀层结与梅雨暴雨的关系

本文通过1991年梅雨后期，7月1—13日200hPa非地转场分析得出，梅雨期200hPa非地转风主要由积云对流动量转换形成。其次是风速平流，并通过数值试验得以证明。非地转风形成的高层辐散，可进一步使暴雨增强与维持。另外梅雨期非均匀层结对非地转风及暴雨也有较大的影响。     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

Multi-scale modeling of three-phase–three-component processes in heterogeneous porous media

Flow and transport processes in porous media occur on different spatial and temporal scales and may also be locally different. Additionally, the structure of the porous medium itself generally shows a high dependence on the spatial scale.     

Estimation Of Sensible Heat Flux From The Surface Temperature Wave And One-Time-Of-Day Air Temperature Observation, Part Ii – Heterogeneous Surfaces

By means of the algorithm presented in Part I of this study, the temporal course H(t) and the daily mean H of the sensible heat flux H can be determined from measurements of the thermodynamic surface temperature (as a function of time) and from a one-time-of-day air temperature observation. Inaddition to these temperatures, one needs estimates of daily mean wind speed,of the roughness lengths of momentum and heat transfer, and of the displacementheight. In Part I, the algorithm was derived for areas with homogeneous surfaces,i.e., with uniform surface temperature, and the method was verified with measurements taken during several field campaigns. The root mean square error for the temperature difference between surface and air, in the comparison between measurement and model, amounted to one or two kelvin, and the error of H was 10 to 25 per cent. The method can be used to determine the sensible heat flux from measurements of surface temperatures by satellites, but can also be applied to ground based measurements.In Part II, the procedure is generalized for areas that consist of various surface types (sub-regions) with different surface temperatures, and can be usedwhen only a few (at least one) air temperature measurements per day are available over only one of the different sub-regions. This generalization should allow improvements to the estimates for H(t) by means of temperature measurements from, e.g., NOAA/AVHRR or LANDSAT/TM, taking into account the heterogeneity of the area contained in one METEOSAT pixel. Criteria are given as to whether effective (areal mean) surface temperatures and roughness lengths may be used for the computation of H or if the above mentioned generalized procedure has to be applied. The new algorithm is verified by measurements sampled during the field campaigns EFEDA 91 (Spain) and HIBE 89 (Hildesheimer Börde in Germany), and by using synthetic data (due to the lack of measured data) for one further combined surfacetype [soil and water (lakes)].     

On the Parameterisation of the Effective Roughness Length for Momentum Transfer over Heterogeneous Terrain

We develop a parameterisation for the effective roughness length of terrain that consists of a repeating sequence of patches, in which each patch is composed of strips of two roughness types. A numerical model with second-order closure in the turbulent stress is developed and used to show that: (i) the normalised Reynolds stress develops as a self-similar profile; (ii) the mixing-length parameterisation is a good first-order approximation to the Reynolds stress. These findings are used to characterise the blending layer, where the stress adjusts smoothly from its local surface value to its effective value aloft. Previous studies have assumed that this adjustment occurs abruptly at a single level, often called the blending height. The blending layer is shown to be characterised by height scales that arise naturally in linear models of surface layer flow over roughness changes, and calculations with the numerical model show that these height scales remain appropriate in the nonlinear regime. This concept of the blending layer allows the development of a new parameterisation of the effective roughness length, which gives values for the effective roughness length that are shown to compare well with both atmospheric measurements and values determined from the second-order model.     

Evaporation over a Heterogeneous Land Surface: EVA_GRIPS and the LITFASS-2003 Experiment—An Overview

The Evaporation at Grid/Pixel Scale (EVA_GRIPS) project was realised in order to determine the area-averaged evaporation over a heterogeneous land surface at the scale of a grid box of a regional numerical weather prediction or climate model, and at the scale of a pixel of a satellite image. EVA_GRIPS combined surface-based and airborne measurements, satellite data analysis, and numerical modelling activities. A mesoscale field experiment, LITFASS-2003, was carried out in the heterogeneous landscape around the Meteorological Observatory Lindenberg (MOL) of the German Meteorological Service in May and June, 2003. The experiment was embedded in the comprehensive, operational measurement program of the MOL. Experimental determination of surface fluxes on a variety of spatial scales was achieved by employing micrometeorological flux stations, scintillometers, a combination of ground-based remote sensing instruments, and the Helipod, a turbulence probe carried by a helicopter. Surface energy fluxes were also derived from satellite data. Modelling work included the use of different Soil–Vegetation–Atmosphere Transfer schemes, a large-eddy simulation model and three mesoscale atmospheric models. The paper gives an overview on the background of EVA_GRIPS, and on the measurements and meteorological conditions during LITFASS-2003. A few general results are discussed.