一种基于MATLAB的碎屑岩粒度分析方法

利用MATLAB图像处理与统计计算功能对碎屑岩进行粒度分析,可以较好地解决传统粒度分析方法中存在的测量结果不精确、费时费力等问题。首先对碎屑岩镜下图像进行灰度化、二值化、图像增强等处理,准确测量出碎屑岩中各类粒度数据,然后将粒度数据转换为粒度参数运用到沉积环境分析中,更为精确地判别沉积环境。分析结果表明,该方法不但与传统粒度分析方法结果一致,而且测量精度高,方法简单,操作便捷,适用性好。     

利用一维正则化反演进行大地电磁测深数据拟二维反演解释

大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。     

Shearlet域稀疏约束地震数据重建

在地震数据处理流程中,通常对不规则的、稀疏的或者缺失的地震数据进行插值处理,通过插值方法来避免多次波的预测错误和成像假频等现象,使地震数据处理更加精准。Shearlet变换是一种多尺度变换,具有最佳的稀疏性、方向性以及局部化特性。将Shearlet变换与基于Landweber加速下降迭代方法结合起来对地震数据进行插值,在保证求解精度的同时提高了计算效率。信号和噪声在Shearlet域具有不同的分布特点,通过阈值法压制随机噪声,可提高算法的抗噪性。此外,采用jitter采样的方式,更好地压制了假频信息。理论和实际地震数据验证了该方法的有效性。     

空间信息网格服务工作流模型及其应用

现有的空间信息网格应用系统缺乏对计算资源的有效组织及管理方法,无法满足复杂空间业务流程化建模及计算任务协同化处理的需求.从数据与功能相互分离的GIS软件架构角度出发,采用活动网络图的方法提出了一种以服务功能建模为基础的新型空间信息服务工作流模型,该模型通过服务功能的组装拆卸及动态聚合能够实现空间信息服务资源的灵活构建及动态迁移,较好地解决了网格环境下空间计算资源的统一组织与按需分配、计算任务的并行化处理等问题.在此基础上,建立了空间信息网格工作流集群的总体框架并讨论了系统的自动负载均衡及流程资源的动态调度机制.研究成果成功地应用于地质调查信息网格环境下的矿产资源储量测算系统,实现了网格环境下复杂空间应用的工作流处理模式,并有效地提高了业务系统的处理效率.     

化学物相分析在物料成分全分析数据处理中的应用

讨论了化学物相分析在物料成分全分析数据加合平衡计算中的作用。根据对不同类型、不同共生组合且具有代表性的物料成分全分析数据平衡计算过程中所推荐的计算程序和计算结果,介绍了化学物相分析在物料成分全分析数据处理中的实际应用。同时阐述了具体应用对象、应用规则、应用模式以及应用过程中需要注意的问题等。     

一种基于GCV法改进的频域内位场正则化向下延拓的参数选取方法

这里给出了一种在频域内位场正则化向下延拓方法中,用于选取最优正则化参数的极小化泛函。通过不同的正则化迭代方法进行仿真数据和实测航磁数据的向下延拓实验,来验证此泛函用于选取参数的有效性,实验结果证明了此方法用于位场向下延拓是可行和有效的。     

露天矿矿坑涌水量预测计算实例

矿坑涌水量预测计算,是矿山水工环工作的重点。不同的矿山具有不同的边界条件,条件不同计算方法不同,按照《矿坑涌水量预测计算规程》(DZ/T0342—2020),有6种基本方法。以西藏某露天矿为例,通过野外收集的有关水文地质参数,进一步作水文地质条件概化、建立模型,选择与之相适应的水均衡法,进行露天矿矿坑涌水量预测计算,为详查阶段的地质报告提供了详实的基础数据,计算结果比较切合实际。为同类型的露天矿矿坑涌水量预测计算,提供一种案例参考。     

一种求地球化学异常下限的新方法——含量排列法

地球化学异常下限的确定,在矿产勘查和资源预测中占有十分重要的地位。传统计算地球化学异常下限的方法,要求数据服从正态分布或对数正态分布,而且存在计算依据不充分的问题。这里在研究背景和异常关系的基础上,提出了一种求地球化学异常下限的新方法,即含量排列法。该方法不要求数据服从正态分布或对数正态分布,直接用原始数据求异常下限,方法简便。用此方法对广东始兴地区化探数据进行处理,取得了较好的效果。     

视磁化率(SUSC)填图的理论与实践

视磁化率(SUSC)填图是对总磁场强度数据进行高分辨率增强滤波的一种方法.虽然也可以对地面磁测或一维数据进行视磁化率的计算,但用于二维航磁图的分析,更能发挥其巨大的潜力.计算视磁化率最简单直接的方法,是在波数域中进行分支线性滤波(component linear filter),需要作以下四种滤波处理:低通稳定滤波、向下延拓、化到磁极、消去几何因子.方法所依据的统计前提是:假定地     

测量数据中粗差的探测与处理

利用巴尔达传统粗差检验方法检验剔除出数据中的粗差,并通过稳健估计方法处理含有粗差的数据以及传统平差方法处理剔除过粗差的数据,与加入粗差前的数据处理结果作对比,比较不同粗差处理方法的优缺点。经比较得出：在保证一定迭代计算次数的条件下,稳健估计能够达到抵抗粗差的目的;巴尔达粗差探测在数据中只含有一个粗差的情况下,与稳健估计相比,在满足一定精度要求的前提下,更能保证计算效率;在数据含有多个粗差的情况下,可利用稳健估计方法进行处理。     

高层住宅短桩基础的沉降分析

对高层住宅（１９层）采用短桩基础处理,是一种有效的地期处理方法。本文运用四种不同的方法估算桩基础的沉降量,得到了短桩基础的沉降数据范围。通过四种计算结果与实测值的对比分析,提出了更为接近工程实际的计算模式。     

地磁异常数据处理技术在低纬度地区的应用研究——以柬埔寨王国Mesam金矿为例北大核心CSCD

柬埔寨王国位于低纬度地区,地质勘查程度较低,地面高精度磁测技术能够快速获取磁异常信息,进行找矿评价。在低纬度地区磁化方向主要以水平磁化为主,磁性体产生的磁异常以负磁异常为主。本文通过对柬埔寨王国Mesam金矿磁测数据进行倒相180°、化赤和低纬度化极处理对比,认为三种处理方法得到的结果都能够完整地反映磁场分布特点,在不能忽略剩磁影响的低纬度地区对实测数据进行低纬度化极后再进行解释更加贴近真实且细节更加丰富。通过对低纬度化极后的磁异常进行水平总梯度模和小波多尺度分解计算等位场异常处理,认为磁异常水平总梯度模计算能够较好地反映深部地质体的边界范围,小波多尺度分解计算能够提取特定深度地质体的空间信息,在该区选取的磁异常处理技术是有效的,为进一步扩大该金矿区远景储量提供了新的思路。     

基于低频约束的高分辨率Radon变换多次波压制方法研究

抛物线Radon变换是地震资料处理中常用的一种多次波压制方法,常规Radon变换由于采用最小二乘算法的原因,对不同的曲率值采用相同的加权值进行求解,分辨率较低。为提高Radon域数据的分辨率,应采用对不同曲率值应用不同的加权值的方法。一般的思路是利用前一次迭代的结果得到加权矩阵,该算法需要迭代进行,运算量较大。这里拟采用低频约束的方式提高Radon变换的分辨率,即利用前一个频率的运算结果对下一个频率的计算进行约束,该算法分辨率高,计算速度快。模拟数据和实际数据的测试表明,本文方法在多次波压制处理中,可以快速有效地完成多次波去除。     

自动识别多期断层擦痕的一种应力反演算法

由于地质历史上构造应力场的演变，多期断层擦痕数据的存在是应力反演所面临的普遍性问题。以往提出处理多期断层擦痕的应力反演算法都基于硬划分，忽视了数据自身的不确定性，并且一些只是传统的、处理一期断层擦痕的算法的简单延拓。在Fry(1999)的sigma空间里，同期断层擦痕向量具有统一的线性分布趋势，多期断层擦痕向量具有不同的线性分布趋势。在此基础上，本文提出利用模糊线性聚类法来识别多期断层擦痕向量的线性结构。这种算法不仅可以弥补以往算法的上述缺陷，还具有自动、直接、有效，且计算量也较小的优点。     

大孔径闪烁仪观测数据的处理方法研究

大孔径闪烁仪可以获取公里级尺度的显热/潜热通量,其观测数据在农林业、水文、气象等领域的研究中具有越来越重要的作用.以海河流域2008年密云站和馆陶站大孔径闪烁仪的观测数据为例,讨论了数据筛选与处理方法,以及大气不稳定状态下不同处理方法对显热通量造成的影响.结果表明:空气折射指数结构参数应结合电压信号的方差进行计算、采用逐日的日平均波文比系数进行湿度订正、结合空间权重函数计算有效高度、选取Andreas(1988)稳定度函数计算显热通量等是可信的.通过对不同插补方法的比较可知,采用非线性回归方法和以零值替代的方法对不稳定状态和稳定状态下30 min缺失的显热通量进行插补,采用动态线性回归方法对日显热通量的缺失进行插补是可行的.为了解决由于冬天净辐射和显热通量较小,造成能量平衡方程余项法计算的蒸散量产生较大误差的问题,可采用根据大孔径闪烁仪与涡动相关仪观测的日蒸散量建立的关系式进行估算.基于上述方法的研究,形成了一套较为完整的大孔径闪烁仪观测数据的处理流程,保证了在不同下垫面、不同的天气状况条件下都能获取质量可靠、数据连续的大孔径闪烁仪观测的显热/潜热通量.     

滑坡位移分解预测中的平滑先验分析方法

目前在滑坡位移预测研究中,先将滑坡位移数据分解为趋势项及周期项后再分别进行预测已成为普遍做法。平滑先验分析法(Smoothness Priors Approach,SPA)是一种计算过程简单、计算量极小,且能快速分离原始数据趋势项和周期项的数据处理方法。在介绍SPA基本原理基础上,以三峡库区白家包滑坡典型监测点位移数据为例,对通过调节SPA正则化参数而获得的不同趋势项及周期项进行特征分析;进而结合对滑坡变形演化机制过程的先验分析,根据位移分解特征确定合理的参数取值;最后针对不同参数SPA位移分解数据,采用支持向量机进行位移预测对比分析。结果表明,SPA是一种适用于滑坡位移预测的位移分解方法,通过调节正则化参数并结合滑坡变形机制先验分析,能够获得较为合理的位移分解结果,进而提高滑坡位移预测精度。     

位场垂向高阶导数的Tikhonov正则化迭代法

在位场数据处理中,垂向导数具有重要的物理意义。其在一定程度上可以划分不同深度和大小异常源产生的叠加异常,且导数的阶次越高,这种分辨能力就越强,但通常认为高阶导数的换算是不稳定的。本文在Tikhonov正则化求位场垂向高阶导数的基础上,结合迭代法进行逐次逼近,提出了位场高阶导数的Tikhonov正则化迭代法,并且得到Tikhonov正则化迭代法的递推公式。通过对该方法的滤波特性分析可以看出,该方法计算的位场垂向高阶导数具有一定的稳定性及保幅性。模型试验和实际数据的处理表明,该方法计算结果较常规FFT求导法有更高的稳定性和实用价值。     

基于地震同相轴预测的剩余静校正方法研究

为解决野外地表参数测定不准确等因素引起的地震数据剩余静校正问题,通过平面波分解滤波器计算时空变地震同相轴局部倾角,利用局部倾角与同相轴的时空预测量之间的关系计算剩余静校正量,在剩余静校正量计算所需局部平面波方程求解过程中,加入时间方向倾角为常数的整形正则化约束条件,以保证剩余静校正方法的准确性。结果表明,基于地震同相轴预测的剩余静校正方法可以有效地处理共中心点道集数据的剩余静校正问题,同时本方法具有一定的抗噪能力,为实际资料处理提供了较好的适用性。     

应用网函数插值方法进行重力区域改正计算

应用网函数插值方法进行重力区域改正计算是内蒙古物探队和内蒙古大学数学系合作进行的。是应用电子计算机快速处理重力测量数据一种新方法。有关网函数插值方法的实质、原理、推导及计算重力区域场改正的理论讨论、模型试验都已在不同文章中给予介绍。本文着重介绍网函数插值法用来进行区域改正计算的具体方法和地质效果。     

超大数据量的快速排序法

在超大数据量情况下,传统的排序方法及分级优化化法由于其自身的算法限制,难以更有效地提高排序效率。为此,作者在本文中提出了一种快速排序法,该方法的基本思路是用特殊设定的不同分界值把整个数据区域成尽可能多的小单元,每个数据通过与不同分界值的对比迅速归位于相应的单元内,然后对各单元内的数据进行排序,按序边接各单元,最终达到排序的目的,与分组优化法的计算工作量大小对比及试算结果,充分表明了该方法的算法可行,计算量少,排序速度更快等优点,特别适用于超大数据量情况的数据排序。