Optimizing a layered and laterally constrained 2D inversion of resistivity data using Broyden's update and 1D derivatives

Dense profile-oriented resistivity data allows for 2D and 3D inversions. However, huge amounts of data make it practically impossible to do full 2D or 3D inversions on a routine basis. Therefore, a number of approximations have been suggested over the years to speed up computations. We suggest using a combination of Broyden's update on the Jacobian matrix with derivatives calculated using a 1D formulation on a parameterized 2D model of locally 1D layered models. The approximations bring down the effective number of 2D forward responses to a minimum, which again gives us the ability to invert very large sections. Broyden's update is not as useful with a parameterized problem as is the case with a smooth minimum structure problem that has been the usual application. 1D derivatives, however, seem to be very effective when initiating a full 2D solution with Broyden's update. We compare the different methods using two different kinds of data on two synthetic models and on two field examples. The most effective and reliable optimization combines 1D derivatives with a full 2D solution and Broyden's update. When using Broyden's update the Jacobian matrix needs to be reset every once in a while. We do this whenever the difference in data residual from the previous iteration is less than 5%. This combined inversion method reduces the computation time approximately a factor of 3 without losing model resolution.     

3-D finite element modeling for evolution of stress field and interaction among strong earthquakes in Sichuan-Yunnan region

Introduction Sichuan-Yunnan region is a major area with frequent strong earthquakes in Chinese mainland, especially the middle-southern segment of South-North Seismic Zone, where many strong earth-quakes occurred in history. In the past 30 years, Sichuan-Yunnan region has two seismically active periods: one is from Tonghai earthquake in 1970 to Longling-Songpan earthquake in 1976, the other is from Lancang earthquake in 1988 to now. During this two periods, the M=7.7 Tonghai, M=7.1 Dagua…     

遥感影像三维可视化研究

利用DEM高程资料实现遥感影像的三维显示可大大提高遥感监测产品的使用价值，而商用软件成本高操作繁琐普及困难。利用计算机图形学知识和Delphi编程语言，开发了一套遥感影像三维可视化软件。对在开发过程中遇到的有关DEM数据的裁剪与拼接、图像配准、地表模拟、透视变换、消隐处理、光照和阴影效果以及颜色变换等问题的解决方案和实现方法作了简要的介绍。     

青海对流云数值模拟分析

利用中国气象科学研究院三维对流云模式和2002年青海省河南县秋季外场试验取得的资料，进行了数值模拟试验。该地区秋季对流云降水主要为冷云降水，暖雨过程不易启动。降雨主要是由于霰落入暖层融化，雨水的蒸发是雨水减少的主要机制。霰在降水的产生中发挥了重要作用。霰的生成又与冰晶密切相关。冰晶是霰的主要来源，而且也是霰生长的主要因素。初始的霰粒主要由冰霰自动转化生成，而较少由雨滴冻结生成。霰胚通过收集过冷云水和冰晶与霰的碰并又促进了霰的进一步生长。冰晶的生成主要是由于自然冰核的核化，因此，自然冰核的数浓度对整个降水过程都有影响。霰是云中过冷水消耗的主要因素。     

经历剪切变形历史的岩石节理表面粗糙特性分析

基于两类岩石节理剪切试验，采用非接触式激光测距仪，实测了剪切前后节理表面粗糙数据，统计分析了岩石节理表面粗糙特性的三维粗糙参数的物理含义，并将这些参数与相应的二维参数进行了充分比较。结果表明，三维参数与二维参数间存在的近似线性关系，而多个三维参数相对比单一的二维参数更全面有效地描述了节理表面粗糙特性。     

基于GIS的孔隙水文地质层三维空间离散方法

针对自然界中孔隙水文地质层空间分布的不连续性与厚度分布的不均匀性,研究基于GIS的孔隙水文地质层三维空间离散实现的技术路线,提出基于GIS的孔隙水文地质层不规则六面体元的三维空间离散方法.该法不仅能最大限度地保证不规则六面体元中水文地质层类型的一元性,而且可充分利用GIS的空间分析与数据的自动提取功能,快速提取各个计算结点上空间位置坐标与各类计算参数,大大缩短水文地质模型空间离散与相关数据文件组织所需的时间,提高地下水三维有限差分数值模拟的时效性,具有较高的实用价值.     

吉林台水电站厂房高边坡稳定性及加固措施研究

吉林台水电站厂房后边坡开挖形成140～340m的高边坡,保证其长期运行稳定性至关重要。在结合地质资料对边坡的失稳模式进行初步判断的基础上,采用极限平衡分析方法对正常蓄水、暴雨影响、地震荷载等工况条件下边坡的整体抗滑稳定性进行分析,得出暴雨对边坡稳定影响较大的结论,据此提出以排水结合预应力锚索为主的边坡加固增稳措施,实践证明该措施是有效的。     

Three-dimensional forward and backward modelling of diapirism: numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin

Numerical studies of ductile deformations induced by salt movements have, until now, been restricted to two-dimensional (2D) modelling of diapirism. This paper suggests a numerical approach to model the evolution of three-dimensional (3D) salt structures toward increasing maturity. This approach is also used here to restore the evolution of salt structures through successive earlier stages. The numerical methodology is applied to study several model examples. We analyse a model of salt diapirs that develop from an initial random perturbation of the interface between salt and its overburden and restore the evolved salt diapirs to their initial stages. We show that the average restoration errors are less than 1%. An evolutionary model of a 2D salt wall loaded by a 2D pile of sediments predicts a decomposition of the salt wall into 3D diapiric structures when the overburden of salt is supplied by 3D synkinematic wedge of sediments. We model salt extrusion feeding a gravity current over the depositional surface and estimate an average rate of extrusion and horizontal velocity of salt spreading. Faulting of the overburden to salt overhangs initiates new secondary diapirs, and we analyse the growth of these secondary diapirs. We also study how lateral flow effects the evolution of salt diapirs. The shape of a salt diapir can be very different if the rate of horizontal flow is much greater than the initial rate of diapiric growth solely due to gravity. We discuss the applicability of the results of the models to the evolution of Late Permian salt structures in the Pricaspian basin (Russia and Kazakhstan). These structures are distinguishable into a variety of styles representing different stages of growth: salt pillows, diapirs, giant salt massifs, 2D diapiric walls and 3D stocks complicated by large overhangs. The different sizes, shapes and maturities of salt structures in different parts of the Pricaspian basin reflect areal differences in salt thickness and loading history. Our results suggest that the numerical methodology can be employed to analyse the evolution of all salt structures that have upbuilt through younger ductile overburdens.     

用GPS秒信号锁定高频振荡器的方法研究

对用远距离传输的GPS秒信号锁定本地高频晶体振荡器的方法进行了探讨。介绍了高精度的时间数字转换器TDC(Time to Digital Converter)和对滤波的方法进行了探讨。为了降低成本又能够满足设计要求的精度，还考虑了双D／A转换的方法。有关的硬件系统已经通过了调试和一系列的测试和实验，在测量和控制精度方面，能达到设计的要求。此外，提出了这个系统有待解决的问题，如双D／A转换的非线性问题，Kalman滤波抗野值及继续提高精度的问题。