原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
高振山等：基于二阶统计广义最小二乘的地震盲反褶积方法研究
其中 为子波褶积矩阵， 为反射序列 2
Cw q C q −s diag (q).*q
J w 2 =+λ ln(1+ 2 （5 ）
σ δ
向量，s 为地震记录向量。
假设反射系数先验信息服从正态分布，根 2 2
其中 σ 2σ ， 为折中参数，
λ
据Bayesian 公式，则反射系数的后验概率分布 n
为 2 2
δ 2σ 。规则化得到
q
1 T 1 ( )T ( )
nq / 2 −q q n / 2 − C q −s C q −s
P (q s) ∝( 2 ) exp( 2 )( 2 ) s exp( w 2 w )
2πσ 2σ 2πσ 2σ
q q n n
(CT C +Q)q CT s （6 ）
（2 ） w w w
其中σ 、σ 为反射系数和噪声的方差； 其中Q 为约束对角矩阵。观察到（6 ）式
q n
n 、n 为反射系数和地震记录的个数。使 （2 ） 为非线性隐格式，应采用迭代方法求解。在每
q s
式最大时得到反射系数称为最大后验概率 步迭代中，线性系统（6 ）可用矩阵求逆、矩阵
（MAP ）Bayesian 反演估计。两边取自然对数， 分解或共轭梯度法(CG) 方法求解，称为直接
略去常数项，重新定义Bayesian 反演的目标函 法。
数，并进一步改进广义最小二乘为： 直接对（6 ）式求解，但相关矩阵的条件数
仍可能很大，进行预优处理将改善计算性能，
εk qk
J ∑ρ1( ) =+∑ρ3 ( ) （3） 提高反演精度。
k σ k σ
n q 稀疏约束反演问题（6 ）对应的原始问题可
1 2 写作：
ρ (u) u
1 2
文档评论（0）