利用计算机断层成像方法由观测图像重建等离子体层全球密度分布——投影数据不完备问题

在用计算机断层成像方法由EUV观测图像重建等离子体层全球密度分布时,地球的遮挡和有限角度都会导致投影数据不完备,从而无法精确重建出等离子体层的密度分布.本文针对该问题,提出一种基于图像总变差极小化的代数迭代算法.通过重建等离子体层投影数据缺失最为严重的中心子午面,证明该算法能够显著提高重建图像的质量. 并且在IMAGE卫星仅能达到90°的有限投影角度下,此算法重建图像的相关系数可达0.760,而代数迭代算法的相关系数仅为0.696.     

Underground water quality model inversion of genetic algorithm

The underground water qualitymodel with non-linear inversion problem is illposed, and boils down to solving the minimum of nonlinear function. Genetic algorithms are adopted in a number of individuals of groups by iterative search to find the optimal solution of the problem, the encoding strings as its operational objective, and achieving the iterative calculations by the genetic operators. It is an effective method of inverse problems of groundwater, with incomparable advantages and practical significances.     

GIS数据误差处理的阻尼最小二乘算法

应用非线性最小二乘平差方法--阻尼最小二乘法研究了数字化数据误差处理,并将该方法与最小二乘条件平差进行了比较.理论和实验表明,在处理较大数字化数据误差时采用阻尼最小二乘平差方法,减弱了在条件方程线性化过程中产生的模型误差,提高了计算精度.     

频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

海上倾斜缆采集技术具有多样的陷波特征,通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题,本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法,该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明,频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波,能达到拓宽地震记录频带目的.     

病态总体最小二乘问题的谱修正迭代法

针对病态问题岭估计法及正则化方法存在破坏方程的等量关系、解是有偏估计的问题,在最小二乘谱修正迭代法的基础上,提出病态总体最小二乘问题的谱修正迭代法,并推导了总体最小二乘谱修正迭代法及其改进算法的具体公式。通过算例验证和阐明了TLS谱修正迭代法在病态问题中的有效性和易受迭代初值影响的缺点。     

地震随机噪声压缩感知迭代压制方法

复杂地表和复杂介质条件下,随机噪声往往严重影响着复杂地震信号的信噪比,同时深层地球物理目标探查中弱地震信号总是被随机噪声所掩盖,如何有效地压制随机噪声干扰、恢复有效地震信号仍然是高精度地震勘探中的关键问题.压缩感知理论突破了奈奎斯特采样定理的限制,利用有效地震信号的可压缩性和稀疏性,提供了从不可压缩随机噪声中进行有效信号分离的数据原理.本文系统分析压缩感知框架下地震随机噪声压制的稀疏优化反问题,提出了基于迭代软阈值算法的"采集-重建-修复"方案对该问题进行求解.在实现高度稀疏表征的基础上进行地震数据的压缩感知随机观测,通过迭代反演对有效地震信号进行重构,有效提高复杂地震数据的信噪比,同时,当求解稀疏优化问题时,如果出现正则化项引起重构信号衰减现象,可以匹配除偏对衰减的有效信号进行修复.通过与工业标准f-x预测滤波方法进行比较,理论模型和实际数据处理的结果表明,压缩感知迭代噪声压制方法对复杂地震数据中的随机噪声有较好的压制效果,可以有效恢复出被较强非平稳随机噪声干扰的时空变同相轴信息.     

Ensemble retrieval of atmospheric temperature profiles from AIRS

ABSTRACT Satellite-based observations provide great opportunities for improving weather forecasting. Physical retrieval of atmo spheric profiles from satellite observations is sensitive to the uncertainty of the first guess and other factors. In order to improve the accuracy of the physical retrieval, an ensemble methodology was developed with an emphasis on perturbing the first guess. In the methodology, a normal probability density function （PDF） is used to select the optimal profile from the ensemble retrievals. The ensemble retrieval algorithm contains four steps： （1） regression retrieval for original first guess; （2） perturbation of the original first guess to generate new first guesses （ensemble first guesses）; （3） using the ensemble first guesses and nonlinear iterative physical retrieval to generate ensemble physical results; and （4） the final optimal profile is selected from the ensemble physical results by using PDE Temperature eigenvectors （EVs） were used to generate the pertur- bation and generate the ensemble first guess. Compared with the regular temperature profile retrievals from the Atmospheric InfraRed Sounder （AIRS）, the ensemble retrievals RMSE of temperature profiles selected by the PDF was reduced between 150 and 320 hPa and below 400 hPa, with a maximum improvement of 0.3 K at 400 hPa. The bias was also reduced in many layers, with a maximum improvement of 0.69 K at 460 hPa. The combined optimal （CombOpt） profile and a mean optimal （MeanOpt） profile of all ensemble physical results were improved below 150 hPa. The MeanOpt profile was better than the CombOpt profile, and was regarded as the final optimal （FinOpt） profile. This study lays the foundation for improving temperature retrievals from hyper-spectral infrared radiance measurements.     

湍流通量参数化方案的非迭代方法研究

基于Högström (1996) 和Beljaars et al.(1991) 的研究工作, 沿用Louis et al.(1982) 和Launiainen (1995) 的思路, 本文采用多元回归分析方法, 研发了一种采用非迭代方法的湍流通量参数化方案。该方案直接用整体理查森数、 空气动力学粗糙度长度和热力学粗糙度长度对稳定度参数进行参数化, 从而避免了通过循环迭代计算Monin－Obukhov长度。该方案不仅有效地节省了CPU计算时间, 而且其计算结果与迭代方案 (BHH方案) 的计算结果非常接近。     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.     

主动源与被动源地震数据插值及联合数据成像

本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.