基于SP-LSQR的综合孔径辐射成像方法

综合孔径微波辐射计是被动微波遥感领域中一种重要的探测器。综合孔径辐射计成像反演过程在数学上是病态的反问题,传统正则化方法虽然能够有效克服其病态性,但重构误差较大、对干扰噪声鲁棒性不强。与传统正则化方法相比,最小二乘Q-R矩阵(Least Square Q-R matrix,LSQR)分解算法具有计算精度高、数值稳定性好等优点,但收敛速度不够快。对此,文中提出利用子空间预处理LSQR(Subspace Preconditioned LSQR,SP-LSQR)算法进行综合孔径辐射计反演成像,提高收敛速度。仿真结果表明:与Tikhonov正则化和带限正则化相比,SP-LSQR方法不仅可以降低重构误差,而且对干扰噪声的鲁棒性更强。此外,与LSQR方法相比,SP-LSQR方法在不降低成像反演精度的情况下,有效提高了计算效率。

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解).在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组.由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈.针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现.导出了LSQR算法的并行效率估算公式.对两个地区的实际地震层析成像数据进行了效率测试.

基于小波多分辨分析和LSQR的快速EIT图像重建算法

针对LSQR(least square QR-factorization)算法在求解电阻抗层析成像(electrical impedance tomography,EIT)逆问题时,由于矩阵维数高、计算量大而导致重建速度较慢的问题,提出基于小波多分辨分析的LSQR算法(wavelet multi-resolution based least square QR-factorization,WALSQR)。该方法将EIT的图像重建过程投影到低维的尺度空间进行,通过提取有效信息减少数据计算量,明显提高了图像重建速度。同时由于去除了噪声和冗余信息,保证了成像质量。本文将SALSQR方法分别应用于二维、三维EIT成像实验,证明其有效性,并为三维动态EIT图像重建算法的研究奠定了基础。

跨孔层析成像LSQR算法研究

目前广泛应用于工程、环境、水文等领域的探地雷达跨孔层析成像技术具有分辨率高,探测范围广,得到的图像结果直观,便于工作者分析地下地层结构等特点。层析成像的关键在于解系数矩阵为大型稀疏矩阵的方程组,笔者采用LSQR算法。此方法具有速度快,对不适定问题数值稳定以及从迭代过程很容易求得数值分析的数值的优点,用其解大型稀疏矩阵方程组得到的的结果清晰可靠,在各种层析成像中LSQR算法得到了广泛的应用。

基于阻尼LSQR算法层析成像技术在铸件缺陷检测中的应用

在普通超声波检测的基础上引入层析成像技术,形成新的超声波层析成像探伤系统,可以有效克服普通超声波检测定位困难的缺点,能准确检测铸体中缺陷的大小、形状和位置。本文重点对超声波层析成像探伤系统的工作原理和阻尼LSQR算法(DLSQR算法)进行研究,经过数据采集、正演模拟、建立走时方程、反演求解等过程,实现了基于DLSQR算法的图像重建方法。实验表明,DLSQR算法具有良好的稳定性、较小的存储容量和计算量,适合病态问题求解。

复杂岩溶地区跨孔地震CT反演成像研究——以贵阳在建地铁3号线为例

本文以贵阳在建地铁3号线为例,详细研究岩溶复杂地区勘察过程中跨孔地震CT技术的应用,揭露研究区岩溶发育特征。通过对不同算法的对比分析表明,相比于CGLS算法,LSQR算法对于低速异常的识别更灵敏和精确,具有更好的收敛性和可靠性。钻探验证表明,采用LSQR算法和弯曲射线追踪法进行地震层析成像反演,跨孔地震CT探测可达米级精度,获得满意的速度图像,能够查明钻孔之间岩溶埋深、规模、延伸等情况,弥补钻探的局限性,说明基于LSQR算法的地震CT技术探测岩溶是切实可行的。

图像去噪的混合正则化LSQR算法

本文研究了图像去模糊去噪问题.利用正则化技术结合Krylov子空间方法,提出了混合正则化LSQR算法.实验结果表明该算法有效改善了问题的不适定性,获得了逼真度较高的复原图像.

基于LSQR算法的二维声波方程频率域正演模拟与数值实现

频率域正演是频率域全波形反演的基础,高效准确地完成频率域正演计算是目前该领域的一个热门问题。这里利用优化17点差分格式和PML完全匹配层构建波场,并使用LSQR算法求解波场,通过在简单层状模型和Overthrust模型上计算验证可知,LSQR算法可大幅提高频率域正演的计算速度和模拟精度,从而为此方向的研究提供一个可行的技术手段。

基于阻尼LSQR-LMBC的火焰三维温度场重建

光场相机可以解决辐射测温多相机系统光路复杂、同步触发难等问题,在辐射成像三维温度重建时有其独特优势.LSQR是求解基于大型稀疏矩阵最小二乘问题的经典算法,该算法用于重建三维温度场时对温度初值依赖较大,在信噪比较低的情况下重建精度不理想.本文提出阻尼LSQR-LMBC重建算法,通过在LSQR方法中添加阻尼正则化项,提高火焰三维温度场重建的抗噪性能,并结合LMBC算法,实现吸收系数和三维温度场同时求解.在数值模拟部分,随着信噪比逐渐降低,阻尼LSQR的重建效果比LSQR更加稳定,在信噪比达到13.86 dB时,重建精度大约提高30%.阻尼LSQR-LMBC的平均重建误差为6.63%.用丁烷火焰进行了实验,重建的丁烷火焰三维温度场分布符合辐射火焰燃烧的特征,和热电偶的测温数据结果进行对比,相对误差在6.8%左右.

基于LSQR与系统距离聚类算法的专变用户非介入式负荷辨识

专变用户非介入式负荷辨识是实现以专变用户为主体的需求侧响应的重要环节。通过对专变用户负荷稳态过程与暂态过程的特征分析,选取有功功率、无功功率、电流有效值等作为稳态过程特征量,并选取暂态发生前、中、后三个阶段的特征量,如电流有效值、有功功率均值、无功功率均值、持续时间、电流有效值最大值等,构建全面的非介入式负荷辨识稳态过程与暂态过程负荷特征空间。以此为基础,利用最小二乘QR分解算法(Least Square QR,LSQR)进行稳态过程负荷分解获得各种负荷的运行情况。并基于系统距离聚类算法将常见专变用户负荷暂态事件进行分类,进一步辨识出哪一种或者哪一类负荷发生投切动作。最后,采集包括排水泵、搅拌机、水泥螺旋、污水泵、除尘器等工业负荷现场数据,对所提到的方法进行了仿真分析,结果表明所提的非介入式负荷辨识方法可获得较高的准确率。

双选信道下OCDM系统低复杂度均衡

正交Chirp复用(OCDM)是近年来提出的一种新的多载波体系,通过菲涅尔变换,获得一组正交Chirp信号,实现了CSS的最大频谱效率。该文介绍了OCDM系统的基本原理,重点研究了OCDM系统的低复杂度均衡算法。双选信道下,传统的MMSE均衡算法性能下降,提出一种基于近似带状矩阵的阻尼LSQR算法,作为求解稀疏矩阵的最小二乘迭代算法。为了缓解快速时变信道中的ICI,提出一种基于近似带状矩阵的LSQR-BDFE算法,结合判决反馈均衡,通过LSQR算法迭代计算。仿真结果表明,双选信道下,OCDM系统比OFDM系统有着更好的BER性能,所提出的LSQR-BDFE算法和带状阻尼LSQR算法,比MMSE均衡算法有着性能优势。

用改进的光滑约束最小二乘正交分解法实现电阻率三维反演

对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解(LSQR)法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法.

井间地震层析成像技术的应用

阐述了井间地震层析成像的正反演技术——基于最短路径法的射线追踪正演技术和基于LSQR算法的反演技术的原理,并且通过数值模型对正反演算法进行了较细致的研究,重建了地下构造原像,效果较好。结合工程应用实例,进一步表明井间地震层析成像技术能够更清晰、更客观地反映出地下构造情况,为地质结构研究和后期施工提供更详实的资料。最后,对层析成像的分辨率进行了探讨,指出了层析成像质量的影响因素。

井间地震层析成象高精度算法及其在环境工程勘探中的应用

井间地震层析成象在环境工程勘探中的应用前景广阔，走时层析成象的射线追踪方法和高精度的解大型稀疏矩阵的方法是其中的关键。文中讨论了ＬＳＱＲ算法系列，数值模拟证明了它的优越性，最后应用于鹏飞大厦场地。

利用地震和爆破资料反演宁夏及邻区浅层三维速度结构

将有限差分正演及LSQR反演方法应用到宁夏及邻区,利用该区2003—2009年7年间的地震和爆破记录,经过重新细致分析,得到了具有4个以上台站清楚记录到初至波Pg和Pn到时的地震426次、爆破52次,剔除到时离散较大的数据后共有2788条记录.获得了不同深度的速度平面分布图像以及不同位置的深度速度剖面图像.得出深度0—5km的平均速度约为5.1km/s,深度5—25km的平均速度约为5.9km/s,深度25—30km的平均速度约为6.4km/s;浅层低速层主要分布在银川地堑、中卫清水河盆地,高速层主要分布在银川地堑西部的贺兰山区.此结论与该地区地质结构相符.深层高速层在中卫—同心断裂带一带最为显著,低速层主要分布在银川地堑及以西部分地区、固原清水河断裂带等.

互相关系数自约束的重力三维反演与高效求解

本文将拟合残差计算所得的互相关系数作为先验信息,与深度加权函数同时引入到重力正则化反演的模型约束中,以提升反演结果的可靠性。针对三维反演中的大型线性方程组问题,引入阻尼LSQR(最小二乘QR分解)算法,结合等效几何格架技术,将大矩阵按照模型单元划分为若干个子矩阵进行存储与运算。理论模型计算结果表明:同时利用互相关系数和深度加权的模型自约束反演,能较清晰地反映真实异常体;基于分块矩阵的阻尼LSQR算法求解线性方程组较直接法节省了几千甚至上万倍的存储量,且计算速率提高了数倍,可在普通计算机上实现较大规模的反演计算。将其应用于云南芦子园铁铅锌铜多金属矿床隐伏花岗岩体定位,取得了良好的效果。

高分辨率过井地震层析成像的级联算法及其应用

将过井地震层析成像的级连算法应用到垦71区井间地震资料处理中,得到了高分辨率、精确的井间纵波速度结构。首先根据初至出射角,旋转合成得到清晰的纵波波场,再精确拾取初至到时。初始速度模型是参考测井资料的分层模型,第一级反演算法为分辨率较低但稳定的BPT反演算法,为第二级高分辨率反演算法——阻尼LSQR迭代反演算法提供模型。层析结果显示出清晰的泥、砂互层,与测井参数构造解释对比,高、低速带与砂/泥岩层基本吻合,最高能分辨1/4地震波波长的薄层。研究证明级连算法能很好地兼顾高分辨率和高精度,大大提高了井间地震走时层析成像的质量,可以广泛地应用于油田的勘探和检测。

核磁共振双TW测井数据联合反演与流体识别

针对核磁共振测井双TW观测数据分析和流体识别的需要,研究了基于全局搜索的遗传算法和局部搜索的最小二乘法的联合反演算法,实现了核磁共振双TW观测数据处理.首先,研究了饱和油气水岩石物理模型的核磁共振双TW观测模式的测井响应机理;然后,利用全局搜索性能优良的遗传算法,对核磁共振回波差数据进行了反演,计算出了流体的核磁弛豫性质及其体积;最后,以遗传算法的反演结果为初值,利用阻尼最小二乘方法对双TW回波串进行更精细的反演,计算出了双TW的T2分布、孔隙度和流体饱和度.理想模型的合成数据和实际测井资料应用表明,遗传算法与最小二乘方法相结合是一种行之有效的联合反演方法,能较好地实现核磁共振测井双TW观测数据的处理和流体评价.

Determination of 3-D Velocity Anomalies of the Nanbei Tectonic Zone of China Based on Local Earthquakes

: In this paper, 3-D velocity images of the crust and upper mantle beneath the Nanbei tectonic zone of China are constructed using P-wave travel time residuals of earthquakes, with the data supplied by China's seismic networks and the International Seismic Centre.

反演近地表物性参数的地震层析成像方法综述

针对反演近地表物性参数的地震层析成像方法的发展历程开展了论述。首先,阐述了传统的联合代数重建法(SIRT)、最小二乘QR分解(LSQR)等求解方程组的方法以及最新发展的局部最优化算法(包括伴随状态法与改进的散射积分法)的模型更新方式;其次,从简单到复杂论述了层析成像方法在射线理论、有限频理论、波动理论下的发展历程;再次,从初至波到时频率依赖的初至波到时初至波波形反演方面阐述了反演精度逐步提升的技术关键;然后,介绍了目前对地下物性参数要求提升的情况下,地震层析成像方法从单参数速度到各向异性等多参数反演的最新研究动态;最后,分析了地震层析成像这种高度非线性的反演过程的模型参数化、反演策略、Hessian矩阵利用等问题。综上所述表明:地震层析成像在实际应用中仍具有不可替代的作用;通过不断挖掘数据中更多的信息,其理论方法与实现方式正在向全波形反演靠拢。