金属矿地震散射波场的数值模拟研究

地震方法用于寻找石油、煤田等沉积层控矿床已取得许多成功经验和大量的成果；对于金属矿可依据微扰理论，通过地震波散射响应的强弱推断矿体。模拟实验结果表明，散射波相干性的好坏与杂乱散射体的不均匀性有关；不均匀性越强，产生的散射相干性就越好。     

边坡稳定可靠性的随机有限元分析

随机有限元法可以处理土性参数的变异性和空间相关性。对二阶摄动随机有限元法的摄动理论和程序进行了研究，提出了偏导矩阵的组集方法，采用正态分布随机变量的正交变换法来提高计算效率。考虑土性参数随机场作用和土性指标之间的互相关性，建立边坡局部抗剪和总体稳定性可靠度的随机有限元分析模型，对某实际土坡进行了可靠度计算，计算结果较为符合实际。     

Estimating groundwater recharge rates in the southeastern Holstein region, northern Germany

In the southeastern Holstein region, located to the east of the metropolitan zone of Hamburg, northern Germany, a groundwater investigation program was conducted from 1984 to 2000 by the State Agency for Nature and Environment (Landesamt für Natur und Umwelt, LANU) of Schleswig-Holstein, Germany, with the aim of providing long-term, ecologically acceptable groundwater management plans for the region. The focal point of the investigation comprised the determination of groundwater recharge rates. The investigation method was based on the transfer of available lysimeter results from other regions to comparable regions within the area studied. With the help of lysimeter equations, potential amounts of percolation water were calculated. The groundwater recharge rate was then determined after subtraction of the surface runoff which was calculated for the entire area. All computations were performed with a spreadsheet program. Groundwater recharge rates were calculated for two areas. One consisted of roughly determining groundwater recharge rates for the total region (1,392 km²) of southeastern Holstein. The overall goal of these investigations was to identify potential areas of water exploitation. Areas in which groundwater recharge rates are high and groundwater outflow is low are particularly suited to water exploitation, since inflow rates into deeper aquifers are high. These areas are located on the flanks of the Elbe and Stecknitz River valleys. Subsurface groundwater runoff to these lowlands would be reduced through groundwater withdrawal. However, the resulting decline in shallow groundwater tables would be so small that it would have no detrimental ecological effects. Groundwater recharge rates were also calculated for a 110-km² area in the outskirts of Hamburg (Grosshansdorf model area) which is intensively developed for water supply. These investigations showed that the amount of groundwater recharge is already being withdrawn to a large extent. Approximately 65% of the recharge rate is currently withdrawn by the waterworks in this area, thus making further increases in exploitation rates unjustifiable from an ecological point of view. Electronic Publication     

Upscaling Uncertain Permeability Using Small Cell Renormalization

Sedimentary rocks have structures on all length scales from the millimeter to the kilometer. These structures are generally associated with variations in rock permeability. These need to be modeled if we are to make predictions about fluid flow through the rock. However, existing computers are not powerful enough for us to be able to represent all scales of heterogeneity explicitly in our fluid flow models—hence, we need to upscale. Small cell renormalization is a fast method for upscaling permeability, derived from an analogue circuit of resistors. However, it assumes that the small scale permeability distribution is known. In practice, this is unlikely. The only information available about small scale properties is either qualitative, derived from the depositional setting of the reservoir, or local to the wells as a result of coring or logging. The influence of small scale uncertainty on large scale properties is usually modelled by the Monte Carlo method. This is time-consuming and inaccurate if not enough realisations are used. This paper describes a new implementation of renormalization, which enables the direct upscaling of uncertain small-scale permeabilities to produce the statistical properties of the equivalent coarse grid. This is achieved by using a perturbation expansion of the resistor-derived equation. The method is verified by comparison with numerical simulations using the Monte Carlo method. The prediction of expected large-scale permeability and its standard deviation are shown to be accurate for small cell standard deviations of up to 40% of the mean cell value, using just the first nonzero term of the perturbation expansion. Inclusion of higher order terms allows larger standard deviations to be modeled accurately. Evaluation of cross-terms allows correlations of actual cell values, over and above the background structure of mean cell values. The perturbation method is significantly faster than conventional Monte Carlo simulation. It needs just two calculations whereas the Monte Carlo method needs many thousands of realisations to be generated and renormalized to converge. This results in significant savings in computer time.     

用切贝谢夫配点法消除地球自振常微分方程组在地心处的奇异性

用切贝谢夫配点法求解地球自振常微分方程组，无需进一步改化即可消除这组方程在地心处的奇异性，并能获得高精度的结果     

地面总辐射的异常减少对气候影响的数值试验

在与实际台风场近似一致的环境流场条件下，应用柱会标中的正压无辐散模式，讨论了台风中扰动的稳定性问题。研究表明，台风暖心结构的台风外围适当强度的冷空气活动有利限台风中扰动发展，且台风愈强，台风中扰动愈易发展。     

不同方案求解非线性化学动力学方程组的比较

在一个简化的气相化学模式的基础上，比较了Ｈｙｂｒｉｄ、ＱＳＳＡ、Ｓｋｌａｒｅｗ三种计算方案求解非线性化学动力学方程组的差异，并引入误差分配和线性组合两种质量守恒技术，分析它们对模拟结果的影响．研究结果表明：不同方案对有机烃类浓度的计算几乎没有影响，但对其它气态物浓度有一定的影响．用ＱＳＳＡ和Ｓｌａｒｅｗ方案预测的结果更为相近，而采用Ｈｙｂｒｉｄ方案得到的气态物浓度的平衡时间略为提前．从计算效率来看，ＱＳＳＡ方案优于其它两种方案，两种质量守恒技术均能改善模拟结果．     

A SIMPLIFIED APPROACH TO MODELING 3D SEDIMENT-LADEN TURBULENT FLOWS

ＡＳＩＭＰＬＩＦＩＥＤＡＰＰＲＯＡＣＨＴＯＭＯＤＥＬＩＮＧ３ＤＳＥＤＩＭＥＮＴ－ＬＡＤＥＮＴＵＲＢＵＬＥＮＴＦＬＯＷＳＤｏｎｇｈｕｏＺＨＯＵａｎｄＳａｍＳ．Ｙ．ＷＡＮＧＡｂｓｔｒａｃｔ：Ａ３－ｄｎｕｍｅｒｉｃａｌｍｏｄｅｌｆｏｒｓｉｍｕｌａｔｉｎｇｓ...     

Stochastic perturbation analysis of groundwater flow. Spatially variable soils,semi-infinite domains and large fluctuations

As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.     

反演二维弹性波动方程组多参数的渐近展开法

采用摄动法对二维弹性波动方程组的系数反问题进行渐近近似,给出了反演系数的具体计算方法。用摄动法解微分方程组的系数反问题还是首次。此外,还采用了地球物理分层模型,这种分层模型便于更精细地反映地球不同圈层的物理特征,更接近于实际情况。同时对进一步提高反演分辨率、增加反演参数、对病态方程组如何克服不适定性等问题,作了具体阐述,具有重要的理论意义与应用价值。