A 2D finite volume model for bebris flow and its application to events occurred in the Eastern Pyrenees

FLATModel is a 2D finite volume code that contains several original approaches to improve debris-flow simulation. Firstly, FLATModel incorporates a ＂stop-and-go＂ technique in each cell to allow continuous collapses and remobilizations of the debris-flow mass. Secondly, flow velocity and consequently yield stress is directly associated with the type of rheology to improve boundary accuracy. Thirdly, a simple approach for entrainment is also included in the model to analyse the effect of basal erosion of debris flows. FLATMODEL was tested at several events that occurred in the Eastern Pyrenees and simulation results indicated that the model can represent rather well the different characteristics observed in the field.     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

某电厂氯化物在地下水中迁移规律研究

为了研究污染物在地下水中的迁移规律,本文以某电厂地下水中氯化物作为预测因子进行数值模拟。对研究区含水层物质组成、边界条件、初始条件及源汇项等做简化,采用GMS软件建立水文地质模型,在简化的水文地质模型的基础上开展地下水中污染物迁移规律研究。研究结果表明：在污染物排放规律简化为连续恒定排放的点源工况下,（1）污染物平面迁移范围以泄漏点为中心呈椭圆形扩散,污染晕长轴为地下水流动方向;（2）污染物最大迁移距离及污染影响范围均随时间不断增加,但是变化速率不断降低;（3）距离污染源较近的位置,污染物浓度随时间不断增加,但是浓度变化速率逐渐降低;距离污染源较远的位置,污染物浓度随时间不断增加,初期浓度变化速率随时间不断增加,后期浓度变化速率随时间逐渐降低。     

Footings under seismic loading: Analysis and design issues with emphasis on bridge foundations

The paper provides state-of-the-art information on the following aspects of seismic analysis and design of spread footings supporting bridge piers: (1) obtaining the dynamic stiffness (“springs” and “dashpots”) of the foundation; (2) computing the kinematic response; (3) determining the conditions under which foundation–soil compliance must be incorporated in dynamic structural analysis; (4) assessing the importance of properly modeling the effect of embedment; (5) elucidating the conditions under which the effect of radiation damping is significant; (6) comparing the relative importance between kinematic and inertial response. The paper compiles an extensive set of graphs and tables for stiffness and damping in all modes of vibration (swaying, rocking, torsion), for a variety of soil conditions and foundation geometries. Simplified expressions for computing kinematic response (both in translation and rotation) are provided. Special issues such as presence of rock at shallow depths, the contribution of foundation sidewalls, soil inhomogeneity and inelasticity, are also discussed. The paper concludes with parametric studies on the seismic response of bridge bents on embedded footings in layered soil. Results are presented (in frequency and time domains) for accelerations and displacements of bridge and footing, while potential errors from some frequently employed simplifications are illustrated.     

下地壳流变与造山带同挤压期地壳伸展的动力学关系

应用不可压缩非牛顿粘性流体的本构关系,对造山带同挤压期下地壳流变及其与上地壳构造伸展的动力学关系进行了二维有限元数值模拟. 结果表明,在板块侧向挤压下,当造山带山根下陷和地表隆起达一定程度后,地壳不同层圈岩石将发生复杂的粘性流变. 流变的运动学方式和分布范围不仅与时间有关、同时还受地壳厚度转变带形态的制约. 在构造挤压和山体荷载达到弹性平衡状态后,地壳流变首先发生在造山带下地壳山根,但经一定的Maxwell时间后,流变将不断局限于造山带前缘的厚度转变带. 这一流变方式的变化是导致造山带浅部地壳动力学转变的主要原因. 它造成造山带内上地壳最小主应力从近水平挤压不断转化为近水平拉张,由此使造山带前陆发生挤压冲断的同时,山体的核部发生上地壳的拉张伸展. 最后,应用这一结果讨论了青藏高原南缘南北向地壳伸展的动力学性质.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

Prediction of Water Temperature as Affected by a Pre‐Constructed Reservoir Project Based on MIKE11

Reservoir construction can lead to much more water stored in front of the dam and significantly increase heat storage capacity of the reservoir waters, thus resulting in different distribution pattern of water temperature in reservoir area compared to river. Especially for large reservoir, the obvious stratification of water temperature will appear in the reservoir with deeper water levels. Meanwhile, the low water temperature will be observed in the downstream river due to the operation of the reservoir. The vertical numerical simulation model for reservoir from MIKE 11 was used to predict the changes of water temperature of Wuxikou Reservoir to check the effects of the reservoir construction on water temperature. The water temperature prediction model was developed to simulate the water temperature of the reservoir and the discharged outflow water. The predicted results can contribute to assessing the feasibility of the pre‐constructed project based on the environmental influence of water temperature.     

IMU/DGPS辅助空三下几种布点方案的精度比较

主要介绍了在IMU/DGPS辅助空三的情况下,通过对不同的布点方案进行空三加密,对这些空三成果进行精度比较与分析,分析POS辅助空三在无地面控制点或少量的地面控制点的情况下能获得怎样的精度,能满足多大比例尺的成图要求。     

Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration

Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J₂ term. Excluding the J₂ term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum.