你主要是想用谱模拟的方法得到地震子波吧。它是根据地震道的振幅谱,模拟出子波的振幅谱,主要是通过地震记录x(t)的自相关函数和地震子波 b(t)的自相关函数相同来求取的。 其中假设反射系数序列是白化的,并且与噪声不相关的情况下,地震数据的振幅谱是地震子波的振幅谱倍数。这个方法只能得到振幅谱,而没法的到相位谱,且实际地震子波是混合相位的,所以只能求出最小相位子波(这个方法的局限是要求地震子波是最小相位的),反正不是很符合实际。
可参考一下这篇文章 ,《一种新的混合相位反褶积方法》文章作者:伊振林 王润秋 中国石油大学资源与信息学院CNPC物探重点实验室
这篇文章中介绍的可以用Z变换得到混合相位,由于对最小相位子波做Z变换,那么他们的根都在单位圆外Z
i, 取几个根如Z
1Z
2Z
3Z
4,,把他们变成1/Z
i,再反变换回去,这样他们的根既有园外的又有园内的,所以就成混合相位了,至于合适不合适只能用其他判别准则来区分了。文章里用P范数与熵准则来联合判定的。
有了合适的子波,反褶积应该就不成问题了吧。
[ 此帖被glr2008在2009-11-26 14:10重新编辑 ]