Method for calculating occurrence of underground structural plane based on in-hole television

Abstract

The invention belongs to the field of geological exploration, and particularly relates to a method for calculating the occurrence of an underground structural plane based on an in-hole television. After drilling work in any drilling direction and drilling angle is completed and an in-hole television result is obtained, calculation work is carried out on the basis of obtaining drilling hole information and structural plane information; according to the method, the occurrence of the underground structural plane can be calculated based on the in-hole television in any drilling direction and drilling angle, and the true tendency and the true dip angle of the structural plane in the geographic orientation can be obtained.

Description

Technical Field

[0001] The invention belongs to the field of geological exploration, and in particular relates to a method for estimating the occurrence of underground structural surfaces based on in-hole television. Background Art

[0002] The underground environmental conditions in urban built-up areas are complex, with various pipelines, integrated pipeline corridors, underground rail transit, buried ditches, caverns, underground shopping malls and other buildings (structures) densely distributed. When it is necessary to ascertain the geological conditions under existing buildings (structures) or the geological conditions in front of the tunnel excavation working surface, inclined drilling or horizontal drilling at different angles can be used for survey work.

[0003] The structural surface dip, inclination and other occurrence information calculated by in-hole television in the survey results are the interpretation results under the assumption of vertical drilling and cannot represent the occurrence characteristics of the actual structural surface.

[0004] Therefore, in order to obtain statistical information on the occurrence of structural surfaces under any drilling direction and drilling angle, a method for estimating the occurrence of underground structural surfaces is urgently needed. Summary of the invention

[0005] In order to solve the above technical problems, the present invention provides a method for estimating the occurrence of underground structural surfaces based on in-hole television. The method can realize the estimation of the occurrence of underground structural surfaces based on in-hole television under any drilling direction and drilling angle.

[0006] The present invention is achieved through the following technical solutions: A method for estimating the occurrence of underground structural surfaces based on in-hole television, the method comprising: (1) Conduct drilling work and obtain drilling information and structural surface information based on in-hole television; (2) constructing a three-dimensional coordinate system: in the three-dimensional coordinate system, the east-west direction of the geographic direction is the X-axis and the east is the positive direction, the north-south direction of the geographic direction is the Y-axis and the north is the positive direction, and the vertical direction is the Z-axis and the upward direction is the positive direction; (3) converting the two geographic azimuths of the drilling direction and the apparent inclination into coordinate angles of the three-dimensional coordinate system; (4) Construct the trajectory of the borehole mouth circle in the three-dimensional coordinate system; (5) Using the coordinate points of the hole mouth circle and the structural surface information interpreted by the in-hole television, calculate the elliptical trajectory at the structural surface; (6) Screening the coordinate point A1 at the maximum value of the Z-axis coordinate, the coordinate point A2 at the minimum value of the Z-axis coordinate, and the coordinate points A3 and A4 at the middle value of the Z-axis coordinate in the elliptical trajectory of the structural surface; (7) Calculate the direction of the structural surface using the coordinates of coordinate points A3 and A4, determine the inclination of the structural surface using the coordinates of A1 and A2, and convert the inclination of the structural surface into a geographic coordinate angle, which is the true inclination of the structural surface under the geographic orientation; (8) Use the coordinates of points A3 and A4 to construct the spatial line L, calculate the distance from point A1 to line L and the difference in Z value, and obtain the true inclination angle.

[0007] Furthermore, in step (1), the obtained drilling information includes: coordinates and elevation of the center of the drilling hole, drilling direction, inclination and drilling diameter; the obtained structural surface information includes: footage, apparent inclination and apparent inclination of the structural surface center point interpreted by in-hole television.

[0008] Furthermore, step (4) is specifically as follows: using the coordinates and elevation of the center of the borehole, the drilling direction, the inclination angle and the borehole diameter, the coordinate points (x1, y1, z1), (x2, y2, z2), (x3, y3, z3) ... (x n , y n , z n ), that is, the trajectory of the drilling hole circle in the three-dimensional coordinate system is obtained.

[0009] Further, step (5) is specifically as follows: using the coordinate point of the trajectory circle O at the borehole mouth in the three-dimensional coordinate system and the footage, apparent inclination, and apparent inclination angle of the center point of the structural surface, the coordinate point (x1 ’ , y1 ’ , z1 ’ ), (x2 ’ , y2 ’ , z2 ’ ), (x3 ’ , y3 ’ , z3 ’ )……(x n ’ , y n ’ , z n ’ ), that is, the elliptical trajectory on the structural surface is obtained.

[0010] Furthermore, in step (6), the coordinate point at the maximum Z-axis coordinate in the elliptical trajectory of the structure surface obtained by screening is recorded as A1(x max ’ , y max ’ , z max ’ ), the coordinate point of the minimum Z-axis coordinate in the elliptical trajectory is recorded as A2(x min’ , y min ’ ,z min ’ ), the coordinate point at the middle value of the Z axis coordinate in the elliptical trajectory is recorded as A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and A4(x avg2 ’ ,y avg2 ’ , z avg2 ’ ).

[0011] Furthermore, step (7) specifically includes: using coordinate point A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4(x avg2 ’ , y avg2 ’ , z avg2 ’ ) Calculate the structural surface strike α1, and then calculate the structural surface inclination or the opposite direction of the inclination α2: The calculation formula of structural surface trend α1 is: ; The calculation formula of structural surface inclination or the opposite direction of inclination α2 is: α2=α1+π / 2; Through the coordinate point A1(x max ’ , y max ’ , z max ’ ) and coordinate point A2(x min ’ , y min ’ , z min ’ ) to determine the quadrant on the XOY plane where the structural surface inclination is located; if the quadrant on the XOY plane where the structural surface inclination is located is the same as the quadrant where α2 is located, then α2 is the inclination of the structural surface in the three-dimensional coordinate system; if the quadrant on the XOY plane where the structural surface inclination is located is opposite to the quadrant where α2 is located, then α2±π is the inclination α3 of the structural surface in the three-dimensional coordinate system; that is: ; Convert the inclination of the structural surface in the three-dimensional coordinate system into the true inclination α in the geographic coordinate angle: the conversion formula is as follows: .

[0012] Furthermore, step (8) specifically includes: using coordinate point A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4(x avg2 ’ , y avg2 ’ , z avg2 ’ ) construct the space line L, find the unit direction vector (m,n,p) of the space line L, the vector (M,N,P) formed by the coordinate point A3 and the coordinate point A1, and calculate A1(x max ’ , y max ’ , z max ’ ) to the space line L: ; Calculate the coordinates of point A1(x max ’ , y max ’ , z max ’ ) and the coordinate point A3(x avg1 ’ , y avg1 ’ , z avg1 ’ ) ’ = z max ’ -z avg1 ’ , and then calculate the true inclination angle β of the structural surface: β = arcsin(z ’ / d).

[0013] Furthermore, in step (1), the drilling work refers to drilling work in any drilling direction and drilling angle.

[0014] Beneficial technical effects of the present invention: (1) The method for estimating the occurrence of underground structural surfaces based on in-hole television provided by the present invention gradually estimates the actual coordinate point positions of the structural surface according to the drilling process in a single coordinate system without the need for coordinate system transformation.

[0015] (2) The method for estimating the attitude of underground structural surfaces based on in-hole television provided by the present invention has a simple and easy-to-understand calculation process and is widely applicable.

[0016] (3) The method provided by the present invention for estimating the occurrence of underground structural surfaces based on in-hole television can be calculated in any drilling direction and drilling angle, such as vertical, inclined, and horizontal, without any restrictions on usage conditions.

[0017] (4) The method provided by the present invention for estimating the attitude of underground structural surfaces based on in-hole television screens the coordinate characteristic values ​​for estimation, and the calculation error is smaller. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 It is the apparent inclination and inclination angle of the structural surface in the hole TV; Figure 2 This is a flow chart of a method for estimating the occurrence of underground structural surfaces based on in-hole television in an embodiment of the present invention; Figure 3 Schematic diagram of a method for estimating the attitude of underground structural surfaces based on in-hole television in an embodiment of the present invention. DETAILED DESCRIPTION

[0019] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0020] On the contrary, the present invention covers any substitution, modification, equivalent method and scheme made on the essence and scope of the present invention as defined by the claims. Further, in order to make the public have a better understanding of the present invention, some specific details are described in detail in the detailed description of the present invention below. Those skilled in the art can fully understand the present invention without the description of these details.

[0021] A method for estimating the occurrence of underground structural surfaces based on in-hole television. The method provided by the present invention is to perform the estimation work based on the coordinates and elevation (m) of the center of the borehole, the drilling direction (geographic azimuth °), the inclination (°) and the borehole diameter (mm) of the borehole opening, and the structural surface information such as the footage (m) of the center point of the structural surface, the apparent inclination (geographic azimuth °), and the apparent inclination (°) interpreted by the in-hole television after the drilling work of any drilling direction and drilling angle has been completed and the in-hole television results have been obtained. By establishing a three-dimensional rectangular coordinate system, the coordinates of all trajectory points on the structural surface are gradually estimated according to the drilling process using the drilling information and the structural surface information, and after screening the characteristic values, the true inclination and true inclination of the structural surface in the geographical orientation are calculated.

[0022] like Figure 2-Figure 3As shown, the method includes: (1) Performing drilling work and obtaining borehole information and structural surface information based on in-hole television; the borehole information obtained includes: coordinates and elevation of the borehole center, drilling direction, inclination and borehole diameter; the structural surface information obtained includes: footage, apparent inclination and apparent inclination of the structural surface center point interpreted by in-hole television. The drilling work refers to drilling work in any drilling direction and drilling angle; (2) constructing a three-dimensional coordinate system: in the three-dimensional coordinate system, the east-west direction of the geographic direction is the X-axis and the east is the positive direction, the north-south direction of the geographic direction is the Y-axis and the north is the positive direction, and the vertical direction is the Z-axis and the upward direction is the positive direction; (3) converting the two geographic azimuths of the drilling direction and the apparent inclination into coordinate angles of the three-dimensional coordinate system; (4) Construct the trajectory of the borehole mouth circle in the three-dimensional coordinate system; Specifically, the coordinate points (x1, y1, z1), (x2, y2, z2), (x3, y3, z3)…(x n , y n , z n ), that is, the trajectory of the borehole opening circle in the three-dimensional coordinate system is obtained. This step can be achieved by using the space vector method, a conventional technology.

[0023] (5) Using the coordinate points of the hole mouth circle and the structural surface information interpreted by the in-hole television, calculate the elliptical trajectory at the structural surface; specifically: using the coordinate points of the trajectory circle O at the hole mouth in the three-dimensional coordinate system and the footage, apparent inclination, and apparent inclination of the center point of the structural surface, calculate the coordinate points of the trajectory ellipse O' at the structural surface (x1 ’ , y1 ’ , z1 ’ ), (x2 ’ , y2 ’ , z2 ’ ), (x3 ’ , y3 ’ ,z3 ’ )……(x n ’ , y n ’ , z n ’ ), that is, the elliptical trajectory on the structural surface is obtained.

[0024] Where: x n ’ =x1+△x1+△x2,y n ’=y1+△y1+△y2,z n ’ =z1+△z1+△z2; Where: △x1, △y1, △z1 are the x value, y value, and z value of the footage increase; △x2, △y2, △z2 are the x value, y value, and z value of the structural surface compared to the circular surface perpendicular to the drilling direction at the center of the structural surface; (6) Screening the coordinate point A1 at the maximum value of the Z-axis coordinate in the elliptical trajectory of the structural surface, the coordinate point A2 at the minimum value of the Z-axis coordinate, and the coordinate points A3 and A4 at the intermediate value of the Z-axis coordinate; Specifically, the coordinate point at the maximum value of the Z-axis coordinate in the elliptical trajectory of the structural surface obtained by screening is recorded as A1 (x max ’ , y max ’ , z max ’ ), the coordinate point of the minimum Z-axis coordinate in the elliptical trajectory is recorded as A2(x min ’ , y min ’ , z min ’ ), the coordinate point at the middle value of the Z axis coordinate in the elliptical trajectory is recorded as A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and A4(x avg2 ’ , y avg2 ’ , z avg2 ’ ).

[0025] (7) Calculate the direction of the structural surface using the coordinates of coordinate points A3 and A4, determine the inclination of the structural surface using the coordinates of A1 and A2, and convert the inclination of the structural surface into a geographic coordinate angle, which is the geographic orientation of the true inclination of the structural surface; Specifically, using the coordinate point A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4 (x avg2 ’ , y avg2 ’ , z avg2 ’ ) Calculate the structural surface strike α1, and then calculate the structural surface inclination or the opposite direction of the inclination α2: The calculation formula of structural surface trend α1 is: ; The calculation formula of structural surface inclination or the opposite direction of inclination α2 is: α2=α1+π / 2; Through the coordinate point A1(x max ’ , y max ’ , z max ’ ) and coordinate point A2(x min ’ , y min ’ , z min ’ ) to determine the quadrant on the XOY plane where the structural surface inclination is located; that is, if x min ’ -x max ’ >0 and y min ’ -y max ’ >0, it is in the first quadrant; if x min ’ -x max ’ <0 and y min ’ -y max ’ >0, it is in the second quadrant; if x min ’ -x max ’ <0 and y min ’ -y max ’ <0, it is in the third quadrant; if x min ’ -x max ’ >0 and y min ’ -y max ’ <0, it is in the fourth quadrant; If the quadrant on the XOY plane where the structural surface inclination is located is the same as the quadrant where α2 is located, then α2 is the inclination of the structural surface in the three-dimensional coordinate system; if the quadrant on the XOY plane where the structural surface inclination is located is opposite to the quadrant where α2 is located, then α2±π is the inclination α3 of the structural surface in the three-dimensional coordinate system; that is: ; Convert the inclination of the structural surface in the three-dimensional coordinate system into the true inclination α in the geographic coordinate angle: the conversion formula is as follows: .

[0026] (8) Use the coordinates of coordinate points A3 and A4 to construct the spatial straight line L, calculate the distance from coordinate point A1 to straight line L and the difference in Z value, and calculate the true inclination angle. avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4(x avg2 ’ ,y avg2 ’ , z avg2 ’ ) construct the space line L, find the unit direction vector (m,n,p) of the space line L, the vector (M,N,P) formed by the coordinate point A3 and the coordinate point A1, and calculate A1(x max ’ , y max ’ , z max ’ ) to the space line L: ; Calculate the coordinates of point A1(x max ’ , y max ’ , z max ’ ) and the coordinate point A3(x avg1 ’ , y avg1 ’ , z avg1 ’ )’s Z value difference z’= z max ’ -z avg1 ’ , and then calculate the true inclination angle β of the structural surface: β=arcsin(z' / d).

[0027] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for estimating the occurrence of underground structural surfaces based on in-hole television, characterized in that: The method comprises: (1) Conduct drilling work and obtain drilling information and structural surface information based on in-hole television; (2) constructing a three-dimensional coordinate system: in the three-dimensional coordinate system, the east-west direction of the geographic direction is the X-axis and the east is the positive direction, the north-south direction of the geographic direction is the Y-axis and the north is the positive direction, and the vertical direction is the Z-axis and the upward direction is the positive direction; (3) converting the two geographic azimuths of the drilling direction and the apparent inclination into coordinate angles of the three-dimensional coordinate system; (4) Construct the trajectory of the borehole mouth circle in the three-dimensional coordinate system; (5) Using the coordinate points of the hole mouth circle and the structural surface information interpreted by the in-hole television, calculate the elliptical trajectory at the structural surface; (6) Screening the coordinate point A1 at the maximum value of the Z-axis coordinate, the coordinate point A2 at the minimum value of the Z-axis coordinate, and the coordinate points A3 and A4 at the middle value of the Z-axis coordinate in the elliptical trajectory of the structural surface; (7) Calculate the direction of the structural surface using the coordinates of coordinate points A3 and A4, determine the inclination of the structural surface using the coordinates of A1 and A2, and convert the inclination of the structural surface into a geographic coordinate angle, which is the true inclination of the structural surface under the geographic orientation; (8) Use the coordinates of points A3 and A4 to construct the spatial line L, calculate the distance from point A1 to line L and the difference in Z value, and obtain the true inclination angle.

2. According to claim 1, a method for estimating the occurrence of underground structural surfaces based on in-hole television, characterized in that: In step (1), the obtained drilling information includes: the coordinates and elevation of the center of the drilling hole, the drilling direction, the inclination angle and the drilling diameter; the obtained structural surface information includes: the footage, apparent inclination and apparent inclination angle of the structural surface center point interpreted by in-hole television.

3. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 1, characterized in that: Step (4) is as follows: using the coordinates and elevation of the center of the borehole, the drilling direction, the inclination angle and the borehole diameter, the coordinate points (x1, y1, z1), (x2, y2, z2), (x3, y3, z3) ... (x n , y n ,z n ), that is, the trajectory of the drilling hole circle in the three-dimensional coordinate system is obtained.

4. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 3 is characterized in that: Step (5) is as follows: using the coordinate point of the trajectory circle O at the borehole mouth in the three-dimensional coordinate system and the footage, apparent inclination, and apparent inclination angle of the center point of the structural surface, calculate the coordinate point (x1 ’ , y1 ’ , z1 ’ ), (x2 ’ , y2 ’ , z2 ’ ), (x3 ’ ,y3 ’ , z3 ’ )……(x n ’ , y n ’ , z n ’ ), that is, the elliptical trajectory on the structural surface is obtained.

5. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 4, characterized in that: In step (6), the coordinate point of the maximum Z-axis coordinate in the elliptical trajectory of the structure surface obtained by screening is recorded as A1(x max ’ , y max ’ ,z max ’ ), the coordinate point of the minimum Z-axis coordinate in the elliptical trajectory is recorded as A2(x min ’ , y min ’ , z min ’ ), the coordinate point at the middle value of the Z axis coordinate in the elliptical trajectory is recorded as A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and A4(x avg2 ’ , y avg2 ’ , z avg2 ’ ).

6. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 5, characterized in that: Step (7) specifically includes: using coordinate point A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4 (x avg2 ’ , y avg2 ’ , z avg2 ’ ) Calculate the structural surface strike α1, and then calculate the structural surface inclination or the opposite direction of the inclination α2: The calculation formula of structural surface trend α1 is: ； The calculation formula of structural surface inclination or the opposite direction of inclination α2 is: α2=α1+π / 2； Through the coordinate point A1(x max ’ , y max ’ , z max ’ ) and coordinate point A2(x min ’ , y min ’ , z min ’ ) to determine the quadrant on the XOY plane where the structural surface inclination is located; If the quadrant on the XOY plane where the structural surface inclination is located is the same as the quadrant where α2 is located, then α2 is the inclination of the structural surface in the three-dimensional coordinate system; if the quadrant on the XOY plane where the structural surface inclination is located is opposite to the quadrant where α2 is located, then α2±π is the inclination α3 of the structural surface in the three-dimensional coordinate system; that is: ； Convert the inclination of the structural surface in the three-dimensional coordinate system into the true inclination α in the geographic coordinate angle: the conversion formula is as follows: 。 7. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 5, characterized in that: Step (8) specifically includes: using coordinate point A3 (x avg1 ’ , y avg1 ’ , z avg1 ’ ) and coordinate point A4(x avg2 ’ , y avg2 ’ , z avg2 ’ ) construct the space line L, find the unit direction vector (m,n,p) of the space line L, the vector (M,N,P) formed by the coordinate point A3 and the coordinate point A1, and calculate A1(x max ’ , y max ’ , z max ’ ) to the space line L: ； Calculate the coordinates of point A1(x max ’ , y max ’ , z max ’ ) and the coordinate point A3(x avg1 ’ , y avg1 ’ , z avg1 ’ ) ’ =z max ’ -z avg1 ’ , and then calculate the true inclination angle β of the structural surface: β=arcsin(z ’ / d)。 8. The method for estimating the occurrence of underground structural surfaces based on in-hole television according to claim 1, characterized in that: In step (1), the drilling work refers to drilling work in any drilling direction and drilling angle.

Images