华北东部地壳三维模型建模方法.DOCVIP

华北东部地区三维地壳模型建模方法 段永红 中国地震局地球物理勘探中心 华北克拉通地区，特别是首都圈地区是我国DSS剖面探测工作时间最长, 测线最密集的区域。从1976 年的乐亭-张家口剖面到2011年完成的盐城-包头剖面，不同时期施工的测线共有42条,在首都圈地区测线间隔大约为30-50km，外围测线密度逐步降低。20 世纪80 年代使用的地震仪多是模拟地震仪,20世纪90 年后, 随着电子技术和计算机的发展，数字地震仪逐渐取代了模拟地震仪，原来的模拟记录也得到数值化，数字化地震记录的使用大大提高了资料处理解释的效率和精度。随着地震仪器的增多，DSS剖面的观测系统越来越完善, 使我们对地壳的解释精度也大大提高。 本研究收集了42条剖面的原始资料和处理结果，在对这些研究成果进行综合分析的基础上，提出建模准则，构建了网格为0.25°×0.25°×10 km（10 km以上为0.25°×0.25°×1 km）的三维地壳结构模型（HBCrust.1），模型范围为经度112°-120°，纬度36°-42°的华北克拉通中、东部地区。 1坐标转换 传统上，深地震测深二维探测剖面的成果采用的坐标系统是桩号系统（在测线定义一个零桩号，其余测点位置均以距零桩号的距离来定义），因此需要将所有剖面的桩号转换为经纬度地理坐标。转换的办法是找出数据采集时的原始的记录班报，根据原始记录炮点和测点的经纬度坐标，将剖面上的桩号转换成经纬度坐标。 2.剖面资料网格化 DSS剖面的速度结构有两种资料，一种是保存有Seis83射线追踪得到的数字化原始数据，这样的资料是主要部分。对有数字化资料的剖面，首先沿剖面对界面（G、C和M）和速度按横向5.0 km，垂向10km以上1.0km、10km以下2.0km的间隔进行网格化。另外由于探测剖面完成时间跨度太大，资料处理人员的更迭，少量的剖面研究成果仅保留图形化的成果。对少量只有图像资料的剖面，首先采用图形数字化将图形资料进行数字化，然后再按照同样的网格间距进行网格化。 3.剖面交叉点的处理 由于获取宽角反射/折射剖面探测的资料的时间跨越近四十年，资料质量、资料处理的人员和方法均不同，因此在剖面相交的区域,同一点在不同剖面上的速度值或界面深度值会有一些差别,在建立三维模型时就必须进行适当的处理来保证模型的连续性。本文的处理的原则是：根据剖面的质量进行加权平均，即观测系统较为完善,观测资料良好的剖面给予最大权重为1，上世纪90年代以后的剖面基本属于这一类；对观测系统基本满足需要，解释合理的剖面给予0.8的权重；对一些观测系统不完善，解释不太合理的剖面的权重为0，最后取加权平均值作为最后的结果。 4. 三维网格化插值 在得到各条剖面的网格化参数值后，下一步就是要将剖面上网格化参数值进行三维网格化。水平网格间距的大小反映了所建地壳模型的分辨尺度，所以网格间距的大小要与测线的密度相匹配，一般情况下，要保证每个网格内都要有采样点。从图可以看到，在首都圈地区测线较为密集，测线间距大约30-50 km，而在外围测线密度变小，因此本文水平网格间距选取0.25?。这样的网格间距在研究区的大部分区域都能保证地壳模型具有与测线密度相匹配的分辨，在模型边缘特别是四个角分辨稍差一些。 克里金数据网格化技术利用随机函数理论，把空间各点的采样值看成随机变量的现实，且把网格节点上数值的估计归结为随机函数的最佳无偏估值问题，因而能很好地适应采样数据与估计结果的随机性（魏文博等，2007），所以本文采用普通克里金插值方法，首先分层插值，最后拼接成地壳模型，这样既保证了插值精度，同时也保证了各层之间速度的独立性。 5. 质量检测 在得到三维地壳模型后，攫取不同方向的剖面进行检查，对一些奇异点进行检查校正，最后得到研究区三维模型。