含地表上长波长静校正量计算方法研究.docxVIP

含地表上长波长静校正量计算方法研究 hd4块位于塔里木盆地。该区域的表面条件复杂，覆盖着不规则的蜂窝沙丘和浮土。地震数据的高频信号吸收衰减相对严重。同时，环境相对较差，导致附近表面的干扰和散射能量相对发育。因此，原始数据的质量差、信噪比低。由于表面波的缓慢和起伏的变化，面波发育，频率范围宽。原始单次激发记录的上一破裂曲线扭曲特征明显，静态校正问题相对突出。由于目标层埋深（约3.5s），结构的横向变化和振幅小于（小于6ms），因此静校正是该区域低信噪比地震数据处理的重要因素。 1 静校正量的计算 实际资料观测显示, 该沙漠地区的潜水面 (高速层顶界面) 分布比较稳定.在炮点井深设计时可以获得该潜水面的深度或高程信息, 而在其上覆盖的沙丘和浮土可视为低速带 (风化层) 区域, 疏松地表速度通常为350 m/s, 同时地表高程可有GPS等手段精确测量, 再根据替换速度和固定基准面高程等信息可以用于计算沙丘曲线静校正量.现已知的测量资料 (图1) 有:炮点井深为 由于震源位置通常选择在潜水面以下某一深度处, 使得地震波能量能够更大程度地向地下传播, 因此, 炮点静校正量的计算可以不考虑潜水面位置和近地表低速带厚度等因素的影响.采用公式 ( =[ 与炮点静校正量的计算过程不同, 检波点静校正量则需要考虑近地表和潜水面位置的影响, 因此需要先将检波点校正到潜水面, 再将其校正到固定基准面.在这个计算过程中, 不能忽略风化层 (如厚度、压实系数等因素) 对静校正量计算的影响, 该影响表现为速度随深度的有规律变化.因此, 采用公式 ( = ( 这里压实系数求解的主要采用最小二乘曲线拟合方法: 构建目标函数: =∑[ ( 其中, 对于给定的样点数 =∑[ ( 根据塔里木河南区块大量的试验和总结, 考虑到其与HD4区块相邻, 因此选择压实系数为 在垂直方向上, 由于采用炮点和检波点高程等信息计算出的单程时间作为静校正量, 没有完全考虑近地表风化层速度的横向变化, 因此本文方法只能解决长波长静校正问题, 短波长静校正量则需要通过剩余静校正方法来解决. 2 应用效果分析 图2为本工区的炮点和检波点以及潜水面的高程图, 这些量均可通过测量得到.从图2 (a) 和图2 (b) 可知, 沙漠地表存在的沙丘数量众多, 其最大起伏高差为40 m, 而潜水面则较为平缓 (图2 (c) ) , 且始终位于地表以下某一深度处, 无高速层出露地表, 这为沙丘曲线静校正量的计算提供了有利条件.综合炮点井深、炮检点高程、替换速度、固定基准面和压实系数等信息, 经式 (1) 和式 (2) 计算, 最终得到了全工区的沙丘曲线静校正量, 并与野外测量的静校正量进行对比 (表1和图3) .从表1 (随机选择了工区的炮点和检波点桩号及其对应的静校正量) 可以看出, 由沙丘曲线公式计算得到的炮检点静校正量与野外测量得到的静校正量数值除个别点差异较大外, 基本符合的较好, 其中前者计算的精度较高.从图3可以看出, 野外测量的炮检点静校正量直方图表现为稀疏特征, 其原因在于未测量点的数值通常是通过邻点线性插值得到, 同时舍弃了小数位从而丢失了精度, 而沙丘曲线炮检点静校正量表现在接近正态分布特征, 保留的数值精度较高. 分析图4可知, 由于沙丘地表的起伏变化, 造成了原始单炮同相轴扭曲和不连续变化特征, 经过野外测量静校正和沙丘曲线静校正后及其对应的一次剩余静校正后 (剩余静校正所采用的时窗和参数均相同) , 这种扭曲特征得到了有效压制, 同相轴变得更加连续, 其中, 经沙丘曲线静校正后同相轴的连续性更好.从叠加剖面 (图5) 分析可知, 在未应用静校正时, 剖面信噪比较低, 且连续的同相轴相对较粗 (非同相叠加所致) , 而应用野外测量静校正和沙丘曲线静校正量及其对应的一次剩余静校正量后, 信噪比得到了有效改善, 同相轴连续性增强, 同时也可以看出, 经沙丘曲线静校正后叠加剖面的效果与经野外测量静校正后的效果相当或更好些. 3 静校正计算精度与野外测量静校正量关系 1) 针对沙漠地表的特点, 并根据已有的测量资料, 可以构建固定基准面在地表之上的沙丘曲线静校正和最小二乘曲线拟合求解压实系数的计算公式. 2) 与野外测量静校正量进行对比可知, 沙丘曲线静校正的计算精度更高;同时在单炮和叠加剖面方面的应用效果相当或更好. 3) 该方法的关键是测量的物理量需准确, 同时压实系数需要通过统计并利用最小二乘法拟合得到, 它的优点在于计算过程简单、快捷, 不需要过多的人力资源等投入. 4) 在潜水面或高速层顶界面出露地表情况, 文中计算公式并不适用, 此时可以通过进一步合理修改式 (2) 中的对数项, 从而避免对数项为0或负值的问题.