Conditions for instabilities in collapsible solids including volume implosion and compaction banding

We review conditions for material instabilities in porous solids induced by a bifurcation of solution into non-unique strain rate fields. Bifurcation modes considered include jumps in the strain rate tensor of ranks one and higher representing deformation band and diffuse instability modes, respectively. Eigenmodes (e-modes) are extracted for each type of instability to fully characterize various frameworks of deformation in collapsible solids. For diffuse instability these e-modes are determined from a homogeneous system of linear equations emanating from the condition of zero jump in the stress rate tensor, which in turn demands that the tangent constitutive tensor be singular for the existence of nontrivial solutions. For isotropic materials we describe two types of singularity of the constitutive tensor: (a) singularity of the constitutive matrix in principal axes, and (b) singularity of spin. Accordingly, we derive the e-modes for each type of singularity. We utilize the singularity of the constitutive matrix in principal axes as a precursor to volume implosion in collapsible solids such as loose sands undergoing liquefaction instability and high-porosity rocks undergoing cataclastic flow. Finally, we compare conditions and e-modes for volume implosion and compaction banding, two similar failure modes ubiquitous in granular soils and rocks.Supported by U.S. Department of Energy, Grant DE-FG02-03ER15454, and U.S. National Science Foundation, Grants CMS-0201317 and CMS-0324674.     

体元结构模型在三维地震模型正演模拟研究中的应用

本文采用双三次康氏（Coons）曲面模拟地震界面,提出了用体元结构来表示复杂构造地区的三维地震模型,同时以此为基础建立了灵活的动态射线追踪过程,实现了快速的、更具实用性的三维地震模型的正演模拟．     

频域激发极化法中体极化与面极化内在联系的探讨

矿化的体极化介质,其激发极化效应来自其中所含电子导电矿物颗粒的表面极化效应。基于上述认识,本文从理论上推导出含球状电子导体的体极化介质复电阻率的“微观”数学表达式。将它与体极化介质复电阻率的“宏观”表达式--Cole-Cole模型相比较,得出了:1.Cole-Cole模型中的特征参量与球状电子导体的表面极化系数、含量以及导体半径等的综合关系式;2.电子导体表面极化系数在频域中的一般表达式。这表达式同目前在实验室中观测到的结果是一致的。     

地热资源评价的蒙特卡罗法

目前国内外评价地热资源的方法 ,均将热储参数作为定值而未考虑其随机特性。蒙特卡罗法是以热储参数的随机特征为基础建立的 ,并且通过计算可以给出结果的可信度。本文在论述蒙特卡罗法基本原理的基础上 ,结合开封市地热系统 ,分别利用蒙特卡罗法和单元热储体积法对地热资源进行了计算 ,经对比分析得出的地热资源量为 133 0 95 5×10 1 3 Kcal,相当于 1 5 476× 10 1 2 度电或 7 7381× 10 8t煤所产生的热量 ,为开封市地热资源的合理开发提供了科学依据。     

三维医学图象处理与分析系统

本文介绍了一个我们自主开发的基于微机的三维医学图象处理与分析系统。系统运行于Ｗｉｎｄｏｗｓ ９５／Ｗｉｎｄｏｗｓ ＮＴ环境,主要包括图象预处理,切片重组,三维显示,三维体显示,手术模拟等功能。可广泛用于ＣＴ,ＭＲＩ等医学图象的处理与分析,具有很好的应用前景。     

Optimizing the dataset size of a topo-bathymetric survey for Hammam Debagh Dam,Algeria

The capacity of water reservoirs may be considerably reduced by sediment transport and accumulation. This phenomenon may occur quickly in regions prone to erosion, such as Hammam Debagh Dam region in Algeria. For better management of water resources, it is important to estimate the volume of a reservoir by regular and frequent topo-bathymetric surveys. Topo-bathymetric surveys with a high point measurement density are expensive and time-consuming. To reduce this expense and to increase the frequency of topo-bathymetric surveys, it is required to optimize their size and density. The complete survey of the reservoir created by the Hammam Debagh Dam in Algeria was examined in comparison with 19 reduction subsets using geostatistics. The study showed that the complete dataset of 32,164 points, corresponding to a density of 42 points ha^-1 could be reduced to 10,000 points, i.e. 13 points ha^-1, without affecting the estimation of the water volume in the reservoir. This reduction could considerably simplify the work load, reduce the expenses to do these surveys, and, therefore, allow for an increased frequency of topo-bathymetric surveys.     

Experimental and numerical study of current-induced artificial upwelling

The aim of the current paper is to investigate hydrodynamic characteristics of the artificial upwelling induced by ocean currents. Experiments were performed in a flume at different density difference heads, horizontal current velocities and upwelling pipe diameters. A three-dimensional computational fluid dynamics (CFD) model was employed on wider range of parameters for further analysis. The performance of the numerical model has been confirmed by the experimental findings. The present results show that the volume flow rate of current-induced artificial upwelling is influenced by geometrical parameters and inclination angle of the pipe, the horizontal current velocity and vertical distribution of water density. In ideal two-layer density stratified water, the critical current velocity to generate upwelling linearly increases with the increase of the density difference, and the maximum rising height for upwelling is inversely proportional to the density difference. Feasibility analysis was taken by using current and density profiles of the East China Sea near Dongji Islands, which provides an useful reference for engineering practice.     

Water masses and currents of deep circulation southwest of the Shatsky Rise in the western North Pacific

We conducted full-depth hydrographic observations in the southwestern region of the Northwest Pacific Basin in September 2004 and November 2005. Deep-circulation currents crossed the observation line between the East Mariana Ridge and the Shatsky Rise, carrying Lower Circumpolar Deep Water westward in the lower deep layer (θ<1.2 °C) and Upper Circumpolar Deep Water (UCDW) and North Pacific Deep Water (NPDW) eastward in the upper deep layer (1.3–2.2 °C). In the lower deep layer at depths greater than approximately 3500 m, the eastern branch current of the deep circulation was located south of the Shatsky Rise at 30°24′–30°59′N with volume transport of 3.9 Sv (1 Sv=10⁶ m³ s⁻¹) in 2004 and at 30°06′–31°15′N with 1.6 Sv in 2005. The western branch current of the deep circulation was located north of the Ogasawara Plateau at 26°27′–27°03′N with almost 2.1 Sv in 2004 and at 26°27′–26°45′N with 2.7 Sv in 2005. Integrating past and present results, volume transport southwest of the Shatsky Rise is concluded to be a little less than 4 Sv for the eastern branch current and a little more than 2 Sv for the western branch current. In the upper deep layer at depths of approximately 2000–3500 m, UCDW and NPDW, characterized by high and low dissolved oxygen, respectively, were carried eastward at the observation line by the return flow of the deep circulation composing meridional overturning circulation. UCDW was confined between the East Mariana Ridge and the Ogasawara Plateau (22°03′–25°33′N) in 2004, whereas it extended to 26°45′N north of the Ogasawara Plateau in 2005. NPDW existed over the foot and slope of the Shatsky Rise from 29°48′N in 2004 and 30°06′N in 2005 to at least 32°30′N at the top of the Shatsky Rise. Volume transport of UCDW was estimated to be 4.6 Sv in 2004, whereas that of NPDW was 1.4 Sv in 2004 and 2.6 Sv in 2005, although the values for NPDW may be slightly underestimated, because they do not include the component north of the top of the Shatsky Rise. Volume transport of UCDW and NPDW southwest of the Shatsky Rise is concluded to be approximately 5 and 3 Sv, respectively. The pathways of UCDW and NPDW are new findings and suggest a correction for the past view of the deep circulation in the Pacific Ocean.     

A GIS-based methodology for hazard mapping of small volume pyroclastic density currents

We present here a methodology implemented within a geographical information system (GIS) for hazard mapping of small volume pyroclastic density currents (PDCs). This technique is implemented as a set of macros written in Visual Basic for Applications (VBA) that run within GIS-software (i.e. ArcGIS). Based on the energy line concept, we calibrated an equation that relates the volume (V) and the mobility (ΔH/L) of single PDCs using data from Soufrière Hills volcano (Montserrat) and Arenal volcano (Costa Rica). Maximum potential run-outs can be predicted with an associated uncertainty of about 30%. Also based on the energy line concept and with data from Soufrière Hills volcano and Mt. St. Helens (USA), we were able to calibrate an equation that predicts the flow velocity as a function of the vertical distance between the energy line and the ground surface (Δh). Velocities derived in this way have an associated uncertainty of 3 m s⁻¹. We wrote code to implement these equations and allow the automatic mapping of run-out and velocity with the inputs being (i) the height and location of the vent (ii) the flow volume and (iii) a digital elevation model (DEM) of the volcano. Dynamic pressure can also be estimated and mapped by incorporating the density of the pyroclastic density current (PDC). This computer application allows the incorporation of uncertainties in the location of the vent and of statistical uncertainties expressed by the 95% confidence limits of the regression model. We were able to verify predictions by the proposed methodology with data from Unzen volcano (Japan) and Mayon volcano (The Philippines). The consistencies observed highlight the applicability of this approach for hazard mitigation and real-time emergency management.     

现代测量数据绘制技术研究

介绍完成几何数据至图像数据的转变的有关算法 ,并讨论了体绘制技术的理论和方法 ,为最终实现测量数据可视化提供了技术基础。