共查询到20条相似文献，搜索用时 46 毫秒
1.
使用偏移垂直地震剖面(VSp)的旅行时间,应用迭代最小平方反演方法可以估算地下的二维倾角。利用合成数据进行试验表明,用直达波或一次反射的旅行时间,由这种反演技术能够估算一个钻井附近的倾角。这种反演方法实质上是一种“剥层”方法,即先估算浅层的倾角,然后估算深部的倾角。几个实例表明,当震源偏移距与层深度的比率很小时,一次反射是主要的波。旅行时间反演还需要精心算层速度和正确的校正静。如果不考虑这些困难和普遍存在的非单值性问题,旅行时间反演可以产生一个正确的地层模型,以用于实际数据的检验。 相似文献
2.
广义共轭梯度算法 总被引：11，自引：1，他引：11
在地球物理反问题的求解过程中,共轭梯度(CG)法是一种经典的、很有价值的主要算法之一。本文在经典的共轭梯度标准基础上,将进一步发展,推导出了求解阻尼最小二乘(LS)解和最小绝对值偏差(LAD)意义下的迭代再加权最小平方(IRLS)解的标准算法形式,从而使得CG法的应用更具一般性。为了更好地理解GCG法的性能,文中还给出了两个例子,并将计算结果与公认的、好的求解病态问题的奇异值分解(SVD)算法的计算结果进行了比较,结果表明:GCG法亦具有很强的求解病态问题的能力,精度高,且运算速度快。此外,GCG法还具有两个显着的特点:①算法简单、编程灵活;②可以保持系数矩阵的稀疏特征。 相似文献
3.
有限单元和边界单元交连耦合法在水文地质参数反演中的应用 总被引：1，自引：0，他引：1
这里提出一个新的水文地质参数的数值反演方法。这个反演方法是以有限单元和边界单元的交连耦合法(CFEBEM)以及最小二乘法(LSM)为基础构成的。CFEBEM方法同有限单元法(FEM)一样灵活,同边界单元法(BEM)一样精确。应用该反演方法不但可以确定水文地质参数的分布,同时也可获得地下水流场的计算结果。文中还研究了三种反演迭代技术。在太原地下水系统研究中的应用说明这种方法是十分有效的。 相似文献
4.
5.
6.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献
7.
反演问题一直是地球物理的重点研究问题之一。本文以长导线源半航空瞬变电磁为例,对时间域半航空数据进行一维反演研究。将自适应正则化(ARIA)反演方法和阻尼最小二乘(LM)反演方法相结合,提出一种改进的反演方法:自适应正则化-阻尼最小二乘(ARIA-LM)反演算法。本文选用三层理论模型对比阻尼最小二乘算法、自适应正则化反演算法与本文方法,并使用四层理论地层模型来验证算法的有效性和可靠性。结果表明:该算法在一定程度上对高阻分辨率增强,能有效提高反演的精度及稳定性。 相似文献
8.
基于广义逆和奇异值分解理论,研究核磁共振T2谱反演的截断法和阻尼法。首先给出反演问题的M-P广义逆解;然后对系数矩阵进行奇异值分解求其M-P广义逆,对奇异值进行截断或者加阻尼,保证解的稳定性兼顾其分辨率;最后利用迭代技术实现解的非负约束。数值仿真实验和岩心T2谱反演表明:对于双峰谱,在信噪比SNR≥10时,截断法和阻尼法反演T2谱的分辨率和稳定性都很高;截断法截断因子选择范围小,在SNR<10时解的分辨率和稳定性差;阻尼法阻尼因子选择范围大,在SNR<10时解的稳定性好,但短组分分辨率低,长组分分辨率高。对于三峰谱,当SNR≥20时,截断法和阻尼法反演的T2谱光滑连续,具有三峰结构;当SNR<20时两者均不能反演出与构造谱接近的三峰谱。因此,2种方法可应用于SNR≥10的双峰谱和SNR≥20的三峰谱反演;从分辨率和稳定性、平滑噪声及截断和阻尼因子的选取上,阻尼法略优于截断法。 相似文献
9.
构造应力场边界载荷反演的有限元逆逼近法 总被引：2，自引：0，他引：2
构造应力场问题的边界条件、加载方式是地应力场分析的难题,由于构造运动未知和地质构造复杂,无法直接求解。边界力反演是进行构造应力场研究的重要手段。从弹性力学基本方程出发,针对油气储层地应力场反问题,结合对区域构造应力场认识,利用阻尼最小二乘法建立了应力场反演的约束模型,提出了一种用于反演构造应力场边界力的反分析方法,该方法根据部分实测点的应力值进行边界力反演,并与有限元方法和有限差分法相结合,以边界力为参量逼近目标函数。根据研究区少量地应力实测资料采用该方法进行最优化计算,使得计算应力场与实测应力达到最优拟合,以弥补目前的回归反演方法和边界载荷调整法的不足来提高计算精度。算例表明了该方法的有效性,可用于解决构造应力场边界力的反演问题。 相似文献
10.
基于积分方程法的大地电磁三维反演 总被引：1，自引：0，他引：1
采用积分方程法计算均匀半空间中三维异常体的大地电磁响应,然后利用阻尼最小二乘法反演得到地下三维异常体电阻率分布。在计算灵敏度矩阵时,将灵敏度矩阵分为线性项和体现反演过程非线性的偏微分项进行求解,提高了计算精度。理论模型反演结果表明,反演较准确,对初值依赖性小,且只针对异常体剖分单元进行反演,反演效率较高。 相似文献
11.
地球物理反演是指根据一组观测结果来估计假设的地下模型的参数。因为模型所伴隋的响应可能是模型参数的非线性函数,因此在反演中非线性最小二乘方法是有用的。反演的一种普通类型是利用马奎特—列文伯格法实现迭代的阻尼最小二乘。通常,这种方法是通过常规的方法解有关的法方程来实现的。然而对于同样的法方程采用奇异值分解法(SVD)在计算精度上就会有明显的提高。迭代的最小二乘模拟在地球物理的多种问题中得到了应用,用两个实例来说明:(1)地震子波的反褶积;(2)根据地表重力资料推断一个隐伏台阶的位置。再广而言之,非线性最小二乘反演可以用于估计任意一组能得出适当数学描述的地球物理观测结果所对应的地下模型。 相似文献
12.
岩石复电阻率Dias模型及其参数求取方法 总被引：4，自引：0，他引：4
Dias模型是众多描述岩石复电阻率的频散特性模型中的一种,它不但参数少,且可由这些参数导出具有特定物理意义的参数。针对Dias模型参数的求取及反演解的不确定性评价,这里提出了一种模拟退火阻尼最小二乘联合反演方法。该方法首先将参数求解空间离散化,利用模拟退火找到全局最优解区间,然后利用阻尼最小二乘法在解空间进行求解。该方法有效地克服了原来两种方法的反演速度慢、反演结果依赖于初值的缺点,并实现了对反演结果多解性的评价。研究结果表明,对于岩石复电阻率数据的处理解释,Dias模型的使用比传统的Cole模型及其变形更具有潜在的优势。 相似文献
13.
14.
对提高核磁共振测井差分谱识刖油气的方法进行了研究.介绍了核磁共振测井求解差分谱的原理,井用NUMAR公司的P型核磁处理软件DPP时P型核磁的原始数据进行了差分谱处理;根据限制最小平方误差法,反演得到T2谱,并对P型核磁测井原始数据进行了处理,与DPP处理的结果进行了对比,效果比较理想.通过比较分析,限制最小平方误差法对于提高低信噪比的油气显示很有意义,由于MAP-Ⅱ算法在低信噪比条件下反演的差分谱偏离真实较远.给实际应用带来较大的不确定性.文章提出的方法为低信噪比的数据处理问题提出了一种可能的解决手段. 相似文献
15.
有人给出了一种旨在提高物探实测数据“信噪比”的频率域“最佳线性数字滤波器”.而最小平方滤波器则是上述滤波器的参变量λ恒取1时的特例.经验表明,最小平方滤波器用于重、磁数据处理,缺乏灵活性.因此,为适应实际资料中信号和干扰成份在组合上的千差万别,研究最佳线性滤波器参变量λ如何具体取值问题是很有必要的.本文讨论了λ在滤波器中的作用及其数值求解方法:编制了计算λ值的专用程序(为节约篇幅,未予刊出 相似文献
16.
固定点源激电测深(FPS)是一种有效的激发极化测深勘探方法,但在数据解释中,常规的定量解释法只适应较为简单的地电条件,在起伏地形条件下其解释效果会受到较为明显的影响。针对常规的定量解释方法所存在的不足与缺陷,这里从有限单元法正演模拟入手,采用最小二乘法进行反演计算,实现了在起伏地形条件下的对地电模型进行模拟计算,并将圆弧交汇-相对强度计算结果与最小二乘反演计算结果进行对比。结果表明:最小二乘反演方法能适应起伏地形对模型进行反演计算,其反演效果优于圆弧交汇-相对强度法计算结果。 相似文献
17.
18.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
19.
《岩土力学》2019,(11):4486-4496
深部岩体应力量值确定一直是区域地应力场研究中的热点和难点问题,常用的深部岩体应力量值预测方法多以拟合经验公式外推确定为主,缺乏理论依据且可靠性低。研究工作基于应力多边形理论,结合侧压力系数K与震源机制解应力形因子R,在大范围应力预测结果基础之上,对深部应力量值进行了精细限定。以河北易县紫荆关地区的原地应力测量为参考,结合震源机制解资料,估算获得6、11、19km深度处最大侧压力系数K_(max)分别为1.07±0.07、1.14±0.14、1.09±0.09,最小侧压力系数K_(min)分别为0.85±0.15、0.88±0.12、0.86±0.14,并以此计算最大、最小水平主应力,结果的相对偏差在6%～17%之间;利用震源机制解和水压致裂法联合确定的研究区最大水平主应力方向为N44.4°E。研究提出并建立的联合震源机制参数R和应力多边形限定深部应力状态为深部岩体应力状态估算提供了新的途径。 相似文献
20.
地球物理反问题存在多解性和不稳定性,正则化是得到稳定解的重要手段之一,正则化方法中稳定因子的设计与正则化因子的选择是两项重要研究内容。稳定因子的主要功能是对模型的解空间进行限制,以减少多解性,求得稳定解;为更好的定位地下目标体,研究工作采用了最小支持稳定因子,并给出了统一表达形式。正则化因子在模型误差与数据误差之间起到平衡作用,研究工作在＂L曲线法＂自动选取正则化因子的基础上,提出一种修正的计算正则化因子算法。以2.5维直流电阻率为例,运用非线性共轭梯度法,进行了试算研究,结果表明了修正的＂L曲线法＂算法快速、稳定,最小支持稳定因子具有较好的聚焦特性,有利于实现陡变边界反演。 相似文献