基于arcgis的地震采集质量分析评价方法.docxVIP

基于arcgis的地震采集质量分析评价方法 0 工区观测系统设计 随着地震勘探的发展，地震资料采集工作逐渐扩展到低信噪比区域，显示出不同局外地质条件复杂的特点，地震勘探数据随着减少（接收通道数量增加、覆盖时间增加、面元变小、多重等）而增加。 ESRI公司开发的ArcGIS框架具有很强的海量数据空间分析处理和建模功能,利用它可对采集现场炮点和检波点的分布情况进行建模,完成工区观测系统的设计。将经过坐标配准后的卫星遥感影像及数字高程模型(DEM)数据加载到ArcGIS中,可以为地震勘探施工设计和野外采集准确地提供地貌特征、地表岩性、人文、交通、水文、构造特征等直观性信息 现有的地震采集质量分析评价方法大多是通过计算炮点的属性,然后对它们进行简单的数值分析,如KLSeis地震采集软件是将炮点属性的分析结果以柱状图的形式显示,GeoMountain软件系统也只是将炮点属性的分析结果归类为一级、二级和废品三个等级,对炮点的具体属性值在观测系统中不进行投影分析。这里通过对工区内各炮的能量、分辨率和信噪比属性在观测系统中的对应位置进行散点渲染分析和插值网格成图分析,并结合高分辨率的卫星影像和构造图等各种地质图件,可以对采集质量进行综合分析评价。因此,该技术将会在地震勘探野外采集质量定性科学化评价中,发挥比较重要的作用。 1 基于apefity文件的观测系统设计 基于ArcGIS的地震采集质量分析评价原理是:①根据地震采集的SPS文件中炮点和检波点的大地坐标和高程信息生成炮点和检波点的shapefile文件,shapefile文件是ArcGIS通用的一种数据格式,将两个shapefile文件和卫星遥感影像及数字高程模型(DEM)数据加载到ArcGIS中叠合显示,利用它们叠加的直观信息,可在室内精心做好各种复杂地表如高陡山区、丘陵、水网及城镇等地区的观测系统设计,图1是对四川某工区观测系统的设计;②加载工区内所有地震数据,利用每一炮的均方根振幅值计算该炮的能量,使用李庆忠 2 射击点属性分析的结果是投影的 2.1 炮点分辨率散点成图 利用ArcGIS的分级渲染功能对炮点shapefile中的信噪比、能量和分辨率属性值进行渲染,用某一种颜色表示某一范围的属性值,就可以得到该工区内所有炮点三个属性的散点成图。如图2为大庆某工区部分炮点能量的散点成图。从图2 中可以看出,呈现黄色、橙色或红色的炮点能量值比较大,炮点采集质量较好。 2.2 基于克氏原螯虾的d东北部插值算法 利用ArcGIS的插值功能对炮点shapefile中的信噪比、能量和分辨率属性值进行插值计算,得到该工区内所有炮点三个属性的网格图。所用到的插值方式有反向距离权重插值(IDW) 式(1)中: 式(2)中,d 克里金插值方法最早是由南非金矿工程师克里金(D.G.Krige)于20世纪50年代提出来的,根据不同的条件又分为普通克里金插值、泛克里金插值、协同克里金插值等。 这里采用的是普通克里金法插值算法 式(4)中:Cov(x 图3是利用反向距离权重插值法和克里金插值法对大庆某工区部分炮点的能量进行插值生成的网格图。图4是利用两种插值方式生成的网格图和上面散点渲染图的叠合显示。因为散点渲染图反映的是真实值,从图4以看出,利用反向距离权重插值法生成的网格图与散点渲染图的吻合度更高,更符合实际情况,图4(a)中呈现黄色、橙色或红色的区域的炮点能量值较高,采集质量较好。 3 射击点属性分析的结果是综合投影 为了分析结果更能综合地反映炮点的采集质量,将炮点的能量、信噪比和分辨率基于一定的规则进行融合,具体分两种融合模式。 3.1 基于校正系数的融合模式 某一炮点能量、信噪比和分辨率的融合值QC按式(5)求得: 其中:SN是该炮的信噪比; 3.2 模式的融合值qc 我们可根据工区和采集条件的实际情况,更改炮点能量、信噪比和分辨率所占的权重比例,该融合模式的融合值QC可按式(6)得: 由于炮点的信噪比、能量及分辨率的值的范围不在同一数量级上,在同一权重的情况下,容易造成只凸显大值的效果。为了消除这一影响,这里首先将炮点的能量、信噪比和分辨率的值按公式(7)做了归一化,将它们的值归一化到[0,255]区域: 其中:V 3.3 kse,kps 图5为大庆某工区部分炮基于校准系数属性融合模式和基于变权重属性融合模式的散点渲染成图,其中,k=0.4,Ksn =1/3,Kse =1/2,Kres =1/6。图6为该工区两种属性融合模式的散点渲染成图和插值网格化成图的叠合显示。从图6(a)中可以看出,在三个属性权重相等的情况下,某个属性凸显大值的情况比较明显,在图6(b)中经过归一化之后的属性融合比较平滑,没有凸显大值的情况,比较符合实际情况,其中呈现黄色、橙色或红色的区域的炮点采集质量较高。 4 地震采集构