建立涌水量曲线方程的最小残差法

本文在综述了目前稳定流涌水量经验方程建立方法的基础上,根据最小二乘法原理,提出了一种新的方法—最小残差法,并给出了计算实例。计算表明,最小残差法是建立稳定法涌水量经验方程的一种较好方法,计算过程在微机上易于实现。     

Error propagation in equations for geochemical modeling of radiogenic isotopes in two-component mixing

This paper presents error propagation equations for modeling of radiogenic isotopes during mixing of two components or end-members. These equations can be used to estimate errors on an isotopic ratio in the mixture of two components, as a function of the analytical errors or the total errors of geological field sampling and analytical errors. Two typical cases (“Small errors” and “Large errors”) are illustrated for mixing of Sr isotopes. Similar examples can be formulated for the other radiogenic isotopic ratios. Actual isotopic data for sediment and basalt samples from the Cocos plate are also included to further illustrate the use of these equations. The isotopic compositions of the predicted mixtures can be used to constrain the origin of magmas in the central part of the Mexican Volcanic Belt. These examples show the need of high quality experimental data for them to be useful in geochemical modeling of magmatic processes.     

数据采集中特征点粗差定位方法的探讨

本文分析了特征点应满足的几种条件，给出了一次范围最小（Ｌ１）法，最小二乘法以及粗差探测法进行特征点粗差定位方法，并以模拟实例为例对它们进行了比较，得出了一些有意义的结论。     

确定多重现期暴雨公式参数的二次优化算法

四参数非线性多重现期暴雨公式在城市排水规划设计中有着广泛的应用。搜索算法与最小二乘法是优化计算的两个简单有效的方法,不过单独直接用于求解四参数非线性多重现期暴雨公式的参数比较困难。提出耦合最小二乘法及搜索算法确定多重现期暴雨公式参数的二次优化算法,该方法可以一次得到多重现期暴雨强度公式的参数,参数优化过程不需要图解。研究表明,计算结果比较客观;成果精度高,例题的平均绝对均方差为0.035mm/min。与遗传算法、蚁群算法等比较,该方法计算原理容易理解,计算简便,可以用Excel进行参数优化计算。     

岩土工程中位移监测方法及精度分析

在岩土工程位移监测中,选择哪种方法取决于监测时的条件;本文对监测中常用的测小角法、角度前方交会方法、距离前方交会方法、全站仪后方交会及极坐标综合方法等各种方法进行了研究。在目前各种文献中,一般都只是给出了监测点坐标的计算方法,充其量还给出监测点沿坐标轴方向的中误差以及点位中误差,而监测点其它方向的误差情况则无从计算。基于此,本文对以上各种监测方法的误差椭圆进行了推导,从而可以计算监测点沿各个方向误差并估算精度。     

原地生成宇宙成因核素测年技术思考（五）： 单核素计算地表最大抬升速率方法

在以往研究中,原地生成宇宙成因核素测年技术虽然理论上可以利用两种核素（如10Be和26Al）联立方程组同时解出暴露年龄和侵蚀速率,但在实践中却常采用单种核素（常用10Be）分别解出最小暴露年龄和最大侵蚀速率,26Al数据只用于辅助判断。这主要是因为受目前误差水平限制,10Be和26Al的测试数据结果,常常并未能达到联立两个方程解出暴露年龄和侵蚀速率两个参数的精度（当然也有其他因素影响）,这在投影图上常表现为被认为只有一次性暴露历史的样品,投影结果在稳态侵蚀岛之外,只有加上误差,才能进入或靠近稳态侵蚀岛。同理,在推导出包含生成速率加速率的计算等式后,虽然理论上抬升区长期暴露的样品也能够利用两种核素（如10Be和26Al）同时解出侵蚀速率和抬升速率,但实践上10Be和26Al浓度投影结果可能并未在抬升的正常投影区（稳态侵蚀岛上曲线附近）,从而解出的侵蚀速率和抬升速率合理性存疑。对此,参考计算最大侵蚀速率的原理,提出利用长期暴露样品的10Be浓度来解出抬升区地表最大抬升速率的计算方法（,从而可以对碰撞造山抬升区的抬升速率上限进行限定。     

Accuracy and efficiency of time integration methods for 1D diffusive wave equation

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.     

福建数字遥测台网单台震级偏差研究

使用福建数字遥测台网产出的观测报告,计算各台站单台震级与台网平均震级偏差和标准离差,分析了产生震级偏差的原因,并给出了台基校正值。结果显示,震级偏差主要来自于各台站的综合影响,进行台基校正是有意义的。     

Modelling in-situ upgrading of heavy oil using operator splitting method

The in-situ upgrading (ISU) of bitumen and oil shale is a very challenging process to model numerically because of the large number of components that need to be modelled using a system of equations that are both highly non-linear and strongly coupled. Operator splitting methods are one way of potentially improving computational performance. Each numerical operator in a process is modelled separately, allowing the best solution method to be used for the given numerical operator. A significant drawback to the approach is that decoupling the governing equations introduces an additional source of numerical error, known as the splitting error. The best splitting method for modelling a given process minimises the splitting error whilst improving computational performance compared to a fully implicit approach. Although operator splitting has been widely used for the modelling of reactive-transport problems, it has not yet been applied to the modelling of ISU. One reason is that it is not clear which operator splitting technique to use. Numerous such techniques are described in the literature and each leads to a different splitting error. While this error has been extensively analysed for linear operators for a wide range of methods, the results cannot be extended to general non-linear systems. It is therefore not clear which of these techniques is most appropriate for the modelling of ISU. In this paper, we investigate the application of various operator splitting techniques to the modelling of the ISU of bitumen and oil shale. The techniques were tested on a simplified model of the physical system in which a solid or heavy liquid component is decomposed by pyrolysis into lighter liquid and gas components. The operator splitting techniques examined include the sequential split operator (SSO), the Strang-Marchuk split operator (SMSO) and the iterative split operator (ISO). They were evaluated on various test cases by considering the evolution of the discretization error as a function of the time-step size compared with the results obtained from a fully implicit simulation. We observed that the error was least for a splitting scheme where the thermal conduction was performed first, followed by the chemical reaction step and finally the heat and mass convection operator (SSO-CKA). This method was then applied to a more realistic model of the ISU of bitumen with multiple components, and we were able to obtain a speed-up of between 3 and 5.