一种基于控制振速下确定城区隧道爆破药量的安全施工区间的方法
【专利摘要】本发明公开了一种基于已知控制振速下进行城市浅埋隧道爆破设计时，被保护建筑位于隧道轮廓线外，确定不同设计药量所对应的隧道施工位置的方法。本发明根据隧道及建筑物的形状特征，采用合理方法建立隧道曲线、建筑物轮廓线的坐标系和轨迹方程；以建筑物轮廓延长线将隧道曲线部分分成若干段，推导出与设计药量相对应的隧道段安全距离的分段函数表达式。建筑物到爆源最近距离对应的最小药量与进尺决定后，根据求得的函数表达式，可以精确获取各爆破药量在隧道的分界位置/桩号，确定隧道设计进尺的安全范围。为复杂建筑环境下确定隧道进尺的安全爆破范围提供了一种精确的方法。
【专利说明】
-种基于控制振速下确定城区隧道爆破药量的安全施工区间 的方法
技术领域
[0001] 本发明设及一种爆破中安全施工区间确定的方法，特别适于保证低振速要求的隧 道轮廓线外被保护建筑物存在的情况下爆破时安全区间的确定。
【背景技术】
[0002] 近年来城市隧道爆破成为研究热点，但对一定振速下确定安全药量对应的施工位 置研究较少。由于爆源与被保护物间距离变动导致设计的分段药量/进尺不断变化，实际施 工通常是在距被保护物最近点先按最小药量掘进一段距离，此段内距离增加药量仍保持不 变；当距离增大到可按正常药量施工的位置时再增加药量。但现在城市浅埋隧道面临的普 遍问题是，隧道埋深与建筑物尺寸相近，不能将建筑作为一个点处理，故一定振速、药量下 得到对应施工区段变得非常困难，若已知被保护建筑所需最小药量与进尺后，最小进尺与 正常进尺的分界粧号是必须精确确定的，否则振速极易超标。
【发明内容】
[0003] 本发明方法是针对隧道爆破振速控制的要求，克服传统的将被保护建筑物简化为 一点，而提出的一种确定隧道爆破安全区间的方法。本发明方法尤其适用于复杂环境下的 城市浅埋隧道爆破安全药量对应区间的确定。
[0004] 本发明方法的原理是：
[0005] 建立隧道曲线的坐标系，得出隧道曲线方程及建筑物墙根线方程;把墙根线延长 与隧道曲线相交，交点和隧道圆弧段与缓和曲线段的分界点可把隧道划分为不同类型的区 段;根据所划分的区段的不同类型分别推导出隧道的粧号对应距被保护建筑物的距离公式 S-R;由已知交点和分界点粧号可知各点所对应的安全距离W及最小距离;根据正常施工进 尺可得对应安全药量，即可反算出其对应安全距离，将其由交点和分界点计算得出的距离 对比可知正常进尺所对应的区段，将正常进尺对应的安全距离带入S-R关系式就可得到正 常进尺对应的隧道粧号。
[0006] 本发明的基于控制振速下确定城区隧道爆破药量的安全施工区间的方法，该方法 包括如下步骤：
[0007] 1)隧道轨迹坐标系的建立与方程表述：
[000引隧道轨迹曲线采用圆弧曲线、缓和曲线和直线段组合，W隧道掘进首爆点轨迹连 线作隧道曲线，按圆弧曲线、缓和曲线两段分别建立坐标系，建立坐标系后可得隧道圆弧段 轨迹方程和缓和曲线段轨迹方程；
[0009] 2)建筑物墙根线方程表述：
[0010] 建筑物各边线轨迹方程称为墙根线方程，药量计算中最短距离产生于爆源与建筑 物在地面轮廓线中的某点联线，计算时W建筑物底层墙根处振速作为控制基准，找到墙根 线拐角点的坐标，从而得出建筑物墙根线方程；
[0011] 3)最短距离计算分段方法：
[0012] 对建筑物墙根线作延长线，将建筑物平面划分为两个部分，一部分为建筑物边线 为宽度的延长带，此区域内隧道到建筑物的最短距离即为隧道上点到对应墙根线的垂直距 离；另一部分为正对建筑物墙角的部分，此区域内隧道到建筑物的最短距离即为隧道上点 到对应墙角的水平距离；
[0013] 建筑物在爆屯、曲线圆弧外，延长墙根线分别与隧道曲线交于4个控制点，把隧道曲 线划分为5段，3种类型:直线段、缓和曲线段、圆弧曲线段，3种类型中又可分为隧道的缓和 曲线段上的点到建筑物边线或角点的距离，隧道圆弧段上的点到建筑物边线或角点的距离 四种情况；
[0014] 4)爆源到建筑物距离与隧道相应位置关系的推导：
[0015] 首先确定隧道爆屯、曲线圆弧段和缓和曲线段方程，又由分段方法对隧道曲线的划 分，分别将对应的隧道曲线方程与空间中点到点、点到直线的距离公式联立就可导出步骤 (3)中所述四种情况下隧道爆屯、曲线上的任一点到建筑物的距离公式S-R，即，隧道曲线上 任意粧号S所对应的位置到建筑物的最短距离R;
[0016] 5)已知S-R关系式后的具体求解：
[0017] ①建筑物墙根线延长至隧道曲线得到的四个控制点及隧道缓和曲线与圆曲线交 点HY的粧号，代入式S-R关系式求得四个控制点相应距离，并找到S-R关系式的最小值Rmin;
[0018] ②将设计药量/进尺带入萨氏公式得到建筑物距爆源空间距离R，判断得R在步骤 ①中的对应区段；
[0019] ③在该区段内选用S-R关系式中相应公式，用MATLAB的solve函数计算程序，筛除 不在本段位置解的if判值函数，编写程序自动逆解并筛选出符合条件S，得到两个粧号S值 分别对应隧道至建筑物距离为R两处位置，从而确定城区隧道爆破药量的安全施工区间。
[0020] 其中，在所述步骤4中，所述S-R公式中，已知S即可求得R;反之，已知R也可反算出 S。
【附图说明】
[0021 ]图1为隧道轨迹曲线建系示意图；
[0022] 图2为墙根线延长线的平面划分；
[0023] 图3为墙根线切割隧道示意图；
[0024] 图4为坐标系的建立及曲线分段；
[0025] 图5为渝中连接隧道周边建筑环境示意图；
[0026] 图6为渝中连接隧道走向及建筑物位置图；
[0027] 图7为隧道粧号一空间距离变化曲线；
[0028] 图8为隧道粧号一最大单段药量变化曲线；
[0029]
【具体实施方式】
[0030] 本发明公开了一种基于已知控制振速下进行城市浅埋隧道爆破设计时，被保护建 筑位于隧道轮廓线外，确定设计药量所对应的隧道施工位置的方法，
[0031] 本发明方法的步骤是：
[0032] 1)隧道轨迹坐标系的建立与方程表述:隧道轨迹曲线采用实际工程中普遍使用的 圆弧曲线、缓和曲线和直线段组合，W隧道掘进首爆点轨迹连线作隧道曲线，按圆曲线、缓 和曲线两段分别建立坐标系；
[0033] ①圆弧段坐标系的建立:设原点0与隧道变坡点同高，且原点0的水平投影与圆弧 圆屯、重合，Z轴指向竖直正上方，Y轴指向隧道变坡点，X轴按右手螺旋法则由Y轴叉乘Z轴得 到，建立坐标系后可得隧道圆弧段轨迹方程：
[0034]
[0035] 式中：r为圆弧段半径，Sl为Y轴与隧道曲线交点粧号，S为隧道曲线任意点(x，y，z) 粧号，i为隧道坡度，r为圆曲线半径，S为缓和曲线上任一点(x，y，z)到直缓点ZH弧长；
[0036] ②缓和曲线段坐标系的建立:缓和曲线用于连接直线段和圆曲线，节点称直缓点 (ZH)和缓圆点化Y)。建立坐标系为：原点0'位于直缓点处，Z轴指向正上方，Y轴背向缓和曲 线曲率中屯、一侧，垂直于直线段，X轴按右手系垂直于Y轴，建立坐标系后可得隧道缓和曲线 段轨迹方程：
[0037]
[003引式中:r为圆弧段半径，i为隧道坡度，So为缓和曲线总长度，ZZH为直缓点Z坐标；
[0039] 2)建筑物墙根线方程表述:一般地，因建筑物底层与其它层相比距地下爆源最近， 分段药量计算中最短距离只可能产生于爆源与建筑物在地面轮廓线中的某点联线，因此计 算时W建筑物底层墙根处振速作为控制基准。建筑物地面轮廓线为长方形，为方便将各边 线轨迹方程称为墙根线方程。
[0040] 设建筑物墙根拐角处坐标A(xo，yo，zo)，找到墙根线拐角点的坐标，即可知道墙根 线方程为：
[0041]
[0042] 式中:kl-墙根线斜率，ZO-建筑物相对高程最低点坐标；
[0043] 3)最短距离计算分段方法:对建筑物墙根线作延长线，可W将建筑物平面划分为 两个部分（图2)，一部分为建筑物边线为宽度的延长带，此区域内隧道到建筑物的最短距离 即为隧道上点到对应墙根线的垂直距离；另一部分为正对建筑物墙角的部分，此区域内隧 道到建筑物的最短距离即为隧道上点到对应墙角的水平距离；
[0044] 建筑物在爆屯、曲线圆弧外，4条延长线分别与隧道曲线交于4个控制点DEFG，把隧 道曲线划分为5段（图3);直缓点为0'（或ZH)，缓圆点HY位于EF之间（图4)。对于隧道DG段之 间的部分，计算隧道上各点到建筑物的最短距离就分为G-F段、F-HY段、HY为E段、E-D段，分 别为:求缓和曲线上点到直线距离、缓和曲线上点到点的距离、圆弧曲线上点到点的距离、 圆弧曲线上点点到直线的距离四种。运四种情况是最为常见的情形，具有一般代表性。
[0045] 4)爆源到建筑物距离与隧道相应位置关系的推导：隧道曲线上任意粧号S所对应 的位置到建筑物的最短距离R为：
[0046]
[0047] 式中：kAM为建筑物的AM边线在坐标系XOY下的斜率;k'AN为建筑物的边线AN在坐标 系X'Y'Z'下的斜率，（XA，yA，ZA)为点A在坐标系XOY下的坐标；（x'A，y'A，z'A)分别为点A在坐 标系X'Y'Z'下的坐标，Sd,泣，Shy，Sf，Sg，Szh为D、E、HY、F、G、ZH点的粧号；
[004引S-R式中仅有粧号S-个变量，其他参数在工程中均为常数，故已知粧号的情况下， 即可算出对应的隧道上点到建筑物的最短距离;反之在给定安全距离R(即萨道夫斯基公式 中的距离参数R)时，可求得对应粧号。
[0049] 5)已知S-R关系式后的具体求解：
[0050] ①建筑物墙根线AM、AN延长至隧道曲线交于D、F、E、G及隧道缓圆交点HY的粧号，将 其代入式S-R对应方程求得相应距离Rd、Rf、Re、Rg、姑Y，并找到S-R式最小值Rmin;
[0051] ②将设计药量/进尺带入萨氏公式得到建筑物距爆源空间距离R，判断得R在步骤1 中的对应区段；
[0052] ③在该区段内选用式S-R中相应公式，用MATLAB的solve函数计算程序，筛除不在 本段位置解的if判值函数，编写程序自动逆解并筛选出符合条件S，得到两个粧号S值分别 对应隧道至建筑物距离为R两处位置。
[0053] 现根据具体的实施方式对其进行说明。
[0054] 本发明所依托的工程为重庆市渝中区渝中连接隧道，隧道上方为渝中区主要商业 中屯、，建构筑物分布密集。隧道曲线由圆曲线+缓和曲线+直线组成，隧道坡度i = 4%，圆曲 线半径r = 260.25m，缓和曲线参数s〇 = 200m。隧道轮廓线外有重点文物罗汉寺大雄宝殿，控 制振速0.5cm/s W下，其余所有建筑物要求振速不超过1. Ocm/s。因文物罗汉寺大雄宝殿距 离隧道较近且控制振速低于正上方建筑物，不能把其简化为一点考虑，在实际施工过程中 就需要准确确定由大雄宝殿距离隧道的最小距离确定的小进尺变化为正常施工进尺时的 隧道粧号。
【具体实施方式】 [0055] 包括W下步骤：
[0056] 1)隧道轨迹坐标系的建立与方程表述:根据隧道左线爆屯、曲线作图6，Y轴与曲线 交点（变坡点)位置粧号ZK14+072.268(粧号S只取V'后面部分），设为高程起始点，即Y轴与 隧道曲线交点粧号为Si = 72.268。从粧号ZK14+072.268~400,圆曲线、缓和曲线两部分。粧 号ZK14+200处为缓圆点化Y)，粧号ZK14+400处为直缓点（ZH)，即Shy = 200，Szh = 400dW直缓 点为原点建系O' X' Y' Z '，并保证坐标系O ' X' Y' Z '与OXYZ各自的XY平面重合。建立坐标系后 分别推导出隧道爆屯、曲线在两个坐标系下的缓和段和圆弧段曲线方程。
[0057] 2)建筑物墙根线方程表述：找出建筑物的角点A、M、N分别在坐标系0' X' Y' Z '与 OXYZ下的坐标，确定角点坐标后就可导出两个坐标系下墙根线的方程。
[0058] 3)最短距离计算分段方法:对墙根线作延长线，可将建筑物平面划分为①②两个 部分，如图2所示，其中区域①内隧道到建筑物的最短距离即为隧道上点到对应墙根线AM、 AN的垂直距离；其中区域②内隧道到建筑物的最短距离即为隧道上点到对应角点A、M、N的 水平距离;延长线与隧道爆屯、曲线交于〇、6立、6四点，可将爆屯、曲线分为5段，结合已知缓圆 点HY可W划分为:GF段、F点到缓圆点HY、缓圆点HY到E、邸段四段。
[0059] 4)爆源到建筑物距离与隧道相应位置关系的推导:GF段完全处于缓和曲线部分。G 点粧号ZK14+337.2，F点粧号ZK14+269.7，A点在坐标系O'X'Y'Z'下坐标为（203.6835，- 47.0942，27.0363)，直线AN斜率为1.75。带入参数后，参照式S-R中第一式计算;F点到缓圆 点HY为缓和曲线段，F点粧号ZK14+269.7，缓圆点HY粧号ZK14+200，参照S-R式中第二式计 算;缓圆点到E点为圆弧段，参照S-R式中第S式计算。E点粧号ZK14+169.46，缓圆点HY粧号 ZK14+200，A点在OXYZ坐标系下的坐标为(95.8053,262.9929,27.0363) ;ED段及DW下全部 是圆曲线。E点粧号ZK14+169.46，D点粧号ZK14+156.70，建筑物AM边在0乂町坐标系下直线斜 率为0.7892，带入参数，参照S-R式中第四式计算。
[0060] 5)已知S-R关系式后的具体求解:在上述各区段明确后，采用MTLAB 2014a程序对 式S-R编程求解，就可根据隧道爆源轨迹上每点到建筑物的空间距离R，很方便获取隧道上 相应的位置粧号。通过式S-R求得H点粧号ZK14+189.5，最小值Rmin = 30.41m(如图6)，现需求 出可采用1.0kg药量施工的隧道临界位置并确保振速在0.5cm/sW下。根据萨氏公式算出 1. Okg药量对应的爆源距离36.7m，带入式S-R第二式和第四式，均仅有粧号S-个参数，采用 Matlab中的solve函数编程求解方程，即得到对应粧号ZK14+210.13、ZK14+165.06,如图6 中B、C所示。最近点H位于BC之间。其工程意义是从H点Wo.7m进尺/0.5kg药量施工，当向前 (向B)施工到ZK14+210.13、或向后（向C)施工到ZK14+165.06粧号时，装药量增加为1.0kg (或进尺增为1.8m)，可满足大雄宝殿安全振速控制要求。表1为左线各区段参考药量表。
[0061]
[0062] 表 1
[0063] 如果渝中隧道设计时需考虑不同分段药量时的施工安全距离，利用S-R式和 Matlab软件编程计算隧道爆屯、曲线上的各点到建筑空间最短距离，绘制出对应的粧号一最 短距离变化曲线（S-R)，如图7所示。又因为在隧道爆破设计中，主要考虑装药量-隧道对应 区段的施工。故根据萨道夫斯基公式，将建筑物到爆源距离带入即可得粧号一药量变化曲 线（图8)。在W最小装药量进尺施工后，W多大药量或进尺施工并满足振速不超标就在图中 一目了然表现出来，可W非常方便地进行不同地段的爆破药量设计，例如欲设计I.化g的分 段药量，只能在K14+158.29向减小粧号方向、或K14+213.93向增加粧号方向施工才能满足 振速要求。运对隧道安全振速下的药量设计有重要应用价值。 W上为本发明方法的一种具体实施案例，本发明并不仅局限于特定的实例，只要采用 相同方法来求解当隧道附近被保护建筑物不可简化为一点时，爆破由小进尺变化到正常进 尺时的区间范围，均应视为本发明方法的同一类型。 W上所述仅是本发明优选实施方式，应当指出，对于本技术领域的普通技术员来说，在 不脱离本发明技术原理的前提下，还可W做出若干改进和变型，运些改进和变型也应该视 为本发明的保护范围。
【主权项】
1. 一种基于控制振速下确定城区隧道爆破药量的安全施工区间的方法，其特征在于， 该方法包括如下步骤： 1) 隧道轨迹坐标系的建立与方程表述： 隧道轨迹曲线采用圆弧曲线、缓和曲线和直线段组合，以隧道掘进首爆点轨迹连线作 隧道曲线，按圆弧曲线、缓和曲线两段分别建立坐标系，建立坐标系后可得隧道圆弧段轨迹 方程和缓和曲线段轨迹方程； 2) 建筑物墙根线方程表述： 建筑物各边线轨迹方程称为墙根线方程，药量计算中最短距离产生于爆源与建筑物在 地面轮廓线中的某点联线，计算时以建筑物底层墙根处振速作为控制基准，找到墙根线拐 角点的坐标，从而得出建筑物墙根线方程； 3) 最短距离计算分段方法： 对建筑物墙根线作延长线，将建筑物平面划分为两个部分，一部分为建筑物边线为宽 度的延长带，此区域内隧道到建筑物的最短距离即为隧道上点到对应墙根线的垂直距离； 另一部分为正对建筑物墙角的部分，此区域内隧道到建筑物的最短距离即为隧道上点到对 应墙角的水平距离； 建筑物在爆心曲线圆弧外，延长墙根线分别与隧道曲线交于4个控制点，把隧道曲线划 分为5段，3种类型:直线段、缓和曲线段、圆弧曲线段，3种类型中又可分为隧道的缓和曲线 段上的点到建筑物边线或角点的距离，隧道圆弧段上的点到建筑物边线或角点的距离四种 情况； 4) 爆源到建筑物距离与隧道相应位置关系的推导： 首先确定隧道爆心曲线圆弧段和缓和曲线段方程，又由分段方法对隧道曲线的划分， 分别将对应的隧道曲线方程与空间中点到点、点到直线的距离公式联立就可导出步骤(3) 中所述四种情况下隧道爆心曲线上的任一点到建筑物的距离公式S-R，即，隧道曲线上任意 粧号S所对应的位置到建筑物的最短距离R; 5) 已知S-R关系式后的具体求解： ① 建筑物墙根线延长至隧道曲线得到的四个控制点及隧道缓和曲线与圆曲线交点HY 的粧号，代入式S-R关系式求得四个控制点相应距离，并找到S-R关系式的最小值心^ ② 将设计药量/进尺带入萨氏公式得到建筑物距爆源空间距离R，判断得R在步骤①中 的对应区段； ③ 在该区段内选用S-R关系式中相应公式，用MATLAB的solve函数计算程序，筛除不在 本段位置解的if判值函数，编写程序自动逆解并筛选出符合条件S，得到两个粧号S值分别 对应隧道至建筑物距离为R两处位置，从而确定城区隧道爆破药量的安全施工区间。2. 根据权利要求2所述的方法，其特征在于，在所述步骤4中，所述S-R公式中，已知S即 可求得R;反之，已知R也可反算出S。
【文档编号】E21D9/14GK106014422SQ201610350097
【公开日】2016年10月12日
【申请日】2016年5月24日
【发明人】龚敏, 李永强, 陈华森, 孟祥栋, 吴昊骏, 聂永进, 胡广风, 宋书, 宋书一, 管学铭, 余鸿飞, 刘洪忆, 张明高, 赵振振, 王从刚
【申请人】北京科技大学, 重庆巨能建设(集团)有限公司