无结构网格上二维浅水流动的数值模拟

提出在无结构网格上建立有限体积高性能格式族的统一框架,通过引入跨单元界面法向数值通量的逆风分解,将一维Osher、TVD两种通量分裂格式自然地推广至二维浅水方程组。给出了各种情况下浅水方程组有限体积法边界处理的计算公式。最后,利用该格式对穹包溢流、陆地动边界和河口潮流三个问题进行计算。结果表明,这类守恒格式具有高精度、无振荡性以及处理复杂流态过渡、自动捕俘间断和模拟陆地动边界的功能。     

The SOMS model and its application to Lake Neuchâtel

A three-dimensional numerical circulation model (SOMS) based on primitive equations is described. The algorithm, by which Coriolis and vertical diffusion terms are treated implicitly while mass is still conserved exactly (algebraically), is discussed in detail. The model is applied to Lake Neuchâtel (Switzerland), to determine the general circulation under influence of the most prevailing wind.     

Development of the instantaneous unit hydrograph using stochastic differential equations

Recognizing that simple watershed conceptual models such as the Nash cascade ofn equal linear reservoirs continue to be reasonable means to approximate the Instantaneous Unit Hydrograph (IUH), it is natural to accept that random errors generated by climatological variability of data used in fitting an imprecise conceptual model will produce an IUH which is random itself. It is desirable to define the random properties of the IUH in a watershed in order to have a more realistic hydrologic application of this important function. Since in this case the IUH results from a series of differential equations where one or more of the uncertain parameters is treated in stochastic terms, then the statistical properties of the IUH are best described by the solution of the corresponding Stochastic Differential Equations (SDE's). This article attempts to present a methodology to derive the IUH in a small watershed by combining a classical conceptual model with the theory of SDE's. The procedure is illustrated with the application to the Middle Thames River, Ontario, Canada, and the model is verified by the comparison of the simulated statistical measures of the IUH with the corresponding observed ones with good agreement.     

Conditioning of coefficient matrices of Ordinary Kriging

The solution of a set of linear equations is central to Ordinary Kriging. Computers are commonly applied because of the amount of data and work involved. There has, until recently, been little attention devoted toward the conditioning of kriging matrices. This article considers implications of conditioning upon numerical stability, instead of on robustness which has been the main focus of past work. The effect of properties of the stationary covariance matrix on the conditioning of the kriging matrix is discussed. The relationship between the covariance and autocorrelation functions allows some conclusions about the conditioning of covariance matrices, based on past work in deconvolution. The conditioning of some coefficient matrices of stationary kriging, defined in terms of either the semivariogram or the covariance, is examined.     

有限元素法全倾角波动方程偏移

本文试图从两个方面提高波动方程偏移的效果:1.使用全方向波动方程进行偏移,使之适用于任意倾角界面的反射;2.使用有限元素法提高微分方程数值解法的精度。     

A review on stochastic differential equations for applications in hydrology

Fundamentals of the theory of stochastic calculus and stochastic differential equations (SDE's) which are finding increasing application in water resources engineering are reviewed. The basics of probability theory, mean square calculus and the Wiener, white Gaussian and compound Poisson processes are given in preparation for a discussion of the general Itô SDE with drift, diffusion and jump discontinuity terms driven by Gaussian white noise and compound Poissionian impulses. Also discussed are stochastic integration and the derivation of moment equations via the Itô differential rule. The lierature of SDE's is reviewed with an emphasis on the more accessible sources.     

利用二维DLT及光束法平差进行数字摄像机标定

利用共线方程和二维DLT之间的对应关系导出了由二维DLT的8个参数表达的主纵线方程式，讨论了主点初值和求解方法，给出了单张像片摄像机参数分解不惟一性及临界序列的证明，详细推导了利用二维DLT参数分解摄像机外方位元素初值的实用算法，论述了利用光束法平差进行摄像机标定的数学模型，实际图像数据实验取得了很好的结果，验证了本文所提出的摄像机标定算法的可行性。     

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     

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

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

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

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