1. 引言
高难山地的地震勘探野外采集施工面临着前所未有的巨大挑战。伴随地震勘探项目规模的持续扩张以及采集参数的不断优化,在装备、人力和物力方面的投入呈几何倍数增长。相应地,地震勘探施工过程中的安全风险也急剧攀升。采取何种科学有效的措施来降低施工风险,消除生产过程中的安全隐患,已然成为地震勘探领域的关键课题[1]-[4]。值得关注的是,新兴技术在地震勘探中的应用,有望为降低施工风险发挥重要作用。
2. 高难山地激发点优选技术
2.1. 研究背景
在山地地震勘探采集施工领域,钻井工序始终占据着至关重要的地位。在这一环节中,常常需要数名施工人员合力抬动重达百斤的钻机,在复杂的地形中艰难前行并开展打井作业。整个施工过程中,设备种类繁多且数量庞大,使得施工负担极为沉重。尤其当施工场景切换至高难山地时,设备与施工人员所面临的风险呈指数增长。出于对“以人为本”安全理念的坚守,在施工设计阶段,常常不得不采用舍弃高难山体区域激发点的方式。但这种处理方式却会引发大片激发点空白区的出现,影响地震勘探成像效果,严重降低了勘探质量,进而导致质量与安全之间陷入难以协调的矛盾困境[5] [6]。在一定程度上,这无疑阻碍了高难山地地震勘探工作的顺利推进。
2.2. 研究目的
在高难山地地震勘探中,施工面临着重重挑战,而精准确定激发点位置以最大程度地规避施工风险,无疑成为达成项目成功的关键突破点。行之有效的激发点优选策略能够大幅降低项目在高难山地施工中的不确定性及风险,有力保障项目的顺利开展,为项目的高品质发展与效益提升奠定坚实基础。历经多年的不懈探索与实践积累,形成了高难山体区激发点优选技术。该技术在有效保障设备和人员安全的基础上,能够顺利完成施工任务,为高难山地地震勘探事业提供了强有力的技术支撑。此技术的目的在于提升高难山体区激发点的准确性和合理性,显著降低施工风险,确保山地地震勘探项目在地形环境极其恶劣的情况下能够高效、安全地实施。
2.3. 技术原理
地震勘探从部署到采集施工过程中,首先依据地质目标进行技术设计,均匀布置物理点,即理论点。随后,在施工设计中遵循满足目的层覆盖次数等原则,结合地面障碍、河流以及禁止施工区域等因素重新布设物理点,即设计点。该激发点优选技术主要应用于施工设计阶段,在地形极为陡峭、施工风险极大的山体区域,基于安全考虑,若按传统方法会出现激发点全部舍弃或难以完成激发点部署的情况。首先计算施工区域起伏度和坡度[7]-[10],按照表1对地形难度进行分级。
Table 1. Terrain risk classification table
表1. 地形风险分级表
起伏度(米) |
坡度(度) |
风险分级 |
P6 > 20 |
S6 > 50 |
极高风险 |
20 ≥ P5 > 15 |
50 ≥ S5 > 40 |
高风险 |
15 ≥ P4 > 10 |
40 ≥ S4 > 30 |
中高风险 |
10 ≥ P3 > 5 |
30 ≥ S3 > 20 |
中风险 |
5 ≥ P2 > 2 |
20 ≥ S2 > 10 |
中低风险 |
2 ≥ P1 ≥ 0 |
10 ≥ S1 ≥ 0 |
低风险 |
起伏度计算方法:本文选取三阶不带权差分法计算地形坡度[11] [12]。若定义3 × 3的DEM网格编号(见图1),则有计算公式:
(1)
式(1)中,S为坡度,d为DEM的单元尺寸,d分别改为对应方向的dx和dy;Zi (i = 1,2,3…,9)为中心点Z5周围各网格点的高程。
Figure1. DEM 3 × 3 partial moving window
图1. DEM 3 × 3局部移动窗口
坡度计算方法:本文选取标准差法来衡量地形的起伏程度[13]-[16],其计算公式为:
(2)
式(2)中,P是起伏度,以矩阵格式进行采样,采样间隔大于DEM数据原始采样间隔,Z1、Z2、…、Zn是计算起伏度矩阵各网格高程,m是n个采样点高程算数平均值。
然后在理论激发点部署基础上,结合地形风险分级,根据测量施工标准,三维地震激发点沿激发线方向最大< ±1/2激发点距,垂直激发线方向最大< ±1/2激发线距,对激发点位置进行优选。尽量将施工设计激发点部署在中、低风险区内。
最后,运用相关软件计算偏移后的覆盖次数,得到覆盖次数图,确保符合设计要求。覆盖次数计算公式为:
N = M × R/(2 × S) (3)
式(3)中,N是覆盖次数,M是接收道数,R为道间距,S为激发点距。
这样不仅能解决地形特别险峻的山体区激发点部署问题,同时能够提高采集施工设计和实测炮点的符合率。
2.4. 技术实施方法
Figure 2. Comparison chart before and after optimization design
图2. 优化设计前后对比图
Figure 3. Comparison chart of fold numbers before and after the optimized design of shooting points
图3. 激发点优化设计前后覆盖次数对比图
极高难山体区,沿炮排方向和连续的相对低风险区,对激发点进行优化部署。某工区设计激发点距100米,激发线间距300米。如图2所示,理论点25385990.5位于地形风险极大的区域,在它旁边的偏移允许范围内有一条连续的中、低风险区,可保障人员和设备安全。同时,此区域沿垂直激发线方向并距离理论点105米,小于激发线间距的一半,符合选点要求。该线路还经过相邻理论点25365990.5,能够有效减少人员在施工中的折返里程。因此将理论点25385990.5移到25385988.5点位处,作为设计施工点位,能够降低作业点风险等级,缩短施工时长,减少人员劳累度,从而保障人员及设备安全。
按照上诉方法对高难山体区的激发点进行偏移,完成施工设计,计算偏移后覆盖次数。如图3所示,是某三维工区的覆盖次数图(横纵坐标为工区经纬度位置,出于保密原因隐去),根据色标和覆盖次数值可以看出优化设计后最低覆盖次数612次,为部署满覆盖672次的91.07%,完全达到设计要求。能在保障安全的前提下,完成地质任务。
3. 技术应用效果
3.1. 工区概况
某三维勘探项目位于库车坳陷北部构造带,南部进入阳霞凹陷,施工面积达971.02 km2。如图4所示,工区北部山体高大、沟壑纵横。属于典型的山高难山地,海拔高程变化大,复杂山体区切割剧烈。给地震勘探物理点部署和采集施工带来了巨大的挑战。
Figure 4. General map of the work area
图4. 工区概况图
3.2. 应用效果
Table 2. Table of changes in shooting point risk classification before and after optimization design
表2. 优化设计前后激发点风险分级变化表
地形风险 |
理论设计 |
优化设计 |
变化率(%) |
极高风险 |
3387 |
383 |
−88.69 |
高风险 |
4910 |
2117 |
−56.88 |
中高风险 |
5197 |
4549 |
12.47 |
中风险 |
4744 |
2211 |
53.39 |
中低风险 |
4391 |
2564 |
41.61 |
低风险 |
4742 |
4218 |
11.05 |
本次三维项目高难山地地震勘探激发点优选技术进行施工设计,在高难山体区优选激发点27,371个,大幅降低了高风险区激发点,变化率见表2。施工班组尽量沿着风险相对较低的区域施工(如图5),保障了200多名施工人员的人身安全。
Figure 5. Partial map of the optimized shooting points
图5. 优选激发点局部图
Figure 6. Distribution map of shooting points in high-difficulty mountainous areas
图6. 高难山体区激发点分布图
Figure 7. Field processing profile
图7. 现场处理剖面图
该项目最终完钻符合率高达92.18% (如图6),相较以往同区域项目提升了12%。与此同时,从现场处理剖面图能够清晰看出,构造顶部成像效果良好,如图7所示。由此可见,该激发点优选技术不但在极大程度上降低了高难山地区施工人员及设备的安全风险,减少了成本投入,而且还能够有力保障地震采集项目质量,具备极大的推广应用价值。
4. 结论
本文以某高难山地施工设计激发点部署为例,基于长期的山地地震勘探生产实践,对高难山地进行了风险分级,总结提出了针对高难山地的激发点优选技术。通过沿着相对低风险区进行激发点优选,在满足覆盖次数设计要求的同时大幅降低了高难山体区施工作业风险等级。最大限度的解决了多年来西北无人区施工风险巨大,激发点部署难,实测点到位率低的问题。实现了保障企业员工生命健康安全和满足地质任务的双赢。广泛适用于高难山地地震采集施工,具有极大的推广应用价值。
NOTES
*通讯作者。