岩爆破坏过程能量释放的数值模拟

岩爆是地下岩石开挖中的一种工程灾害现象,是岩体结构发生破坏时,由于内部储存的弹性能释放并转换为动能而造成动力形式的破坏。岩爆破坏过程中的能量释放与岩体在应力峰值前后的应力-应变特性紧密相关。另外,施工中开挖速度引起的加载速率的变化也会对岩爆的产生有明显影响。以岩体全过程应力-应变曲线试验为基础分析岩爆破坏过程。分析中采用的模型考虑了岩石峰值后应力-应变特性及加载速率的影响。运用数值方法对岩石洞室的开挖过程进行了模拟,在模拟中对岩体破坏的发生及弹性能释放过程进行了分析。数值分析结果显示,岩体洞室开挖过程中岩石破坏由岩体表面向岩体内部发展,岩石的弹性能释放率也随着破坏的发展而不断增加。分析结果还显示,岩体破坏时的弹性能释放速率会随着开挖速率的提高而明显增加。     

地形对大地电磁测深(MTS)资料的影响

地形影响是山区ＭＴＳ工作的主要问题之一，对此做了理论分析并用均匀半空间表面二维地形模型进行了数值模拟。结果表明，地形影响不仅与ＭＴＳ资料的种类、极化模式、场源电磁波周期及测点位置有关，而且还依赖于不同的地形形态和尺度。此外，还根据悬崖模型计算，给出了关于地形影响范围的估     

一种导航卫星中长期轨道预报方法

基于神经网络混合建模的思想提出一种针对导航卫星的中长期轨道预报方法,在原动力学模型的基础上引入神经网络模型作为补偿,从而获得新的预报模型。在训练过程中神经网络通过学习动力学模型轨道预报误差来掌握其变化规律,并在预报过程中为动力学模型预报提供补偿,以提高预报精度。对GPS卫星动力学模型中长期预报误差的特点进行分析,然后根据所得结论提出混合模型的中长期(15 d以上)预报方案,最后通过对GPS卫星的仿真试验证明混合模型的改进效果,结果表明新方法在15~40 d的预报上表现出很好的改进效果。     

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性, 文中选用二次场算法, 将背景电阻率设置为水平层状且各向异性, 场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组, 通过Galerkin加权余量法形成大型稀疏有限元方程, 采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位; 进而, 利用滑动平均方法得到二次矢量位和标量位在空间的导数, 由此得到二次电磁场; 通过一维模型对算法的可靠性进行验证, 与此同时, 针对实际复杂海洋电磁模型, 比较有限元模拟结果与积分方程模拟结果, 进一步验证算法精度.若干计算结果均表明, 文中算法具有良好的通用性, 适用于井中电磁、航空电磁, 环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

基于正演模拟的复合曲流带砂体内部构型解剖研究及应用

Q油田是渤海大型河流相沉积的稠油油田,其北区NmⅣ1主力含油砂体为复合曲流带沉积,经过16年的注水开发,油水关系复杂,平面注水受效不均的问题日益严重,如何优化注釆结构来改善开发效果已成为急需解决的问题。优化现有注釆井网的关键是认识复合曲流带砂体内部不同单河道砂体的切叠关系及展布特征,即对目标砂体开展内部不同单河道的构型单元分布特征研究。首先,依据岩电标定的测井曲线特征、局部发育的薄夹层以及曲流河纵向下切模式,对复合砂体内部不同单河道砂进行划分与对比;然后基于单井沉积微相划分结果,开展井间剖面相对比分析,表明复合砂体内部存在5种不同单河道砂体切叠模式;建立5种正演概念模型,利用地震正演模拟合成地震记录,总结了不同单河道砂体间顶面高程差异、点坝边界、废弃河道等曲流河砂体构型边界的波形特征;利用地震属性提取及多种地质构型方法,综合分析,最终实现了海上大井距条件下复合曲流带砂体内部不同单一曲流河砂体构型边界的精细识别;最后由注采动态资料的响应关系评估不同的切叠模式对注水效果的侧向遮挡程度,并以此为依据,优化现有的注釆井网并取得了非常好的生产效果。     

自动建模客观分型方法在非汛期人工增雨中的应用

根据安徽省人工影响天气业务系统需要,利用NCEP 500 hPa高度场和风场格点资料及安徽省81站10年非汛期逐日降水资料,在客观分析天气形势基础上,结合预报员实践经验给出的不同系统影响安徽省时的关键区,计算机自动建模,形成逐日天气型日历,实现环流分型客观化。此方法识别条件实用,模型物理意义明确,自动化程度高,自动生成的逐日天气型日历,为安徽省非汛期人工增雨作业指挥系统提供了很好的天气预报参考依据。     

ACTIVE system for monitoring volcanic activity: A case study of the Izu-Oshima Volcano,Central Japan

A system is proposed for the monitoring of changes in the underground structure of an active volcano over time by applying a transient electromagnetic method. The monitoring system is named ACTIVE, which stands for Array of Controlled Transient-electromagnetics for Imaging Volcano Edifice. The system consists of a transmitter dipole used to generate a controlled transient electromagnetic (EM) field and an array of receivers used to measure the vertical component of the transient magnetic field at various distances, with automatic operation of both units. In order to verify the performance of the proposed system, numerical and field experiments were carried out by application of the system to the Izu-Oshima volcano, where a remarkable change in the apparent DC resistivity over time had been detected in association with the eruption in 1986.     

凹坑地形风流结构对污染物散布的模拟研究

利用建立的三维非静力高分辨率高阶湍流闭合模式与随机游动扩散模式研究了一个深凹露天矿区污染物散布的规律,同时在风洞中进行了示踪实验。结果表明,由于凹坑内复环流结构的存在,使得坑内污染物浓度较大,且浓度最大值出现在源的上风侧。数值试验与风洞试验结果吻合较好。     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.