CN114662257A - Trenchless pipeline path curvature calculation method and planning method thereof - Google Patents

Trenchless pipeline path curvature calculation method and planning method thereof Download PDF

Info

Publication number
CN114662257A
CN114662257A CN202210469258.0A CN202210469258A CN114662257A CN 114662257 A CN114662257 A CN 114662257A CN 202210469258 A CN202210469258 A CN 202210469258A CN 114662257 A CN114662257 A CN 114662257A
Authority
CN
China
Prior art keywords
drill
inclination angle
drill rod
drill bit
data point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210469258.0A
Other languages
Chinese (zh)
Inventor
张志余
金宇阳
金宇鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HUANGSHAN GOLDEN LAND ELECTRONICS Inc
Original Assignee
HUANGSHAN GOLDEN LAND ELECTRONICS Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by HUANGSHAN GOLDEN LAND ELECTRONICS Inc filed Critical HUANGSHAN GOLDEN LAND ELECTRONICS Inc
Priority to CN202210469258.0A priority Critical patent/CN114662257A/en
Publication of CN114662257A publication Critical patent/CN114662257A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Earth Drilling (AREA)

Abstract

The invention discloses a method for calculating the curvature of a non-excavation pipeline path and a planning method thereof, wherein the method for calculating the curvature of the non-excavation pipeline path comprises the following steps: s1) estimating the position of the drill bit each time a drill rod is added: 1.1) recording the ground position O directly above the drill2The original position (J, W, G, Q); 1.2) then converting the four parameters into three-dimensional coordinate values: 1.2.1) the drill generates a data point for each advance of a drill rod, the relative height of the drill at the current data point being Hi=Hi‑1+d*sin(Qi) Wherein i is the number of the drill rod, the advancing length of the drill rod is d, H0Setting the first drill rod of the drill bit as i-1; 1.2) converting the three-dimensional coordinates of the drill bit: set the GPS coordinate height of the origin O to G0Converting the GPS standard coordinate and the relative height of the drill bit into a new GPS coordinate to obtain the coordinate (J) of any data pointi,Wi,G0+Hi) (ii) a 1.3) convert the new GPS coordinates for each data point to geocentric coordinates (x)i,yi,zi) (ii) a S2) calculating the inclination angle change of the geocentric coordinates of each drill rod.

Description

一种非开挖管道路径弯曲度计算方法及其规划方法A method for calculating the curvature of trenchless pipeline path and its planning method

技术领域technical field

本发明涉及非开挖技术领域,尤其是一种非开挖管道路径弯曲度计算方法及非开挖管道路径规划方法。The invention relates to the technical field of trenchless, in particular to a method for calculating the curvature of a trenchless pipeline path and a method for planning a trenchless pipeline path.

背景技术Background technique

随着城市建设的大规模发展,需要在城市中铺设截污管或能源(液化气、天然气等)供应管,较常用的方法是开挖槽来埋管埋线,这会造成环境污染,引起交通堵塞,并且存在施工安全隐患。With the large-scale development of urban construction, it is necessary to lay sewage interception pipes or energy supply pipes (liquefied gas, natural gas, etc.) Traffic jams and construction safety hazards exist.

因此,目前也已开发使用了非开挖铺管技术,即一种利用岩土钻掘手段,在路面不挖沟、不破坏大面积地表层的情况下,铺设、修复和更换地下管线的施工技术。使用非开挖技术具有周期短、成本低、污染少、安全性能好等优点,而且不会影响正常的交通秩序。Therefore, trenchless pipe laying technology has also been developed and used, that is, a construction method of laying, repairing and replacing underground pipelines without digging trenches on the road surface and without damaging the surface layer of a large area by means of geotechnical drilling and excavation. technology. The use of trenchless technology has the advantages of short cycle, low cost, less pollution, good safety performance, etc., and will not affect the normal traffic order.

非开挖铺管技术应用较广的为水平导向前进法,其是利用非开挖导向仪引导装有钻头的钻杆进行定向前进来实现。非开挖导向仪包括提供钻头实时的工况-深度、倾角以及钟点方向,让地面的操作人员实时掌握钻孔轨迹以便对后续的操作进行及时的修正,以保证按既定的路线轨迹精确定向,完成非开挖铺管。由此可见,非开挖铺管技术对于非开挖导向仪的精确测量有着很高的要求。The widely used trenchless pipe laying technology is the horizontal steered advance method, which is realized by using a trenchless steerer to guide the drill pipe equipped with the drill bit for directional advancement. The trenchless steering instrument includes providing real-time working conditions of the drill bit - depth, inclination and clock direction, allowing operators on the ground to grasp the drilling trajectory in real time so as to make timely corrections to subsequent operations to ensure accurate orientation according to the established route and trajectory. Complete trenchless pipe laying. It can be seen that the trenchless pipe laying technology has high requirements for the precise measurement of the trenchless guide.

非开挖施工过程中一项衡量施工质量的重要指标是要看管道铺设的路径是否在工程设计允许的弯曲度之内。弯曲度若不在设计允许的范围内,可能会影响管道承受负载后的正常工作、使用安全、及使用寿命,严重时可能会造成重大事故。因此,施工过程中需要有必要的手段和措施来检查并监督施工路线的弯曲度。In the process of trenchless construction, an important indicator to measure the construction quality is to see whether the path of the pipeline laying is within the allowable curvature of the engineering design. If the bending degree is not within the allowable range of the design, it may affect the normal operation, safety and service life of the pipeline after being subjected to load, and may cause major accidents in serious cases. Therefore, necessary means and measures are needed to check and supervise the curvature of the construction route during the construction process.

目前所使用的方法是检查路径沿管道的倾角变化,而倾角是通过重力加速传感器检测的。值得注意的是这里的倾角是管道在垂直面与地平面的角度。也就是说,倾角只反映管道在垂直面的弯度变化;而管道的左右弯曲变化并没有体现出来。尽管多数施工是选择转大弯的设计方法,以避免管道承受急弯的窘境,但因地下情况复杂,施工过程没法保证管道始终符合“大弯”的要求;最令人担心的是施工过程中没有任何工具和方法及时检测并发现弯度的变化情况。因此,施工质量的好坏完全寄希望于施工团队的施工经验和行业操守。The method currently used is to examine the change in the inclination of the path along the pipe, which is detected by a gravitational acceleration sensor. It is worth noting that the inclination here is the angle between the vertical plane of the pipe and the ground plane. That is to say, the inclination angle only reflects the camber change of the pipeline in the vertical plane; the left and right bending changes of the pipeline are not reflected. Although most of the construction chooses the design method of turning a big bend to avoid the dilemma of the pipeline being subjected to a sharp bend, but due to the complicated underground conditions, the construction process cannot guarantee that the pipeline will always meet the requirements of "big bend"; the most worrying thing is that during the construction process There are no tools and methods to detect and detect changes in camber in time. Therefore, the quality of construction depends entirely on the construction experience and industry ethics of the construction team.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的第一个技术问题是针对上述现有技术存在的不足,提供一种非开挖管道路径弯曲度计算方法,确保施工质量。The first technical problem to be solved by the present invention is to provide a method for calculating the curvature of a trenchless pipeline path to ensure the construction quality, aiming at the shortcomings of the above-mentioned prior art.

本发明所要解决的第二个技术问题是提供一种应用上述计算方法的非开挖管道路径规划方法。The second technical problem to be solved by the present invention is to provide a trenchless pipeline path planning method using the above calculation method.

本发明解决上述第一个技术问题所采用的技术方案为:一种非开挖管道路径弯曲度计算方法,所述方法应用的非开挖钻探系统包括具备GPS数据采集功能的导向仪、钻杆和钻头,其特征在于:所述方法包括如下步骤:The technical solution adopted by the present invention to solve the above-mentioned first technical problem is as follows: a method for calculating the curvature of a trenchless pipeline path. And drill bit, it is characterized in that: described method comprises the steps:

S1)推算每次增加钻杆时钻头的位置:S1) Calculate the position of the drill bit each time the drill pipe is added:

1.1)记录钻头正上方地面位置O2的原始位置(J,W,G,Q),上述四个参数依次为经度J、纬度W、高度G以及倾角Q,倾角Q是钻头和地平面的夹角;1.1) Record the original position (J, W, G, Q) of the ground position O 2 directly above the drill bit. The above four parameters are the longitude J, latitude W, height G and inclination angle Q in turn. The inclination angle Q is the clamp between the drill bit and the ground plane. horn;

1.2)然后将四个参数转换成三维坐标值,三维坐标以地心为坐标原点:1.2) Then convert the four parameters into three-dimensional coordinate values, and the three-dimensional coordinates take the center of the earth as the coordinate origin:

1.2.1)设钻头的起点O的相对高度H0为0,钻杆的推进长度为d,每推进一个钻杆则钻头产生一个数据点,当前数据点的钻头的相对高度为Hi=Hi-1+d*sin(Qi),其中i为钻杆的编号,设置钻头的第一个钻杆记为i=1,每增加一个钻杆,i+1;1.2.1) Set the relative height H 0 of the starting point O of the drill bit to be 0, the advancing length of the drill pipe to be d, and the drill bit will generate a data point every time a drill pipe is advanced, and the relative height of the drill bit of the current data point is H i =H i-1 +d*sin(Q i ), where i is the number of the drill pipe, the first drill pipe of the drill bit is marked as i=1, and i+1 is added for each additional drill pipe;

1.2)换算钻头的三维坐标:1.2) Convert the three-dimensional coordinates of the drill bit:

设起点O的GPS坐标高度为G0,将GPS标准坐标和任意一个钻杆的相对高度换算成新的GPS坐标,得到任意一个数据点的坐标为(Ji,Wi,G0+Hi);Let the GPS coordinate height of the starting point O be G 0 , convert the GPS standard coordinates and the relative height of any drill pipe into new GPS coordinates, and obtain the coordinates of any data point as (J i , Wi , G 0 +H i ) );

1.3)将每一个数据点的新的GPS坐标换算成地心坐标(xi,yi,zi);1.3) Convert the new GPS coordinates of each data point into geocentric coordinates (x i , y i , z i );

S2)计算每一个数据点的地心坐标的倾角变化。S2) Calculate the inclination change of the geocentric coordinates of each data point.

根据本发明的一个方面,在步骤S2)中,计算倾角变化包括如下步骤:According to one aspect of the present invention, in step S2), calculating the change of inclination includes the following steps:

2.1)计算每一个数据点在地心坐标内三个平面的倾角:2.1) Calculate the inclination of each data point in three planes in geocentric coordinates:

X平面倾角Qxi=atan(|(yi-yi-1)/(zi-zi-1)|)X-plane inclination angle Qx i =atan(|(y i -y i-1 )/(z i -z i-1 )|)

Y平面倾角Qyi=atan(|(xi-xi-1)/(zi-zi-1)|)Y plane inclination angle Qy i =atan(|(x i -x i-1 )/(z i -z i-1 )|)

Z平面倾角Qzi=atan(|(yi-yi-1)/(xi-xi-1)|)Z plane inclination angle Qz i =atan(|(y i -y i-1 )/(x i -x i-1 )|)

其中,起点O的Qx0,Qy,Qz0为0;Among them, Qx 0 , Qy and Qz 0 of the starting point O are 0;

2.2)计算每一个数据点相对于前一个数据点在三个平面的倾角变化:2.2) Calculate the change in inclination of each data point relative to the previous data point in three planes:

X平面倾角变化qxi=|Qxi-Qxi-1|X-plane inclination change qx i =|Qx i -Qx i-1 |

Y平面倾角变化qyi=|Qyi-Qyi-1|Y plane inclination change qy i =|Qy i -Qy i-1 |

Z平面倾角变化qzi=|Qzi-Qzi-1|Z-plane inclination change qz i =|Qz i -Qz i-1 |

通过倾角变化得到表征管道弯曲度的钻杆的弯曲度。The bending degree of the drill pipe, which characterizes the bending degree of the pipeline, is obtained through the change of the inclination angle.

为排除异常数据,当倾角计算出现无穷大异常时,将倾角按90°记录。In order to exclude abnormal data, when there is an infinite anomaly in the calculation of the inclination angle, the inclination angle is recorded as 90°.

根据本发明的另一个方面,在步骤S2)中,计算倾角变化包括如下步骤:According to another aspect of the present invention, in step S2), calculating the inclination angle change comprises the following steps:

2.1’)根据每相连的三个数据点的地心坐标解析出一个三维空间的圆,并求出该圆的半径R,如果R不为无穷大,则相邻两个钻杆之间的倾角变化的弧度是2d/R,其中d是每一根钻杆的推进长度;如果R为无穷大,则倾角变化为0;2.1') According to the geocentric coordinates of each connected three data points, a circle in three-dimensional space is analyzed, and the radius R of the circle is obtained. If R is not infinite, the inclination angle between two adjacent drill pipes changes. The radian of is 2d/R, where d is the advancing length of each drill pipe; if R is infinite, the change in inclination is 0;

2.2’)将弧度换算成角度。2.2') Convert radians to degrees.

本发明解决上述第二个技术问题所采用的技术方案为:一种非开挖管道路径规划方法,其特征在于:所述规划方法包括如下步骤:The technical solution adopted by the present invention to solve the above-mentioned second technical problem is: a method for planning a trenchless pipeline path, characterized in that: the planning method includes the following steps:

首先,通过如上所述的非开挖管道路径弯曲度计算方法得到管道路径弯曲度;First, obtain the curvature of the pipeline path through the method for calculating the curvature of the trenchless pipeline path as described above;

然后判断是否符合工程设计标准:把每个钻杆的倾角变化和工程要求的指标比较,如果倾角变化小于工程指标,则施工满足工程要求,否则不满足工程要求。Then judge whether it meets the engineering design standards: compare the change of inclination angle of each drill pipe with the index required by the project. If the change of the inclination angle is less than the project index, the construction meets the project requirements, otherwise it does not meet the project requirements.

为确保路径满足要求,由此,如果不满足工程要求,则将当前推进的钻杆回退。To ensure that the path meets the requirements, therefore, if the engineering requirements are not met, the currently advancing drill pipe is retracted.

与现有技术相比,本发明的优点在于:通过使用GPS设备,精准计算各个方向弯度变化,从而确保施工质量。Compared with the prior art, the present invention has the advantages that: by using GPS equipment, the curvature changes in all directions can be accurately calculated, thereby ensuring the construction quality.

附图说明Description of drawings

图1为本发明实施例所使用的钻探系统钻探路径示意图。FIG. 1 is a schematic diagram of a drilling path of a drilling system used in an embodiment of the present invention.

具体实施方式Detailed ways

下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。The following describes in detail the embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein the same or similar reference numerals refer to the same or similar elements or elements having the same or similar functions.

在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,由于本发明所公开的实施例可以按照不同的方向设置,所以这些表示方向的术语只是作为说明而不应视作为限制,比如“上”、“下”并不一定被限定为与重力方向相反或一致的方向。此外,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", " Rear, Left, Right, Vertical, Horizontal, Top, Bottom, Inner, Outer, Clockwise, Counterclockwise, Axial, The orientations or positional relationships indicated by "radial direction", "circumferential direction", etc. are based on the orientations or positional relationships shown in the accompanying drawings, which are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying the indicated devices or elements. must have a specific orientation, be constructed and operated in a specific orientation, and since the disclosed embodiments of the present invention can be arranged in different orientations, these directional terms are only for illustration and should not be regarded as limiting, such as "up", "Down" is not necessarily defined as a direction opposite or in line with the direction of gravity. Furthermore, features delimited with "first", "second" may expressly or implicitly include one or more of that feature.

实施例一Example 1

本发明的工作基础是精准GPS位置信息(其所应用的非开挖钻探系统包括具备GPS数据采集功能的导向仪),施工步骤和方法可参考本申请人的在先申请,申请号为202010609895.4的中国专利所公开的航向计算方法,以及申请号为202010251440.X的中国专利所公开的导向装置及路线规划方法。The working basis of the present invention is precise GPS location information (the trenchless drilling system to which it is applied includes a directional instrument with GPS data acquisition function), and the construction steps and methods can refer to the applicant's previous application with the application number of 202010609895.4 The heading calculation method disclosed in the Chinese patent, the guiding device and the route planning method disclosed in the Chinese patent with the application number of 202010251440.X.

设想钻探系统的钻头行进的轨迹在一个三维坐标系内,如果知道每次增加钻杆时钻头的坐标位置(钻头位于行进轨迹最前方的钻杆的前端部,为现有技术),便可知道每根钻杆在三维空间内的姿态以及钻杆之间的姿态变化。姿态变化所造成的线路弯曲即是本发明所关心的内容。It is assumed that the trajectory of the drill bit of the drilling system is in a three-dimensional coordinate system. If the coordinate position of the drill bit is known each time the drill pipe is added (the drill bit is located at the front end of the drill pipe at the front of the travel trajectory, which is the prior art), we can know The attitude of each drill pipe in three-dimensional space and the attitude change between drill pipes. The line bending caused by the attitude change is what the present invention is concerned with.

参见图1,显示了一种非开挖钻探系统的钻探路径(即钻头经过的路径,铺入管道时,即为管道路径)示意图,其中,O为钻探的起始点,O1为钻探系统的钻头的当前位置,O2为钻头正上方地面的位置。在本实施例中,钻探系统的钻杆的弯度通过如下步骤得到:Referring to FIG. 1, a schematic diagram of the drilling path of a trenchless drilling system (that is, the path passed by the drill bit, when the pipe is laid, is the pipeline path) is shown, wherein O is the starting point of drilling, and O 1 is the drilling system. The current position of the drill, O2 is the position of the ground directly above the drill. In this embodiment, the camber of the drill pipe of the drilling system is obtained through the following steps:

S1)根据GPS位置及重力加速传感器数据推算每次增加钻杆时钻头的位置:S1) Calculate the position of the drill bit every time the drill pipe is added according to the GPS position and the data of the gravity acceleration sensor:

1.1)钻头正上方地面位置O2原始记录的位置为(J,W,G,Q),上述四个参数分别为经度(J)、纬度(w)、高度(G)以及倾角(Q),这是来自非开挖的钻探系统的导向仪的原始数据(操作非开挖钻探系统的导向手所在的位置),其中,倾角Q是钻头和地平面的夹角;1.1) The original recorded position of the ground position O 2 directly above the drill bit is (J, W, G, Q), and the above four parameters are longitude (J), latitude (w), altitude (G) and inclination (Q), respectively. This is the raw data from the steerer of the trenchless drilling system (where the pilot hand operating the trenchless drilling system is located), where the dip angle Q is the angle between the drill bit and the ground plane;

1.2)然后是将这四个参数变成三维坐标(地心为坐标原点)的步骤:1.2) Then there are the steps of turning these four parameters into three-dimensional coordinates (the center of the earth is the origin of the coordinates):

1.2.1)按倾角Q计算每一个钻杆相对起点O的高度H:1.2.1) Calculate the height H of each drill pipe relative to the starting point O according to the inclination angle Q:

设钻头的起点O(钻头入土之前的位置)的相对高度H0为0,钻杆的推进长度为d(如钻杆完整推进,d为钻杆的长度,此时d为固定常数),每推进一个钻杆后钻头产生一个数据点,当前数据点的钻头的相对高度则为Hi=Hi-1+d*sin(Qi),其中i为钻杆的编号,以设置钻头的第一个钻杆开始计数,记为1,每增加一个钻杆,i+1;The relative height H 0 of the starting point O of the drill bit (the position before the drill bit goes into the soil) is set to 0, and the advancing length of the drill pipe is d (if the drill pipe is fully advanced, d is the length of the drill pipe, at this time d is a fixed constant), each time After advancing a drill pipe, the drill bit generates a data point, and the relative height of the drill bit at the current data point is H i =H i-1 +d*sin(Q i ), where i is the number of the drill pipe to set the number of the drill bit. A drill pipe starts to count, recorded as 1, and i+1 for each additional drill pipe;

这里假定钻头入土(朝下)的倾角为负,出土(朝上)的角度为正;Here, it is assumed that the inclination angle of the drill bit into the soil (downward) is negative, and the angle of the unearthed (upward) is positive;

1.2.2)换算钻头的三维坐标:1.2.2) Convert the three-dimensional coordinates of the drill bit:

设起点O的GPS坐标高度为G0,将GPS标准坐标和上一步计算的相对高度变换成新的GPS坐标,得到钻头所经路线上任意一个数据点的坐标为(Ji,Wi,G0+Hi),注意:此处原始GPS高度Gi被G0+Hi替代;Set the GPS coordinate height of the starting point O as G 0 , transform the GPS standard coordinates and the relative height calculated in the previous step into new GPS coordinates, and obtain the coordinates of any data point on the route passed by the drill bit as (J i , W i , G 0 +H i ), note: the original GPS altitude G i is replaced by G 0 +H i here;

1.2.3)将每一个数据点的新的GPS坐标换算成地心坐标,即将(Ji,Wi,G0+Hi)变换成(xi,yi,zi),具体转换算法及公式可参考公开的大地物理知识;1.2.3) Convert the new GPS coordinates of each data point into geocentric coordinates, that is, transform (J i , Wi , G 0 +H i ) into ( xi , y i , z i ) , the specific conversion algorithm And the formula can refer to the public geophysics knowledge;

S2)按下面步骤计算相邻两个数据点的倾角变化,即相邻两个钻杆的倾角变化,从而得到钻杆的弯曲度:S2) Calculate the inclination angle change of two adjacent data points according to the following steps, that is, the inclination angle change of two adjacent drill pipes, so as to obtain the curvature of the drill pipe:

2.1)使用下面的公式计算钻头的每一个数据点在地心坐标内三个平面(两两互相垂直的X、Y和Z三个平面)的倾角:2.1) Use the following formula to calculate the inclination of each data point of the drill bit in three planes (X, Y, and Z planes that are perpendicular to each other) in geocentric coordinates:

X平面倾角Qxi=atan(|(yi-yi-1)/(zi-zi-1)|)X-plane inclination angle Qx i =atan(|(y i -y i-1 )/(z i -z i-1 )|)

Y平面倾角Qyi=atan(|(xi-xi-1)/(zi-zi-1)|)Y plane inclination angle Qy i =atan(|(x i -x i-1 )/(z i -z i-1 )|)

Z平面倾角Qzi=atan(|(yi-yi-1)/(xi-xi-1)|)Z plane inclination angle Qz i =atan(|(y i -y i-1 )/(x i -x i-1 )|)

其中,起点O的Qx0,Qy,Qz0初始化为0;Among them, Qx 0 , Qy, Qz 0 of the starting point O are initialized to 0;

2.2)计算钻头的每一个数据点相对于前一个数据点在三个平面的倾角变化:2.2) Calculate the change in inclination of each data point of the drill bit relative to the previous data point in three planes:

X平面倾角变化qxi=|Qxi-Qxi-1|X-plane inclination change qx i =|Qx i -Qx i-1 |

Y平面倾角变化qyi=|Qyi-Qyi-1|Y plane inclination change qy i =|Qy i -Qy i-1 |

Z平面倾角变化qzi=|Qzi-Qzi-1|Z-plane inclination change qz i =|Qz i -Qz i-1 |

在上述步骤2.1)中,计算倾角时可能会出现异常现象,比如,zi-zi-1等于零,导致出现零除异常,即倾角计算出现无穷大异常时,将倾角按90°记录。In the above step 2.1), anomalies may occur when calculating the inclination angle, for example, zi -z i -1 is equal to zero, resulting in a zero division anomaly, that is, when an infinite anomaly occurs in the inclination angle calculation, the inclination angle is recorded as 90°.

施工时,钻杆的弯曲度(即管道路径弯曲度)可能就是每一根钻杆倾角的变化,因此,只要取倾角的正弦并换算成百分比即可。比如,施工的钻杆长度为10米,工程要求每一钻杆倾角的变化百分比不得超过10%,这便意味着倾角不得超过6度。During construction, the curvature of the drill pipe (that is, the curvature of the pipeline path) may be the change of the inclination of each drill pipe. Therefore, it is only necessary to take the sine of the inclination and convert it into a percentage. For example, the length of the drill pipe for construction is 10 meters, and the engineering requires that the percentage change of the inclination angle of each drill pipe shall not exceed 10%, which means that the inclination angle shall not exceed 6 degrees.

弯曲度定义为单位长度内的弦高,其计算方法是:(d/q)*(1-COS(q/2))。其中d是钻杆的长度,q是倾角变化(单位是弧度,取qx,qy,qz之中最大的即可)。The curvature is defined as the chord height per unit length and is calculated as: (d/q)*(1-COS(q/2)). Among them, d is the length of the drill pipe, and q is the change of the inclination angle (the unit is radian, and the largest of qx, qy, and qz can be taken).

具体实施时,管道路径的规划方法,在弯曲度计算结果的基础上,判断是否符合工程设计标准:把上面的计算得来的每一个钻杆的倾角变化和工程要求的指标比较,如果计算得来的倾角变化小于工程指标,则施工满足工程要求,否则不满足工程要求。In the specific implementation, the planning method of the pipeline path, on the basis of the calculation result of the bending degree, judges whether it meets the engineering design standard: compare the inclination angle change of each drill pipe calculated above with the index required by the project, if the calculated result is If the change of the inclination angle is less than the engineering index, the construction meets the engineering requirements, otherwise it does not meet the engineering requirements.

钻头轨迹的弯曲度计算及监控是在导向记录的基础上完成的,因此弯曲度具体计算要求有连续完整的导向记录,即至少需要在每一次增加钻杆时记录一次当前钻头所处的位置及倾角,即钻头正上方地面的GPS位置和地下钻头的倾角。在工程开始之前需要记录起始点的位置及钻头的初始倾角。每完成一次导向记录需要激发一次弯曲度计算。如果弯曲度不在工程允许的范围内,需要警示导向手采取必要的措施,如回退,以保证钻头轨迹的弯曲度在施工允许的范围之内。The curvature calculation and monitoring of the drill bit trajectory is completed on the basis of the steering record, so the specific calculation of the curvature requires a continuous and complete steering record, that is, at least one record of the current drill bit position and The inclination angle is the GPS position of the ground directly above the drill bit and the inclination angle of the underground drill bit. The position of the starting point and the initial inclination of the drill should be recorded before the project starts. Every time a guide recording is completed, a curvature calculation is required to be activated. If the curvature is not within the allowable range of the project, it is necessary to warn the guide hand to take necessary measures, such as retreating, to ensure that the curvature of the drill bit track is within the allowable range of the construction.

实施例二Embodiment 2

在本实施例中,与上述实施例一的不同之处在于,步骤S2)中,包括如下步骤:In this embodiment, the difference from the above-mentioned first embodiment is that in step S2), the following steps are included:

2.1’)根据每相连的三个数据点解析出一个三维空间的圆来,并求出该圆的半径R,如果R不是无穷大,则相邻两个钻杆之间的倾角变化应是2d/R(弧度),其中d是每一根钻杆的钻进长度(通常为钻杆的长度,如果推进的总长度不是两个完整的钻杆,使用两次推进的总长度替代即可);如果R是无穷大,倾角变化则为0;从三点解析圆的半径可参考相关线性代数教课书;2.1') According to each connected three data points, a circle in three-dimensional space is analyzed, and the radius R of the circle is obtained. If R is not infinite, the change of inclination angle between two adjacent drill pipes should be 2d/ R (radian), where d is the drilling length of each drill pipe (usually the length of the drill pipe, if the total length of the advancing is not two complete drill pipes, the total length of the two advancing can be used instead); If R is infinite, the change of inclination angle is 0; the radius of the circle from three points can be analyzed by referring to the relevant linear algebra textbooks;

2.2’)将弧度换算成角度,弧度和角度之间的变换算法可参考公开的三角几何资料。2.2') Convert radian into angle, and the conversion algorithm between radian and angle can refer to the public triangular geometry data.

Claims (6)

1. A trenchless drilling system applied to a trenchless pipeline path curvature calculation method comprises a guide instrument with a GPS data acquisition function, a drill rod and a drill bit, and is characterized in that: the method comprises the following steps:
s1) estimating the position of the drill bit each time a drill rod is added:
1.1) recording the ground position O directly above the drill2The four parameters are longitude J, latitude W, height G and inclination angle Q in turn, the inclination angle Q is the included angle between the drill bit and the ground plane;
1.2) then converting the four parameters into three-dimensional coordinate values, wherein the three-dimensional coordinate takes the geocenter as a coordinate origin:
1.2.1) setting the relative height H of the starting point O of the drill00, the advancing length of the drill rod is d, the drill bit generates a data point every time one drill rod is advanced, and the relative height of the drill bit of the current data point is Hi=Hi-1+d*sin(Qi) Wherein i is the number of the drill rods, the first drill rod of the drill bit is set to be 1, and i +1 is added every time one drill rod is added;
1.2) converting three-dimensional coordinates of the drill bit:
set the GPS coordinate height of the origin O to G0Converting the GPS standard coordinate and the relative height of any drill rod into a new GPS coordinate to obtain the coordinate (J) of any drill rodi,Wi,G0+Hi);
1.3) convert the new GPS coordinates for each data point to geocentric coordinates (x)i,yi,zi);
S2) calculating a change in inclination of the geocentric coordinates of each data point.
2. The trenchless pipe path tortuosity calculation method of claim 1, wherein: in step S2), calculating the inclination change includes the steps of:
2.1) calculating the inclination angle of each data point in three planes in the geocentric coordinate:
x plane inclination angle Qxi=atan(|(yi-yi-1)/(zi-zi-1)|)
Y plane dip angle Qyi=atan(|(xi-xi-1)/(zi-zi-1)|)
Z plane inclination angle Qzi=atan(|(yi-yi-1)/(xi-xi-1)|)
Wherein Qx of the starting point O0,Qy,Qz0Is 0;
2.2) calculating the change of inclination angle of each data point relative to the previous data point in three planes:
variation of X plane inclination qxi=|Qxi-Qxi-1|
Variation of tilt angle in Y plane qyi=|Qyi-Qyi-1|
Variation of Z plane inclination qzi=|Qzi-Qzi-1|
And obtaining the bending degree of the drill rod for representing the bending degree of the pipeline through the inclination angle change.
3. The method for calculating tortuosity of a trenchless pipeline path of claim 2 wherein: when infinite abnormality occurs in the calculation of the inclination angle, the inclination angle is recorded according to 90 degrees.
4. The trenchless pipe path tortuosity calculation method of claim 1, wherein: in step S2), calculating the inclination change includes the steps of:
2.1') resolving a three-dimensional circle according to the geocentric coordinates of each connected three data points, and calculating the radius R of the circle, wherein if the R is not infinite, the radian of the change of the inclination angle between two adjacent drill rods is 2d/R, wherein d is the propelling length of each drill rod; if R is infinite, the inclination angle is changed to 0;
2.2') convert radians to angles.
5. A trenchless pipeline path planning method is characterized in that: the planning method comprises the following steps:
firstly, obtaining the curvature of the pipeline path by the calculation method of the curvature of the trenchless pipeline path according to any claim 1 to 4;
then judging whether the engineering design standard is met: and comparing the inclination angle change of each drill rod with the indexes of engineering requirements, wherein if the inclination angle change is smaller than the engineering indexes, the construction meets the engineering requirements, and otherwise, the construction does not meet the engineering requirements.
6. The trenchless pipeline path planning method of claim 5, wherein: and if the engineering requirements are not met, retracting the currently pushed drill rod.
CN202210469258.0A 2022-04-28 2022-04-28 Trenchless pipeline path curvature calculation method and planning method thereof Pending CN114662257A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210469258.0A CN114662257A (en) 2022-04-28 2022-04-28 Trenchless pipeline path curvature calculation method and planning method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210469258.0A CN114662257A (en) 2022-04-28 2022-04-28 Trenchless pipeline path curvature calculation method and planning method thereof

Publications (1)

Publication Number Publication Date
CN114662257A true CN114662257A (en) 2022-06-24

Family

ID=82037222

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210469258.0A Pending CN114662257A (en) 2022-04-28 2022-04-28 Trenchless pipeline path curvature calculation method and planning method thereof

Country Status (1)

Country Link
CN (1) CN114662257A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6389360B1 (en) * 1999-01-13 2002-05-14 Vermeer Manufacturing Company Automated bore planning method and apparatus for horizontal directional drilling
CN111411893A (en) * 2020-04-01 2020-07-14 黄山金地电子有限公司 Guiding device of drilling system and drill rod route planning method applying same
CN111411939A (en) * 2020-04-01 2020-07-14 黄山金地电子有限公司 Method for calculating drill bit depth of trenchless drilling system
CN113756721A (en) * 2020-05-29 2021-12-07 宁波金地电子有限公司 Method for eliminating inclination angle accumulation error of drilling system
CN113863856A (en) * 2020-06-29 2021-12-31 黄山金地电子有限公司 Drilling route ground course calculation method of drilling system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6389360B1 (en) * 1999-01-13 2002-05-14 Vermeer Manufacturing Company Automated bore planning method and apparatus for horizontal directional drilling
CN111411893A (en) * 2020-04-01 2020-07-14 黄山金地电子有限公司 Guiding device of drilling system and drill rod route planning method applying same
CN111411939A (en) * 2020-04-01 2020-07-14 黄山金地电子有限公司 Method for calculating drill bit depth of trenchless drilling system
CN113756721A (en) * 2020-05-29 2021-12-07 宁波金地电子有限公司 Method for eliminating inclination angle accumulation error of drilling system
CN113863856A (en) * 2020-06-29 2021-12-31 黄山金地电子有限公司 Drilling route ground course calculation method of drilling system

Similar Documents

Publication Publication Date Title
CN107654852B (en) A positioning device and positioning method in a pipeline based on pipeline segment length and pipeline connector detection
US9581018B2 (en) System and method for determining soil characteristics and drilling instability during horizontal directional drilling
CN101387198A (en) Monitoring method for drilled wellbore trajectories
CN103322989A (en) Measuring device and method of dynamic real-time measurement of position and posture of pipe pushing jack
AU2021103747A4 (en) Deep Displacement Monitoring Device and Method
CN203422091U (en) Device for dynamically measuring positions and attitudes of push bench in real time
US11556680B1 (en) Test method based on test system for five-dimensional space effects of ground surface settlement caused by dual-mode shield construction
CN105737751A (en) Vertical storage tank deformation monitoring system and method
CN104251688A (en) Method for linear directional measurement by utilizing laser in pipe jacking engineering
CN111429575A (en) A three-dimensional visual monitoring method, system, equipment and storage medium
CN107607086A (en) A kind of deep foundation underground engineering combined type tilt measurement of complicated narrow space
CN103967479B (en) A kind of rotary steerable drilling enters target prediction of situation method
CN116128290A (en) Intelligent evaluation digital twin system for urban built-up area shield method construction environment risk
CN116718142A (en) Landslide deep position data fusion monitoring method based on inclinometer and GNSS
CN111414705B (en) Three-dimensional space traversable evaluation method for bump-proof design of borehole trajectory of cluster well
CN206862331U (en) Geology internal displacement three-dimension monitor system
CN106096796A (en) A kind of soft soil base sedimentation course prediction method mapped based on coordinate complex indexes
CN107169224A (en) Great diameter and long pile tube pile drivability analysis method is carried out based on CPTU tests
CN116164735A (en) A three-dimensional attitude measurement system and method for water-saving immersed tubes
CN104655191B (en) Multi-parameter and three-dimensional monitoring method and monitoring probe for reservoir bank of water-level-fluctuating zone
CN111411939B (en) Method for calculating drill bit depth of non-excavation drilling system
CN114662257A (en) Trenchless pipeline path curvature calculation method and planning method thereof
CN107219335B (en) Pipeline connector detection method based on complex continuous wavelet transform
CN104060983A (en) Wired terrestrial magnetism while drilling guide instrument and measurement method
CN209310685U (en) A high-fill soil layered settlement monitoring device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination