Finite cube elements method for calculating gravity anomaly and structural index of solid and fractal bodies with defined boundaries

In this paper, we present the finite cube elements method (FCEM); a novel numerical tool for calculating the gravity anomaly g and structural index SI of solid models with defined boundaries and variable density distributions, tilted or in normal position (e.g. blocks, faulted blocks, cylinders, spheres, hemispheres, triaxial ellipsoids). Extending the calculation to fractal objects, such as Menger sponges of different orders and bodies defined by polyhedrons, demonstrates the robustness of FCEM. In addition, approximating the cube element by a sphere of equal volume makes the calculation of gravitation and related derivatives much simpler. In gravity modelling of a sphere, cubes with edges of 100 m and 200 m achieve a good compromise between running time and overall error.
Displaying the distribution of SI of the studied models on contour maps and profiles will have a strong impact on the forward and inverse modelling of potential field data, especially for Euler deconvolution.
For Menger sponges, plots of gravity elements g and its derivatives show similar patterns independent of fractal order. Moreover, both the pattern and magnitude of SI are independent of fractal order, allowing the use of SI as a new invariant measure for fractal objects. However, SI pattern and magnitude strongly depend on the depth to the buried bodies as do other elements
In this study, we also present a new type of plot; the structural index against distance variation diagrams from which we extract the three critical SI ( CSI ) values, one per axis. The inversion of gravity anomaly data at CSI values gives the optimal mean location of the buried body.     

Instantaneous deformation from continuous GPS: contributions from quasi-periodic loads

Continuous GPS (CGPS) coordinate time-series are known to experience repeating deformation signals with seasonal and other periods. It is unlikely that these signals represent perfect sinusoids with temporally constant amplitude. We develop an analysis method that accommodates temporal variations in the amplitudes of sinusoidal signals. We apply the method to simulated coordinate time-series to numerically explore the potential consequences of neglecting decadal variation in amplitude of annual motions on the residual-error spectra of CGPS measurements, as well as potential bias in estimates for secular site velocity. We find that secular velocity bias can be appreciable for shorter time-series, and that residual-error time-series of longer duration may contain significant power in a broad band centred on semi-annual frequency if temporal variation in the amplitude of annual motions is not accounted for in the model used to reduce the observations to residuals. It may be difficult to differentiate the bandpass filtered signature of mismodelled loading signals from power-law noise, using residual-error spectra for shorter time-series. We provide an example application to a ∼9-yr coordinate time-series for a CGPS station located in southern California at Carbon Creek Control Structure (CCCS), which is known to experience large amplitude seasonal motions associated with the Santa Ana aquifer system.     

一种机载雷达点云数据的快速分类方法

机载激光雷达扫描技术可以以点云形式快速获取地形表面高精度三维信息。基于激光雷达扫描数据及建筑物本身的拓扑信息就可以对建筑物进行精确的重建,而重建中最关键的技术是对点云数据进行分类,进而进行地物识别。对大规模三维点云数据进行快速分类,提出一种采用区域分割结合基于最小二乘平差的多项式拟合方法,将大量离散的三维点云分割后进行多项式拟合,并将二维数据分类转化为一维数据分类。在分类的基础上,将建筑物几何规则作为约束条件提取了房屋边缘。实验分析表明,该方法既能去除多余噪声,又能有效保留特征点,分类的总误差率低于3%。     

基于尺度律的裂隙网络生成及其对地下厂房稳定性的影响分析

以锦屏二级为工程背景,通过对猫猫滩闸址区域内断层的统计,得到断层数量与尺度的幂律关系,并通过分段线性拟合的方法估算了断层系统的尺度-数量分维发现,当断层迹长 0.1 km时,分维 0.487 5;当断层迹长 0.1 km时,分维 1.496 1。通过该尺度律关系,结合蒙特卡洛法,推演了更小规模的节理的空间分布。运用DDA数值计算程序,分析并模拟了大水沟地下厂房的开挖过程,获得了高度节理化的由开挖引起的应力、位移变化图。     

基于间断位移场的隧道围岩承载能力研究

岩体是岩石块体和分割这些块体的边界的组合体。沿边界面(结构面)的局部剪切变形使得位移场产生间断,对岩体的力学性能产生巨大的影响。采用光滑函数组合来描述滑移面上的间断,得到反映岩体沿着边界局部剪切的数学模型。利用此局部剪切模型得出了应力–滑移关系具有强化阶段时圆形坑道围岩自身承载能力的解析表达式,研究了岩体滑移面对坑道围岩自身承载能力的影响。通过分析得到,自承能力在深部坑道围岩的稳定中发挥了重要作用,应力–滑移关系的强化段对围岩承载力有着很大的影响;围岩压力与坑道的埋深、坑道的尺寸以及围岩的物理力学性质等具有直接关系;围岩出现不稳定的动力现象(岩爆)与此岩体参数m关系极大,而岩体参数与接触面的软化程度、能量转移、耗散结构的几何尺寸以及弹性模量有关。     

浙江省妙西多金属矿床成矿模型及找矿标志

笔者通过大量地质、地球化学、球物理和矿床（点）资料研究,提出浙江省妙西多金属矿床成矿模型,该模型正在为普查项目所证实。妙西多金属矿矿床是典型的岩浆期后热液矿床,矿体主要受本区的三级构造即岩体接触带构造和地层不连续界面构造控制。地球物理测量提供的深部信息和地质构造标志是本区多金属矿直接找矿标志,地表物化探资料在普查工作中可以起辅助作用。     

节理峰值后归一化位移软化模型

切向荷载作用下贯通节理岩体峰值强度后的剪切应力、位移变形规律一般呈双曲线形式发展,即峰值后强度衰减速率逐步变小而趋于0,到达残余强度。直剪试验表明,非贯通节理岩体峰值后剪切应力、位移呈“S”型趋势发展,即强度衰减速率先是逐步增大,到达一定程度后逐步减少,最后趋于稳定达到残余强度值。贯通节理岩体峰值强度后的本构关系并不适合于描述该特征。提出一个新的非线性归一化位移软化本构模型,通过归一化剪切应力R与归一化剪切位移D的指数函数关系体现峰值后剪切应力、位移的发展趋势,即模型采用无量纲表达式。归一化剪切应力R是峰值后强度降? p-? 与峰值强度、残余强度差值? p-? r的比值;归一化位移D是峰值后剪切位移、峰值位移差值? -? p与残余位移、峰值位移差值? r -? p的比值。对含不同起伏角的非贯通节理岩体在常法向荷载条件下的直剪试验数据进行拟合分析,模型预测结果与实测值具有相当高的吻合度,验证了模型的正确性     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

塔中地区构造应力场数值模拟研究

构造应力场研究对于油气运聚分析及储层特征描述具有重要的理论和实际意义。本文根据塔中隆起断裂特征及其演化史,结合钻井资料,选择中下奥陶统碳酸盐岩为主要研究对象,应用大型有限元分析软件ANSYS,在处理断层问题时运用非连续的接触算法,模拟了中奥陶世末、中泥盆世末、三叠纪末、新近纪末的构造应力场,研究了该区不同时期的古应力分布特征。研究表明,最大主压应力受东部车尔臣-星星峡走滑断裂影响显著,在隆起区沿断裂带呈条带状低值分布,在不同时期塔中Ⅰ号断裂带无一例外地位于最大主压应力的最低区,反映了塔中断裂带尤其是塔中Ⅰ号断裂带是油气运移的有利指向区。最小主压应力沿断裂带呈条带状高值分布,在断裂的上盘、断裂走向发生变化和断裂的倾末端是张应力高值区,为张裂缝发育的重要构造部位。最小主压应力和最大剪应力受车尔臣-星星峡走滑断裂的影响明显减弱。