InSAR与概率积分法联合的矿区地表沉降精细化监测方法

结合InSAR与概率积分法的优势,提出一种InSAR和概率积分法联合进行矿区地表沉降的精细化监测方法.该方法首先计算InSAR时序累积沉降盆地,进而建立判别大梯度形变的约束条件,区分沉降边缘与沉降中心.对于形变较小的沉降边缘,保留InSAR结果,而对于大梯度沉降中心,则结合InSAR与概率积分模型建立矿区工作面的沉降盆地,并通过空间插值,获取地理坐标系下连续的地表沉降信息,最终得到完整的矿区地表沉降结果.论文以山东某矿区为研究区域,采用2016年10月16日 2018年3月4日期间的21景SAR影像和工作面水准实测数据对该方法的可行性和精度进行了验证.结果 表明,该方法能够在减少水准监测工作量的前提下,获得与实际情况相吻合的沉降结果,其监测能力明显优于常规InSAR和概率积分法,可有效弥补两种技术单独在矿区地表沉降监测中的不足,获取更为准确、可靠的矿区地表沉降信息.     

基于双精度与四精度的重力场解算精度分析

基于动力学方法比较分析了双精度与四精度模式下重力场模型的解算精度,主要包括缔合勒让德函数计算、数值积分器及重力场反演结果。结果显示,在勒让德函数计算方面,部分角度在双精度模式下计算至1 900阶以后会出现溢出问题,而在四精度模式下任何角度都满足精度要求,并且计算结果比双精度模式高8个量级。数值积分器Adams预测校正法积分1 d的位置和速度误差,在四精度模式下比在双精度模式下高4个量级。在精密轨道反演重力场计算方面,动力学方法在双精度及四精度模式下反演结果一致,统计其计算至60阶的累计大地水准面误差为1.29×10~(-5 )m,这是因为动力学方法的线性误差相对计算误差而言是主要误差;非线性动力学方法在四精度模式下比在双精度模式下高7个量级,其大地水准面误差分别为8.92×10~(-15) m和8.16×10~(-8) m。     

新疆艾比湖流域枯、丰水期三维荧光光谱特性及其与水质的关系

以艾比湖流域主要入湖河流为研究对象,在5月（丰水期）和8月（枯水期）分别沿博尔塔拉河（博河）和精河进行采样,采用平行因子模型（PARAFAC）和三维荧光区域积分法对水体三维荧光特性进行研究并对其与水质的关系在枯、丰水期下的变化进行探讨.结果表明①河流DOM在枯水期与丰水期都含有C1（240、425 nm） UVC类腐殖质,C2（225、290 nm）紫外区内络氨酸类有机物,C3（230/280、330 nm）蛋白类有机物,C4（265、260 nm）腐殖质类共4种组分.通过对水体三维荧光进行区域积分可以看出DOM荧光成分的占比在不同时期的变化.博河在枯水期时EEM光谱中的区域Ⅲ富里酸含量低于丰水期,枯水期时区域Ⅱ芳香类蛋白质、区域Ⅳ可溶性微生物代谢物以及区域Ⅴ类腐殖质酸高于丰水期;对于精河来说,区域Ⅱ芳香类蛋白质和区域Ⅳ可溶性微生物代谢物在枯水期的含量高于丰水期,区域Ⅲ富里酸和区域Ⅴ类腐殖质酸的含量枯水期低于丰水期,这表明水体腐质化程度较高.②本研究选取了一些常规的荧光指数来描述枯、丰水期水体的荧光指数特性.经研究发现,精河的荧光指数、自生源指数和腐殖化指数在不同时期的变化幅度较小,而博河的变化幅度较大.③将荧光指数与水质参数进行相关性分析并建模,结果表明枯水期自生源指数（BIX）与化学需氧量呈显著正相关,相关系数R=0.688;丰水期时BIX与铵态氮浓度呈显著负相关,相关系数R=-0.493.通过对比分析艾比湖主要入湖河流的三维荧光光谱特性与水质在枯、丰水期时的关系进一步表明水体中DOM的特性以及在枯、丰水期下的差异,为艾比湖流域的治理改善提供一定的理论支持和参考依据.     

机器学习辅助下的概率积分法参数预计模型寻优

收集整理了多组地表移动观测站资料作为训练样本和检验样本,以工作面地质采矿条件为输入集,概率积分法预计参数为输出集,利用机器学习方法对概率积分法预测参数进行了预测。选取支持向量机、BP神经网络和偏最小二乘法3种机器学习方法对训练样本进行训练,利用训练所得模型预测检验样本中的概率积分法预测参数,并将预测结果与观测站实测值进行对比。结果表明,利用支持向量机预测下沉系数、主要影响角正切值及水平移动系数的精度最高,其平均相对误差分别达到7.46%、4.00%、13.17%;拐点偏距及开采影响传播角利用偏最小二乘法预计精度最高,平均相对误差分别为10.83%、0.88%;总体而言支持向量机的预测精度最为稳定。     

基于改进的散射积分算法的初至波走时层析

初至波走时层析是获取近地表速度结构的一种常用方法.随着采集技术的不断发展,可使用的数据量迅速增多,传统的基于射线追踪和解方程组的地震走时层析成像方法面临着内存占用大、方程求解不稳定等问题.为了解决这些问题,本文基于前人在波形反演研究中提出的一种改进的散射积分算法,提出了一种预条件最速下降法初至波走时层析.该方法无需存储核函数矩阵与Hessian矩阵即可方便地实现目标函数梯度的计算与预条件,且该方法计算效率高、求解稳定、易于并行.数值实验结果表明,该方法可以获得与传统方法精度相当的反演结果,但所占用的内存大幅减小.     

顾及多方向观测值权比反演地球重力场的动力积分法

考虑到不同坐标系下各个方向观测值对反演地球重力场的频谱贡献不同,建立了顾及多方向观测值权比的动力积分法,并利用该方法反演了高精度的GOCE HL-SST卫星重力场模型.首先,分析了不同坐标系下各个方向观测值与地球重力场信息的响应关系,其中惯性系(IRF)下X、Z方向的观测值分别对扇谐系数、带谐系数最为敏感,Z方向的解算精度在全频段均略高于X、Y方向;地固系(EFRF)下各个方向的独立解算精度均与能量守恒法的解算精度相当;局部指北坐标系(LNOF)下X、Z和Y三个方向的解算精度依次递减,且Y方向在47阶附近有明显"驼峰"现象.其次,比较了不同坐标系下顾及三个方向观测值权比的加权解算模型,其中加权联合解算模型精度在20至70阶次均明显优于等权解算模型,在带谐项和共振阶次精度提升明显,且LNOF下的加权联合解算精度要优于IRF和EFRF.最后,比较了GOCE和CHAMP卫星的模型解算精度,采用本文计算方法,仅利用2个月GOCE轨道观测值解算的模型精度优于包含更长观测时段信息的AIUB-CHAMP01S和EIGEN-CHAMP03S模型,且略优于ASU-GOCE-2months模型.     

模拟三维裂纹问题的多尺度扩展有限元法

扩展有限元法模拟裂纹时独立于网格,因此该方法是目前求解裂纹问题最有效的数值方法。为了在计算代价不大的情况,实现大型结构分析中考虑小裂纹或提高裂纹附近精度,在裂纹附近一般采用小尺度单元,其他区域采用大尺度单元。提出了分析三维裂纹问题的多尺度扩展有限元法,在需要的地方采用小尺度单元。基于点插值构造了六面体任意节点单元。所有尺度单元都采用8节点六面体单元,这样六面体任意节点单元可方便有效地连接不同尺度单元。采用互作用积分法计算三维应力强度因子。边裂纹和中心圆裂纹算例分析结果表明,该方法是正确和有效的。     

Kirchhoff积分法OBS数据地震波场延拓

针对海底地震仪(OBS)数据基准面校正的特殊性,采用Kirchhoff积分法对中国南海实际OBS资料进行了基准面校正处理,将炮点延拓到OBS所在处的海底面上,取得了预期的效果。分析表明该延拓法对地震波场的分辨率及信噪比影响不大,且与静校正方法相比较,该方法能够得到地下地层较真实的速度场,这对于后续的动校正、叠加等常规处理十分重要。     

土石坝动力分析的改进二次加速度积分法

土石坝动力分析常用的是Wilson-? 积分方法,但该方法用于分析高频问题和大步长问题时精度较低。针对该缺点提出一种改进的二次加速度积分法,该算法与一步二次加速度积分法相比,计算精度不变、稳定性要求容易满足,计算用时大大减少。对算法的基本原理、稳定性及精度进行了分析,并给出一个土石坝动力分析的算例,验证了该算法的可行性。与常用的Wilson-? 法、Newmark法和中心差分法进行了数值比较。计算结果表明：该算法在分析高频问题和大步长问题时,计算精度较高。     

一元线性类储量计算积分法应用简介

本文在《储量计算积分法概要》基础上,对若干新术语和常用的一元线性类基本用法用了简要介绍。并结合示例重点说明了相邻工程和断面间各地质变量呈线性变化时的储量计算过程及其应用效果（计算精度一般可提高５％￣３０％或更高）。本法简便易行,且有利于提高储量报告质量和矿产的合理开发及充分利用对一般矿产储量计算及其他行业中类似问题均有较高的实和价值。