基于多网格的频率域全波形反演(英文)

频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有效改善问题极小值缺陷,加快反演的收敛速度。以局部非均匀的三孔模型和Marmousi模型的数值模拟结果验证了该算法的有效性。     

大地电磁全信息资料三维共轭梯度反演研究(英文)

在对张量阻抗数据、倾子数据和共轭梯度算法深入分析的基础上,我们实现了大地电磁全信息资料三维共轭梯度反演算法。基于全信息资料的三维共轭梯度反演研究,探讨了同时利用五个电磁场分量整理得到的大地电磁资料进行三维反演定量解释的方法以及全信息数据在三维反演中的作用。理论模型合成数据的反演结果表明,在三维反演中使用张量阻抗和倾子数据结合的全信息数据的反演结果优于只使用张量阻抗数据(或只使用倾子数据)的反演结果,提高了反演结果的分辨率和可信度。合成数据的反演算例也验证了所实现的大地电磁全信息资料三维共轭梯度反演算法的正确性和稳定性。     

基于G.729的固定码书搜索的改进性研究

在详细阐述G.729协议的CS-ACELP算法工作原理的基础上,分析了该算法各功能模块的算法复杂度,针对其总体复杂度较高,对硬件实现条件要求较高的问题,对复杂度相对较高的固定码书搜索部分采用减少脉冲搜索位置和补充搜索相结合的方式进行了改进,加快了码书的搜索速度,降低了算法复杂度.通过对改进后的算法进行实验仿真,验证了其可行性和有效性.     

频域航空电磁法地形影响及校正研究

为研究起伏地形对航空电磁响应的影响,利用二维频域交错网格有限差分方法( FDFD) 计算模型中电磁场分布情况,并且用不完全乔里斯基预处理的双共轭梯度法求解大型稀疏复线性代数方程组。设计山脊和山谷两种起伏地形,同时改变坡度大小,比较不同坡度对电磁响应的影响程度。在此基础之上,提出线性校正方法,该方法具有简单快速的优点,弥补了耗时巨大的带地形的反演方法的不足。任意地形下的地形校正后的结果表明,该方法可以准确指示地下异常的存在。     

三维频率域航空电磁反演研究

航空电磁数据的三维解释由于数据量大需要有高效的反演算法作为支撑.本文利用两种目前主流的数值优化技术（非线性共轭梯度和有限内存的BFGS法）实现了三维频率域航空电磁反演,并进一步比较了两种方法的有效性和运算效率.在反演过程中,为了更好地反演异常体的空间位置,模型方差矩阵中的光滑系数在反演起始阶段取值较大;当数据拟合差下降趋于平缓时,再利用较小的光滑因子约束反演过程来实现聚焦和获得精确的反演结果.理论数据反演表明这两种优化策略具有相似的内存需求,但是有限内存的BFGS技术比非线性共轭梯度法在计算时间和模型反演分辨率上具有一定的优越性,因此有限内存BFGS法更适合于求解大规模三维反演问题. 模型试验进一步表明目前主流的迭代法求解技术不适合大规模航空电磁数据反演,未来移动平台多源电磁数据快速正反演可通过引入矩阵分解技术来实现.     

大地电磁资料精细处理和二维反演解释技术研究(一)——阻抗张量分解与构造维性分析

本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术.     

缅甸Ms7.2地震前电离层TEC的异常变化

基于IGS提供的电离层TEC资料和中国地壳运动观测网络提供的GPS资料解算的单站VTEC数据,采用滑动平均处理方法,对2011年3月24日缅甸7.2级地震震前的电离层TEC资料进行了处理和分析,结果表明：在震前6～8天,出现电离层TEC异常减小现象,震前3天出现电离层TEC异常增加现象,排除太阳活动和地磁扰动的影响后,3月16日、17日和18日出现的电离层TEC异常减小,特别是在18日出现的异常减小及呈现出的共轭结构,可能与缅甸7.2级地震有关。     

基于稀疏约束的自适应正则化迭代重建算法

本文针对ASD-POCS算法中约束项权重对不同应用的多变性引起的算法鲁棒性差等问题,提出了一种基于稀疏约束的自适应正则化迭代重建算法,该算法采用一种Lagendijk型的正则化策略构造最优化问题,分别采用局部方差、图像能量估计自适应地求取加权对角矩阵和全局正则化参数。最优化问题的求解过程中,采用SART算法和共轭梯度法求解保真项和约束项最优化问题。实验结果表明,AR-SART-CG算法能更好地权衡恢复图像边缘和平滑噪声的关系,更好地调节保真项和约束项的权重,得到更高质量的重建图像。     

可控源音频大地电磁三维共轭梯度反演研究

可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.     

离轴二次曲面面型参数测量归算方法的研究

为了在制造过程中保证离轴非球面的光学特性,需要对其主要的面型参数进行严格控制,除了常规形式的顶点曲率半径和非球面系数外,还需要控制镜面的离轴量和离轴角.本文基于普通的二次曲面方程,推导了在新坐标系下的二次曲面表达式,并在这个坐标系下,利用Powell共轭方向法求取镜面的顶点曲率半径、二次曲面系数、离轴量和离轴角.通过对抛物面和双曲面两种形式不同外形的镜面数据的计算,拟合的RMS值可以达到10nm的精度,从而表明该测量计算方法可以有效地控制离轴面加工时的参数误差.