实现稀疏反褶积的预条件双共轭梯度法

地震勘探稀疏反褶积计算一般要导出一个Toeplitz矩阵的线性系统，通常可以用矩阵求逆、Levison递推及共轭梯度等方法直接求解。当Toeplitz矩阵的条件数很大时，数值稳定性差，甚至无法求解。使用共轭梯度法，在矩阵的对角元素上加入规则化因子，可以改善这种情况，但不能彻底解决数值稳定性和精度问题。若求解最小二乘问题的原始问题，结果会好些。线性系统形式的细微改变，将导致不同的数值计算特性。在规则化策略基础上，可巧妙地构造稀疏反褶积的问题原型，引入预条件，采用双共轭梯度法求解，从而实现稀疏反褶积，获得较好结果。数值算例表明，预条件双共轭梯度法比直接稀疏反褶积方法收敛快、精度高。     

基于预优共轭梯度的磁化强度反演方法及应用

基于预优共轭梯度算法的反演方法在求解磁化强度大规模欠定线性方程组时可以改善方程组的条件数,并提高算法的纵向分辨率。该方法对磁性体形状无要求,且收敛速度快、稳定性好,能较直观地反映磁性体特征。本文通过理论模型试验结果分析证明了基于预优共轭梯度的磁化强度反演方法的稳定可靠性,北衙铁金矿区万硐山矿段L32线磁异常反演结果表明了该方法应用于实际资料处理与解释的有效性和可靠性。     

使用共轭梯度的三维电阻率正演模拟与反演

研制了利用共轭梯度松驰法的三维直流电阻率正演模拟与反演算法。在正演网络的模拟计算中，将一种用于预处理的不完全Ｃｈｏｌｅｓｋｙ分解法与稀疏矩阵程序结合起来，便可得到一种快速有效的算法。侧面和底部的边界条件用阻抗条件来标度，而阻抗条件考虑了由于电源流的一些配置而在边界上产生的电流。     

基于共轭梯度法的全波形反演

全波形反演是一种全新的地震成像方法,主要利用全波形信息反演地下介质参数,通过非线性优化波场理论值和观测值的残差实现波形反演。基于时间域声波方程,建立了时间域波场残差目标函数,以分层模型为例,分别从波场精度、目标函数收敛性和运行时间3个方面,比较了共轭梯度（Conjugate Gradient,CG）算法和拟牛顿算法（Broyden-Fletcher-Goldfarb-Shanno,BFGS）反演的效果。同时,应用共轭梯度法对正、逆断层模型和Marmousi模型进行了速度结构反演。反演结果表明：共轭梯度法计算效率较高,反演得到的速度模型精度更高,反演效果较好,是一种有效的波形反演方法。     

基于共轭梯度法的垂直有限线源三维电阻率反演

利用垂直有限线源研究油田注水分布和剩余油分布成为一种新型电测方法。本文对该方法实际应用的有效性进行了分析和说明,并利用共轭梯度迭代技术实现了垂直有限线源三维电阻率反演。结果表明,反演稳定可靠,计算速度快。     

利用共轭梯度方法的激发极化三维快速反演

利用共轭梯度方法实现了激发极化(IP)三维快速反演。首先,利用共轭梯度方法反演电位数据,得到地下的三维电阻率模型,(由于避免了直接求偏导数矩阵,因此反演计算速度很快。)然后,以此电阻率模型为地下电导率分布,再反演IP数据得到三维极化率分布理论模型。试算结果表明其效果较好。     

用改进的共轭梯度算法解波阻抗反演问题

研究快速有效的波阻抗反演算法,是地震反演中一直追求的目标。这里从目标函数出发,推导的反演方程,在常规共轭梯度算法的基础上,对算法搜索的共轭梯度方向提出了具体的改进措施,并在阻尼项及松弛因子的给定上实现了自适应化。另外,在整个算法过程中实现了向量运算,无需额外的内存。通过合成记录及实测资料试算对比发现,该算法具有精度高,运算速度快及抗病态能力强的特点。     

井地大地电磁非线性共轭梯度二维反演研究

本文提出并验证了地下大地电磁测深数据的正、反演算法和程序,为井地大地电磁测深和井地可控源音频大地电磁测深这一前沿研究课题提供了地球物理解释方法.使其更加完善和有效.在Rodi(2001)发表的大地电磁场二维NLCG反演理论和算法的基础上,推导了井地大地电磁场二维NLCG反演算法,改编反演的代码,开发出二维NLCG反演程...     

ZTEM二维非线性共轭梯度反演研究

ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。     

网格剖分对大地电磁测深反演精度的影响

网格剖分的合适与否,影响着反演结果的精度。用不同模型、不同反演方法,对比验证粗细两种不同的网格对反演精度的影响。研究结果表明,从整体上看两种不同网格剖分方式,粗网格在深部比细网格模拟精度高,但是近地表开始阶段,粗网格与细网格波动幅度都较大,都脱离正常值,粗网格模拟精度不如细网格;粗网格的整体变化幅度比细网格缓和,其中对低阻异常体反演精度较好,但对高阻异常体反演精度较差,与真实值相去甚远。     

电阻率三维反演稳定性和可靠性研究

电阻率三维反演参数太多,将导致反演不稳定性和高度非唯一性等诸多方面的问题。比较传统的最小二乘反演方法,利用共轭梯度的电阻率三维最小构造反演显示了更好的效果。讨论了初始模型选择及数据误差对三维反演结果的影响,说明了反演方法的稳定性和可靠性。     

隧道勘探AMT数据二维非线性共轭梯度反演的关键参数探讨

研究了二维非线性共轭梯度反演对于铁路隧道音频大地电磁(AMT)勘探的适应性。通过模型反演实验,确定控制反演进程的关键参数选取方法,以便于在不损失分辨率的前提下避免假异常。实验结果表明,在类似的地电条件下:①采用TE+TM模式联合反演;②使用浅部细分网格以更好地拟合地形;③不使用1000 Hz以上易受干扰的数据;④使用较小的正则化因子;⑤选用中、浅层电阻率的均匀空间作为初始模型;以上这些手段能提升反演结果的准确性,更适应铁路隧道勘探解释需要。     

基于深度加权的重力梯度张量数据的3D聚焦反演

针对重力梯度数据聚焦反演结果中存在的"上漂"现象,在经典Tikhonov正则化理论框架下,引入最小支撑泛函数对反演模型进行约束以避免反问题解的不稳定,并针对重力梯度数据聚焦反演中存在的趋肤效应,在模型目标函数中引入指数深度加权函数。通过理论模型,对部分重力梯度张量分量进行了单独以及联合聚焦反演,验证了基于深度加权的聚焦反演方法的有效性,并将该反演方法运用到涩北一号气田区的实际数据中,反演结果较好地反映出气田位置。     

基于遗传算法的叠前地震波形反演构建虚拟井曲线

叠前地震波形反演能够提供详细的地下地层特征，但由于其计算量大、数据和模型之间的非线性、目标函数的多极值和反演结果的多解性使叠前地震反演的实施成为一大难点。本文通过叠前和叠后混合反演技术采用遗传算法实现了深海无井地震反演。遗传算法思想简单、易于实现和使用、具有隐含并行性和全局搜索能力等优点，基于遗传算法的叠前地震波形反演得到了与井中实际数据基本吻合的速度和密度数据。     

基于地震叠后反演与属性分析的页岩气有利储层预测

中国南方页岩气勘查区块多为勘探新区,仅部署少量二维地震和地质钻井,且资料品质普遍较差,页岩气储层预测很难像成熟勘探区一样开展系统性研究工作,如何对现有数据进行充分挖掘是解决勘探新区页岩气储层预测的关键。为此,综合利用不依赖初始模型的叠后约束稀疏脉冲反演方法和多种频率类属性分析技术,开展页岩气有利储层预测。地震反演能够获得多种参数参与储层预测,属性分析能够从多种角度识别储层,通过多种信息综合应用和叠合分析以提高储层预测精确性和有效性,为新区页岩气勘探提供了一种新的思路。该项技术应用于南方某区块,圈定了页岩气的有利目标区,在该区部署的两口水平井均获得高产工业气流。