The identifiability of parameters in a water quality model of the Biebrza River, Poland

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.     

Numerical modeling of 3-D terrain effect on MT field

Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain. Project supported by the National Natural Science Foundation of China.     

τ-P域地震资料三维叠前深度偏移及其实用算法

本文提出了在τ-P域实现三维叠前深度偏移的计算方法，该方法把地震道集由时-空域变换到τ-P域，然后对共P数据体进行相移或相移加插值偏移，将偏移后的各共P数据体叠加得到叠前深度偏移结果．该算法能够实现陡倾地层的正确归位，降低偏移处理维数，减少计算机内存需求，易于并行处理，但是计算量和所花费的计算机时间仍然很大．为此我们又提出了针对目的层的分时（层）偏移方法，由于只考虑目的层，可节省大量的计算时间．文中以二维实际地震资料及三维数值模型资料的处理说明了上述方法的可行性与有效性．     

普通克里格法在矿产储量计算中的应用

普通克里格法能充分利用样品信息，有效地提高块段平均品位及矿产储量精度。以刁泉银铜矿为实例，介绍了该方法的理论和方法应用以及储量计算流程图。     

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

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

高密度电阻率法在地下涵洞勘测中的应用

本文介绍的高密度电阻率法是一种新兴的电阻率方法，与常规电阻率方法相比，高密度电阻率法具有成本低，效率高，测点密集，反映地电信息更丰富，更直观以及资料解释简单方便等特点。本文运用高密度电阻率法勘测地下涵洞，实践证明，高密度电阻率法获得的图形，质量高，解释误差小，是一种方便实用，效果良好，能广泛应用于工程勘测的物探方法。     

Environmental Properties of Minerals and Contaminants Purified by the Mineralogical Method

The investigation of the environmental properties of minerals, i.e., environmental mineralogy, is a branch of science dealing with interactions between natural minerals and spheres of the Earth surface as well as a reflection of global change, prevention of ecological destruction, participation in biomineralogy, and remediation of environmental pollution. Pollutant treatment by natural minerals is based on the natural law and reflects natural self-purification functions in the inorganic world, similar to that of the organic world - a biological treatment. A series of case studies related to natural self-purification, which were mostly completed by our group, are discussed in this paper. In natural cryptomelane there is a larger pseudotetragonal tunnel than that formed by [MnO6] octahedral double chains, with an aperture of 0.462-0.466 nm2, filled with K cations. Cryptomelane might be a real naturally-occurring mineral of the active octahedral molecular sieve (OMS-2). CrⅥ-bearing wastewater can be treate     

NHC硅化法及其应用

NHC硅化加固方法 ,已在数十项工程中推广应用 ,均取得满意效果 ,介绍了这一方法的作用原理及其特性 ,并列举了在几项工程中的应用实例。     

A Least-squares Approach to Depth Determination from Numerical Horizontal Self-potential Gradients

We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.