密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

EPR investigation of the methyl radical,the hydrogen atom and carbon oxide radicals in Maxixe-type beryl

The EPR spectra of Maxixe-type beryl contain a large number of overlapping signals. The angular dependence of the 1:3:3:1 signal typical for the CH₃ radical shows that this radical is located at the center of the channel cavity with its symmetry axis parallel to the crystal c-axis and is rotating around this axis. Its EPR spectrum is axially symmetric with g _// = 2.00263, g _⊥ = 2.00249 and A_// = 2.288 mT, A_⊥ = 2.256 mT. These anisotropies have the opposite signs of those found for surface-adsorbed methyl radicals. Hydrogen atoms are located at position 2a at the center of the beryl cavity and the EPR parameters of the narrow doublet signal are A ₀ = 1,407 MHz and g = 2.00230. Another doublet signal, which is broader and has axial symmetry with g _// = 2.00265, g _⊥ = 2.00625 and A_// = 0.895 mT, A_⊥ = 0.885 mT, could come from a HCO₃ radical. One narrow and easily saturated signal with g _// = 2.00227 and g _⊥ = 2.00386 is interpreted to arise from a carbon monoxide radical in the beryl channel, oriented with its axis parallel to the crystal c-axis. Additional weak doublet lines, which have similar g values as the carbon monoxide radical, are created by nearby hydrogens. A powder spectrum with g _// = 2.0017 and g _⊥ = 2.0004 appears upon UV irradiation of the single crystal and is easily saturated. This spectrum is interpreted to arise from a carbon dioxide radical, which rotates around its symmetry axis.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

滑坡体三维地质建模与可视化分析

针对滑坡地质研究的自身特点,提出了面向滑坡地质体三维建模的NURBS-BRep混合数据结构和地质结构单元实体构造技术.通过对滑坡基础地质数据的预处理、滑坡数字地形和滑动面的三维建模、地质结构面的系统构造、地质结构体的生成和显示,形成了一套完整的滑坡三维地质模型的构建方法.将该方法应用于某水库滑坡,建立了相应的三维滑坡地质模型,并基于此模型进行了三维剖切分析、数字钻孔、等值线生成、滑块自动剖分、滑坡失稳可视化动态模拟和滑坡体方量精确计算等一系列实用的可视化分析,为滑坡稳定性的准确计算和客观评价提供了有力的支持.     

一种新的波浪变形三维数值模式──0-1混合型边界元

为提高边界元法模拟三维波场波浪变形的数值计算精度,借鉴常数元和线性元剖分方式、波势函数及波势函数法向导数对单元节点设定的各自适应性,提出了一种新的单元剖分模式──0-1混合型边界元,以控制和减缓由于计算误差累计而造成的波浪数值计算上的“横向振动”,借此结合边界元法的分区模式可实现较大范围的波场线性波浪变形计算,并为时域内的波浪非线性变形计算提供时间步长的数值保证.     

渤海三维海洋温度和海流数值预报

"十·五"期间,我国开展了三维海洋温度和海流数值预报的业务化研究工作.经过3年的努力,渤海三维海洋温度和海流数值预报系统研制完成,并于2003年10月,开始试预报.本文对该预报系统以及运行情况进行了介绍,并分析了所存在的问题和发展方向.