C波段雷达在新疆区域数值预报系统中的应用初探

基于新疆区域数值预报系统（Desert Oasis Gobi Rapid Analysis Forecast System,下称DOGRAFS）,开展了同化C波段雷达资料对2010年10月6日发生在新疆的一次强降水过程预报结果影响的试验分析。其中设计包括不同化任何资料、同化常规资料、同化雷达反射率因子、同化雷达径向风、同时同化反射率因子和径向风五组试验,重点分析了雷达资料同化对此次天气过程降水、温度以及风速模拟效果的影响。结果表明：（1）同化雷达径向风和同时同化径向风和反射率因子相对于其他三组试验,对降水预报的TS和ETs评分更高;（2）相对于其他三组试验,同化雷达径向风和同时同化径向风和反射率因子对模式垂直方向上的温度、风速预报偏差具有一定的改善效果;（3）对于地面2 m温度和10 m风速而言,同化常规观测资料比其他四组试验预报的平均偏差和均方根误差更小;其它四组试验误差相当,差别不明显,表明同化雷达对近地面层温度和风的影响不明显。本研究旨在探索C波段雷达观测资料在新疆区域数值预报系统中的适用性,为今后雷达观测资料在业务系统中的应用提供参考。     

BIM技术在深基坑工程中的应用探讨

BIM技术是一种新型的设计方法,它实现了设计从二维空间向三维空间变革,从静态的信息向动态的信息发展。在大量查阅国内外有关基坑支护设计研究文献的基础上,综述了基坑支护设计方法,分析了不同设计软件优缺点和设计方法需解决的关键性理论问题,提出了基于BIM技术在深基坑工程中的应用设想,分析了BIM技术应用在深基坑工程中需解决的关键性技术问题,为相关研究提供参考借鉴。     

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

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

九瑞矿集区三维电性结构研究及找矿意义

九瑞矿集区作为全国重点整装勘查区之一,深部找矿工作一直是研究重点和热点.本文在综合分析研究区区域地质资料、对典型地层和岩石标本进行电阻率参数测试和统计分析的基础上,利用研究区面积性AMT数据进行三维电性结构研究.为获得AMT最佳反演参数,选择一条典型AMT剖面开展二维非线性共轭梯度(NLCG)反演,进而总结出利用TM极化模式数据、拉格朗日乘子等于3的二维反演方案.综合利用地质、钻探和电性特征对全区23条AMT剖面进行二维约束反演,并将所有测线的电阻率反演结果三维网格化,得到了研究区的三维电性结构.最后,我们通过对三维电性结构的综合分析,确定研究区地下地层、岩体展布、断裂构造分布特征以及基底起伏情况.与已知矿床对比,我们发现浅部北西断裂与北东向区域性深部断裂交汇处、基底隆起凹槽边缘以及岩体边缘为成矿有利区.     

九瑞矿集区重磁三维约束反演及深部找矿意义

研究九瑞矿集区成矿地层和岩体的三维分布特征可为区域成矿背景和成矿规律研究提供新的信息,实现研究区深部及外围找矿突破.本文首先综合分析研究区区域地质及地层和岩石密度和磁化率特征,然后将区域地质和12条地质剖面信息利用自行开发的GIF2UBCmodel程序构建剩余密度和磁化率参考模型和上下边界约束模型,实现了带约束重磁数据三维反演.反演结果揭示九瑞矿集区剩余密度和磁化率三维分布特征,解译了矿集区主要地层和岩体分布特征.该区存在邓家山—东雷湾以及武山—丁家山两个大型岩体;盖层褶皱与结晶基底隆起方向一致,是深部基地隆起的继承.通过与矿集区已知矿床对比发现,高磁性岩体边缘和基底隆起(高密度体)两翼为成矿有利区.     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

基于谱元法的广域电磁法三维正演模拟

为了提高广域电磁法三维正演精度和效率, 本文提出采用基于Gauss-Lobatto-Legendre(GLL)多项式的谱元法进行广域电磁法三维正演模拟.首先从麦克斯韦方程组出发, 推导了二次场满足的电场双旋度控制方程, 利用伽辽金加权残差法将微分形式的边值问题转换为积分弱形式, 再通过单元剖分和高阶正交基函数插值对全域问题进行离散, 最后通过求解大型线性方程组得到全局的数值解.层状介质模型验证了算法的正确性和精度, 三维地电模型分析了算法稳定性以及广域电磁法响应特征.研究表明谱元法是进行广域电磁法三维正演模拟的有效方法, 具有高精度、低网格依赖性等优势.

     

基于自适应有限元正演的大地电磁法三维反演算法研究

目前大地电磁法三维反演大多基于有限差分正演.反演使用规则六面体网格无法有效模拟复杂地形, 同时正反演网格同套存在严重影响反演可靠性的问题.针对上述两个问题, 本文利用自适应有限元算法, 开发了基于非规则六面体的高精度大地电磁三维正演方法; 反演中, 使用独立的反演网格和正演网格来提高正演计算精度和反演可靠性.这一思路既保证了正演响应和灵敏度矩阵计算的精确性, 又降低了因反演参数过多造成的不唯一性.最后, 通过算例验证了正演算法的精确性和反演算法对地形处理的有效性.

     

利用极化率参数监测煤矿滞后突水的可行性

煤矿工作面底板下方隐伏断层、陷落柱等小尺度地质构造不仅难以识别,而且在采动影响下,还将进一步发育,形成潜在的隐伏导水通道,容易发生滞后突水.本文提出通过实时监测极化率参数的变化实现对由小规模地质构造引起的煤矿滞后突水预警.采用三维激发极化法装置形式,利用理论分析、数值模拟和物理模型实验相结合的方法,发现采空区底板下陷落柱构造突水过程中极化率参数变化规律.随着水位上升,异常区极化率逐渐增大,异常位置更加明显,可根据极化率参数的极大值点的位置,通过交汇作图方法确定构造异常体中心的平面位置,可根据极化率参数值大小变化监测其水位变化.结果表明通过对极化率参数的监测,对隐伏构造突水过程进行及时预警具有可行性.

     

渤海MN油田南区相控地质建模

在精细地层对比、沉积相分析和隔夹层分布等研究基础上,笔者采用相控方法进行渤海MN油田南区储层三维地质建模。在沉积相模拟中,以沉积模式为蓝本,将沉积亚相展布、井震结合预测的泥岩夹层分布等确定性认识和随机模拟技术相结合,同时依据各沉积亚相内储层发育的特征和差异,选择不同的“砂包泥”和“泥包砂”模拟思路,实现逐级相控、定量表征储层的空间分布。研究结果表明,以精细油田地质特征认识为基础,采用多条件地质约束的相控建模方法,可以更好地反映沉积认识,提高储层建模的精确度。