山西义兴寨金矿床金矿物颗粒的产出及其成矿动力学意义

山西义兴寨金矿床的金主要呈包体金和裂隙金两种产出状态.文章根据成矿地球化学动力学理论,论证了金元素在含矿断裂发生脆性破裂和脆-韧性扩张活动中析出的化学反应类型,以及受反应动力学条件的制约,形成矿石中金不同产出状态的规律.指出热液矿床中普遍存在的包体金为热液成矿早期阶段断裂发生脆性破裂和成矿流体沸腾,导致快速反应结晶的产物;而裂隙金和晶隙金为成矿晚期金等矿质在流体中残余富集,同时成矿断裂再次发生脉动破裂,即成矿热力学演化与有利的构造动力学因素耦合作用的结果.得出金矿物颗粒的产出状态具有成矿动力学意义的结论.     

端头扩大型土锚的新发展 ——船锚式注浆张开型土锚

作为提高锚固性能的一种有效手段，端头扩大型土锚被广泛应用于实际工程。介绍了端头扩大型锚杆的各种扩孔方法，重点阐述一种新型的船锚式注浆张开型土锚。带有锚抓的新型土锚，不需要专门的扩孔机械，在注浆作用下能连续完成端头扩张与灌浆，具有施工方便、锚固面大、效果好等优点。通过对土锚进行张开试验，证明了在上海等软土地区通过灌浆压力使锚抓张开并形成扩大头是能够实现的。     

Ground water contaminant transport by nondivergence-free, unsteady and nonstationary velocity fields

Pore flow velocity is assumed to be a nondivergence-free, unsteady, and nonstationary random function of space and time for ground water contaminant transport in a heterogeneous medium. The laboratory-scale stochastic contaminant transport equation is up scaled to field scale by taking the ensemble average of the equation by using the cumulant expansion method. A new velocity correction, which is a function of mean pore flow velocity divergence, is obtained due to strict second order cumulant expansion (without omitting any term after the expansion). The field scale transport equations under the divergence-free pore flow velocity field assumption are also derived by simplifying the nondivergence-free field scale equation. The significance of the new velocity correction term is investigated on a two dimensional transport problem driven by a density dependent flow.     

Perspectives on present-day sea level change: a tribute to Christian le Provost

In this paper, we first discuss the controversial result of the work by Cabanes et al. (Science 294:840–842, 2001), who suggested that the rate of past century sea level rise may have been overestimated, considering the limited and heterogeneous location of historical tide gauges and the high regional variability of thermal expansion which was supposed to dominate the observed sea level. If correct, this conclusion would have solved the problem raised by the IPCC third assessment report [Church et al, Cambridge University Press, Cambridge, pp 881, 2001], namely, the factor two difference between the 20th century observed sea level rise and the computed climatic contributions. However, recent investigations based on new ocean temperature data sets indicate that thermal expansion only explains part (about 0.4 mm/year) of the 1.8 mm/year observed sea level rise of the past few decades. In fact, the Cabanes et al.’s conclusion was incorrect due to a contamination of abnormally high ocean temperature data in the Gulf Stream area that led to an overestimate of thermal expansion in this region. In this paper, we also estimate thermal expansion over the last decade (1993–2003), using a new ocean temperature and salinity database. We compare our result with three other estimates, two being based on global gridded data sets, and one based on an approach similar to that developed here. It is found that the mean rate of thermosteric sea level rise over the past decade is 1.5±0.3 mm/year, i.e. 50% of the observed 3 mm/year by satellite altimetry. For both time spans, past few decades and last decade, a contribution of 1.4 mm/year is not explained by thermal expansion, thus needs to be of water mass origin. Direct estimates of land ice melt for the recent years account for about 1 mm/year sea level rise. Thus, at least for the last decade, we have moved closer to explaining the observed rate of sea level rise than the IPCC third assessment report.     

Effect of fluid pressure on rock compressive failure in a nearly impermeable crystalline rock: Implication on mechanism of borehole breakouts

We conducted laboratory true triaxial experiments in the nearly impermeable Pohang rhyolite to investigate failure mechanisms under ‘dry’ and ‘wet’ rock conditions. Under ‘dry’ conditions prismatic specimens were jacketed all around to prevent confining fluid penetration. Under ‘wet’ conditions one pair of the specimen faces was left unjacketed and in direct contact with the confining fluid (kerosene) applying the least principal stress in an attempt to simulate the case of an unlined borehole wall. In both testing setups the true triaxial compressive strength for a given least principal stress increases significantly as the intermediate principal stress rises. The unjacketed rhyolite strength is, however, only 60 to 85% of the strength under dry conditions, depending on the magnitude of the intermediate principal stress. In dry rhyolite the failure process begins upon dilatancy onset, followed by microcrack localization, and ending in a steeply dipping shear fracture. On the other hand, brittle fracture in wet specimens occurs almost immediately after the onset of dilatancy by the development of one or more through-going extensile fractures subparallel and adjacent to one of the unjacketed faces, resembling the extensile cracks leading to borehole breakouts in crystalline rocks. We infer that upon dilatancy the confining fluid intrudes and quickly propagates newly opened stress-induced microcracks subparallel to the unjacketed faces, leading to ‘early’ failure.     

Modeling thermal deformation of VLBI antennas with a new temperature model

Temperature variations at very long baseline interferometry (VLBI) sites cause thermal deformations of the VLBI antennas and corresponding displacements of the VLBI reference points. The thermal deformation effects typically contain seasonal and daily signatures. The amplitudes of the annual vertical motion of the antenna reference point can reach several millimeters, depending on the design of the antenna structure, on the material, and on the environmental effects such as global station position, station height and climatology effects. Simple methods to correct this effect use the difference of the environmental temperature with respect to a defined reference temperature, the antenna dimensions, the elevation of the antenna, the material of antenna structure. Applying these simple models for thermal deformation in the VLBI data analysis improves the baseline length repeatability by 3.5%. A comparison of these simple models with local thermal deformation measurements at the antennas in Onsala and Wettzell show that the local measurements and the modeled corrections agree well when the temperature of the antenna structure is used, but agree less good when the surrounding air temperatures are used. To overcome this problem we present a method to model temperature penetration into the antenna structures, that allows to model thermal deformation effects that agree with the observed vertical deformation of the Onsala and Wettzell radio telescopes with a root mean square deviation of 0.07 and 0.13 mm, respectively. Possible implementations in the VLBI analysis are presented, and the definition of an adequate reference temperature is discussed.     

硫酸盐渍土在多次冻融循环时的盐胀累加规律

通过室内模拟试验,研究了硫酸盐渍土在历经多次冻融循环时的盐胀累加规律,得出了盐胀累加的不同类型,给出了累加盐胀率与循环次烽之间的关系式,进一步揭示了盐胀累加的 和影响因素。     

Simulation of seismicity due to fluid migration in a fault zone

Spatio-temporal variation of rupture activity is modelled assuming fluid migration in a narrow, porous fault zone formed along a vertical strike-slip fault in a semi-infinite elastic medium. The principle of the effective stress coupled to the Coulomb failure criterion introduces mechanical coupling between fault slip and the pore fluid. The fluid is assumed to flow out of a localized high-pressure fluid compartment in the fault at the onset of earthquake rupture. The duration of the earthquake sequence is assumed to be much shorter than the recurrence period of characteristic events on the fault. Both an earthquake swarm and a foreshock–main-shock sequence can be simulated by changing the relative magnitudes of the initial tectonic stress, pore fluid pressure, fracture strength and so on. When an inhomogeneity is introduced into the spatial distribution of fracture strength, high complexity is observed in the spatio-temporal variation of rupture activity. For example, the time interval between two successive events is highly irregular, and a relatively long quiescence of activity is sometimes observed in a foreshock–main-shock sequence. The quiescence is caused by the temporary arresting of rupture extension, due to an encounter with fault segments having locally high strengths. The frequency–magnitude statistics of intermediate-size events obey the Gutenberg–Richter relation. The calculations show the temporal variation of the b value during some foreshock sequences, and the degree of the change seems to depend on the statistical distribution of the fracture strength.     

Understanding of hydraulic properties from configurations of stochastically distributed fracture networks

A systematic investigation of the effect of configurations of stochastically distributed fracture networks on hydraulic behaviour for fractured rock masses could provide either quantitative or qualitative correlation between the structural configuration of the fracture network and its corresponding hydraulic behaviour, and enhance our understanding of appropriate application of groundwater flow and contaminant transport modelling in fractured rock masses. In this study, the effect of block sizes, intersection angles of fracture sets, standard deviations of fracture orientation, and fracture densities on directional block hydraulic conductivity and representative elementary volume is systematically investigated in two dimensions by implementing a numerical discrete fracture fluid flow model and incorporating stochastically distributed fracture configurations. It is shown from this investigation that the configuration of a stochastically distributed fracture network has a significant quantitative or qualitative effect on the hydraulic behaviour of fractured rock masses. Compared with the deterministic fracture configurations that have been extensively dealt with in a previous study, this investigation is expected to be more practical and adequate, since fracture geometry parameters are inherently stochastically distributed in the field. Moreover, the methodology and approach presented in this study may be generally applied to any fracture system in investigating the hydraulic behaviours from configurations of the fracture system while establishing a ‘bridge’ from the discrete fracture network flow modelling to equivalent continuum modelling in fractured rock masses. Copyright © 2007 John Wiley & Sons, Ltd.     

Approximate analysis for relative motion of satellite formation flying in elliptical orbits

This paper studies the relative motion of satellite formation flying in arbitrary elliptical orbits with no perturbation. The trajectories of the leader and follower satellites are projected onto the celestial sphere. These two projections and celestial equator intersect each other to form a spherical triangle, in which the vertex angles and arc-distances are used to describe the relative motion equations. This method is entitled the reference orbital element approach. Here the dimensionless distance is defined as the ratio of the maximal distance between the leader and follower satellites to the semi-major axis of the leader satellite. In close formations, this dimensionless distance, as well as some vertex angles and arc-distances of this spherical triangle, and the orbital element differences are small quantities. A series of order-of-magnitude analyses about these quantities are conducted. Consequently, the relative motion equations are approximated by expansions truncated to the second order, i.e. square of the dimensionless distance. In order to study the problem of periodicity of relative motion, the semi-major axis of the follower is expanded as Taylor series around that of the leader, by regarding relative position and velocity as small quantities. Using this expansion, it is proved that the periodicity condition derived from Lawden’s equations is equivalent to the condition that the Taylor series of order one is zero. The first-order relative motion equations, simplified from the second-order ones, possess the same forms as the periodic solutions of Lawden’s equations. It is presented that the latter are further first-order approximations to the former; and moreover, compared with the latter more suitable to research spacecraft rendezvous and docking, the former are more suitable to research relative orbit configurations. The first-order relative motion equations are expanded as trigonometric series with eccentric anomaly as the angle variable. Except the terms of order one, the trigonometric series’ amplitudes are geometric series, and corresponding phases are constant both in the radial and in-track directions. When the trajectory of the in-plane relative motion is similar to an ellipse, a method to seek this ellipse is presented. The advantage of this method is shown by an example.