发明名称 一种基于炮孔钻进信息的隧道超前地质预报方法
摘要 本发明公开了一种基于炮孔钻进信息的隧道超前地质预报方法,其步骤:A.在预报实施前先建立可钻性指数的标准数据库[K<sub>i</sub>];B.超前地质预报实施的具体步骤是:1)对掌子面的超前炮孔进行编号;2)将采集得到的钻进速率v、钻进压力p、钻杆转速w及岩渣、回水量及回水颜色逐一记录;3)根据可钻性指数公式(7)<img file="DDA00003184459700011.GIF" wi="178" he="131" />计算每组参数的可钻性指数K;4)将计算的可钻性指数k与[K<sub>i</sub>]进行对比,若K∈(min[K<sub>i</sub>],max[K<sub>i</sub>]),则本钻进段的基本围岩级别为i级;5)在结合岩渣成分、回水情况和掌子面围岩情况,根据表2对围岩级别进行综合判定;6)综合预报的围岩级别与隧道开挖后的确定的实际围岩级别进行对比、验证;7)将修正后的数据扩充到[K<sub>i</sub>]中。本方法操作方便、费用低廉,准确高效,可以进行实时的超前地质预报。
申请公布号 CN103308946B 申请公布日期 2016.08.31
申请号 CN201310176656.4 申请日期 2013.05.14
申请人 中国科学院武汉岩土力学研究所 发明人 王浩;田湖南;焦玉勇;覃卫民;张国华;张秀丽
分类号 G01V1/40(2006.01)I;G01V1/28(2006.01)I 主分类号 G01V1/40(2006.01)I
代理机构 武汉宇晨专利事务所 42001 代理人 王敏锋
主权项 一种基于炮孔钻进信息的隧道超前地质预报方法,包括如下步骤:在炮孔的钻进中,对钻进信息进行采集,所述的钻进信息包括:钻进速率v、钻进压力p、钻杆转速N、岩渣、回水量及回水颜色;基于能量守恒原理,由采集到的钻进速率v、钻进压力p、钻杆转速N三个参数计算出围岩的可钻性指数Kˊ可钻性指数Kˊ的计算公式是<img file="FDA0001023599060000011.GIF" wi="211" he="118" />建立典型围岩级别的可钻性指数标准数据库[Ki],其步骤是:通过在i级围岩中进行n次炮孔钻进并获得钻进信息,由式<img file="FDA0001023599060000012.GIF" wi="185" he="117" />计算得到可钻性指数Ki1'、Ki2'、Ki3'……Kin',得到i级围岩的可钻性指数数据库[Ki]',i=Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ,分别表示Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级围岩;对i级围岩的可钻性指数数据库[Ki]'内的数据Ki1'、Ki2'、Ki3'……Kin',根据概率统计方法进行数据筛选,筛选后的数据作为可钻性指数标准数据库[Ki];数据筛选步骤包括:根据概率统计方法,i级围岩的可钻性指数数据库[Ki]'内的数据Ki1'、Ki2'、Ki3'……Kin′服从正态分布:可钻性指数的均值<img file="FDA0001023599060000013.GIF" wi="63" he="71" />为:<maths num="0001"><math><![CDATA[<mrow><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><msub><mi>K</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&prime;</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023599060000014.GIF" wi="1013" he="134" /></maths>标准差Si′为:<maths num="0002"><math><![CDATA[<mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msup><msub><mi>K</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&prime;</mo></msup><mo>-</mo><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023599060000015.GIF" wi="1006" he="159" /></maths>可钻性指数均值的置信水平为1‑α的置信区间为:<maths num="0003"><math><![CDATA[<mrow><mo>(</mo><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>-</mo><mfrac><mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mo>&prime;</mo></msup></mrow><msqrt><mi>n</mi></msqrt></mfrac><msub><mi>t</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>,</mo><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>+</mo><mfrac><mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mo>&prime;</mo></msup></mrow><msqrt><mi>n</mi></msqrt></mfrac><msub><mi>t</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023599060000016.GIF" wi="1030" he="158" /></maths>式中:α为分位点,α取值范围是0.01~0.40;t<sub>α</sub>/<sub>2</sub>(n‑1)为t分布中,样本数为n时,t分布中α/2分位点对应的数值,通过查阅概率论手册中的t分布表来确定;将符合:<maths num="0004"><math><![CDATA[<mrow><msup><msub><mi>K</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&prime;</mo></msup><mo>&Element;</mo><mrow><mo>(</mo><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>-</mo><mfrac><mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mo>&prime;</mo></msup></mrow><msqrt><mi>n</mi></msqrt></mfrac><msub><mi>t</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>,</mo><msup><msub><mover><mi>K</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>&prime;</mo></msup><mo>+</mo><mfrac><mrow><msup><msub><mi>S</mi><mi>i</mi></msub><mo>&prime;</mo></msup></mrow><msqrt><mi>n</mi></msqrt></mfrac><msub><mi>t</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>...</mo><mo>...</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023599060000021.GIF" wi="1541" he="155" /></maths>的数据纳入可钻性指数标准数据库中,形成i级围岩的可钻性指数标准数据库[Ki];通过计算出的可钻性指数与可钻性指数标准数据库[K<sub>i</sub>]进行对比,计算的可钻性指数K∈(min[K<sub>i</sub>],max[K<sub>i</sub>]),钻进段的围岩级别为i级;结合岩渣成分、掌子面围岩情况,对岩石性状进行预判;回水情况结合基本围岩级别,对围岩级别做出综合预报,至此完成对隧道掌子面前方围岩级别的预报;按照实际围岩级别对计算的可钻性指数K分类;按照可钻性指数标准数据库[K<sub>i</sub>]建立的步骤,将计算的可钻性指数K扩充到可钻性指数标准数据库[K<sub>i</sub>]中,完成数据反馈和数据库修正。
地址 430071 湖北省武汉市武昌区水果湖街小洪山2号