叠前反演参数概率分布函数描述.pdfVIP

·594· 中国地球物理2011 (17)油气田与煤田地球物理勘探 叠前反演参数概率分布函数的描述 Prestackinversion ofthe parameterdescription function distribution probability 田玉昆‘周辉张玲玲张红静 TianYukunZhouHuiZhangLingling ZhangHongjing 北京昌平102249 中国石油大学(北京)CNPc物探重点实验室 1．引言 贝叶斯参数估计理论已经广泛地应用于地球物理反演领域，能够有效地利用地球物理各方面的资料作 为先验信息，加入到反演当中，可以很好地对反演进行约束，提高解的稳定性。贝叶斯估计中涉及到弹性 参数的概率密度分布问题，通常的方法是假设弹性参数符合某种分布，一般情况下常用的是高斯分布。然 后在某种贝叶斯估计的准则下，对参数进行估计。该假设常常存在这种已知分布类型的概率分布函数，能 否准确描述地球物理参数分布的问题。事实上，当采用高斯假定时，数据集中单一的离群点会导致难以接 受的反演结果。通常，采用“长尾型”密度函数来表示这类数据集的不确定性。另外，长尾性表示地层的 细节信息，能够更好地体现薄层信息和保护微弱反射，对于提高反演分辨率来说更有意义。本文采用一种 多参数控制的正态反高斯分布函数描述反演参数的分布，合理的控制长尾巴的程度，从而提高反演的性能。 2．基本原理 从贝叶斯理论可知，模型的后验概率函数可以表示为模型的似然函数与先验信息的乘积，即： p(聊Id)。!旦型掣Ⅸp(dl脚)P(m) (1) p～口J 模型的先验函数，如何更精确的表示p(m)是本文着重研究的问题。这里用正态反高斯分布来描述模型参数的分布： faq(x)] (2) 厂(石)=兰％exp[p(工)]K。 。 一 刀g～X， ．—————一广———————————= 不同的长尾型函数对于反演结果求取的准确程度，也有着差异。例如，Laplace函数在对参数进行建模时， 有可能夸大了参数分布的长尾程度，从而得到了糟糕的反演结果。正态反高斯分布的形状由四个参数表征， 由于四个参数取值的灵活性，使得正态反高斯分布能够对任意长尾程度的参数进行建模。所以采用多参数控 制的正态反高斯分布函数来表示模型的先验函数，能够充分利用参数控制分布的衰减速度口、偏斜因子∥、 尺度参数6和平移参数∥等因素，使得所建立模型更接近实际资料的分布。有助于提高反演分辨率。 3．结论与讨论 正态反高斯概率分布函数能够更准确地描述地球物理参数的先验信息，对比常用的高斯分布，其解的 精确性和稳定性都有很大提高。并且长尾信息表示地层的细节信息，能够保护弱反射和薄层信息，换言之 就是能够获得高分辨率的反演结果。同时，对多种长尾型概率密度函数之间解的对比发现，多参数控制的 正态反高斯分布函数能够更好地表示模型的先验信息，求得的结果也更加接近真实解。 参考文献 ProblemandMethodsforModelParameter forIndustri