离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

Crack damage stress as a composite function of porosity and elastic matrix stiffness in dolomites and limestones

Twenty-five uniaxial compression tests were performed to determine stress at onset of dilation, referred to herein as “the crack damage stress,” in heterogeneous dolomites and limestones. A simplified model for crack damage stress (σ_cd) is developed here using porosity, elastic modulus, Poisson's ratio and three empirical coefficients. The model shows that when porosity decreases and elastic modulus increases, σ_cd rapidly increases and approaches its maximum value. On the other hand, when porosity increases and elastic modulus decreases, σ_cd rapidly decreases and approaches its minimum value. The proposed model is validated for six heterogeneous limestone and dolomite formations which are widely distributed in Israel.     

Timing of Fracture Formation and Duration of Mineralization at the Hishikari Deposit, Southern Kyushu, Japan

Abstract. The Hishikari epithermal gold deposit consists of the Honko-Sanjin zone and Yamada zone. On the basis of observed crosscutting relationships of veins, vein system in the Hishikari deposit was classified into the early and late veins. They differ from orientation, width, control by unconformity, the amount of displacement and structural features, but have similarity in gold grade and Au/Ag ratios. K-Ar ages are presented for the classified veins to clarify timing of fracturing and duration of mineralization.
Age determination revealed that the early veins in the Honko-Sanjin zone range from 0.86 to 1.11 Ma and that the late veins range from 0.73 to 0. 84 Ma. The economically most important mineralization in the Hishikari deposit occurred around 0.90 Ma in a very short period (about 0.05 million years).
The distribution of mineralization age in the whole Hishikari deposit, occurrences of slickensides and ground water temperatures imply that mineralization lasted longer in the southern part of the deposit. As compared the early veins with late veins, structural changes including width, strike, dip and vein characteristics occurred at around 0.85 Ma.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

X线及CT检查在外伤性骨质疏松中的应用及讨论

目的 讨论CT及X线检查在外伤性骨质疏松中的应用。方法 1998－2000年200例外伤复查病人行CT及X光检查。结果 2－3个月后出现骨疏松的病人113例，3－4个月后出现骨疏松的66例，未见骨疏松的21例。结论 骨折2－3个月后复查病人，应适量降低CT的窗宽，窗位和X线的照射剂量。     

岩石圈流变机制的确定及影响岩石圈流变强度的因素

讨论了确定岩石圈流变强度中存在的问题，根据Byerlee定律，给出了在Anderson断层系统下三种断层的摩擦滑动强度公式，利用小标本实验结果外推，用大标本实验结果作约束，得到岩石圈几种典型岩石的脆性破裂规律，利用上述结果和传统的方法，分别得到了鄂尔多斯和山西裂谷两个典型地区的流变强度随深度的变化。结果表明，以往的计算对岩石圈流变强度的估计过高，对脆性形变区估计不足，流变机制估计不对；岩石圈有力学分层的特性，但各地分层的深度范围不同，使得成层和力学作用变得复杂，讨论了水、应变率以及多相矿物对流变强度的影响。     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

河北平原水系密度与隐伏活动构造的关系

以1926年出版的比例尺为1：5万的顺直地形图作为基础图件，用地理信息系统技术对河北平原作水系密度计量分析，揭示隐伏在第四纪巨厚的沉积物之下活动构造的信息，研究结果证实水系线密度统计图与使用其它地球科学方法获得的成果，例如基底活动断块图，地貌类型图和第四纪厚度图都有很好的对应，特别与布格重力图对应更好，研究结果发现水系密度图与布格重力异常图在总体一致性较好的前提下，多处反向异常区均为近代强震发生的区域，提出这种地表水系密度统计研究的成果不仅可以揭示深部隐伏的活动构造，而且具有较明确的时间含义。