一种基于双重距离的空间聚类方法

传统聚类方法大都是基于空间位置或非空间属性的相似性来进行聚类,分裂了空间要素固有的二重特性,从而导致了许多实际应用中空间聚类结果难以同时满足空间位置毗邻和非空间属性相近。然而,兼顾两者特性的空间聚类方法又存在算法复杂、结果不确定以及不易扩展等问题。为此,本文通过引入直接可达和相连概念,提出了一种基于双重距离的空间聚类方法,并给出了基于双重距离空间聚类的算法,分析了算法的复杂度。通过实验进一步验证了基于双重距离空间聚类算法不仅能发现任意形状的类簇,而且具有很好的抗噪性。     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Adsorption of fluoride on goethite surfaces-implications on dental epidemiology

Fluoride ion interaction with synthetically prepared goethite has been investigated over a range of pH values (4–9) and F^– concentrations (10^–3–10^–5 M). The amount of F^– retained by goethite suspensions was found to be a function of pH, media ionic strength, F^– concentration, and goethite concentration. The lowest ionic strength (0.001 M KNO₃) gave the highest adsorption medium. Uptake was minimal at pH >7 and increased with decreasing pH. Thermodynamic properties for fluoride adsorption at 298 K and 323 K were investigated. The isosteric heat of adsorption (H _r) was calculated and the heterogeneity and homogeneity of the surface examined for goethite. In view of the importance of fluoride in dental health, the interaction of fluoride on goethite in the physical environment has important implications on dental epidemiology.     

江西省安福地区地热水化学特征研究

为研究江西安福地区水文地球化学特征及控制因素, 本文采集了15组样品, 采用水化学、同位素分析等方法进行研究。结果表明: 研究区地热水以Na-HCO3型水为主, 地下水以HCO3-Na·Ca型为主, 地表水由西南向东北从HCO3-Na·Ca向HCO3-Ca型水演化。水化学组分演化过程主要受岩石风化作用控制, 地层封闭性较差, 水中的Na+、HCO– 3、Sr2+来源于硅酸岩风化溶解。由稳定同位素特征可知, 研究区地热水补给来源为大气降水。研究区热储温度为49.8~101.4 ℃, 地热水循环深度为1 502.6~1 513.6 m。热水在沿断裂带上升过程中与浅层冷水发生混合, 其混合比例为76.1%~87.5%。研究成果为安福地区水循环演化提供依据, 有利于地热资源的合理开采与保护。     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

An examination of directional discontinuities and magnetic polarity changes around interplanetary sector boundaries usingE > 2 keV electrons

Past studies of interplanetary magnetic sector boundaries have been based on the assumption that one can determine the field polarities by comparing the field directions with those of the nominal Parker spiral angles. Previous investigators have found evidence for decreases of B, the magnitude of the magnetic fieldB, and increases of , the angle betweenB and the ecliptic plane, at sector boundaries. Others have argued that the characteristic thickness of sector boundaries exceeds that of tangential discontinuities, making sector boundaries a separate class of structures.We use a simple technique for inferring the polarities of interplanetary magnetic fields based on the assumption thatE > 2 keV electrons are always flowing along the magnetic field away from the Sun. Electron data from the UC Berkeley experiment on the ISEE-3 spacecraft are used to examine periods around several apparent sector boundaries in 1978 and 1979. We compare properties of (a) boundaries with field polarity changes and (b) large-angle ( > 60°) directional discontinuities with no field polarity changes. We find no significant differences between the sector boundaries and the directional discontinuities in terms of associated decreases in B or of values of . These results suggest no significant difference between sector boundaries and directional discontinuities other than the change in field polarities. Within limited statistics we find that about half the polarity changes would not have been identified using a requirement that > 90° and that half of the > 120° discontinuities would have been misidentified as polarity changes.     

Doppler velocity measurements made with a scanning photoelectric magnetograph

We discuss the problems connected with the measurements and evaluation of line-of-sight velocities, obtained with a scanning photoelectric magnetograph using a line-shifter with enhanced sensitivity. We bring arguments for the validity of the results of our photoelectric Doppler velocity recordings. We have found a network of cellularly shaped patterns in the distribution of photo-electrically measured line-of-sight motions, upflowing in the magnetically quiet (blue-shifted) and downflowing in magnetically active (red-shifted) areas of the photosphere, if the mean velocity level is estimated for a sufficiently large measured area. The features of both directions are mutually complementary. We demonstrate the effect of the shift of the reference zero velocity level on the topology of the line-of-sight velocity maps, and the dependence of this level on the size of the area from which it is estimated.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Comparison of structural and least-squares lines for estimating geologic relations

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.