多自由度块行压缩存储技术及大型稀疏方程组的求解

在地球物理电磁勘探领域有限元数值模拟中,最后都会得到一个大型稀疏的复系数线性方程组,受计算机内存空间的限制,必须根据有限元刚度矩阵的稀疏性对其进行压缩存储。由于电磁场有限元计算的自由度大都在三个以上,因而提出了适合多自由度的块按行压缩稀疏存储方案,并通过存储格式的转换,把块按行压缩方式转换成流行的,大型稀疏矩阵的行压缩存储格式,以便于求解。用求解大型稀疏方程组的Krylov子空间方法中的稳定双共轭梯度（Bicgstab）方法,收敛速度快,精度高,而且稳定性好,结合ilu预处理技术,可以大大提高求解大型稀疏方程组的效率。     

浅析《详图集》编制过程中的检查与验收

介绍了《详图集》在编制过程中检查与验收工作的必要性以及原则,并对《详图集》检查与验收工作的主要内容和实施方法进行了详细地阐述。     

音频大地电磁法强干扰压制试验研究

压制非相关噪声已经有多种比较成熟的技术方法,如多次叠加、远参考、Robust阻抗估计等,但这些方法对压制矿集区存在的强电磁噪声(通常属于相关噪声)基本没有效果.本文在前期提出的针对矿集区大地电磁强噪声压制的数学形态滤波基础上,开展了音频大地电磁法强干扰压制的试验研究.在四川西昌某AMT测点附近,布置了大功率电磁发送站,通过接地导线向地下供20A的方波电流,同时在供电期和间歇期用凤凰公司MTU-5A仪器观测4个水平分量的电磁场.在室内采用数学形态滤波和阈值法相结合的手段在时间域对含有人工源强干扰的电磁场数据进行处理,以视电阻率对比、频谱和极化方向等参数对处理效果进行评价.结果表明:结构元素恰当的数学形态滤波可以对人工源强干扰进行识别、去除,阈值法则可以消除形态滤波后的脉冲干扰,二者结合可以有效地压制AMT中的强干扰噪声,提高数据可靠性.本文结果对AMT实际资料处理具有重要的实用和参考价值.     

复杂地下异常体的可控源电磁法积分方程正演

可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力.

     

全信息矢量相干度在MT数据质量评价中的应用效果分析

采用最小二乘法计算张量阻抗要素时,可以利用各种计算结果之间的差异来衡量噪声影响的程度,而全信息矢量相干度正是基于这种思想提出来的,用以评价数据质量.以庐枞矿集区MT数据为例,通过对部份具有典型特征的数据进行质量评价,并与常见的数据质量评价方法进行对比,验证了全信息矢量相干度在MT数据质量评价中的作用,并分析了各种评价方法的应用效果,为MT数据的进一步处理提供了有力的保证.     

浅谈常规地图的审校方法

通过对常规地图产品审校方法的研究探讨,总结了一些有实用价值的经验,以满足地图生产全过程的技术指导与检查,确保地图产品的质量。     

黄河三角洲埕北海域底质沉积物粒度特征及泥沙输运分析

通过黄河三角洲埕北海域表层沉积物的粒度资料、实测潮流数据以及50 a一遇的波浪要素,分析了埕北海域表层沉积物的分布类型以及粒度参数特征,并计算了表层沉积物的临界起动应力、日平均单宽输沙通量以及波流共同作用下的单宽推移质输沙率。结果表明:研究区表层沉积物主要是粉砂质砂、砂质粉砂、粉砂和泥四种,分选中等偏差,近对称到正偏,中等尖锐到尖锐峰态;对大潮期间埕北海域潮流底应力和表层沉积物的临界起动应力的计算说明潮流在大部分区域具有起动和输运泥沙的作用;秋季潮流作用下的日平均单宽输沙通量均值为372.32 kg/（m·d）,方向为涨潮流向。利用贝克尔（Bijker）和经典的Van Rijn公式计算出的50 a一遇的波流共同作用下单宽推移质日平均净输沙率结果近似相等。     

基于字典学习的音频大地电磁数据处理

音频大地电磁（Audio Magnetotelluric,AMT）信号常常受到持续性人文噪声影响,这类噪声使用远参考法和Robust阻抗估计等常规方法往往效果不佳.为此,本文从噪声的规律与特征出发,提出一种新的AMT数据处理方法.首先通过字典学习方法从观测数据中自主学习到人文噪声的特征结构,构建冗余字典,然后利用学习到的冗余字典,分离出AMT数据中的人文噪声.为验证方法的有效性,首先进行了合成数据的仿真试验,然后在四川凉山进行了针对性的野外试验研究,最后将本文方法应用于庐枞矿集区实测数据的处理.结果表明,本文方法能够快速、准确地分离出AMT信号中的人文干扰,保留有用信号,显著改善AMT数据质量.

     

音频大地电磁法“死频带”畸变数据的Rhoplus校正

为研究AMT数据的"死频带"畸变特征及其校正方法,在长江中下游地区,利用凤凰MTU-5A仪器,系统观测了音频电磁场的"死频带"数据.对不同季节、不同时段的AMT"死频带"数据进行了分析.提出AMT"死频带"频率域畸变数据的Rhoplus校正方法,给出了该方法的适用条件、关键技术与评价方案,并提供了大量实测数据论证其应用效果.观测实验表明:在"死频带"内天然音频电磁场信号强度极低,正交电磁场分量相干度极低;造成了阻抗视电阻率、相位及相位张量等数据的畸变;畸变范围可达10 k~100 Hz,以5~1 kHz频段最为明显;并且数据畸变程度与观测时段密切相关,在长江中下游地区为秋冬季畸变强夏季弱,日间畸变强夜间弱.处理实例表明:Rhoplus校正方法可以有效地处理AMT"死频带"内畸变的视电阻率、相位数据,相干度阈值与人机交互相结合的处理策略快速客观,所得结果曲线光滑连续,数据与实测未畸变数据基本一致.     

基于频点CSAMT-维最小构造反演

传统CSAMT线性反演方法,往往存在对初始模型依赖性强和对某些低阻薄层反映不灵敏的弱点,这里在CSAMT最小二乘反演的基础上,在目标函数加入光滑约束条件,构成最小构造反演。另外,根据频点个数和地表电阻率来构建初始模型,使反演的初始模型更加灵活、简单、合理、更易应用于实际数据处理中。