CN103605152B - 一种基于障碍物安全区的自动炮点偏移方法 - Google Patents

Abstract

本发明实施例提供一种基于障碍物安全区的自动炮点偏移方法，包括：步骤1，获取多个障碍物类型以及范围；步骤2，根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域；步骤3，叠加合并所述多个障碍物的安全区域，生成障碍物安全区；步骤4，计算落在所述障碍物安全区内的炮点，并按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外。本发明实施例的基于障碍物安全区的炮点偏移方法，可以自动计算炮点偏移位置，为三维地震采集施工选择提供炮点偏移参考。

Description

一种基于障碍物安全区的自动炮点偏移方法

技术领域

本发明涉及地球物理勘探技术，尤其涉及一种基于障碍物安全区的炮点偏移方法。

背景技术

随着地震勘探领域逐渐向人口密集的村镇、交通发达的城区和复杂山地转移，各种地物如居民区、工厂设施、道路、管线、高压线、水库、堤坝等，这些影响地震采集的地物称为障碍物。为保证障碍物的安全无损，炮点布设应与障碍物保持一段距离，这种禁止放炮的区域称作障碍物安全区。炮点可布设区在在障碍物安全区外。观测系统设计后，部分炮点落在障碍物中，或落在障碍物安全区中，不能保障安全激发，需要离开建筑物管线等地物一段距离，或者落在存在工农关系问题区域，不能按照设计的观测系统进行施工。如何设计观测系统，避开障碍物，怎样有效利用障碍物间隙可放炮区域，进行炮点设计，成为采集施工人员关心的问题。

目前在建筑物附近进行地震采集施工中，炮点激发时需要符合有关安全距离的要求。工区中有多个障碍物时，靠野外现场施工保障安全距离难度较大。目前依靠采集设计人员凭经验，或依据高清卫片、航片凭借经验实施炮点偏移，不能精确控制安全距离，存在较大的安全风险。

当前国内外尚无实现基于障碍物安全区自动进行炮点偏移的方法。

发明内容

本发明的目的是实现一种基于障碍物安全区的自动炮点偏移方法，在观测系统设计时，通过障碍物的安全距离设定，叠加计算工区障碍物安全区，自动计算炮点偏移位置，为三维地震采集施工选择提供炮点偏移参考。

为了达到上述目的，本发明实施例提供一种基于障碍物安全区的自动炮点偏移方法，包括：步骤1，获取多个障碍物类型以及范围；步骤2，根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域；步骤3，叠加合并所述多个障碍物的安全区域，生成障碍物安全区；步骤4，计算落在所述障碍物安全区内的炮点，并按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外。在所述步骤4中，所述按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外，包括：将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置；若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直测线方向移动，直到移动到合适的位置。

进一步地，在一实施例中，在所述步骤1中，所述获取多个障碍物类型以及范围，包括：通过野外测量组实地实测影响炮点施工地物作为障碍物；根据工区内的卫星照片或高清航片定义障碍物的类型。

进一步地，在一实施例中，在所述步骤2中，所述根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域，包括：根据所述障碍物的建筑年份、建筑构架以及管线走向，确定所述障碍物的安全区域。

进一步地，在一实施例中，在所述步骤3中，叠加合并所述多个障碍物的安全区域，生成障碍物安全区，包括：将安全区重叠的障碍物合并，融合多边形边界，叠加计算工区障碍物安全区。

进一步地，在一实施例中，在所述步骤4中，将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置，包括：计算过障碍物的沿测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；炮点S_i纵向偏移坐标为：

$\begin{array}{c}{S}_{ix}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)cos\left(\mathrm{\θ}\right)\\ S_{iy}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)cos\left(\mathrm{\θ}\right)\end{array};$ 其中，S_ix、S_iy为偏移后炮点S_i的x、y坐标偏移量；i为1、2、3……m，m为所述过障碍物的沿测线的炮点构成线上落在障碍物中心点一侧的点数；O_d为沿测线构成线离中心位置最近的炮点到障碍物边界的距离；R_d为道距；j为1、2、3……j＝i，若炮点S_i偏移后落在已有的炮点上，j＝i+1；θ为测线的方位角。

进一步地，在一实施例中，在所述步骤4中，若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直测线方向移动，直到移动到合适的位置，包括：计算过障碍物的垂直测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；炮点Sx_i横向偏移坐标为： $\begin{array}{c}{\mathrm{Sx}}_{ix}=\left(\mathrm{int}\right(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}){\mathrm{RL}}_d+{\mathrm{jRL}}_d)\sin (\mathrm{\θ}\±\mathrm{\π}/2)\\ {\mathrm{Sx}}_{iy}=\left(\mathrm{int}\right(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}){\mathrm{RL}}_d+{\mathrm{jRL}}_d)\cos (\mathrm{\θ}\±\mathrm{\π}/2)\end{array};$ 其中，Sx_ix、Sx_iy为偏移后炮点Sx_i的x、y坐标偏移量；i为1、2、3……m，m为过障碍物的垂直测线的炮点构成线上落在障碍物中心点一侧的点数；Ox_d为垂向构成线离中心位置最近的炮点到障碍物边界的距离；RL_d为接收线距；j为1、2、3……j＝i，若炮点Sx_i偏移后落在已有的炮点上，j＝i+1；θ为测线的方位角。

进一步地，在一实施例中，还包括：当井炮施工变动后，所述障碍物的安全区域发生变动，按照所述步骤2、3、4进一步将炮点偏移出所述障碍物安全区外，并计算炮点偏移后的观测系统属性覆盖次数、炮检距分布和方位角分布等，保证面元内属性分布均匀，并进行对比，确认炮点偏移后的效果。

本发明实施例的基于障碍物安全区的自动炮点偏移方法，可以自动计算炮点偏移位置，为三维地震采集施工选择提供炮点偏移参考。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的基于障碍物安全区的自动炮点偏移方法的流程示意图；

图2为本发明一具体实施例中的观测系统上的障碍物分布图；

图3为本发明一具体实施例中的障碍物安全区的计算示意图；

图4为本发明一具体实施例中的炮点偏移结果示意图；

图5为本发明一具体实施例中的炮点偏移后的观测系统覆盖次数图；

图6为本发明一具体实施例中的炮点偏移后的观测系统炮检距范围图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例的基于障碍物安全区的自动炮点偏移方法的流程示意图。如图所示，本实施例的基于障碍物安全区的炮点偏移方法包括：步骤S101，获取多个障碍物类型以及范围；步骤S102，根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域；步骤S103，叠加合并所述多个障碍物的安全区域，生成障碍物安全区；步骤S104，计算落在所述障碍物安全区内的炮点，并按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外。

在所述步骤S104中，所述按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外，包括：将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置；若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直方向移动，直到移动到合适的位置。

在本发明实施例中，进行炮点偏移之前，需要按采集设计施工要求设计理论观测系统，例如设置道距、炮距、接收线距以及炮线距等等。

在本实施例的步骤S101中，可以通过野外测量组实地实测影响炮点施工地物作为障碍物，并且根据工区内的卫星照片或高清航片定义障碍物的类型。例如，房屋、桥梁、水堤、输电通讯线路以及输油输气管线等，都属于障碍物，这些障碍物也属于不同的类型。

在本实施例的步骤S102中，根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域。地震勘探采集施工中，井炮作业需要离开障碍物一段安全距离，对施工安全至关重要。实际施工中多种障碍物的形状大小不一，不能完全按照单个障碍物的形状范围进行设计偏移，需要引入障碍物安全区控制参数，建立起多个障碍物之间的安全区域，炮点布设在工区可布设区中，自动进行炮点的选择。

不同的障碍物，如房屋、桥梁、水堤、输电通讯线路和输油输气管线等设施，设定不同的安全距离，以保障布设炮点后，炮点爆炸后，建筑物不至于被地震波冲击损伤。如根据房屋类型，某一井炮施工确定药量的炮点施工，土坯房安全距离50米，砖房安全距离40米，水泥混凝土安全距离25米，管线安全距离75米。不同影响炮点施工的建筑物地物影响施工的距离是不一样的。

在物探小队施工前，根据工区地物类型的不同，施工前甲方有明确规定不同井炮药量，相应房屋、管线等地物的安全施工距离，不至于引起纠纷。

在本实施例的步骤S103中，叠加合并所述多个障碍物的安全区域，生成障碍物安全区。工区内多个地物构成的障碍物都具备不同的安全距离，多个障碍物的安全区会有部分重叠，将安全区域重叠的障碍物合并，融合多边形边界，叠加计算工区障碍物安全区范围。

影响井炮施工的障碍物的安全区域，反应在二维平面图上，就是一个个不同的多边形。障碍物的安全区域，就是多边形按照安全距离膨胀扩展，形成的区域构成新的多边形区域。相邻有重叠的多边形叠加，是利用地理信息系统GIS的缓冲区分析原理，相邻多边形叠加到一起，边界融合形成一个新的多边形。

使用步骤S104中统计的点按照采集设计偏移规则，自动计算偏移出障碍物安全区外炮点。即：按横纵方向进行炮点偏移，最大横、纵向偏移距离通过确定的最浅目的层的深度得到。从离障碍物安全区最近的位置开始，按整道距依次进行偏移，多个炮点均匀分布在障碍物两侧。其中，

在所述步骤S104中，将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置。其计算过程是：

计算过障碍物的沿测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；炮点S_i纵向偏移坐标为：

$\begin{array}{c}{S}_{ix}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)sin\left(\mathrm{\θ}\right)\\ S_{iy}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)\cos \left(\mathrm{\θ}\right)\end{array};$ 其中：

S_ix、S_iy为偏移后炮点S_i的x、y坐标偏移量；

i为1、2、3……m，m为所述过障碍物的沿测线的炮点构成线上落在障碍物中心点一侧的点数；

O_d为沿测线构成线离中心位置最近的炮点到障碍物边界的距离；

R_d为道距；j为1、2、3……j＝i，若炮点S_i偏移后落在已有的炮点上，j＝i+1；

θ为测线的方位角。

若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直测线方向移动，直到移动到合适的位置。其计算过程是：

计算过障碍物的垂直测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；炮点Sx_i横向偏移坐标为：

$\begin{array}{c}{\mathrm{Sx}}_{ix}=\left(\mathrm{int}\right(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}){\mathrm{RL}}_d+{\mathrm{jRL}}_d)\sin (\mathrm{\θ}\±\mathrm{\π}/2)\\ {\mathrm{Sx}}_{iy}=\left(\mathrm{int}\right(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}){\mathrm{RL}}_d+{\mathrm{jRL}}_d)\cos (\mathrm{\θ}\±\mathrm{\π}/2)\end{array};$ 其中：

Sx_ix、Sx_iy为偏移后炮点Sx_i的x、y坐标偏移量；

i为1、2、3……m，m为过障碍物的垂直测线的炮点构成线上落在障碍物中心点一侧的点数；

Ox_d为垂向构成线离中心位置最近的炮点到障碍物边界的距离；

RL_d为接收线距；

j为1、2、3……j＝i，若炮点Sx_i偏移后落在已有的炮点上，j＝i+1；

θ为测线的方位角。

如果炮点纵向偏移超过纵向最大偏移距离，同时横向偏移超过横向最大距离，则作为空点处理。

使用步骤S105的炮点偏移，计算炮点偏移后的观测系统属性，保证面元内属性分布均匀。这一步步骤，是为了最后检测炮点偏移，是否可以保障观测系统面元属性的均匀性，满足观测系统的变观设计需求。

覆盖次数、炮检距分布和方位角分布等直接反应观测系统面元属性均匀性，是衡量观测系统设计优劣的重要指标，最终影响后续资料处理的和资料解释。由于三维观测系统是多方位的，覆盖次数均匀一些，各方位上的远、近炮检距基本均匀，可以能够获取多方位信息，从而提高成像效果，在资料处理中求取的会更准确一些，使解释结果更接近地下的真实情况。

在本实施例中，当井炮施工变动后，所述障碍物的安全区域发生变动，按照所述步骤S102、S103、S104进一步将炮点偏移出所述障碍物安全区外。

通过上述实施例，可以得到，本发明实施例的基于障碍物安全区的炮点偏移方法，可以根据观测系统、障碍物安全区准确地计算出落在障碍物安全区域内炮点偏移。具体实施例：

本发明可以应用到复杂地表观测系统变观设计中，可以实现根据障碍物安全区进行炮点偏移，同时也可以实现同一观测系统不同炮点偏移下的覆盖次数和炮检距的对比，具体实现情况如下：

设置观测系统：32L3S180R，道距40米，炮距40米，接收线距120米，炮线距120米；

(1)要实现观测系统的变观，需要有测量实测、高清卫星照片、航片提供的障碍物位置形状类型分布图，如图2所示。

(2)如图2所示，计算这些障碍物安全区。例如房屋使用安全距离50米，堤坝安全距离100米，管线安全距离75米。其中，图3左图为障碍物分布图，图3中间的图式为每个障碍物的安全区分布图，图3右图为进行叠加以后的障碍物安全区分布图。

(3)进行炮点偏移计算，如图4所示，可以看出炮点都已经偏移出设定的安全区域。其中，图4中圆圈表示炮点。

(4)如图5所示，为炮点偏移后的观测系统覆盖次数图；图6为炮点偏移后的观测系统炮检距范围图。由图5和图6可以看出，炮点偏移满足要求，可以提供给野外井炮施工作为参考。

覆盖次数、炮检距分布和方位角分布等直接反应观测系统面元属性均匀性，是衡量观测系统设计优劣的重要指标，最终影响后续资料处理的和资料解释。覆盖次数、炮检距分布图能展示观测系统面元属性的均匀性，可以说明炮点偏移后，能够满足观测系统面元属性的均匀性。

(5)可以根据采集施工变动，障碍物的安全距离变动，重新进行炮点偏移计算。

例如，在(2)中为了保证浅层资料不缺失，在障碍物附近减小药量，对安全距离进行变动，房屋使用安全距离40米，堤坝安全距离70米，管线安全距离50米，重新计算障碍物的安全区域，然后进行炮点偏移。

综上所述，本发明实施例的基于障碍物安全区的自动炮点偏移方法，可以自动计算炮点偏移位置，为三维地震采集施工选择提供炮点偏移参考。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种基于障碍物安全区的自动炮点偏移方法，其特征在于，所述基于障碍物安全区的炮点偏移方法包括：

步骤1，获取多个障碍物类型以及范围；

步骤2，根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域；

步骤3，叠加合并所述多个障碍物的安全区域，生成障碍物安全区；

步骤4，计算落在所述障碍物安全区内的炮点，并按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外；

在所述步骤4中，所述按照采集设计炮点偏移的原则，将所述落在所述障碍物安全区内的炮点偏移出所述障碍物安全区外，包括：

将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置；

若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直测线方向移动，直到移动到合适的位置。

2.根据权利要求1所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，在所述步骤1中，所述获取多个障碍物类型以及范围，包括：

通过野外测量组实地实测影响炮点施工地物作为障碍物；

根据工区内的卫星照片或高清航片定义障碍物的类型。

3.根据权利要求2所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，在所述步骤2中，所述根据所述多个障碍物类型以及范围，确定所述多个障碍物的安全区域，包括：

根据所述障碍物的建筑年份、建筑构架以及管线走向，确定所述障碍物的安全区域。

4.根据权利要求3所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，在所述步骤3中，叠加合并所述多个障碍物的安全区域，生成障碍物安全区，包括：

将安全区重叠的障碍物合并，融合多边形边界，叠加计算工区障碍物安全区。

5.根据权利要求4所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，在所述步骤4中，将需要偏移的炮点首先沿测线纵向开始偏移，炮点偏移的位置为道距的整数倍，炮点从所述障碍物安全区的边界，按所述道距的整数倍，沿测线纵向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整道距沿测线方向移动，直到移动到合适的位置，包括：

计算过障碍物的沿测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；

炮点S_i纵向偏移坐标为：

$\begin{array}{c}{S}_{ix}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)sin\left(\mathrm{\θ}\right)\\ S_{iy}=\left(\mathrm{int}\right(\frac{O_d}{R_d})R_d+{\mathrm{jR}}_d)\cos \left(\mathrm{\θ}\right)\end{array};$ 其中，

S_ix、S_iy为偏移后炮点S_i的x、y坐标偏移量；

i为1、2、3……m，m为所述过障碍物的沿测线的炮点构成线上落在障碍物中心点一侧的点数；

O_d为沿测线构成线离中心位置最近的炮点到障碍物边界的距离；

R_d为道距；

j为1、2、3……j＝i，若炮点S_i偏移后落在已有的炮点上，j＝i+1；

θ为测线的方位角。

6.根据权利要求5所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，在所述步骤4中，若炮点移动超过纵向最大距离，在原位置垂直测线横向偏移，偏移位置为接收线距的整数倍，炮点从所述障碍物安全区的边界，按所述接收线距的整数倍，垂直测线横向方向偏移出所述障碍物安全区，如果在偏移位置已经存在炮点，跳过已经存在炮点接着按整接收线距垂直测线方向移动，直到移动到合适的位置，包括：

计算过障碍物的垂直测线的炮点构成线与障碍物的交点，取两交点中心位置作为炮点偏移方向控制，炮点在两交点中心位置的同侧偏移；

炮点Sx_i横向偏移坐标为：

$\begin{array}{c}{\mathrm{Sx}}_{ix}=\left(\mathrm{int}\left(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}\right){\mathrm{RL}}_d+{\mathrm{jRL}}_d\right)\sin \left(\mathrm{\θ}\±\mathrm{\π}/2\right)\\ {\mathrm{Sx}}_{iy}=\left(\mathrm{int}\left(\frac{{\mathrm{Ox}}_d}{{\mathrm{RL}}_d}\right){\mathrm{RL}}_d+{\mathrm{jRL}}_d\right)\cos \left(\mathrm{\θ}\±\mathrm{\π}/2\right)\end{array};$ 其中，

Sx_ix、Sx_iy为偏移后炮点Sx_i的x、y坐标偏移量；

i为1、2、3……m，m为过障碍物的垂直测线的炮点构成线上落在障碍物中心点一侧的点数；

Ox_d为垂向构成线离中心位置最近的炮点到障碍物边界的距离；

RL_d为接收线距；

j为1、2、3……j＝i，若炮点Sx_i偏移后落在已有的炮点上，j＝i+1；

θ为测线的方位角。

7.根据权利要求1所述的基于障碍物安全区的自动炮点偏移方法，其特征在于，还包括：

当井炮施工变动后，所述障碍物的安全区域发生变动，按照所述步骤2、3、4进一步将炮点偏移出所述障碍物安全区外。