“数字海底”--海底数据集成管理与可视化研究

针对目前"数字海底"建设中存在海底调查数据集成管理与三维可视化表达不足等问题,本文采用GIS技术,以黄河水下三角洲埕北海域为研究对象,利用地理空间数据集成理论与三维建模方法,建立了埕北海域三维海底空间数据库,实现了研究区域地形数据、地层数据、钻孔数据、表层沉积物数据的有效组织管理与可视化,并在此基础上,采用ArcGIS Engine 10.0,结合.NET平台,在Visual Studio 2010开发环境下,利用C#语言进行了二次开发,实现了基于C/S (Client/Server)架构的三维海底虚拟仿真系统的开发,设计了一套面向埕北海域的三维海底虚拟仿真原型系统。     

Reconstruction and analysis of sub-plinian tephra dispersal during the 1530 A.D. Soufrière (Guadeloupe) eruption: Implications for scenario definition and hazards assessment

The last magmatic eruption of Soufrière of Guadeloupe dated at 1530 A.D. (Soufrière eruption) is characterized by an onset with a partial flank-collapse and emplacement of a debris-avalanche that was followed by a sub-plinian VEI 2–3 explosive short-lived eruption (Phase-1) with a column that reached a height between 9 and 12 km producing about 3.9 × 10⁶ m³ DRE (16.3 × 10⁶ m³ bulk) of juvenile products. The column recurrently collapsed generating scoriaceous pyroclastic flows in radiating valleys up to a distance of 5–6 km with a maximum interpolated bulk deposit volume of 11.7 × 10⁶ m³ (5 × 10⁶ m³ DRE). We have used HAZMAP, a numerical simple first-order model of tephra dispersal [Macedonio, G., Costa, A., Longo, A., 2005. A computer model for volcanic ash fallout and assessment of subsequent hazard. Comput. Geosci. 31, 837–845] to reconstruct to a first approximation the potential dispersal of tephra and associated tephra mass loadings generated by the sub-plinian Phase 1 of the 1530 A.D. eruption. We have tested our model on a deterministic average dry season wind profile that best-fits the available data as well as on a set of randomly selected wind profiles over a 5 year interval that allows the elaboration of probabilistic maps for the exceedance of specific tephra mass load thresholds. Results show that in the hypothesis of a future 1530 A.D. scenario, populated areas to a distance of 3–4 km west–southwest of the vent could be subjected to a static load pressure between 2 and 10 kPa in case of wet tephra, susceptible to cause variable degrees of roof damage. Our results provide volcanological input parameters for scenario and event-tree definition, for assessing volcanic risks and evaluating their impact in case of a future sub-plinian eruption which could affect up to 70 000 people in southern Basse-Terre island and the region. They also provide a framework to aid decision-making concerning land management and development. A sub-plinian eruption is the most likely magmatic scenario in case of a future eruption of this volcano which has shown, since 1992, increasing signs of low-energy seismic, thermal, and acid degassing unrest without significant deformation.     

Na+,K+,Rb+,Cs+//SO2-4-H2O体系固液相平衡热力学模拟

锂云母、铯榴石等矿物的硫酸盐焙烧浸出,是从中提取锂铷铯钾等碱金属元素的重要步骤,所得硫酸盐浸出液可近似看作Li⁺,Na⁺,K⁺,Rb⁺,Cs+ // SO₄2--H₂O六元水盐体系;其在浓缩过程呈现复杂的成盐特征,增加了各碱金属元素分离提取的难度,研究其相关水盐体系相图和相平衡性质对理解浸出液中碱金属的分离具有重要的理论指导意义。本文采用Pitzer-Simonson-Clegg超额Gibbs自由能方程和具有内在热力学一致性的CALPHAD方法,构建了Na⁺,K⁺,Rb⁺,Cs⁺//SO₄2--H₂O五元体系的多温固液相平衡热力学模型。借助该模型:1)实现了对Na₂SO₄₊Rb₂SO₄₊H₂O、Na₂SO₄₊Cs₂SO₄₊H₂O和K₂SO₄₊Cs₂SO₄₊H₂O三元体系的多温溶解度等温线的准确描述,并进一步构建了这些体系在240.15~380.15 K温度范围内的多温相图;2)实现了对K₂SO₄+Rb₂SO₄₊H₂O和Rb₂SO₄₊Cs₂SO₄₊H₂O体系中存在的固溶体相—水溶液平衡关系的准确描述,确认了(K,Rb)₂SO₄连续固溶体、以及(Rb,Cs)₂SO₄(α)和(Rb,Cs)₂SO₄(β)间断固溶体存在的合理性;3)预测了上述五元体系中所包含的全部4个四元子体系(Na₂SO₄₊ K₂SO₄₊ Rb₂SO₄+H₂O、Na₂SO₄ + K₂SO₄₊Cs₂SO₄₊ H₂O、Na₂SO₄₊ Rb₂SO₄₊Cs₂SO₄ + H₂O和K₂SO₄₊Rb₂SO₄₊Cs₂SO₄₊ H₂O)的相图,明确了四元相图中固溶体的类型及形成区域。然而,模型预测结果与现有四元体系实验数据间存在普遍的不一致性,尤其是固溶体形成区域,但模型预测结果从非实验角度给出了重要参考,对现有实验提出了挑战,未来对这些体系开展更为精细的实验研究非常必要。     

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.     

Small-scale ductile shear zones as transposed rheologically weak domains: A numerical modeling investigation and practical application

We test the hypothesis that small ductile shear zones are developed from initial rheologically weak domains. We regard weak domains as ellipsoidal inhomogeneities and apply Eshelby's formalism extended for power-law viscous materials to investigate the kinematics and finite strain evolution of the partitioned flow fields in weak domains. We show that, under an imposed bulk flow field, weak domains, regardless of their initial shapes and orientations, deform into zone-like features at relatively small bulk strains and the deformation paths inside weak domains have all characteristics expected in ductile shear zones. We apply our model to the Cap de Creus area Spain, where abundant small ductile shear zones exist. To relate the fabrics inside these shear zones with the regional deformation, we take a multi-scale approach. We assume that the area contained many weak domains which were randomly shaped and oriented initially and were deformed into shear zones eventually in response to the regional flow. We constrain the regional flow field by the fabric patterns across the area and compute numerically the partitioned flow fields in individual weak domains. The latter are related to fabrics inside shear zones. Our model reproduces first-order features of field-observed fabrics. Although the deformation path of each shear zone was close to simple shearing, the deformation of the whole belt was not. Our approach also resolves the strain compatibility problem for a finite-sized shear zone embedded in a far less deformed country rock.     

Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada

 Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity–concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12%; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations. Received, September 1996 Revised, September 1997, February 1998 Accepted, February 1998     

Tsunami Inundation Modeling and Mapping using ALTM- and CARTOSAT-Derived Coastal Topographic Data

Coastal topography is the principal variable that affects the movement of the tsunami wave on land. Therefore, land surface elevation data are critical to a tsunami model for computing extent of inundation. Elevation data from India's remote sensing satellite CARTOSAT-1 are available for the entire Indian coastline, while elevation data collected using Airborne Laser Terrain Mapper (ALTM) are only available for selected sections of the coastline. This study was carried out to evaluate the suitability of CARTOSAT-1 and ALTM elevation data sets in the tsunami inundation modeling. Two areas of the coastal Tamil Nadu that were severely affected during the December 2004 tsunami and surveyed extensively for mapping the extent of inundation were selected as the study areas. Elevation data sets from ALTM, CARTOSAT-1 and field measurement collected using Real-time Kinematic GPS (RTK-GPS) were compared for these areas. The accuracy of ALTM and CARTOSAT-1 data, the significance of interpolation methods and data used on model outputs were studied. The analysis clearly revealed that the elevation accuracy of CARTOSAT-1 data (+/?2m) was much lower than ALTM data (+/?0.6m). However, it was found that despite the differing elevation accuracy, both ALTM and CARTOSAT-1 can be used to produce tsunami inundation maps for open coasts with an accuracy of 185 m (2 grid cells) at 75% and 50% confidence level, respectively.     

南粤古驿道历史建筑数字化保护与VR场景构建

南粤古驿道遗留了众多文化遗产,具有较高的历史、艺术、科学价值,但历史建筑损坏较为严重.为推进南粤古驿道活化利用,本文采用精细三维建模技术与虚拟现实(VR)技术实现了南粤古驿道的数字化保护.首先对南粤古驿道遗址进行了分析,然后从数据采集、点云数据处理、几何模型加工、纹理贴图等方面阐述对历史遗迹进行精细三维建模和VR场景构...     

Large to intermediate-scale aquifer heterogeneity in fine-grain dominated alluvial fans (Cenozoic As Pontes Basin,northwestern Spain): insight based on three-dimensional geostatistical reconstruction

Facies reconstructions are used in hydrogeology to improve the interpretation of aquifer permeability distribution. In the absence of sufficient data to define the heterogeneity due to geological processes, uncertainties in the distribution of aquifer hydrofacies and characteristics may appear. Geometric and geostatistical methods are used to understand and model aquifer hydrofacies distribution, providing models to improve comprehension and development of aquifers. However, these models require some input statistical parameters that can be difficult to infer from the study site. A three-dimensional reconstruction of a kilometer scale fine-grain dominated Cenozoic alluvial fan derived from more than 200 continuously cored, closely spaced, and regularly distributed wells is presented. The facies distributions were reconstructed using a genetic stratigraphic subdivision and a deterministic geostatistical algorithm. The reconstruction is only slightly affected by variations in the geostatistical input parameters because of the high-density data set. Analysis of the reconstruction allowed identification in the proximal to medial alluvial fan zones of several laterally extensive sand bodies with relatively higher permeability; these sand bodies were quantified in terms of volume, mean thickness, maximum area, and maximum equivalent diameter. These quantifications provide trends and geological scenarios for input statistical parameters to model aquifer systems in similar alluvial fan depositional settings.     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.