多普勒雷达速度场缺测区域填补技术的数值模拟

针对多普勒天气雷达探测风场时,常因大面积无回波区的存在,对提取平均散度和平均形变信息的精度有较大影响。利用Sirmans建议的技术来模拟降水回波信号,得到含有0、1、2阶谐波的平面线性速度场并进行加噪处理,得到近乎真实的多普勒速度场。对模拟的平面线性速度场人为设定无回波区,利用VAD技术和迭代法对平面线性速度场做连续性缺口和非连续性缺口的迭代法填补。数值模拟结果表明:对于模拟的无噪声平面线性速度场,迭代法可以基本无误差填补连续性累积缺口在10°到180°的速度场无回波区;对连续性累积缺口,在σ_v=2m·s~(-1),SNR≥5 dB和SNR=20 dB,σ_v≤4 m·s~(-1)的情况下,连续性累积缺口在120°以内,迭代法能较高精度地填补速度场无回波区,填补后0、1、2阶谐波误差绝对值基本能控制在15%以内,80 km距离圈上缺口处迭代前后速度相对误差均在30%以内;对非连续性累积缺口为0°～180°的平面线性速度场,在附加不同噪声条件下,迭代法均可较好地填补累积缺口在180°以内的平面速度场,且80 km距离圈上缺口处迭代前后速度值误差均能控制在15%以内。说明利用迭代法填补多普勒速度场无回波区,填补效果较好,精度较高,这对改善从多普勒速度场中提取平均散度和平均形变信息的精度有极大的帮助作用。     

利用四元数改进大旋转角坐标变换模型

大旋转角坐标变换模型的迭代解法依赖于初值的确定。用四元数构造旋转矩阵,建立了三维坐标变换的牛顿迭代公式,并提出了一种初值构造算法。利用实测数据和模拟数据对该算法进行了验证,并与其他算法进行比较。结果表明,该初值构造算法使得基于四元数的大旋转角坐标变换模型更加稳健。     

共轭梯度法在图像重建中的应用

图像重建常常被转化为解非线性无约束极值问题,通过范数极小化推导出共轭梯度法的一般算法。通过对模拟数据和实际工件断层扫描数据进行图像重建,估计了算法的有效性,结果表明,与最速下降法相比,此算法更适用于不完全投影数据的图像重建,在保证重建图像拟合度的同时,大大提高了重建速度。     

重磁位场转换计算中迭代法的综合分析与研究

处理转换计算在重磁资料解释中发挥着重要的作用,但一些计算如向下延拓、化极等有时是很不稳定的,在频率域中则表现为其转换因子具有明显的放大作用,所以其FFT理论计算结果是不稳定的.因此,很多研究工作都是围绕增加计算的稳定性、提高计算效果进行的,其中迭代法是近来在研究中受到普遍重视的方法技术,并取得了较好的成果.但也存在对迭代法研究还不够深入,对其存在的缺点认识不够充分、客观等问题,例如,迭代法进行延拓及化极等计算时,对一些具体应用虽能在一定程度上获得较好的计算结果,但却存在计算结果并不会随着迭代次数的增加而得到持续改善的问题,对于原本不稳定的计算,迭代法在迭代次数比较大时,所得的结果依然是不稳定的.为此,本文在对迭代法进行分析研究的基础上,进一步推导了迭代法的通式,并分析了对迭代法收敛性影响的各种因素.分析结果表明:迭代法收敛到FFT理论直接计算结果的决定因素是计算过程中如何选择原始数据到目标数据的映射函数;在选择了合适的映射函数的情况下,迭代次数不仅仅是决定计算成本,而是决定结果好坏的关键因素;增加迭代次数虽然能够使计算收敛到FFT直接计算理论结果,但如果该理论结果本身就是不稳定的,则迭代法计算如果收敛,也是收敛到一个不稳定的结果.所以针对位场处理转换中一些不稳定计算采用迭代法,并没有从根本上解决计算的不稳定性问题.     

位场向下延拓迭代法收敛性分析及稳健向下延拓方法研究

地磁导航作为一种新的无源导航方式,具有重要的国防意义.构建空间地磁数据库是实现地磁导航的基础,位场延拓是解决地磁数据库构建的有效方法.积分-迭代法是一种解决位场大深度向下延拓的实用方法.本文着重对积分-迭代法的收敛性进行了分析,从数学角度证明积分-迭代法能够收敛到直接下延法理论解.同时对积分-迭代法的抗干扰性进行了初步分析,当观测数据含有噪声时,积分-迭代过程中使得噪声得到累加,影响延拓数据的精度.本文利用正则化方法和递增型维纳滤波方法,提出了波数域位场向下延拓新算法.模型检验表明,新算法稳定、抗干扰能力强、计算速度快.     

迭代法建立的低频模型在储层反演中的应用

储层预测的精度直接制约着油气开发的经济性与有效性,而低频信息的补偿是改善储层反演效果的有利途径之一.常规的测井曲线内插低通滤波建模方法难以精细表征复杂地质背景下的低频模型.为降低研究区特殊地质体对储层反演带来的影响,采用迭代法建立低频模型,建模中分步考虑压实作用、特殊岩性等因素影响,并通过反演效果不断迭代更新低频模型,最终准确建立反映不同岩性(包含火山岩)的反演低频模型.研究表明:该方法能够保证储层反演结果的可靠性,为储层精细描述提供了资料;同时基于反演结果,对东海西湖凹陷W气田主力层(P1)有利储层进行重新刻画,优化了开发井水平段的位置,该井实施结果证实钻遇了优质储层.用迭代法建立的低频模型为复杂地质背景地区储层预测提供了一种更有效的建模方法.     

GPS精密星历的轨道内插方法比较

在高精度的GPS测量数据处理过程中,获取高精度卫星轨道是重要的环节,内插是获取任意历元的精密轨道信息的重要手段,对拉格朗日插值法、牛顿插值法、Neville逐次线性内插法和切比雪夫插值法的效果做了细致的比较,并在此基础上进一步探讨了短时间外推的效果。     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

迭代法井间地震衍射层析成像

从各向同性非均匀介质中弹性波波动方程出发，导出了在弱不均匀条件下，频率－波数域衍射场方程，给出了格林函数解析解。在电磁波地球物理层析成像数据合成的基础上，用类“体电流法”，将解析解线性化，得到未作Ｂｏｒｎ近似或Ｒｙｔｏｖ近似的衍射场数值解，便于计算合成数据。采用迭代法进行反演，每步迭代中，采用等比Ｎｅｕｍａｎｎ级数解；反问题用ＳＶＤ法求解，探讨了阻尼因子的取法。数值模型试验表明，这种算法是有效和成     

一种大地坐标变换的快速算法

大地纬度的解算是从地心直角坐标到大地坐标变换的关键，传统上多采用迭代法，效率较低，难以满足实时性较高的应用。给出一种非迭代的，采用有理多项式逼近的方法来计算大地纬度，并且通过与Heikkinen、改进Bowring、Ozone3种算法的比较证明这是一种速度极快的算法，而且在50km以下的范围内，采用这种算法的平均误差不超过1mm。