一种基于面模型的巷道三维自动建模方法研究

研究了一种基于面模型的矿山巷道三维自动建模方法。巷道建模分解为网络建模、断面建模、巷道段建模以及空间交汇几个步骤。对建模过程涉及的各算法进行了分析,基采用C++编程实现整个建模过程,并总结了有待进一步研究的问题。     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation

Unstructured-grid models grounded on semi-implicit, finite-volume, Eulerian–Lagrangian algorithms, such as UnTRIM and ELCIRC, have enjoyed considerable success recently in simulating 3D estuarine and coastal circulation. However, opportunities for improving the accuracy of this type of models were identified during extensive simulations of a tightly coupled estuary–plume–shelf system in the Columbia River system. Efforts to improve numerical accuracy resulted in SELFE, a new finite-element model for cross-scale ocean modeling. SELFE retains key benefits, including computational efficiency of existing semi-implicit Eulerian–Lagrangian finite-volume models, but relaxes restrictions on grids, uses higher-order shape functions for elevation, and enables superior flexibility in representing the bathymetry. Better representation of the bathymetry is enabled by a novel, “localized” vertical grid that resembles unstructured grids. At a particular horizontal location, SELFE uses either S coordinates or SZ coordinates, but the equations are consistently solved in Z space. SELFE also performs well relative to volume conservation and spurious oscillations, two problems that plague some finite-element models. This paper introduces SELFE as an open-source code available for community use and enhancement. The main focus here is on describing the formulation of the model and on showing results for a range of progressively demanding benchmark tests. While leaving details to separate publications, we also briefly illustrate the superior performance of SELFE over ELCIRC in a field application to the Columbia River estuary and plume.     

A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

苏皖下扬子区晚白垩世以来的构造-热历史：浦口组砂岩磷灰石裂变径迹证据

对苏皖下扬子区上白垩统浦口组（K2P）的三个砂岩样品进行了磷灰石裂变径迹（AFT）研究。结果显示：2个样品的AFT合并年龄为（88．8±4．4）Ma（径迹长度为（12．0±0．3）μm）和（82．1±6．8）Ma（径迹长度为（14．4±0．3）μm）与浦口组沉积年龄相近,说明它们沉降的幅度达到但没有超过AFT部分退火区间,1个样品的AFT合并年龄（117．3±5．9）Ma（径迹长度为（13．3±0．3）μm）大于浦口组沉积年龄,代表物源区抬升、剥露的冷却年龄。根据热历史模拟结果,识别出黄桥事件（110～90Ma）、仪征事件（70～60Ma）和三垛事件（35～22Ma）三期重要的构造事件,并将下扬子区晚白垩世以来的盆地演化划分为四个阶段：110～70Ma断坳复合型伸展盆地、70～35Ma拉张断陷盆地、35～22Ma挤压抬升阶段和22Ma至今坳陷盆地。     

The Structure of A Turbulent Jet in A Crossflow-Effect of Jet-Crossflow Velocity

A comprehensive numerical study on the three-dimensional structure of a turbulent jet in crossflow is performed. The jet-to-crossflow velocity ratio (R) varies in the range of 2 - 16; both vertical jets and inclined jets without excess streamwise momentum are considered. The numerical results of the Standard two-equation k-ε model show that the turbulent structure can be broadly categorised according to the jet-to-crossflow velocity ratio. For strong to moderate jet discharges, i.e. R> 4, the jet is characterized by a longitudinal transition through a bent-over phase during which the jet becomes almost parallel with the main freestream, to a sectional vortex-pair flow with double concentration maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless distance of around 20-60. The similarity coefficients are only weakly dependent on R. The cross-section scalar field is kidney-shaped and bifurcated, vvith distinct double concentr     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.

     

A numerical study of the behavior of fully grouted rockbolts under dynamic loading

The widespread preference for rockbolts in providing adequate rock stability in underground mines and man-made cavities make it necessary to obtain a better understanding of the response of these support systems. It is widely accepted that the fully grouted rockbolts provide better roof stability in areas of very poor roof condition which may be caused by high ground stress conditions. There is little information about the in situ behavior of these systems especially under dynamic loading. Hence, to study the behavior of fully grouted rockbolts under dynamic loading, a numerical modeling study was conducted using the FLAC^3D code. In this study the behavior of three types of fully grouted rockbolts were compared with each other including rockbolts with and without head plate and a yielding type one. The results of analyses indicated that under dynamic loading the fully grouted rockbolts without the head plate are incapable of controlling the rock mass movement. Although fully grouted rockbolts with head plate damp a considerable amount of the dynamic energy through friction as these bolts slide within the grout, but the elongation of the rockbolt is not possible due to the rapid breakage of the rockbolt–grout interface. Yielding rockbolts are the best choice for the absorption of the dynamic stress wave and controlling of the rock mass movement. The obtained results show that the optimal design of yielding rockbolt should be in such a way that after dynamic loading, anchoring part of bolt has a limited movement to prevent stress concentration in the shaft of rockbolt and its breakage.     

Predicting Diel,Diurnal and Nocturnal Dynamics of Dissolved Oxygen and Chlorophyll‐a Using Regression Models and Neural Networks

Human‐induced and natural interruptions with continuous streams of observational data necessitate the development of gap‐filling and prediction strategies towards better understanding, monitoring and management of aquatic systems. This study quantified the efficacy of multiple non‐linear regression (MNLR) versus artificial neural network (ANN) models as well as the temporal partitioning of diurnal versus nocturnal data for the predictions of chlorophyll‐a (chl‐a) and dissolved oxygen (DO) dynamics. The temporal partitioning increased the predictive performances of the best MNLR models of diurnal DO by 45% and nocturnal DO by 4%, relative to the best diel MNLR model of diel DO ($r_{{\rm adj}}^{2} = 68.8\%$ ). The ANN‐based predictions had a higher predictive power than the MNLR‐based predictions for both chl‐a and DO except for diurnal DO dynamics. The best ANNs based on independent validations were multilayer perceptron (MLP) for diel chl‐a, generalized feedforward (GFF) for diurnal and nocturnal chl‐a, MLP for diel DO, GFF for diurnal DO, and MLP for nocturnal DO.     

台湾海峡及其邻区现今构造应力场数值模拟

本文将台湾海峡及其邻区划分为多个具有不同弹性参数的等厚区块,并利用基于弹性理论的有限元数值模型,以GPS观测速度和由震源机制解(FMS)得到的主应力轴方向为约束条件,计算了该研究区的现今构造应力场.结果表明:(1)最大主应力轴在台湾海峡中部和台湾岛弧中部大致为NW-SE向,在东北部呈顺时针旋转,而在南部呈逆时针旋转;(2)研究区水平位移场大致以23°N为界,南强北弱,且在北部呈顺时针旋转,南部呈逆时针旋转;(3)东沙-澎湖-北港隆起因强度大于周边盆地区而成为"阻挡带",导致在其南北出现两个反向的构造逃逸区,并且隆起东南区形变特别强烈,而其西北区则相对较弱.