一种资源勘查钻孔轨迹精准恢复的计算方法与流程

1.本技术涉及资源勘查钻孔的技术领域，尤其是涉及一种资源勘查钻孔轨迹精准恢复的计算方法。


背景技术：

2.资源勘查钻孔是利用钻探设备向地下钻进形成的直径小且深度数百米乃至数千米的柱状圆孔。资源勘探施工中钻孔常常发生偏斜现象，虽然这是一种常见现象，但是因孔斜数据计算精度不足、钻孔轨迹不能精准恢复等原因，当偏斜较大时钻孔评级会受到很大影响，甚至是造成废孔。
3.造成偏斜的原因是在钻进过程中，有岩性、钻探装备性能和钻工操作等。一是岩石的各向异性在钻头处产生歪倒力矩使粗径钻具倾斜而导致孔斜；二是钻进软硬夹层时，钻头处产生钻压差或钻头沿岩石节理面滑移；三是钻进中遇到溶洞或大裂隙；四是在卵砾石层中钻进；在钻进工艺及操作方面的原因：换径不带导向或粗径钻具本身就弯曲；冲洗液量大冲刷孔壁使孔壁间隙过大，造成粗径钻具或岩心管与钻孔轴线产生偏斜角；采用的钻压过大与转速过快，钻杆柱受轴心压力及回转的影响而呈挠曲状态影响了钻具的稳定性，促使粗径钻具歪倒。
4.偏斜作用使钻孔轨迹由设计的垂直直线形变成了弯曲的空间曲线形态，表现为钻孔某些深度、大部井段乃至整个井段的天顶角和方位角与设计存在较大偏差，这种现象一般称为孔斜。孔斜较大或高级勘查阶段必须进行校正，而目前勘探单位最常使用的方法是根据测井工作的孔斜测量数据按折线法进行计算加以校正，此方法操作复杂、精度不高，在填表计算过程中还易于发生人为因素导致的错误。有必要提出更高精度、适应性更好的孔斜校正方法。
5.但是要想精确判断勘探钻孔发生孔斜的程度，就需要求取钻孔的准确空间轨迹。以往通常的做法是在钻探过程中，每隔一段距离，如50m、25m，就测量一次钻孔孔斜数据，钻孔孔斜数据包括天顶角、方位角和孔深；根据钻孔各测点的天顶角、方位角和孔深这些基本要素，采用分段折线的方法计算出各测点的空间坐标，而处于两个测点之间的地层界面，则按线性内插的方式进行计算。
6.以往对孔斜校正的研究并不少，但基本都是沿用折线法分段累进计算原理方式实现的。一般仅着眼于勘查任务所关注的几个层点空间坐标的计算目的，尽管有的实现了自动计算，但并没有消除折线法所带来的误差；即使是现场测井时进行了局部(多在资源层段)加密，也无法消除加密段以前的系统误差；这些方法都存在已知数据固定有限，且无通过数学方法合理增加计算数据的思路和做法，所以也就谈不上钻井轨迹的精准恢复。因此，有限数据来源的折线分段累进不可避免地存在由有限测点所决定的折线计算所带来的误差。同时，如果想计算更多的层点，还需要人为判断或指出各层点所在折线段位置等手工工作。
7.资源钻孔存在偏斜是普遍性现象，要完成孔斜校正等目的就需要精准恢复钻孔的
空间轨迹。目前对钻孔孔斜的数学表达是由测井成果中的方位角、天顶角与井深的有限离散测量成果。以往要提高孔斜计算的精度的做法只有加密测点密度，但一般只针对资源层段进行加密，资源层段以上仍然由较稀的测点控制。所以当前有必要提出更高精度、适应性更好的孔斜校正方法。


技术实现要素：

8.为了提出更高精度、适应性更好的孔斜校正方法，本技术提供一种资源勘查钻孔轨迹精准恢复的计算方法。
9.本技术提供的一种资源勘查钻孔轨迹精准恢复的计算方法采用如下技术方案：一种资源勘查钻孔轨迹精准恢复的计算方法，包括：s1、建立原始数据文件：原始数据包括勘查钻孔的孔口坐标与高程、测井工作获得的孔斜测量成果数据；孔斜测量成果数据包括间隔设置的钻井深度对应的天顶角和方位角数据；s2、根据原始数据文件，确定方位角和天顶角的初始边界条件；方位角的初始值按散点曲线趋势延至井口即可，天顶角的初始值设为0；s3、根据测井工作的孔斜测量成果数据和确定的初始边界条件，可选择多种计算方法，包括线性内插、样条曲线内插法、多项式拟合内插和抛物线内插，实现拟合出方位角和天顶角随钻井深度h变化的曲线的目的；s4、将拟合的曲线转换成密点离散数据：在拟合曲线上进行密点再采样；密点采样间隔δh的可选参数有1、2、5m等。选择的原则要求是，测井曲线越复杂、计算精度要求越高，选择的密点采样间隔应越小；2m一般均能达到精准恢复实际钻孔轨迹的目的，可将默认值设定为2m。
10.s5、求取钻井轨迹：自井口算起，第i个井深间隔δh的水平投影长度δli由下式求得：δli＝δh
×
sinθi；第i个井深间隔dh的垂直投影长度δhi由下式求得：δhi＝δh
×
cosθi；第i个井深间隔dh的水平投影所对应的地理坐标增量δxi和δyi由下式求得：由下式求得：设钻孔孔口高程和坐标分别为h0、x0、y0，则钻孔第m段尾点的空间坐标及高程hm、xm、ym分别由下式计算：按上述公式计算m＝1直到最大井深对应的采样数后，即可求得钻孔实际轨迹。
11.针对任何一个已知深度的地层界面，运用本方法均可自动计算出密点采样间隔前提下其所在的位置，即可以计算出地层分界面位于第m和m 1个采样点之间时，在完成m段计算的基础上，将m 1段内的坐标与高程的增量累计进去即可，高程的增量相应的设为δxi、δyi、δhi，由下式计算：：
12.所述步骤s4中密点采样间隔δh的可选参数有1、2、5m等，选定密点采样间隔后进行步骤s5中的计算，密点采样间隔越小，钻孔轨迹恢复的效果越好，即精度越高。通过多个模型计算分析，设定拟合密点采样密间隔为2m，绝大多数情况下均能满足钻孔轨迹的精准恢复。故可将密点采样间隔默认为2m。
13.通过采用上述技术方案，能够由有限的孔斜测井数据求取精度符合要求的连续曲
线，即和f
θ
＝f(h)，由这两条曲线就能够精准表达钻孔轨迹。这个思路解决了前人在进行钻孔轨迹恢复时，无法有更多的数据来源的问题，提高了计算精度和计算效率；本方法可以便于计算任意井深处的空间坐标，为矿山数字建模的钻孔成果数据建库打下了坚实基础；本方法还可为测井规范中有关孔斜测量要求的优化提供支撑，减少目前实际测井工作中为提高精度采用多种采样间隔反复测井带来的效率低、劳动强度大等问题，从而降低有可能带来的孔内事故的隐患。
14.综上所述，本技术包括以下至少一种有益技术效果：1.运用多种数学曲线拟合算法对天顶角、方位角与钻孔深度的关系的离散曲线进行拟合可得到方位角、天顶角随钻井深度的连续曲线，该曲线即可视作真实的方位角、天顶角随钻井深度的连续曲线；对所得的曲线进行密点采样，再运用积分算法进行计算即可实现法钻孔轨迹的精准恢复；2.运用本方法进行编程，有明确的思路、实现简单、既可完成原始数据的交互输入、计算成果表格的自动输出，还可形成钻孔轨迹的三维立体显示和平剖面投影显示。同时，本方法完善了计算方法，提高了工作效率和成果精度，还可避免以往方法人工填表、选择等步骤容易产生错误的问题；3.采用本技术的计算方法，可得到该钻孔内任何指定层位的空间坐标，计算的任意一层的空间坐标都可以直接用于该层平面成果图的编制；在完成钻孔轨迹精准恢复的基础上，还可以实现向任意走向地质剖面进行钻孔投影；精准恢复的钻孔轨迹也是数字化矿山的重要基础资料之一。所有这些都可以由编程来实现。
15.4.本计算方法还可对资源勘查有关测井规范中关于孔斜测量的要求的优化提供支撑，减少目前实际测井工作中为提高精度采用多种采样间隔反复测井带来的效率低、劳动强度大等问题，从而降低有可能带来的孔内事故的隐患。
附图说明
16.图1为两种理论模型的方位角；图2为两种理论模型的天顶角曲线；图3为理论2的钻孔轨迹水平投影；图4为不同测斜点间隔(50、20、10、5m)情况下传统折线法计算的钻孔轨迹水平投影与理论模型钻孔轨迹的对比；图5为测斜点间隔为50m时抛物线内插法与传统折线法计算的钻孔轨迹水平投影及与理论模型钻孔轨迹水平投影的对比；图6为测斜点间隔为20m时抛物线内插法与传统折线法计算的钻孔轨迹水平投影及与理论模型钻孔轨迹水平投影的对比；图7为理论模型进行50m间隔测斜点采样，在资源段加密测斜点至20m和10m的情况下采用传统折线法计算的钻孔轨迹水平投影；图8为理论模型进行20m间隔测斜点采样，在资源段加密测斜点至10m的情况下采用传统折线法计算的钻孔轨迹水平投影；图9为理论模型的孔斜参数曲线与实际测井曲线的孔斜曲线的复杂程度对比示意图；
图10为方位角初始值的确定示意图；图11为第i段的各种投影及增量示意图。
具体实施方式
17.以下结合附图1-11对本技术作进一步详细说明。
18.本技术实施例公开了一种资源勘查钻孔轨迹精准恢复的计算方法。
19.一种资源勘查钻孔轨迹精准恢复的计算方法包括：建立原始数据文件：本计算方法的数据依据为测井工作的孔斜测量成果数据(天顶角和方位角与观测深度的对应关系。观测间隔不超过20m时不需要进行局部加密测量和相应的数据)和勘查钻孔的孔口坐标与高程；工程关注(编制成果图)的界面(或其他目标)对应的钻进深度。
20.首先，根据初始测井的数据离散曲线，确定方位角和天顶角θ的边界条件(初始值)，原则为：天顶角的初始值恒定为0；方位角的初始值可由实测曲线沿趋势合理延伸至井口来确定。
21.根据测井测斜数据，选择数学内差方法，如线性内插、样条曲线内插法、多项式拟合内插、抛物线内插等，这些方法都有成熟的数学公式和算法。拟合出方位角和天顶角θ随钻井深度h的变化曲线，即能够由有限的测井孔斜数据求取精度符合要求的连续曲线，即和f
θ
＝f(h)，由这两条曲线就能够精准表述钻孔轨迹。这个方法思路解决了前人在进行钻孔轨迹恢复时无法有更多的数据来源的问题，也减少了目前实际测井工作时，为了提高精度采用多种采样间隔反复测井带来的低效问题和有可能的孔内事故的隐患。
22.将连续的曲线进行密点采样转换成离散数据。根据精度要求，选择拟合曲线的再采样，设密点采样间隔为δh。拟合曲线的再采样间隔进行选择，如设计成1m、2m、5m、10m等可选择参数。通过理论模型与实际测井资料计算，综合分析，δh不大于10m时一般情况下的精度要求能够满足；δh＝2m可满足更为复杂情况下的精度要求。因此，将δh＝2m设为默认参数。
23.拟合曲线进行密点采样后，则可以计算得到钻井深有多少个采样间隔；也可以求得第i个采样间隔δh的水平投影长度δli和垂直投影长度δhi。
24.精准的钻井轨迹求取(恢复)。自井口算起，第i个井深间隔δh的水平投影长度δli由下式求得：δli＝δh
×
sinθi(1)第i个井深间隔dh的垂直投影长度δhi由下式求得：δhi＝δh
×
cosθi(2)第i个井深间隔dh的水平投影所对应的地理坐标增量δxi和δyi由下式求得：由下式求得：设钻孔孔口高程和坐标分别为h0、x0、y0，则钻孔第m段尾点的空间坐标及高程hm、xm、ym分别由下式计算：
按上述公式计算m＝1直到最大井深对应的采样数后，我们就可以准确求得(或称恢复)钻孔实际轨迹。
25.如果某个地层分界面位于第m和m 1个采样点之间时，我们就可以在完成m段计算的基础上，将m 1段内的坐标与高程的微小增量(相应的设为δxi、δy、δh
ii
)累计进去即可，计算公式见下：算公式见下：算公式见下：对于任何一个界面深度而言，先按选定的计算间隔求出其所在的采样段数，继而根据公式(5)～(10)进行相应的计算，即可得到该钻孔该层位的空间坐标。所有这些都可以由编程来实现。
26.进行编程时，我们可设置了拟合加密间隔δh的选项，有1、2、5、10、20、25、50m等多个选项。10、20、25、50m三个选项主要用于本方法与目前常用方法进行对比分析，测试计算方法与编程效果，实际运用时并不选择。而1、2、5m的密点采样间隔则用于实际的钻孔轨迹精准恢复来选择，孔斜测井给出的天顶角较小且方位角变化较缓慢时，间隔可选择大一些，如5m；否则，情况较复杂且变化剧烈时则应选择小一点；默认情况下我们设定拟合加密间隔为2m，可以满足绝大多数情况。
27.参考图1和图2，展示的两个理论模型均较实际的钻孔方位角的变化更为剧烈，结果在这种情况下能够精准恢复设计模型的钻孔轨迹，方法运用于实际情况精度会更高。
28.参考图3，从钻孔轨迹的水平投影可以看出理论模型的复杂情况。
29.参考图4，随着全井段测井孔斜采样间隔的不断加密，传统的折线法也可以精确地恢复钻孔轨迹；但是实际情况下是不允许如此操作，一方面测井效率低下，另一方面容易出现孔内事故造成很大损失；过密的采样间隔会造成大量冗余数据，给常规计算表格的操作和计算带来麻烦。本技术就是不需要过度加密，也可以取得加密的效果和精度；而且本方法并不要求实际测井时一定要等间隔采样。
30.参考图5，即使是采用50m间隔的测斜数据，运用本方法，求取的钻孔轨迹与实际情况也很接近；如果用于较为简单与大多数实际情况相接近的模型进行计算，精度依旧可以提高。
31.参考图6，20m间隔的测斜数据，运用本方法，求取的钻孔轨迹与实际情况几乎一致，精度完全达到要求；运用其他的拟合曲线进行计算结论相同。这表明，本方法只需要一条采样间隔正常的钻井测斜成果，无需进行加密测量，即可确保钻孔精准恢复。
32.参考图7，在测斜点间隔为50m时，如果a以深为资源层段，在该段进行孔斜加密测量后，采用传统折线法计算的钻孔轨迹则由两段构成。显而易见，资源段的钻孔轨迹形态更接近于真实情况，但无法消除a点以上段的大间隔所带来的误差。如果b点以深为资源层段，对该段进行加密测斜，得出的结论与a点相同。
33.参考图8，在测斜点间隔为20m时，如果a以深为资源层段，在该段进行孔斜加密测量后，即加密至10m，采用传统折线法计算的钻孔轨迹同样由两段构成。资源段的钻孔轨迹形态与真实情况更接近，但并无本质的提高。如果b点以深为资源层段，对该段进行加密测斜，得出的结论与a点相同。此图表明，当采样间隔达到一定程度时加密采集数据并不能给钻孔轨迹恢复提高精度带来明显效果。
34.参考图9，实际资源钻孔的孔斜程度较理论模型明显简单，则本方法的计算精度更容易保证，常规实际孔斜采样间隔50m基本已经满足要求；为了保险起见，建议实际开展钻孔孔斜测量时，采样间隔选择20m或25m即可。
35.参考图10，当实际测井时，孔斜测量数据没有提供井口处的数据，即缺少初始值，方位角和天顶角的初始值根据图6所示方法确定。。
36.参考图11，直观示意出第i段的各种投影与增量的关系。
37.本技术对数据采样间隔是否等间距没有要求，只要平均采样间隔围绕着20～25m即可，且没有特别大的空白即可，特别大的空白指2～3倍平均间隔以上。
38.本具体实施方式的实施例均为本技术的较佳实施例，并非依此限制本技术的保护范围，故：凡依本技术的结构、形状、原理所做的等效变化，均应涵盖于本技术的保护范围之内。