石臼港污染物三维扩散系数的计算     

张法高  王文素  李培泉  朱兰部  朱沛 《海洋科学》1990,14(5):1-7

  相似文献   

底部边界层颗粒态物质迁移特征的放射性核素分析方法     

李炎  G.W.Berger 《海洋科学》1992,16(4):66-68

利用颗粒态放射性核素携带的颗粒物历经过程信息,我们提出用泥沙扩散方程和颗粒态放射性核素扩散方程联解底部边界层颗粒态物质迁移参数的方法。分析实例的样品取自荷兰Waden Sea南部Balgzand潮滩(砂坪)和Mok湾潮滩(泥坪)的两个站位(BG1和Mok2)。示踪核素为~(234)Th,~(210)Pb和~(137)Cs,其放射性比度由r能谱测出。  相似文献   

In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer     

Robel A. Gebrekristos  Allen M. Shapiro  Brent H. Usher 《Hydrogeology Journal》2008,16(4):629-639

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.  相似文献   

Probabilistic modeling of seismically triggered landslides using Monte Carlo simulations     

Huabin Wang  Gonghui Wang  Fawu Wang  Kyoji Sassa  Yuping Chen 《Landslides》2008,5(4):387-395

The 2004 Chuetsu earthquakes of Niigata (Japan) triggered numerous landslides, and the most widespread types of landslides were highly disrupted, relatively shallow slides and soil (debris) flows. This paper presented a method to evaluate slope instability using Newmark displacement on a pixel-by-pixel basis in a given area. The proposed method was able to integrate Newmark displacement modeling and Monte Carlo simulations within geographical information systems. In the modeling, an empirical attenuation relationship was utilized to calculate Arias intensity over this study area, and the variability of geotechnical parameters was taken into account to calculate coseismic landslide displacement. Before deriving the displacement from related inputs, the Monte Carlo simulations ran 1,500 times and generated 1,500 displacement values for each grid cell, and then means and standard deviations of displacement were calculated and probabilistic distributions can be obtained. Finally, given 10 cm as a threshold value of displacement, estimated probabilities of displacement exceeding 10 cm were shown as a map of seismic landslide hazards. The resulting hazard map was classified into four categories from very low to high level.  相似文献   

顾及DEM误差自相关的坡度计算模型精度分析   总被引：11，自引：1，他引：10  

刘学军  卞璐  卢华兴  朱莹 《测绘学报》2008,37(2):0-249

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。  相似文献   

Information structure analysis for quantitative assessment of mineral resources and the discovery of a silver deposit     

Wang  Yutian  Tan  Bingqi  Wang  Yifeng  Wu  Jiangtao 《Natural Resources Research》1994,3(4):284-294

We propose an information-structure-analysis (ISA) method to quantify the correlations between quantitative and qualitative variables as well as within each type of variable. This method is applied to the evaluation of mineral resources in the western Zheijiang Province of China. The district contains a number of silver-bearing Fe–Cu–Pb–Zn mineral deposits near igneous complexes and FeCuPbZn zones away from the complexes. Silver anomalies occur not only in the known Fe–Cu–Zn–Pb deposits, but also in the country rock, suggesting the possible existence of silver deposits far from the igneous complexes.The tonnage distribution of silver is modeled by Monte Carlo simulation. This simulation is conducted on the basis of the correlations between silver (Ag) and lead (Pb), since no known data on silver is available. The known tonnage distribution of lead in 11 control cells was used to approximate the tonnage distribution of silver in the Monte Carlo simulation. With ISA and Monte Carlo methods, the total amount of potential polymetallic resources in 49 cells in the western Zhejiang Provice is predicted. Significantly, a deposit with about 24 tonnes of silver has been found within our exploration target area.  相似文献   

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix     

C. R. Dietrich  G. N. Newsam 《Mathematical Geology》1995,27(2):207-228

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.  相似文献   

Comparison of structural and least-squares lines for estimating geologic relations   总被引：1，自引：0，他引：1  

Garnett P. Williams  Brent M. Troutman 《Mathematical Geology》1990,22(8):1027-1049

Two different goals in fitting straight lines to data are to estimate a true linear relation (physical law) and to predict values of the dependent variable with the smallest possible error. Regarding the first goal, a Monte Carlo study indicated that the structural-analysis (SA) method of fitting straight lines to data is superior to the ordinary least-squares (OLS) method for estimating true straight-line relations. Number of data points, slope and intercept of the true relation, and variances of the errors associated with the independent (X) and dependent (Y) variables influence the degree of agreement. For example, differences between the two line-fitting methods decrease as error in X becomes small relative to error in Y. Regarding the second goal—predicting the dependent variable—OLS is better than SA. Again, the difference diminishes as X takes on less error relative to Y. With respect to estimation of slope and intercept and prediction of Y, agreement between Monte Carlo results and large-sample theory was very good for sample sizes of 100, and fair to good for sample sizes of 20. The procedures and error measures are illustrated with two geologic examples.  相似文献   

Ductile deformation of garnet in mylonitic gneisses from the Münchberg Massif (Germany)   总被引：1，自引：0，他引：1  

Axel Vollbrecht  Jan Pawlowski  Bernd Leiss  Till Heinrichs  Madlen Seidel  Andreas Kronz 《Tectonophysics》2006,427(1-4):153

Mylonitic gneisses from the Münchberg Massif contain single grains (type I) and polycrystalline aggregates (type II) of garnet displaying a distinct elongation parallel to a macroscopic lineation which is interpreted as the result of ductile deformation. Lattice-preferred orientations of quartz (textures) symmetrical to the macroscopic foliation and lineation and the lack of rotational microfabrics indicate that the bulk deformation was pure shear at least during the latest strain increments. Garnet textures measured by EBSD together with microprobe analyses demonstrate that these two structural types of garnet can be related to two different processes of ductile deformation: (1) For the single grains stretching can be attributed to diffusion creep along grain boundary zones (Coble creep). The related mass transfer is indicated by the fact that primary growth zones are cut off at the long faces of the grains while the related strain shadow domains do not show comparable chemical zoning. Pressure solution and precipitation suitable to produce similar structures can be largely ruled out because retrogressive reactions pointing to the presence of free hydrous fluids are missing. (2) For the polycrystalline garnet aggregates consisting of cores grading into fine-grained mantles, dislocation creep and associated rotation recrystallization can be assumed. Continuous lattice rotation from the core to the outer polycrystalline rim allow a determination of the related dominant slip systems which are {100}<010> and equivalent systems according to the cubic lattice symmetry. The same holds for garnets which appear to be completely recrystallized. For this type of fine-grained aggregates an alternative nucleation model is discussed. Due to penetrative dislocation glide in connection with short range diffusion and the resulting lattice rotation, primary growth zones are strongly disturbed.Since for the considered rock unit of the Münchberg Massif peak metamorphic temperatures between 630 and 670 °C can be assumed, this study clearly demonstrates that the inferred processes of ductile garnet deformation can occur not only in HT regimes as often suggested in the literature even if embedded within a matrix of “low-strength” minerals like quartz, feldspars and micas.  相似文献