Kirchhoff偏移法在探地雷达正演图像处理中的应用
|
摘要
本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法(FDTD)的原理,并用它合成了几种常见的雷达正演剖面;最后利用Kirchhoff积分偏移法对正演所得的雷达剖面进行偏移处理,通过对比偏移处理前后的雷达正演剖面,可知Kirchhoff积分偏移法能使雷达正演剖面中的反射波的归位,绕射波收敛,从而大大地提高了雷达正演剖面的分辨率,更好地指导GPR剖面的地质解释和验证偏移方法的有效性.
Firstly,this article analysis the wave equation,and gets the format consistency of the sound wave equation and the radar wave equation.According to this,it illustrates the feasibility of inducting the migration technique commonly used in the seismic data processing into the GPR data processing.Secondly,this article explains the principles of FDTD,and synthesizes some radar forward simulate sections with it.At last,this article utilizes Kirchhoff 's method of integration migration to do the data processing to the radar section coming from the forward simulate.By comparing the radar forward simulate sections before and after the migration process,it is showed that Kirchhoff 's method of integration migration could make the reflection wave return to original position,the diffraction wave convergent.So that,the resolution of radar forward simulate sections could be enhanced highly,and it could instruct the geology explanation of GPR section and validate the validity of migration method better.
Firstly,this article analysis the wave equation,and gets the format consistency of the sound wave equation and the radar wave equation.According to this,it illustrates the feasibility of inducting the migration technique commonly used in the seismic data processing into the GPR data processing.Secondly,this article explains the principles of FDTD,and synthesizes some radar forward simulate sections with it.At last,this article utilizes Kirchhoff 's method of integration migration to do the data processing to the radar section coming from the forward simulate.By comparing the radar forward simulate sections before and after the migration process,it is showed that Kirchhoff 's method of integration migration could make the reflection wave return to original position,the diffraction wave convergent.So that,the resolution of radar forward simulate sections could be enhanced highly,and it could instruct the geology explanation of GPR section and validate the validity of migration method better.
引文
[1]冯德山,戴前伟.用Visual C++开发地质雷达正演模拟软件[J].物探化探计算技术,2004,26(3):231~234.
[2]冯德山.地质雷达二维时域有限差分正演[D].长沙:中南大学,2003.
[3]李成方.偏移处理在探地雷达数据处理中的研究[D].成都:成都理工大学,2002.
[4]薛桂霞,邓世坤,刘秀娟.吸收边界条件在探地雷达偏移处理中的应用[J].物探化探计算技术,2004,26(3):235~237.
[5]陈树文,刘洪,李幼铭.三维叠前深度偏移的准三维算法研究[J].地球物理学进展,2001,16(4):23~28.
[6]耿建华,黄海贵,马在田.Kirchhoff积分波场延拓基准面静校正方法研究[J].地球物理学进展,1996,24(6):665~669.
[7]刘喜武,刘洪.波动方程地震偏移成像方法的现状与进展[J].地球物理学进展,2002,17(4):582~591.
[8]黄德济,贺振华,包吉山.地震勘探资料数字处理[M].北京:地质出版社,1990.
[9]李世华,杨有发.物探数据处理[M].北京:地质出版社,1995.
[10]陆基孟.地震勘探原理(第三版)[M].北京:石油大学出版社,2001.
[11]李大心.探地雷达方法与应用[M].北京:地质出版社,1994.
[12]赵永辉,吴健生,万明浩,谭春.有限差分法探地雷达波动方程偏移成像[J].物探化探计算技术,2001,23(1):47~51.
[13]邓世坤.克希霍夫积分偏移法在探地雷达图像处理中的应用[J].地球科学—中国地质大学学报,1993,18(3):303~308.
[14]底青云,王妙月.雷达波有限元仿真模拟[J].地球物理学报,1999,42(6):818~825.
[15]Jeffrey E.Lucius,Michael H.Powers.GPR data processingcomputer software for the PC[R].USGS,2002.
[2]冯德山.地质雷达二维时域有限差分正演[D].长沙:中南大学,2003.
[3]李成方.偏移处理在探地雷达数据处理中的研究[D].成都:成都理工大学,2002.
[4]薛桂霞,邓世坤,刘秀娟.吸收边界条件在探地雷达偏移处理中的应用[J].物探化探计算技术,2004,26(3):235~237.
[5]陈树文,刘洪,李幼铭.三维叠前深度偏移的准三维算法研究[J].地球物理学进展,2001,16(4):23~28.
[6]耿建华,黄海贵,马在田.Kirchhoff积分波场延拓基准面静校正方法研究[J].地球物理学进展,1996,24(6):665~669.
[7]刘喜武,刘洪.波动方程地震偏移成像方法的现状与进展[J].地球物理学进展,2002,17(4):582~591.
[8]黄德济,贺振华,包吉山.地震勘探资料数字处理[M].北京:地质出版社,1990.
[9]李世华,杨有发.物探数据处理[M].北京:地质出版社,1995.
[10]陆基孟.地震勘探原理(第三版)[M].北京:石油大学出版社,2001.
[11]李大心.探地雷达方法与应用[M].北京:地质出版社,1994.
[12]赵永辉,吴健生,万明浩,谭春.有限差分法探地雷达波动方程偏移成像[J].物探化探计算技术,2001,23(1):47~51.
[13]邓世坤.克希霍夫积分偏移法在探地雷达图像处理中的应用[J].地球科学—中国地质大学学报,1993,18(3):303~308.
[14]底青云,王妙月.雷达波有限元仿真模拟[J].地球物理学报,1999,42(6):818~825.
[15]Jeffrey E.Lucius,Michael H.Powers.GPR data processingcomputer software for the PC[R].USGS,2002.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心