层状矿体等高线离散化后的四面体建模

基于不规则四面体的三维实体建模方法具有很多优点,但在层状矿体的建模中,该方法却不能准确地描述其层状界面。根据层状矿体建模数据的特点,先将层状矿体的等高线进行等距离散化,再建立不规则四面体三维实体模型,并为该方法设计由点、面和四面体组成的逻辑模型结构。该方法用VC++.net编程得以实现,并使用内蒙某煤矿的实际数据,建立了相应的矿床三维实体模型。实践表明,通过离散化等高线后建立的层状矿体四面体模型,能够准确反映出层状矿体的层状界面,实现对层状矿体的三维实体建模。      

二维Delaunay三角网的任意点删除算法研究

针对目前基于影响域多边形剖分的点删除算法缺陷,提出一种二维Delaunay三角网点删除算法。首先利用具有拓扑关系的三角网搜索影响多边形,并以三角形矢量面积为工具三角剖分影响域多边形,最后通过镶嵌优化后的剖分三角网完成点的删除,且满足Delaunay法则。通过测试证明了算法的可靠性和高效性。     

不规则界面介质中首波的传播路径

通常首波被认为是以下层高速度沿着平坦界面传播的折射波。然而,沿着不规则界面首波的传播路径至今仍不明确。在这里,我们将利用重叠单元(SEMO)方法以模拟的方式来研究这个问题。之所以选择重叠单元方法,是因为这种方法可以准确地计算出沿不规则界面传播的界面波。然后我们利用时间窗将首播从界面波中分离出来,并分析首播沿不规则界面传播进程中能量的变化和走时的变化。分析结果表明,与平坦的界面上首播的传播情况不同,沿不规则界面传播的首播由两种机制产生:折射机制和透射机制,也就是说,首播可能是沿界面传播的折射波或高速介质中向低速介质的透射波。关于首播的这种认识将有助于我们建设更加精确的首播反演方法,例如建立更为精确的首波旅行时层析成像,并由此获得更准确的地下结构模型,这对于研究认识一些特殊地区(如青藏高原)的形成机制是非常重要的。     

基于节理网络有限元的柱状节理岩体等效变形模量结构效应研究

因其特殊的地质成因,柱状节理成为一种具有强烈非连续性及各向异性的特殊岩体结构。以白鹤滩水电工程设计、建设中所遇到的柱状节理玄武岩体为背景,针对柱状节理岩体的等效变形模量问题,在探讨了其结构效应表征参数的基础上,利用节理网络有限元为工具,研究各结构效应表征参数对柱状节理岩体等效变形模量的影响。研究结果表明：在柱轴线横向平面上,对于柱体不规则程度,当柱体完全不规则相对于完全规则情况,等效模量提高大约10%左右;在其他条件不变情况下,当柱体平均边长从0.1 m增加至0.5 m时,等效模量从5.36 GPa增加至23.4 GPa;对不同加载方向的研究结果表明,柱状节理在柱轴横向平面上可以视作各向同性;当节理刚度线性递增时,相应的岩体等效模量也基本符合线性递增的规律。在平行柱体轴线平面,节理组2间距越大,柱状节理岩体的等效变形模量越大;随着错距比的增加,柱状节理岩体的等效变形模量呈现先增加后减少的趋势,当错距比为50%时,节理岩体的等效变形模量取最大值;对于节理刚度的影响,与柱轴横向平面类似的,随柱体轴线平面上节理刚度的线性增加,相应的岩体模量也基本符合线性递增的规律。该研究成果与已有的现场研究成果及三维数值模型成果相比,吻合良好。说明了文中提供的研究结论具有较好的可信度及参考价值。     

用GPS观测资料解算地球自转参数的方法及影响因素分析

研究利用GPS观测资料解算ERP参数的方法。首先,利用全球近100个均匀分布、站址稳定的IGS测站的GPS观测数据解算地球自转参数,并将其与IGS相应产品进行比较,结果显示该解算方案可靠。其次,从实验数据和数理表示方面分析ERP参数解算过程中的轨道约束方案、先验ERP信息、站坐标数量对解算精度的影响,给出一些提高ERP参数解算精度的建议：对ECOM光压模型9参数中径、法向（D、Y）4个周期项采取相对宽松约束、随机脉冲参数采取紧约束,解算效果更佳;降低ERP先验信息对解算精度的影响,重点在于提高UT1-UTC的预报值精度;利用100个跟踪站数据解算,可达到解算效率和解算质量均最佳的效果。     

基于弹性薄板的不规则采动地表沉陷预计模型

文章借助弹性薄板理论,建立地表沉陷变形的预计模型,并给出任意点地表变形的预计表达式。针对地下开挖空间形状的不规则性,推导出了水平空间、倾斜面空间以及不规则开挖空间：椭圆抛物面空间、二次柱面空间以及二次曲面空间地表沉陷的预计公式。本文的研究结果对该类开采沉陷问题的研究具有一定指导意义。     

不同尺度低纬电离层不规则体漂移特性的观测研究

利用我国三亚站GPS短基线接收机阵以及VHF电离层相干散射雷达多波束扫描模式进行同时观测,对2011年10月22日至29日夜间两种不同尺度电离层不规则体的纬圈漂移速度进行比较研究.结果表明,两种不同手段观测的不同尺度不规则体漂移速度,在大小、变化趋势以及时间范围上基本一致.不规则体产生在20:00 LT(LT为北京时)左右,以100～250 m/s的速度东向漂移,并且速度幅值起伏较大;在不规则体发展中后期(21:00LT以后),东向漂移速度减低到50～150m/s,并延续至午夜.两种手段观测的漂移速度差异,一个重要原因可能来自于两者观测范围的不同.此外,GPS短基线接收机阵观测的400米尺度电离层不规则体通常比VHF雷达探测到的3米尺度电离层不规则体存在的时间要长.     

Condensed hyperelements method of non-vertical consistent boundaries for wave propagation analysis in irregular media

The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.     

原始机载LiDAR点云中建筑物激光点的自动提取

提出了一种基于原始离散LiDAR点云的多层结构建筑物激光点自动提取方法。首先对原始离散Li-DAR点云进行不规则三角网构建,然后利用分阶段区域生长的原理逐阶段遍历TIN中的各三角面。依据三角面的法向量方位角判定建筑物的各层结构,依次提取出多层建筑物的墙面及各级屋顶面,最后将提取出的各层屋顶面数据进行分类存储,确保提取出的各栋建筑物激光点之间相互独立且包含各层屋顶结构等细部特征信息,便于后续各建筑物的三维重建。     

局部倾角约束最小二乘偏移方法研究

随着石油勘探难度的进一步加大,地震数据往往存在采样不规则、地震道缺失等现象,如果不对其进行处理,会对后续的地震成像产生影响,引入成像噪音.针对这一问题,一般是通过地震道插值或数据规则化对叠前数据进行处理,然后采用常规的偏移方法进行成像,本文则是将地震成像看作最小二乘反演问题,在共成像点道集引入平滑算子,在共偏移距/角度道集引入平面波构造算子(PWC)进行约束,通过预条件共轭梯度法使得反偏移后数据与输入数据之间的误差达到最小,最终得到信噪比更高、振幅属性更为可靠的成像结果.理论模型和实际资料处理表明,本文方法不仅可以有效压制数据不规则对成像产生的噪音,而且具有更高的成像精度.