backus-gilbert反演理论pps - 武汉大学.ppt

BG反演理论在反褶积中的应用 式中： 是观测数据协方差矩阵 之要素。 必须在方差和分辨率之间取折衷，使目标函数变为： （4.106） 这里 是折衷参数，且 将（4.103）式、（4.105）式代入（4.106）式并利用（4.102）式求解。设： BG反演理论在反褶积中的应用 简化后，得： 式中： （4.107） 而 BG反演理论在反褶积中的应用 同样： BG反演理论在反褶积中的应用 这里 是协方差矩阵。下面讨论（4.107）式的解法： 因矩阵 为对称正定矩阵，可令 。定义转轴后的向量 和 为： （4.108） BG反演理论在反褶积中的应用 且令： （4.109） 在 时， 则（4.107）式的解为： 在求得 后，可按（4.108）式直接计算 ，进而求出平均函数 及 。 BG反演理论在反褶积中的应用 根据褶积定理，时域上的褶积等于频率域上的乘积。如以小写字母表示时间域的函数，大写字母表示其谱，则下式成立： （4.110） 设 代表 的周期， 和 分别代表 和 之傅氏变换，则（4.110）式可表示为由于 是一个带限的，所以 必然也是带限。 BG反演理论在反褶积中的应用 下面，研究如何在频率域上实现反褶积。为了便于理解，还必须从时域讲起。在地震勘探中，常以 表示地震波， 表示反射函数， 表示子波， 表示干扰，因此设： （4.111） 假定： ， ， 和 都是周期为 的周期函数。在连续介质情况下，设 对应于 ， 对应于 ， 对应于 ，我们的目的是求一连续函数 ，以使平均函数 BG反演理论在反褶积中的应用 （4.112） 尽可能接近于 函数。 对（4.111）两端同时进行褶积运算，并利用（4.112）式： 定义： BG反演理论在反褶积中的应用 （4.113） 故 定义： 是 的方差； 是真值 设目标函数： （4.114） 假定噪声是平稳的，则方差项可写成： BG反演理论在反褶积中的应用 （4.115） 根据Parseval定理，则有： 式中， 为周期， ， 分别为噪声和速度谱。 把Parseval定理应用于（4.114）式的第一项，并利用（4.115）式则得： （4.116） 式中， 是子波的频谱 BG反演理论在反褶积中的应用 如以 ， 表示 ， 之共轭，并设： 则得： （4.117） 由于 是带限的，当 时， 。当 时，上式变为： （4.118） BG反演理论在反褶积中的应用 又因为 故 （4.119） 利用（4.118）式，则有： 在频率域上，根据（4.119）式求出 之后，再对之作傅氏反变换即可得 。 Thank You BG线性评价（一） 3、BG展伸准则 讨论： 按（4.68）式计算的 ，可以近似地表示平均函数的分辨宽度。当平均函数 集中在 附近时，其值与平均函数的分辨宽度重合很好。记 ，则 的大小和分辨率成反比。而 越小，则分辨率越高。 由BG展伸准则求出的平均函数的宽度，一般比狄里希来准则求出的平均函数的要宽，但其边叶要小，不存在负值，则相应的 的方差较小。 BG线性评价（一） 3、BG展伸准则 讨论： 基于BG展伸准则工作时，由于每一个 的加权系数都不相同，从而大大增加了计算内积矩阵（4.72）式和解方程（4.79）式的工作量，从而增加了计算成本。 定义 时，应加一系数12。 综上所述，BG线性评价的信息包括在三个量中，即平均值 、平均函数 和分辨率 。只有综合分析以上三种信息才可能取得关于 处地球物理模型的真正信息 BG线性评价（二） 上一节讨论的是在观测数据为有限、精确的情况下，BG的线性评价问题。然而，实际观测数据都具有误差。在有误差的观测数据时，如何进行线性评价？即是本节所要讨论的问题。设： （4.82） 式中： ， 分别为第 个观测数据和方差； 为第 个数据的真值。 下面讨论具有误差的 处的模型的平均值，有： BG线性评价（二） 将 记为： （4.83） 式中： 为平均值之真值； 为 处平均值之