反射系数与子波同时迭代反演的预条件共轭梯度法

在反射系数白噪、子波最小相位的假设下,研究基于线性反演的地震反射系数和子波同时估计问题。在Cauchy准则稀疏反演求解中,应用预条件共轭梯度法实现反射系数和子波同时迭代反演。在迭代求解正则化方程时,用共轭梯度法求解相应的原问题,初猜子波求解也使用该策略。模型数据试算与比较,表明了该算法正确而有效。用实际数据检验算法的实用性,经研究表明,预条件共轭梯度法计算的反射系数和子波,要比直接稀疏反演精度高,而且收敛较快,数值稳定,实用性强。     

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

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

梯度反褶积法及其在矿产勘探中的应用

梯度探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用,但现今针对梯度数据的反演方法还比较少。笔者提出一种新的梯度数据解释方法--梯度反褶积法,该方法以梯度数据为基础获得场源体的位置信息,在计算公式中消除了构造指数,避免了因构造指数选取不当引起的误差。通过理论模型试验证明,梯度反褶积法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于5%,且相对于常规欧拉反褶积法更加稳定、准确。将梯度反褶积法应用于实测航磁梯度异常的解释,获得了地下铁矿的分布状态。     

广义共轭梯度算法

在地球物理反问题的求解过程中,共轭梯度(CG)法是一种经典的、很有价值的主要算法之一。本文在经典的共轭梯度标准基础上,将进一步发展,推导出了求解阻尼最小二乘(LS)解和最小绝对值偏差(LAD)意义下的迭代再加权最小平方(IRLS)解的标准算法形式,从而使得CG法的应用更具一般性。为了更好地理解GCG法的性能,文中还给出了两个例子,并将计算结果与公认的、好的求解病态问题的奇异值分解(SVD)算法的计算结果进行了比较,结果表明:GCG法亦具有很强的求解病态问题的能力,精度高,且运算速度快。此外,GCG法还具有两个显着的特点:①算法简单、编程灵活;②可以保持系数矩阵的稀疏特征。     

基于互相关误差泛函的最小二乘逆时偏移方法

最小二乘偏移可以得到成像质量更高的剖面,有利于更好地进行岩性储层成像和储层参数反演。但是,最小二乘偏移在实用化过程中却很难实现其高精度岩性成像的目标。在分析传统的最小二乘偏移存在的问题后,发展了一种基于相位匹配的最小二乘偏移方法,该方法舍弃了传统的二次型泛函,采用互相关泛函,更多地强调相位信息在最小二乘偏移中的作用,在理论上可以更好地适用于实际资料。通过对互相关误差泛函构建、逆时反偏移合成预测数据、梯度预条件处理、共轭梯度法迭代求解等方面进行研究,探索实现了基于互相关误差泛函的最小二乘逆时偏移方法。简单凹陷模型和复杂salt模型测试表明,本文方法具有较好的适用性和有效性。     

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

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

预条件共轭斜量法及其在求解边值问题中的应用

在论证共轭斜量法误差估计式的基础上，为提高敛速，对系数矩阵进行预处理，提供了减少等价问题条件数的方法，完美地建立了预条件共轭斜量法的实用算法。最后，以泊松方程边值问题为例，通过数值实验说明了该方法的有效性。     

基于结构匹配反褶积的同时震源数据直接成像

针对同时震源数据最小二乘偏移(LSM)计算效率低的问题,笔者研究并提出了同时震源混合数据直接成像的结构匹配反褶积(SMD)方法。通过分窗的方式求取局部Hessian逆,加入局部结构约束,可以压制混合数据成像串扰噪音的同时对成像做照明补偿,提高成像质量。并在复杂sigesbee2A模型的数值测试中验证了SMD方法的有效性。结果表明,SMD能够有效地压制成像中的串扰噪音,且对成像振幅有一定的补偿作用,计算效率相比LSM有巨大提升。     

共轭梯度法在新井参数设计中的应用

通常情况下,含水层的参数可分别借助于Dupuit和Theis井流公式,在局部地段利用稳定流和非稳定流抽水试验获得.在新井设计阶段,一般是依据现有参照井数据进行线性或非线性插值计算,但插值受插值表的约束,效率不高.本次研究采用2次优化方法,用网格法确定所求参数最优估计值作为下一步的叠代初值.用共轭梯度法逐步逼近所求参数的精确解,高效地利用参照井的基本数据得到了精度较高的新井参数,为准确地设计新井提供可靠的依据,实地应用效果良好.     

解大型稀疏方程组的ICCG方法及其计算机实现

针对ＩＣＣＧ算法中的关键步骤,提出了快速有效的计算技巧,以利于实际应用。     

地磁场水平梯度的计算与分析

In this paper，based on the Gauss theory of geomagnetism and the spherical harmonic mi}del of the geomagnetic field，the horizontal gradients of the geomagnetic field in north一south and west-;east dlrections were calculated in Chinese and"its surrounding regions, the contour maps of the horizontal gradients of the geomagnetic components were constracted，and its regularities of distributions were analysed.
The horizontal gradients of the total intensity (F)，horizontal intensity (H) and vertical intensity (Z ) vary mainly with lat相似文献   

估计磁场源位置的倾斜角梯度算法改进

倾斜角梯度法在未知场源物性先验知识和构造指数情况下,可快速估计出场源的边界位置和深度分布,对于大面积实测磁异常的快速解释具有应用前景。笔者在Salem的倾斜角梯度法基础上,提出了一些算法改进。在导数换算上采用空间域有限差分法换算磁异常的水平导数,采用ISVD算法换算磁异常的高阶垂直导数。在剔除坏解上,采用磁异常广义倾斜角梯度的总水平梯度滤波法和场源深度有效性筛选法相结合。理论模型数据试验表明本方法简单快速,效果明显。     

短期电力负荷的组合预测方法

电力负荷的影响因素很多,用单一预测模型进行预测时,因自身的局限,使其预测精度和稳定性不高。组合预测可以综合考虑各单一预测模型的结果,提高预测精度和稳定性。在电力预测中常用的指数平滑法、趋势分析法、Box-Jenkins时间序列法以及灰色预测法的单一预测基础上,建立组合预测模型,采用最小二乘法确定组合预测模型权重,并对某水电电力公司的上网负荷进行预测。实例表明,组合预测法能在短期电力预测中改善预测精度,提高预测的稳定性。     

倾斜板体磁异常总梯度模反演方法

对于二度体磁异常，其总梯度模具有不受磁化方向影响、表达式简单和具有较强分辨叠加异常的特点，针对顶面倾斜板状体提出的梯度模反演方法，能反演倾斜顶面中心的深度、顶斜面的倾角和宽度等参数。     

双频梯度测量探查井巷漏水水力联系

介绍一种利用双频仪解决矿区井巷巷道中漏水水力联系的原理和方法。双频梯度测量法有抗干扰能力强、分辨率高、仪器轻便等特点,实践证明应用效果很好。     

水力坡度在层间氧化带砂岩型铀矿区域预测中的意义

本文根据地下水水动力参数之一──水力坡度所反映出的地下水动力特征和地下水径流强度以及对层间氧化带发育的影响，试图找出适合于层间氧化带发育的水力坡度范围，从而达到利用这一地下水水动力参数指导层间氧化带砂岩型铀矿区域预测工作的目的。     

成分数据的因子分析及其在地质样品分类中的应用

地质样品的化学组成都是成分数据。确定一组混合源成分数据的端元数目、端元的化学组成及各样品中含各端元的绝对份额是地质样品成因分类的关键。求解这一问题相当于确定模型X=LB(lij≥ 0 ,∑lij=1)中X矩阵的秩、B矩阵的行向量和L矩阵的列向量。在确定X矩阵的秩和初始B0 矩阵的基础上 ,用限定最小二乘法循环调整矩阵B0 ,使其bij≥ 0和∑bij=1而得矩阵B ,同时也解得了矩阵L。X矩阵的秩和B0 是在传统的Q型因子分析基础上 ,求最优斜交因子解得到的