[发明专利]基于岩石露头雷达探测的地震正演模拟方法

说明书

技术领域

本发明属于地球物理勘探领域，具体地，涉及一种基于岩石露头雷达探测的地层地震反射正演模拟方法，利用雷达探测岩石露头的内部结构及其地球物理反射信息，模拟地层的地震反射结果。

背景技术

地震反射正演模拟揭示不同地质条件下地层的地震反射特征，为开展地震资料处理和地质解释提供理论基础和指导，是复杂地质条件下地震资料处理和地震地质解释的重要依据。

地震反射正演模拟是伴随着地震资料处理和解释而产生和发展起来的，经过了几十年的发展，逐渐形成了地震数值模拟正演和物理模拟正演两大类不同的方法。前者建立地质模型，根据地震波在地下地层中的传播方程，采用不同的数学算法（射线追踪算法或者波动方程算法等），模拟地震波在给定地质模型条件下的反射特征；后者则是根据地下地层规模，建立尺度缩小的等效物理模型，在实验室模拟野外的地震资料采集过程，在物理模型中设置检波器，采集震源激发后的反射波信息，达到地震正演模拟的效果。

上述两类地震正演模拟方法不断发展，前者的发展主要在两方面，一是地质模型不断复杂化，由单一地质现象的模拟到多种地质特征组合的复杂地质模型；二是模拟算法的发展，从基于射线理论的算法模型到基于波动方程的算法模型，通过模拟数学方法的改进，不断向地震波的真实传播特征逼近；后者的发展主要是表现在物理模型不断复杂化，越来越向地质学家认识的地下复杂地层特征逼近，同时，物理实验采用的检波器、激发震源等硬件也不断精密化，采集的地震反射信息更加精确和全面。

对于基于数值模拟的地震正演模拟来说，目前最大的问题有两点，一是所建立的地质模型虽然不断向复杂化发展，但是距离真实地层的复杂性还有很大差距，二是正演结果受到模拟算法的影响，而目前已经使用的模拟算法都只是波在地层中传播的一种近似。

探地雷达是通过高频电磁波利用有耗介质中介电常数差异探测地下或其他物体内部不可见目标或界面分布的一种浅层地球物理方法。1910年，Leimbach等率先将雷达应用于确定地下埋藏物的位置，20世纪70年代开始迅速发展，被应用于地下浅层地层研究、冰层厚度探测、水文地质调查、冰川研究、考古、排雷、矿产资源勘查、工程地质勘查、无损检测及工程建筑物结构调查等领域。

20世纪60年代，探地雷达开始被应用于地质领域，随着国内外学者对野外岩石露头和现代沉积原型模型研究的重视程度增加，上世纪90年代末起，探地雷达开始被应用于野外岩石露头、浅层地层和现代沉积研究，在这些研究中形成了探地雷达与其他技术相结合的研究方法，如，JensTronicke等（2002）在重力流沉积岩石露头中将GPR、井间雷达层析成像与自然伽玛测井资料相结合，KeumsukLee等（2007）将GPR和激光扫描两种露头研究技术相结合在浊积水道和三角洲前缘岩石露头地质描述研究中取得了良好的效果。

虽然探地雷达和人工地震反射资料是对地层不同特征参数的反映，但是二者都是地层界面反射信息，探测原理具有相似性和可类比性，因此，在探地雷达资料的地质解释方法上，地震地层学的解释方法和概念体系被众多学者接受，AdrianNeal（2004）在研究中论证了二者在岩石露头探测中的相似性，并指出地震地层学的解释方法在探地雷达资料解释中应用是可行的，但是受探地雷达研究尺度的制约，对大尺度的概念（层序和层序边界）不建议应用于探地雷达资料解释中。

已开展的探地雷达在地质研究中的应用主要是利用其进行岩石露头或现代沉积的直接探测，得到地层或沉积体的内部结构，从而对无法直接观察的地质体进行地质描述，用于建立各种地层或沉积体的地质模型。

发明内容

为克服现有技术中地震正演模拟方法存在的缺陷，本发明提供一种基于岩石露头雷达探测的地震正演模拟方法，为复杂地层的地震反射正演模拟提供一种快捷而更加准确的方法，为地震资料处理和地质解释提供依据。

为实现上述目的，本发明采用下述方案：

基于岩石露头雷达探测的地层地震反射正演模拟方法，包括如下步骤：

(1)、选择岩石露头剖面，设置雷达探测测线；

(2)、进行岩石露头剖面测量和地质描绘；

(3)、分别选用100MHz和400MHz雷达天线沿测线对岩石露头剖面进行连续探测，得到相应频率天线探测的雷达剖面数据；

(4)、对步骤(3)采集的雷达数据进行数据成像处理；