Submarine pipeline routing area element space analysis method

Abstract

The invention discloses a submarine pipeline routing area element space analysis method, which comprises the following steps of counting the number and types of elements intersected with a submarine pipeline, counting the distribution interval of the elements of the same type along the axis path of the submarine pipeline, counting the cumulative length distribution of linear elements, and counting the area distribution of polygonal elements. And carrying out statistics on element types and element numbers in the submarine pipeline buffer area to complete submarine pipeline routing area element space analysis. According to the regional element spatial analysis method for the submarine pipeline route, the number of elements intersected with the submarine pipeline, element types, intersection positions, the distribution interval of the elements of the same type along the submarine pipeline, the length distribution condition of the elements of the same type, the area distribution condition of the elements of the same type, and the element types in a submarine pipeline buffer area can be displayed; therefore, the risk analysis of the third-party operation of the submarine pipeline is facilitated.

Description

Technical Field

[0001] The present application relates to the field of submarine pipelines, and in particular to a method for spatial analysis of regional elements of submarine pipeline routing. Background Art

[0002] In the submarine pipeline integrity management, third-party damage in the submarine pipeline routing area is one of the factors that need to be considered. At present, in the submarine pipeline integrity management, the third-party operation trace analysis method is relatively simple, and it is difficult to effectively perform spatial analysis of third-party operation traces, which is reflected in the following aspects:

[0003] (1) It is difficult to analyze the intersection between the third-party operation traces in the submarine pipeline routing area and the plane position of the submarine pipeline;

[0004] (2) It is difficult to conduct statistical analysis on the length, area, and distribution interval of third-party operation traces in the submarine pipeline routing area;

[0005] (3) It is difficult to quantify the traces of third-party operations within the submarine pipeline buffer zone. Summary of the invention

[0006] The present application provides a method for spatial analysis of submarine pipeline routing regional elements, which can facilitate risk analysis of third-party operations of submarine pipelines.

[0007] An embodiment of the present application provides a method for spatial analysis of submarine pipeline routing regional elements, comprising the following steps:

[0008] Calculation of the axis equation of the submarine pipeline in geographic coordinates.

[0009] According to the submarine pipeline axis equation, the submarine pipeline axis is divided into multiple nodes P, and the mileage of the submarine pipeline axis nodes P along the submarine pipeline axis path is calculated.

[0010] According to the geographic coordinates of the element node F, the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the geographic coordinates of the foot C of the element node F on the corresponding segment of the submarine pipeline axis are calculated.

[0011] According to the geographical coordinates of the foot of the perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis and the geographical coordinates of the node P of the submarine pipeline axis, the mileage kp of the path of the foot of the perpendicular C along the submarine pipeline axis is calculated. c .

[0012] When the feature type is a linear feature, the distance between adjacent nodes and the cumulative length d of the feature are calculated based on the geographic coordinates of the feature node F. When the feature type is a polygonal feature, the area A of the feature is calculated based on the geographic coordinates of the feature node F.

[0013] The intersection relationship between the element and the submarine pipeline is determined based on the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis. The mileage kp of the path along the submarine pipeline axis is determined based on the foot of the perpendicular C. c Determine the relative position of the element in the submarine pipeline axis path and the length of the element along the submarine pipeline axis path. According to the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer distance h, analyze whether the element is in the submarine pipeline buffer.

[0014] Count the number and type of features that intersect with the submarine pipeline, count the distribution intervals of the same type of features along the axis path of the submarine pipeline, count the cumulative length distribution of linear features, count the area distribution of polygonal features, count the type and number of features in the submarine pipeline buffer zone, and complete the spatial analysis of the elements of the submarine pipeline routing area.

[0015] In some embodiments, the calculation of the axis equation of the submarine pipeline under geographic coordinates includes:

[0016] The plane position of the submarine pipeline is decomposed into straight line segments and arc segments, and the intervals are divided. The equations of the straight line segments and arc segments are calculated respectively, and then they are combined into the axis equation of the submarine pipeline.

[0017] In some embodiments, calculating the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the geographic coordinates of the foot C of the element node F on the corresponding segment of the submarine pipeline axis includes:

[0018] Calculate the vertical distance s from the element node F to the straight line segment L of the submarine pipeline axis, and the geographical coordinates of the foot C of the perpendicular to the element node F on the straight line segment L of the submarine pipeline axis. The calculation method is:

[0019] The geographic coordinates of feature node F are (x f ,y f ), the equation of the straight line segment L of the submarine pipeline axis is y=kx+b1, x∈[x1,x2], and the calculation formula of the vertical distance s from the element node F to the straight line segment L of the submarine pipeline axis is as follows:

[0020]

[0021] The geographical coordinates of the foot C of the perpendicular to the straight line segment L of the submarine pipeline axis of the element node F are (x c ,y c ), the calculation formula is as follows:

[0022]

[0023] In some embodiments, calculating the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the geographic coordinates of the foot C of the element node F on the corresponding segment of the submarine pipeline axis includes:

[0024] Calculate the vertical distance s from the element node F to the arc segment R of the submarine pipeline axis, and the geographic coordinates of the foot C of the perpendicular to the element node F on the arc segment R of the submarine pipeline axis. The calculation method is:

[0025] The geographic coordinates of feature node F are (x f ,y f ), the equation of the arc segment R of the submarine pipeline axis is (xx r ) 2 +(yy r ) 2 =r 2 , The calculation formula of the vertical distance s from the element node F to the arc segment R of the submarine pipeline axis is as follows:

[0026]

[0027] The geographical coordinates of the foot C of the perpendicular to the arc segment R of the submarine pipeline axis are (x c ,y c ), the calculation formula is as follows:

[0028]

[0029] In some embodiments, the mileage kp of the path of the foot of perpendicular C along the axis of the submarine pipeline is calculated based on the geographical coordinates of the foot of perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis and the geographical coordinates of the node P of the submarine pipeline axis. c ,include:

[0030] The geographic coordinates of the foot of the perpendicular C (x c ,y c ) and the geographic coordinates of node P (x p ,y p ), if the foot of the perpendicular C coincides with the node P, then the mileage kp of the path along the axis of the seabed pipeline is c kp is the mileage of the path from node P to the axis of the submarine pipeline p If the foot of the perpendicular C is between two nodes P, the distance Δd from the foot of the perpendicular C to the previous node P is calculated, and the mileage kp of the path of the foot of the perpendicular C along the axis of the submarine pipeline is c kp is the mileage of the previous node P along the axis of the submarine pipeline P-1 +Δd.

[0031] In some embodiments, the area A of the element is calculated according to the geographic coordinates of the element node F, and the calculation formula is as follows:

[0032]

[0033] Where n is the number of nodes of the polygon feature.

[0034] In some embodiments, judging the intersection relationship between the element and the submarine pipeline according to the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis includes:

[0035] When the vertical distances s from the element node F to the corresponding segment of the submarine pipeline axis are all positive or all negative, it indicates that the element is completely on one side of the submarine pipeline and has no intersection with the submarine pipeline.

[0036] When the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis contains positive and negative values, it indicates that the element and the submarine pipeline have an intersection relationship, and the intersection point is between the adjacent nodes F with positive and negative vertical distances s.

[0037] In some embodiments, the mileage kp of the path along the axis of the seabed pipeline according to the foot of the vertical C is c Determine the relative position of the element along the axis path of the submarine pipeline and the length of the element along the axis path of the submarine pipeline, including:

[0038] According to the mileage kp of the path along the axis of the seabed pipeline according to the foot of the vertical C c The maximum and minimum values ​​in determine the relative position of the feature in the submarine pipeline axis path and the length of the feature along the submarine pipeline axis path.

[0039] In some embodiments, analyzing whether a feature is within the submarine pipeline buffer zone according to the vertical distance s from the feature node F to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer zone distance h includes:

[0040] When the absolute value of the vertical distance s from the element node to the corresponding segment of the submarine pipeline axis is less than h, the element is determined to be within the submarine pipeline buffer zone.

[0041] When the absolute value of the vertical distance s from the element node to the corresponding segment of the submarine pipeline axis is greater than h, the element is determined to be outside the submarine pipeline buffer zone.

[0042] When the absolute values ​​of the vertical distances s from the element node to the corresponding segments of the submarine pipeline axis are partially greater than h and partially less than h, it is determined that the element is partially within the submarine pipeline buffer zone.

[0043] In some embodiments, the distribution intervals of the same type of elements along the axis path of the submarine pipeline are counted, including:

[0044] When the mileage intervals of two elements along the axis path of the submarine pipeline overlap, the union of the two elements is calculated.

[0045] According to an embodiment of the present application, a method for spatial analysis of submarine pipeline routing regional elements is provided, comprising the following steps: Calculate the equation of the submarine pipeline axis under geographic coordinates. According to the equation of the submarine pipeline axis, divide the submarine pipeline axis into multiple nodes P, and calculate the mileage of the submarine pipeline axis node P along the submarine pipeline axis path. According to the geographic coordinates of the element node F, calculate the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis, and the geographic coordinates of the foot of the perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis. According to the geographic coordinates of the foot of the perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis, calculate the mileage kp of the foot of the perpendicular C along the submarine pipeline axis path. c When the feature type is a linear feature, the distance between adjacent nodes and the cumulative length d of the feature are calculated based on the geographic coordinates of the feature node F. When the feature type is a polygonal feature, the area A of the feature is calculated based on the geographic coordinates of the feature node F. The intersection relationship between the feature and the submarine pipeline is determined based on the vertical distance s from the feature node F to the corresponding segment of the submarine pipeline axis. The mileage kp of the path along the submarine pipeline axis is calculated based on the foot of the perpendicular C. c The maximum and minimum values ​​in determine the relative position of the element in the submarine pipeline axis path and the length of the element along the submarine pipeline axis path. According to the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer distance h, analyze whether the element is in the submarine pipeline buffer. Count the number of elements and element types that intersect with the submarine pipeline, count the distribution interval of the same type of elements along the submarine pipeline axis path, count the cumulative length distribution of linear elements, count the area distribution of polygonal elements, count the element types and number of elements in the submarine pipeline buffer, and complete the spatial analysis of the elements in the submarine pipeline routing area. The submarine pipeline routing area element spatial analysis method of the present application can display the number of elements, element types, and intersection positions that intersect with the submarine pipeline, the distribution interval of the same type of elements along the submarine pipeline, the length distribution of the same type of elements, the area distribution of the same type of elements, the element types in the submarine pipeline buffer, the number of elements, etc., which is convenient for the third-party operation risk analysis of the submarine pipeline. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without paying any creative work.

[0047] Figure 1 It is a schematic diagram of the overall process of the spatial analysis method of submarine pipeline routing area elements in an embodiment of the present application.

[0048] Figure 2is a spatial distribution map of the elements of the submarine pipeline routing area in the embodiment of the present application;

[0049] Figure 3 is a distribution diagram of the elements in the direction of the path of the axis of the submarine pipeline in the embodiment of the present application;

[0050] Figure 4 is a statistical graph of the cumulative length of linear elements in the embodiment of the present application;

[0051] Figure 5 is a statistical graph of the areas of polygonal elements in the embodiment of the present application;

[0052] Figure 6 is a submarine pipeline attribute table in an embodiment of the present application;

[0053] Figure 7 is an attribute table of the linear element in the embodiment of the present application;

[0054] Figure 8 is an attribute table of polygon elements in the embodiment of the present application;

[0055] Fig. 9 is a statistical table of intersections between elements and submarine pipelines in an embodiment of the present application;

[0056] Fig.10 It is a statistical table of the number and type of features in the submarine pipeline buffer zone. DETAILED DESCRIPTION

[0057] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application 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 application and are not used to limit the present application.

[0058] See also Figure 1-5 The embodiment of the present application provides a method for spatial analysis of submarine pipeline routing regional elements, comprising the following steps:

[0059] Step 1: Calculate the axis equation of the submarine pipeline in geographic coordinates.

[0060] In the above steps, the axis of the submarine pipeline is the path, and the direction of the path is the direction of the submarine pipeline. The plane position of the submarine pipeline is decomposed into straight line segments and arc segments, and the intervals are divided. The equations of the straight line segments and arc segments are calculated respectively, and then combined into the equation of the axis of the submarine pipeline.

[0061] Step 2: According to the submarine pipeline axis equation, the submarine pipeline axis is divided into multiple nodes P, and the mileage of each node P of the submarine pipeline axis along the submarine pipeline axis path is calculated.

[0062] In the above steps, according to the equation of the submarine pipeline axis, the submarine pipeline axis is divided into multiple nodes P, and the distance Δd between adjacent nodes P is set to 1m. According to the geographic coordinates (x p ,y p ) Calculate the distance dn of each node P along the path starting point, and obtain the mileage kp of each node P along the path p , generate the attribute table of the submarine pipeline axis, including ID, feature name, feature type, geographic coordinates of node P (x p ,y p ), along the axis of the submarine pipeline path mileage kp p .

[0063] Step 3: Based on the geographic coordinates of each node F of the element, calculate the vertical distance s (called "offset distance") from each node F of the element to the corresponding segment of the submarine pipeline axis, as well as the geographic coordinates of the foot C of the perpendicular of each node F of the element on the corresponding segment of the submarine pipeline axis.

[0064] In the above steps, the information of the submarine pipeline routing area elements is first generated into an element attribute table, including ID, element name, element type, element code, geographic coordinates of node F (x f ,y f ).

[0065] Submarine pipelines generally consist of straight segments and arc segments. f ,y f ), calculate the vertical distance s from each node F of the element to the corresponding segment of the submarine pipeline axis, and the geographical coordinates of the foot C of each node F of the element on the corresponding segment of the submarine pipeline axis, including: calculating the vertical distance s from each node F of the element to the straight segment L of the submarine pipeline axis, and the geographical coordinates of the foot C of each node F of the element on the straight segment L of the submarine pipeline axis, and calculating the vertical distance s from each node F of the element to the arc segment R of the submarine pipeline axis, and the geographical coordinates of the foot C of each node F of the element on the arc segment R of the submarine pipeline axis.

[0066] Among them, the vertical distance s from each node F of the element to the straight line segment L of the submarine pipeline axis and the geographical coordinates of the foot C of each node F of the element on the straight line segment L of the submarine pipeline axis are calculated as follows:

[0067] The geographic coordinates of each node F of the element are (x f ,y f ), the equation of the straight line segment L of the submarine pipeline axis is y=kx+b1, x∈[x1,x2], and the calculation formula of the vertical distance s from each node F of the element to the straight line segment L of the submarine pipeline axis is as follows:

[0068]

[0069] The geographical coordinates of the foot C of the perpendicular to the straight line segment L of the submarine pipeline axis of each node F are (x c ,y c ), the calculation formula is as follows:

[0070]

[0071] Calculate the vertical distance s from each node F of the element to the arc segment R of the submarine pipeline axis, and the geographical coordinates of the foot C of the perpendicular of each node F of the element on the arc segment R of the submarine pipeline axis. The calculation method is:

[0072] The geographic coordinates of each node F of the element are (x f ,y f ), the equation of the arc segment R of the submarine pipeline axis is (xx r ) 2 +(yy r ) 2 =r 2 , The calculation formula for the vertical distance s from each node F of the element to the arc segment R of the submarine pipeline axis is as follows:

[0073]

[0074] The geographical coordinates of the foot C of the perpendicular to the arc segment R of the submarine pipeline axis of each node F are (x c ,y c ), the calculation formula is as follows:

[0075]

[0076] Step 4: Calculate the mileage kp of the path of the foot of perpendicular C along the axis of the submarine pipeline based on the geographic coordinates of the foot of perpendicular C of each node F on the corresponding segment of the submarine pipeline axis and the geographic coordinates of each node P on the submarine pipeline axis. c .

[0077] In the above steps, the geographic coordinates of the foot of the perpendicular C (x c ,y c ) and the geographic coordinates of node P (x p ,y p ), if the foot of the perpendicular C coincides with the node P, then the mileage kp of the path along the axis of the seabed pipeline is c kp is the mileage of the path from node P to the axis of the submarine pipeline p If the foot of the perpendicular C is between two nodes P, the distance Δd from the foot of the perpendicular C to the previous node P is calculated, and the mileage kp of the path of the foot of the perpendicular C along the axis of the submarine pipeline is c kp is the mileage of the previous node P along the axis of the submarine pipeline P-1 +Δd.

[0078] The mileage between the feet of perpendiculars of adjacent nodes of the element on the corresponding segments of the submarine pipeline axis (called "segmental mileage along the submarine pipeline axis path") is also calculated.

[0079] Step 5: When the feature type is a linear feature, the distance between adjacent nodes and the cumulative length d of the feature are calculated based on the geographic coordinates of each node F of the feature. When the feature type is a polygonal feature, the area A of the feature is calculated based on the geographic coordinates of each node F of the feature;

[0080] In the above steps, for a single linear feature (1ine), according to the geographic coordinates (x f ,y f ) calculates the distance between adjacent nodes and the cumulative length d of the element. For a single polygon element (polygon), according to the geographic coordinates (x f ,y f ) Calculate the area A of the feature using the following formula:

[0081]

[0082] Where n is the number of nodes of the polygon feature.

[0083] Step 6: Compute the offset distance s of each node F of the element, the geographical coordinates (x c ,y c ), the mileage kp of the path along the axis of the seabed pipeline along the foot of the perpendicular C c , vertical foot C along the axis of the seabed pipeline path segment mileage The cumulative length d of linear features and the area A of polygonal features are added to the feature's attribute table.

[0084] In the above steps, when the feature type is a linear feature, the feature attribute table includes 10 attributes: ID, feature name, feature type, feature code, geographic coordinates of each node F of the feature (x f ,y f ), the offset distance s of each node F of the element, the geographical coordinates (x c ,y c ), the path mileage kp of the vertical foot C along the axis of the seabed pipeline c , the segmental mileage of the path along the axis of the seabed pipeline with the foot of perpendicular C The cumulative length d of the linear features.

[0085] When the feature type is a polygon feature, the feature attribute table includes 10 attributes: ID, feature name, feature type, feature code, geographic coordinates of each node F of the feature (xf ,y f ), the offset distance s of each node F of the element, the geographical coordinates (x c ,y c ), the path mileage kp of the vertical foot C along the axis of the seabed pipeline c , the segmental mileage of the path along the axis of the seabed pipeline with the foot of perpendicular C The area A of the polygon feature.

[0086] Step 7: Determine the intersection relationship between the element and the submarine pipeline based on the vertical distance s from each node F of the element to the corresponding segment of the submarine pipeline axis in the element attribute table.

[0087] In the above steps, when the vertical distances s from each node F of the element to the corresponding segment of the submarine pipeline axis are all positive or all negative, it indicates that the element is completely on one side of the submarine pipeline and has no intersection with the submarine pipeline. When the vertical distances s from each node F of the element to the corresponding segment of the submarine pipeline axis include positive and negative values, it indicates that the element and the submarine pipeline have an intersection relationship, and the intersection point is between adjacent nodes F where the vertical distances s are positive and negative.

[0088] Step 8: Count the number and type of features that intersect with the submarine pipeline.

[0089] Step 9: According to the feature attribute table, the mileage kp of the vertical foot C along the axis of the seabed pipeline c Determine the relative position of the feature in the submarine pipeline axis path and the length of the feature along the submarine pipeline axis path.

[0090] In the above steps, according to the mileage kp of the vertical foot C along the axis of the submarine pipeline in the feature attribute table, c The maximum and minimum values ​​in determine the relative position of the feature in the submarine pipeline axis path and the length of the feature along the submarine pipeline axis path.

[0091] According to the mileage kp of the vertical foot C along the axis of the submarine pipeline in the feature attribute table c , the offset distance s is used to draw a scatter plot of the distribution of elements along the axis path of the submarine pipeline.

[0092] Step 10: Count the distribution intervals of the same type of elements along the axis of the submarine pipeline.

[0093] In the above steps, when the mileage intervals of the two elements along the axis path of the submarine pipeline overlap, the union of the two elements is calculated. Assume that the mileage interval of the axis path of the submarine pipeline of element 1 is [kp c1 , kp c2 ], the mileage interval of the path along the axis of the submarine pipeline of element 2 is [kp c3 , kp c4 ], and kpc1 <kp c3 <kp c2 <kp c4 , whose union interval is [kp c1 -kp c4 ], repeat the above steps to complete the calculation of the distribution interval of a certain type of element on the axis path of the submarine pipeline.

[0094] Step 11: Calculate the cumulative length distribution of linear features.

[0095] In the above steps, the cumulative length distribution of linear features is statistically analyzed based on the feature name and the cumulative length d of the feature in the feature attribute table.

[0096] Step 12: Count the area distribution of polygonal elements.

[0097] In the above steps, the area distribution of polygonal features is counted based on the feature name and the feature area A in the feature attribute table.

[0098] Step 13: Analyze whether the feature is within the submarine pipeline buffer zone based on the vertical distance s from each node F of the feature to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer zone distance h in the feature attribute table.

[0099] In the above steps, when the absolute values ​​of the vertical distances s from each node of the element to the corresponding segment of the submarine pipeline axis are all less than h, the element is judged to be within the submarine pipeline buffer zone. When the absolute values ​​of the vertical distances s from each node of the element to the corresponding segment of the submarine pipeline axis are all greater than h, the element is judged to be outside the submarine pipeline buffer zone. When the absolute values ​​of the vertical distances s from each node of the element to the corresponding segment of the submarine pipeline axis are partially greater than h and partially less than h, the element is judged to be partially within the submarine pipeline buffer zone.

[0100] Step 14: Count the types and quantities of features in the submarine pipeline buffer.

[0101] Step 15. Generate a report on the spatial distribution of elements in the submarine pipeline routing area based on the number and type of elements intersecting with the submarine pipeline, the distribution interval of the same type of elements along the submarine pipeline axis path, the cumulative length distribution of linear elements, the area distribution of polygonal elements, and the type and number of elements in the submarine pipeline buffer zone, and complete the spatial analysis of elements in the submarine pipeline routing area.

[0102] In the above steps, the intersection relationship between the element and the plane position of the submarine pipeline can be judged, the length, area and distribution interval of the element can be statistically analyzed, and the submarine pipeline buffer zone can be analyzed, providing a basis for the analysis of third-party operation traces in the integrity management of the submarine pipeline.

[0103] In summary, the principle of this application is: the axis of the submarine pipeline is the path, and the direction of the path is the direction of the submarine pipeline. The submarine pipeline is divided into multiple nodes P according to the axis equation of the submarine pipeline, and the geographic coordinates (x p ,y p ) Calculate the mileage kp of each node P along the axis of the submarine pipeline c , and corresponds to the geographic coordinates of node P. According to the geographic coordinates (x f ,y f ) calculates the offset distance s of the node and the geographic coordinates (x c ,y c ), calculate the mileage kp of the path of the foot of perpendicular C along the axis of the seabed pipeline according to the geographical coordinates of the foot of perpendicular C c , along the segmented mileage of the underwater pipeline axis path According to the geographic coordinates (x f ,y f ) Calculate the cumulative length d of the linear element and the area A of the polygonal element, and generate the attribute table of the element. According to s in the attribute table of the element, determine whether the element intersects with the submarine pipeline. According to kp in the attribute table of the element c The length of the element along the axis of the submarine pipeline and its relative position on the path, as well as the overlapping relationship of the elements can be calculated. c The spatial distribution of elements along the submarine pipeline axis path can be analyzed by using s. The length distribution of elements of the same type can be statistically analyzed by using d in the element attribute table. The area distribution of elements of the same type can be statistically analyzed by using A in the element attribute table. Whether an element is within the submarine pipeline buffer can be determined by using s in the element attribute table and the submarine pipeline buffer distance h.

[0104] In summary, the spatial analysis method of submarine pipeline routing area elements of the present application can be used for the spatial analysis of third-party operation traces in the submarine pipeline routing area, and can display the number of elements intersecting with the submarine pipeline, the element type, the intersection position, the distribution interval of the same type of elements along the submarine pipeline, the length distribution of the same type of elements, the area distribution of the same type of elements, the element type and the number of elements in the submarine pipeline buffer zone, etc., which is convenient for the risk analysis of third-party operations of submarine pipelines.

[0105] The following is a detailed description with reference to specific embodiments:

[0106] The SG axis of a submarine pipeline consists of 4 straight line segments and 3 arc lines, namely L1, L2, L3, L4, C1, C2, and C3. Its mathematical equation is:

[0107]

[0108] The traces of third-party operations in the submarine pipeline routing area are mainly submarine anchor marks, submarine pits, etc.

[0109] (1) The axis of the submarine pipeline is the path, and the direction of the path is the direction of the submarine pipeline. According to the mathematical equation of the submarine pipeline axis, it is divided into different nodes at a spacing of 1m. Each node has an ID, feature name, feature type, and geographic coordinates (x p ,y p )4 attributes to generate the submarine pipeline attribute table.

[0110] (2) According to the geographic coordinates (x p ,y p ) Calculate the distance dn of each node along the path starting point, and obtain the mileage kp of each node along the path p , and added to the submarine pipeline attribute table, such as Figure 6 As shown in the table, the table is the submarine pipeline attribute table (partial).

[0111] (3) The ID, feature name, feature type, feature code, geographic coordinates (x f ,y f ) to generate feature attribute tables.

[0112] (4) Calculate the offset distance s of each node F and the geographical coordinates (x c ,y c ),in:

[0113] ① For a single element, traverse the y-coordinate and x-coordinate of its node F and determine the corresponding segment interval of the submarine pipeline according to formula (6);

[0114] ② According to the segment interval corresponding to the element node F, select the corresponding calculation formula;

[0115] ③ For the straight line interval, determine its parameters k and b1 according to formula (6), and then calculate the offset distance s of the element node F and the geographical coordinates (x c ,y c );

[0116] ④ For the arc segment interval, determine its parameter x according to formula (6): r ,y r , r, and then calculate the offset distance s of the element node F and the geographical coordinates (x c ,y c );

[0117] (5) Calculate the mileage of the element node F on the path. c ,y c ) and the geographic coordinates (x p ,y p ) is compared. If the foot of the perpendicular C coincides with the node P, then the mileage kp corresponding to the node P is p kp is the mileage of the foot of C c ; If the foot of the perpendicular C is between two nodes on the axis of the submarine pipeline, calculate the distance Δd from the foot of the perpendicular C to the previous node, then the mileage kp corresponding to the foot of the perpendicular C is c is the mileage kp of the previous node p-1 +Δd.

[0118] (6) For a single feature, calculate the segment mileage of the feature along the path

[0119] (7) For a single linear element (line), according to the geographic coordinates (x f ,y f ) Calculate the distance d between adjacent nodes i and the cumulative length d of the elements.

[0120] (8) For a single polygon, according to formula (5) and the geographic coordinates (x f ,y f )Calculate the area A of the polygon feature.

[0121] (9) The offset distance s of each node F and the geographical coordinates (x c ,y c ), mileage kp along the path c , Mileage along the route The cumulative length d of the linear features and the area A of the polygonal features are added to the feature attribute table to generate a complete attribute table with 10 attributes. Among them:

[0122] ① For linear features, including: ID, feature name, feature type, feature code, node geographic coordinates (x f ,y f ), offset distance s, vertical geographic coordinates (x c ,y c ), along the path mileage kp c , segment mileage along the path The cumulative length d, such as Figure 7 As shown in the table in , this is the attribute table of the linear features.

[0123] ②For polygonal features, including: ID, feature name, feature type, feature code, node geographic coordinates (x f ,y f ), offset distance s, vertical geographic coordinates (x c ,y c ), along the path mileage kp c , segment mileage along the path Area A, such as Figure 8 As shown in the table in , the table is the attribute table of the polygon features.

[0124] (10) Determine the intersection relationship between the feature and the submarine pipeline. The intersection relationship between the feature and the submarine pipeline can be determined based on the offset distance s in the feature attribute table, where:

[0125] ① When all s are positive or all negative, it means that the element is completely on one side of the submarine pipeline and has no intersection with the pipeline, such as Figure 7 As shown in the table, the feature has no intersection with the submarine pipeline;

[0126] ② When s contains positive and negative values, it indicates that the element and the submarine pipeline have an intersection relationship, and the intersection point is between adjacent nodes where s is one positive and one negative, such as Figure 8 As shown in the table in the offset distance s, the feature intersects the submarine pipeline.

[0127] (11) Count the number and type of elements that intersect with the submarine pipeline, such as Fig. 9 As shown in the table, the table is a statistical table of the intersection between features and submarine pipelines.

[0128] (12) Calculate the length of the feature along the path: According to the mileage kp along the path in the feature attribute table c The maximum and minimum values ​​can be used to calculate the length of the feature along the path and the relative position of the feature on the path.

[0129] (13) According to the mileage kp along the path in the feature attribute table c , the offset distance s is used to draw a scatter plot of the distribution of features along the path, such as Figure 2 As shown, Figure 2 This is the spatial distribution map of the elements of the submarine pipeline routing area.

[0130] (14) Element overlay analysis. When the mileage kp of two elements along the path c When there is overlap, the union of the two elements can be calculated: Assume that the mileage interval of element 1 on the path is [kp c1 , kp c2 ], the mileage interval of element 2 on the path is [kp c3 , kp c4 ], and kp c1 <kp c3 <kp c2<kp c4 , whose union interval is [kp c1 -kp c4 ].

[0131] (15) Repeat step (14) to complete the calculation of the distribution space of a certain type of elements on the path, such as Figure 3 As shown, Figure 3 This is the distribution map of the seabed anchor marks along the path direction.

[0132] (16) Cumulative length statistics of linear features. According to the feature name in the feature attribute table and the cumulative length d, the length distribution of the same type of features is counted, such as Figure 4 As shown, Figure 4 It is a statistical chart of the cumulative length of linear features (seabed anchor marks).

[0133] (17) Area statistics of polygonal features. According to the feature name and area A in the feature attribute table, the area distribution of the same type of features is counted, such as Figure 5 As shown, Figure 5 This is an area statistics map of polygonal features (seabed pits).

[0134] (18) Determine whether the feature is within the buffer zone. According to the offset distance s and the submarine pipeline buffer distance h in the feature attribute table, analyze whether the feature is within the buffer zone, where:

[0135] ① When the absolute value of s is less than h, it can be determined that the element is in the buffer zone;

[0136] ② When the absolute value of s is greater than h, it can be determined that the element is outside the buffer zone;

[0137] ③When the absolute value of s is partially greater than h and partially less than h, it is partially within the buffer zone.

[0138] like Fig.10 As shown in the table, the table is a statistical table of the number and type of features in the submarine pipeline buffer zone.

[0139] (19) The number and type of elements intersecting with the submarine pipeline are counted, the spatial distribution of the elements on the submarine pipeline path is calculated, the length of the linear elements is counted, the area of ​​the polygonal elements is counted, and it is determined whether the elements are within the buffer zone, thus completing the spatial analysis of the elements within the submarine pipeline routing area.

[0140] In summary, the spatial analysis method of submarine pipeline routing area elements of the present application calculates the attribute table of the elements according to the geographic coordinates of the element nodes and the geographic coordinates of the submarine pipeline nodes, performs spatial analysis of the elements according to the attribute table of the elements, and provides a complete spatial analysis report for the third-party operation trace analysis in the submarine pipeline integrity management.

[0141] In general, the spatial analysis method of submarine pipeline routing area elements of the present application can display the number of elements intersecting with the submarine pipeline, the element type, the intersection position, the distribution interval of the same type of elements along the submarine pipeline, the length distribution of the same type of elements, the area distribution of the same type of elements, the element type and the number of elements in the submarine pipeline buffer zone, etc., which is convenient for conducting risk analysis of third-party operations of submarine pipelines.

[0142] The same or similar numbers in the drawings of this embodiment correspond to the same or similar parts; in the description of this application, it should be understood that if the terms "upper", "lower", "left", "right", etc. indicate an orientation or position relationship based on the orientation or position relationship shown in the drawings, it is only for the convenience of describing this application and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation. Therefore, the terms describing the position relationship in the drawings are only used for illustrative purposes and cannot be understood as a limitation on this patent. For ordinary technicians in this field, the specific meanings of the above terms can be understood according to specific circumstances.

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

Claims

1. A method for spatial analysis of regional elements of submarine pipeline routing, characterized in that: The following steps are involved: Calculation of the axis equation of submarine pipelines in geographic coordinates; According to the submarine pipeline axis equation, the submarine pipeline axis is divided into multiple nodes P, and the mileage of the submarine pipeline axis node P along the submarine pipeline axis path is calculated; According to the geographic coordinates of the element node F, calculate the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis, and the geographic coordinates of the foot C of the perpendicular of the element node F on the corresponding segment of the submarine pipeline axis; According to the geographical coordinates of the foot of the perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis and the geographical coordinates of the node P on the submarine pipeline axis, the mileage kp of the path of the foot of the perpendicular C along the submarine pipeline axis is calculated. c ; When the feature type is a linear feature, the distance between adjacent nodes and the cumulative length d of the feature are calculated based on the geographic coordinates of the feature node F; when the feature type is a polygonal feature, the area A of the feature is calculated based on the geographic coordinates of the feature node F; The intersection relationship between the element and the submarine pipeline is determined based on the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis; the mileage kp of the path along the submarine pipeline axis is determined based on the foot of the perpendicular C. c Determine the relative position of the element in the submarine pipeline axis path and the length of the element along the submarine pipeline axis path; analyze whether the element is in the submarine pipeline buffer zone based on the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer zone distance h; Count the number and type of features that intersect with the submarine pipeline, count the distribution intervals of the same type of features along the axis path of the submarine pipeline, count the cumulative length distribution of linear features, count the area distribution of polygonal features, count the type and number of features in the submarine pipeline buffer zone, and complete the spatial analysis of the elements of the submarine pipeline routing area.

2. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: Calculation of the axis equation of the submarine pipeline under geographic coordinates, including: The plane position of the submarine pipeline is decomposed into straight line segments and arc segments, and the intervals are divided. The equations of the straight line segments and arc segments are calculated respectively, and then they are combined into the axis equation of the submarine pipeline.

3. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: Calculate the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis, and the geographic coordinates of the foot C of the element node F on the corresponding segment of the submarine pipeline axis, including: Calculate the vertical distance s from the element node F to the straight line segment L of the submarine pipeline axis, and the geographical coordinates of the foot C of the perpendicular to the element node F on the straight line segment L of the submarine pipeline axis. The calculation method is: The geographic coordinates of feature node F are (x f ,y f ), the equation of the straight line segment L of the submarine pipeline axis is y=kx+b1, x∈[x1,x2], and the calculation formula of the vertical distance s from the element node F to the straight line segment L of the submarine pipeline axis is as follows: The geographical coordinates of the foot C of the perpendicular to the straight line segment L of the submarine pipeline axis of the element node F are (x c ,y c ), the calculation formula is as follows:

4. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: Calculate the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis, and the geographic coordinates of the foot C of the element node F on the corresponding segment of the submarine pipeline axis, including: Calculate the vertical distance s from the element node F to the arc segment R of the submarine pipeline axis, and the geographic coordinates of the foot C of the perpendicular to the element node F on the arc segment R of the submarine pipeline axis. The calculation method is: The geographic coordinates of feature node F are (x f ,y f ), the equation of the arc segment R of the submarine pipeline axis is (xx r ) 2 +(yy r ) 2 =r 2 , x∈[x3,x4], then the calculation formula of the vertical distance s from the element node F to the arc segment R of the submarine pipeline axis is as follows: The geographical coordinates of the foot C of the perpendicular to the arc segment R of the submarine pipeline axis are (x c ,y c ), the calculation formula is as follows:

5. The method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: According to the geographical coordinates of the foot of the perpendicular C of the element node F on the corresponding segment of the submarine pipeline axis and the geographical coordinates of the node P on the submarine pipeline axis, the mileage kp of the path of the foot of the perpendicular C along the submarine pipeline axis is calculated. c ,include: The geographic coordinates of the foot of the perpendicular C (x c ,y c ) and the geographic coordinates of node P (x p ,y p ), if the foot of the perpendicular C coincides with the node P, then the mileage kp of the path along the axis of the seabed pipeline is c kp is the mileage of the path from node P to the axis of the submarine pipeline p ; If the foot of the perpendicular C is between two nodes P, the distance △d from the foot of the perpendicular C to the previous node P is calculated, and the mileage kp of the path of the foot of the perpendicular C along the axis of the seabed pipeline is c kp is the mileage of the previous node P along the axis of the submarine pipeline P-1 +△d.

6. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: The area A of the feature is calculated based on the geographic coordinates of the feature node F. The calculation formula is as follows: Where n is the number of nodes of the polygon feature.

7. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: The intersection relationship between the element and the submarine pipeline is determined based on the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis, including: When the vertical distances s from the element node F to the corresponding segment of the submarine pipeline axis are all positive or all negative, it means that the element is completely on one side of the submarine pipeline and has no intersection with the submarine pipeline; When the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis contains positive and negative values, it indicates that the element and the submarine pipeline have an intersection relationship, and the intersection point is between the adjacent nodes F with positive and negative vertical distances s.

8. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: According to the mileage kp of the path along the axis of the seabed pipeline according to the foot of the vertical C c Determine the relative position of the element along the axis path of the submarine pipeline and the length of the element along the axis path of the submarine pipeline, including: According to the mileage kp of the path along the axis of the seabed pipeline according to the foot of the vertical C c The maximum and minimum values ​​in determine the relative position of the feature in the submarine pipeline axis path and the length of the feature along the submarine pipeline axis path.

9. A method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: According to the vertical distance s from the element node F to the corresponding segment of the submarine pipeline axis and the submarine pipeline buffer distance h, analyze whether the element is within the submarine pipeline buffer, including: When the absolute value of the vertical distance s from the element node to the corresponding segment of the submarine pipeline axis is less than h, the element is judged to be within the submarine pipeline buffer zone; When the absolute value of the vertical distance s from the element node to the corresponding segment of the submarine pipeline axis is greater than h, the element is judged to be outside the submarine pipeline buffer zone; When the absolute values ​​of the vertical distances s from the element node to the corresponding segments of the submarine pipeline axis are partially greater than h and partially less than h, it is determined that the element is partially within the submarine pipeline buffer zone.

10. The method for spatial analysis of submarine pipeline routing regional elements according to claim 1, characterized in that: Statistics on the distribution intervals of the same type of elements along the axis of the submarine pipeline, including: When the mileage intervals of two elements along the axis path of the submarine pipeline overlap, the union of the two elements is calculated.

Images