Three-dimensional modeling and space coordinate conversion method for multi-source geophysical prospecting information

Abstract

The invention relates to the technical field of tunnel construction modeling, in particular to a three-dimensional modeling and space coordinate conversion method for multi-source geophysical prospecting information, which is space conversion from a geophysical prospecting three-dimensional model under a mileage positioning tunnel relative coordinate system to a three-dimensional model of a WGS84 geodetic coordinate system. The problem of the relation between railway tunnel narrow space coordinates and railway length geodetic space coordinates can be well solved. Through the conversion technology, geophysical prospecting data of different formats and types obtained by different geophysical prospecting devices can form a standard modeling file, three-dimensional modeling work of geophysical prospecting attribute data can be rapidly completed in modeling software in the later period, and finally three-dimensional integration visualization in a GIS system is completed. The rock mass characteristics of the rock mass in front of the tunnel excavation face can be visually and rapidly evaluated through the relative spatial position relation between the models, and therefore a basis is provided for comprehensive analysis of multiple forecasting means of tunnel advanced geological forecasting.

Description

Technical Field

[0001] The invention relates to the technical field of tunnel construction modeling, and in particular to a three-dimensional modeling and space coordinate conversion method of multi-source geophysical exploration information. Background Art

[0002] In recent years, with the continuous advancement of highway and railway projects in the mountainous and canyon areas of the southwest, a large number of tunnel projects have emerged. Most of these tunnel projects are located in areas with active geological movements. During the tunnel construction process, due to the large burial depth, they often face complex engineering geological problems such as high ground stress and alteration zones, because the construction risk is extremely high. In order to explore the construction risks ahead as much as possible, advanced geological forecasters will use different geophysical equipment to collect multi-source geophysical information to complement and verify each other. At present, the comprehensive use of different geophysical equipment and solution methods can accurately predict the geological conditions ahead of the tunnel excavation face, which has been widely used.

[0003] At present, there is no standard rapid modeling file and process for geophysical information data of different formats obtained by various geophysical methods due to the different working principles of the acquisition equipment, different equipment parameters and models, and different parsing protocols for parsing the original data. It is difficult to quickly complete the rapid modeling of physical property information obtained by different geophysical equipment.

[0004] At the same time, since a railway line is long, especially a long trunk railway, the length is hundreds or even thousands of kilometers. However, the railway is a linear project, and the width of a double-track railway tunnel is 7-8m and the height is about 6m, which is obviously not the same scale as the hundreds of kilometers of railway lines. How to connect the railway line with the geophysical exploration work done in the narrow tunnel space has become a problem that needs to be solved urgently. Summary of the invention

[0005] The purpose of the present invention is to: in view of the problems existing in the prior art, the present invention provides a three-dimensional modeling and spatial coordinate conversion method of multi-source geophysical exploration information.

[0006] First, a tunnel space coordinate conversion method is provided; after selecting the center line of the line in the two-dimensional engineering geological plan, it is processed into a multi-segment line, and the railway broken chain mileage is read. After completing the structured storage, the conversion relationship between the tunnel's relatively independent coordinate system and the spatial geodetic coordinate system is established, thereby effectively solving the relationship between the narrow spatial coordinates of the railway tunnel and the geodetic spatial coordinates of the railway length. This method can automatically complete the standard conversion technology of the tunnel space independent coordinate system to the spatial geodetic coordinate system for multi-source geophysical data obtained by various different geophysical equipment, and form a standard modeling file, quickly complete the three-dimensional modeling of geophysical attribute data in the modeling software, and finally complete the three-dimensional integrated visualization in the GIS system.

[0007] In order to achieve the above object, the technical solution adopted by the present invention is:

[0008] A tunnel space coordinate conversion method comprises the following steps:

[0009] Convert the horizontal curve of the tunnel line into the horizontal curve break table data; convert the vertical curve of the tunnel line into the vertical curve break table data; wherein the horizontal curve break table data and the vertical curve break table data both include the pre-break mileage crown number and the pre-break mileage, the post-break mileage crown number and the post-break mileage of the broken chain position, and the plane coordinates of the multi-line segment;

[0010] The three-dimensional coordinate data (mileage number, X0, Y0, Z0) under the relatively independent coordinates of the tunnel are respectively matched with the horizontal curve fracture table data and the vertical curve fracture table data using the mileage number;

[0011] Read the plane coordinates of the multi-line segments in the corresponding horizontal curve break table data as a polyline object; disassemble the polyline object into a single-segment line array; initialize the remaining length of the cache variable; calculate the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0);

[0012] The plane coordinates of the polyline segments in the corresponding vertical curve break table data are read as polyline objects; the Y value of the intersection of the longitudinal extreme line geometric object and the longitudinal extreme line on the construction mileage coordinate is the elevation Z1;

[0013] That is, the conversion of the tunnel's relatively independent coordinate system (mileage number, X0, Y0, Z0) to the spatial geodetic coordinate system (X1, Y1, Z1) is realized.

[0014] The invention provides a tunnel space coordinate conversion method, comprising: converting a tunnel line horizontal curve into horizontal curve fracture table data; converting a tunnel line vertical curve into vertical curve fracture table data; respectively making position correspondence between three-dimensional coordinate data (mileage crown number, X0, Y0, Z0) under relatively independent coordinates of the tunnel and horizontal curve fracture table data and vertical curve fracture table data by using the mileage crown number; reading the plane coordinates of multiple line segments in the corresponding horizontal curve fracture table data as a polyline object; disassembling the polyline object into a single-segment line array; initializing the remaining length of a cache variable; calculating the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0); reading the plane coordinates of multiple line segments in the corresponding vertical curve fracture table data as a polyline object; constructing the Y value of the intersection of the longitudinal extreme line geometric object and the longitudinal extreme line on the mileage coordinate as the elevation Z1; realizing the spatial conversion of the geophysical three-dimensional model under the tunnel relative coordinate system positioned by the mileage to the three-dimensional model of the WGS84 geodetic coordinate system, and being able to well solve the relationship problem between the narrow space coordinates of the railway tunnel and the geodetic space coordinates of the railway length. Through this conversion technology, geophysical data of different formats and types obtained by different geophysical equipment can be formed into a standard modeling file. Later, the three-dimensional modeling of the geophysical attribute data can be quickly completed in the modeling software, and finally the three-dimensional integrated visualization in the GIS system can be completed.

[0015] Further, the horizontal curve of the tunnel line is converted into the horizontal curve fracture table data, and / or the vertical curve of the tunnel line is converted into the vertical curve fracture table data, specifically using the following method:

[0016] Use 2D image software to open the geological plan of the tunnel construction site;

[0017] Select the railway centerline position in the 2D image software;

[0018] Load the plug-in and process the selected line into a polyline;

[0019] Use the plug-in to interactively read the start and end mileage and plane position of the broken chain;

[0020] Export plane coordinates and start and end mileage information to complete structured data conversion.

[0021] Furthermore, load the plug-in and process the selected line into a polyline. The specific operation process is as follows:

[0022] When the line position is a straight line segment or a transition curve segment, extract the node coordinate array;

[0023] When the line position is an arc segment, the arc coordinate array is obtained according to the accuracy through the GetPointAtDist() function based on the spacing;

[0024] Merge the node coordinate array and the arc coordinate array.

[0025] Furthermore, the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0) is calculated. The specific operation process is as follows: traverse the single segment line array, find the corresponding segment where X1 falls into the jth interval, and take out the starting point of the jth segment (seglines[j]) in the single segment line array. end

[0026] Calculate line segment vector

[0027] calculate The right normal vector of When X se =0:X n =Y se , Y n =0; when Y se =0:X n =0, Y n =-X se ; When Y se >0:X n =1, When Y se <0:X n =-1,

[0028] calculate The unit vector

[0029] calculate The unit vector Final calculation of plane coordinates

[0030] Furthermore, the relatively independent coordinate system of the tunnel is established based on the starting point of tunnel excavation, at the center of the tunnel entrance rail bottom with the tunnel excavation direction as the positive axis of the X-axis, the horizontal right deviation direction perpendicular to the tunnel excavation direction as the positive direction of the Y-axis, and the vertical direction relative to the center of the tunnel entrance rail bottom as the positive direction of the Z-axis; X0 is the tunnel mileage value; Y0 is the distance of the obtained geophysical parameter space point offset from the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate; Z0 is the altitude of the geophysical parameter space point relative to the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate.

[0031] The invention also provides a three-dimensional modeling method for multi-source geophysical exploration information.

[0032] A three-dimensional modeling method for multi-source geophysical exploration information includes the following steps:

[0033] Collecting rock physical characteristic parameters of the rock mass in front of the tunnel to be excavated; the rock physical characteristic parameters include at least one of shear wave wavelength, longitudinal wave wavelength and resistivity;

[0034] Processing the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel, wherein the three-dimensional modeling data file I includes (mileage number, X0, Y0, Z0) coordinate data;

[0035] The three-dimensional modeling data file I is subjected to coordinate transformation using the above method to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system;

[0036] The three-dimensional modeling data file II is three-dimensionally modeled using three-dimensional modeling software to obtain a three-dimensional geophysical property model.

[0037] The present invention provides a three-dimensional modeling method of multi-source geophysical information, including collecting rock physical characteristic parameters of the rock mass in front of a tunnel to be excavated; processing the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel, wherein the three-dimensional modeling data file I includes (mileage crown number, X0, Y0, Z0) coordinate data; performing coordinate conversion on the three-dimensional modeling data file I using the above method to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system; performing three-dimensional modeling on the three-dimensional modeling data file II using three-dimensional modeling software to obtain a three-dimensional geophysical property model. The method uses BIM technology to establish a standard modeling format for three-dimensional modeling of geophysical data of different formats and types obtained by different geophysical equipment, and uses three-dimensional modeling software to realize three-dimensional geophysical modeling. After the final conversion into the standard 3dtiles format file of the GIS system, the geophysical 3D model, the tunnel structure BIM model, the ultra-deep directional drilling model, etc. can be expressed in the same space on the computer. Through the relative spatial position relationship between the models, the rock mass characteristics in front of the tunnel excavation face can be quickly evaluated intuitively, thereby providing a basis for the comprehensive analysis of various forecasting methods for tunnel advanced geological forecasting. Most staff do not need to enter the tunnel and can repeatedly use the model to make judgments and decisions.

[0038] Furthermore, it also includes converting the three-dimensional geophysical attribute model into 3dtiles format and using the GIS system to achieve import and integrated visualization.

[0039] Another object of the present invention is to provide a three-dimensional modeling system of multi-source geophysical exploration information.

[0040] A system using the above-mentioned three-dimensional modeling method of multi-source geophysical exploration information comprises:

[0041] The first module collects the rock physical property parameters of the rock mass ahead of the tunnel to be excavated;

[0042] The second module is used to process the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel;

[0043] The third module is used to obtain the spatial coordinate transformation relationship between the tunnel's relatively independent coordinate system and the spatial geodetic coordinate system;

[0044] The fourth module is used to perform coordinate transformation on the three-dimensional modeling data file I to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system;

[0045] The fifth module is used to perform three-dimensional modeling on the three-dimensional modeling data file II using three-dimensional modeling software to obtain a three-dimensional geophysical property model.

[0046] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the three-dimensional modeling method of multi-source geophysical exploration information as described above.

[0047] A computer device comprises: a processor and a memory storing a computer program, wherein the processor is configured to execute the above-mentioned three-dimensional modeling method of multi-source geophysical exploration information when running the computer program.

[0048] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0049] 1. The present invention provides a method for converting spatial coordinates of a tunnel, comprising converting a horizontal curve of a tunnel line into horizontal curve fracture table data; converting a vertical curve of a tunnel line into vertical curve fracture table data; correspondingly positioning three-dimensional coordinate data (mileage number, X0, Y0, Z0) of the tunnel under relatively independent coordinates with horizontal curve fracture table data and vertical curve fracture table data using the mileage number; reading the plane coordinates of multiple line segments in the corresponding horizontal curve fracture table data as a polyline object; decomposing the polyline object into a single line array; initializing the remaining length of a cache variable; calculating the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0); reading the plane coordinates of multiple line segments in the corresponding vertical curve fracture table data as a polyline object; constructing the Y value of the intersection of the longitudinal extreme line geometric object and the longitudinal extreme line on the mileage coordinate as the elevation Z1; and realizing the spatial conversion of the geophysical three-dimensional model under the relative coordinate system of the tunnel located at the mileage to the three-dimensional model of the WGS84 geodetic coordinate system, which can well solve the relationship problem between the narrow space coordinates of the railway tunnel and the geodetic space coordinates of the railway length. Through this conversion technology, geophysical data of different formats and types obtained by different geophysical equipment can be formed into a standard modeling file. Later, the three-dimensional modeling of the geophysical attribute data can be quickly completed in the modeling software, and finally the three-dimensional integrated visualization in the GIS system can be completed.

[0050] 2. The present invention provides a three-dimensional modeling method of multi-source geophysical information, including collecting rock physical characteristic parameters of the rock mass in front of the tunnel to be excavated; processing the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel, the three-dimensional modeling data file I including (mileage crown number, X0, Y0, Z0) coordinate data; using the above method to perform coordinate conversion on the three-dimensional modeling data file I to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system; using three-dimensional modeling software to perform three-dimensional modeling on the three-dimensional modeling data file II to obtain a three-dimensional geophysical attribute model. The method uses BIM technology to establish a standard modeling format for three-dimensional modeling of geophysical data of different formats and types obtained by different geophysical equipment, and uses three-dimensional modeling software to realize three-dimensional geophysical modeling. After the final conversion into the standard 3dtiles format file of the GIS system, the geophysical 3D model, the tunnel structure BIM model, the ultra-deep directional drilling model, etc. can be expressed in the same space on the computer. Through the relative spatial position relationship between the models, the rock mass characteristics in front of the tunnel excavation face can be quickly evaluated intuitively, thereby providing a basis for the comprehensive analysis of various forecasting methods for tunnel advanced geological forecasting. Most staff do not need to enter the tunnel and can repeatedly use the model to make judgments and decisions. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a flow chart of the method for visualizing three-dimensional modeling and spatial coordinate transformation of multi-source geophysical exploration information in the embodiment.

[0052] Figure 2 Schematic diagram of the relatively independent coordinate system of the tunnel.

[0053] Figure 3 It is the TSP (seismic wave method) standard four-column multi-source geophysical 3D modeling data at the PDK282+146 tunnel face.

[0054] Figure 4 Flowchart for conversion of structured data of tunnel line curve.

[0055] Figure 5 The process of processing line positions for CAD plug-ins.

[0056] Figure 6 Query the logical diagram for the broken link list.

[0057] Figure 7 Logic diagram for plane coordinate calculation.

[0058] Figure 8 It is the tunnel line horizontal curve information database.

[0059] Fig. 9 It is the vertical curve feature database of tunnel lines.

[0060] Fig.10 It is a spatial integration diagram of the three-dimensional geophysical property model and the tunnel structure BIM model.

[0061] Fig.11 The figure is a flow chart of the tunnel space coordinate conversion method of the present invention.

[0062] Fig.12 The invention discloses a three-dimensional modeling method for multi-source geophysical exploration information. DETAILED DESCRIPTION

[0063] The present invention will be described in detail below in conjunction with the accompanying drawings.

[0064] 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.

[0065] Example

[0066] like Fig.11 As shown, the present invention provides a tunnel space coordinate conversion method, comprising the following steps:

[0067] Convert the horizontal curve of the tunnel line into the horizontal curve break table data; convert the vertical curve of the tunnel line into the vertical curve break table data; wherein the horizontal curve break table data and the vertical curve break table data both include the pre-break mileage crown number and the pre-break mileage, the post-break mileage crown number and the post-break mileage of the broken chain position, and the plane coordinates of the multi-line segment;

[0068] The three-dimensional coordinate data (mileage number, X0, Y0, Z0) under the relatively independent coordinates of the tunnel are respectively matched with the horizontal curve fracture table data and the vertical curve fracture table data using the mileage number;

[0069] Read the plane coordinates of the multi-line segments in the corresponding horizontal curve break table data as a polyline object; disassemble the polyline object into a single-segment line array; initialize the remaining length of the cache variable; calculate the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0);

[0070] The plane coordinates of the polyline segments in the corresponding vertical curve break table data are read as polyline objects; the Y value of the intersection of the longitudinal extreme line geometric object and the longitudinal extreme line on the construction mileage coordinate is the elevation Z1;

[0071] That is, the conversion of the tunnel's relatively independent coordinate system (mileage number, X0, Y0, Z0) to the spatial geodetic coordinate system (X1, Y1, Z1) is realized.

[0072] like Fig.12 As shown, a three-dimensional modeling method of multi-source geophysical exploration information includes the following steps:

[0073] Collecting rock physical characteristic parameters of the rock mass in front of the tunnel to be excavated; the rock physical characteristic parameters include at least one of shear wave wavelength, longitudinal wave wavelength and resistivity;

[0074] Processing the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel, wherein the three-dimensional modeling data file I includes (mileage number, X0, Y0, Z0) coordinate data;

[0075] The three-dimensional modeling data file I is subjected to coordinate transformation using the above method to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system;

[0076] The three-dimensional modeling data file II is three-dimensionally modeled using three-dimensional modeling software to obtain a three-dimensional geophysical property model.

[0077] This conversion method is a spatial conversion of the geophysical three-dimensional model in the relative coordinate system of the tunnel located by mileage to the three-dimensional model of the WGS84 geodetic coordinate system, which can well solve the relationship between the narrow spatial coordinates of the railway tunnel and the geodetic spatial coordinates of the railway length. Through this conversion technology, the geophysical data of different formats and types obtained by different geophysical equipment can be formed into a standard modeling file. Later, the three-dimensional modeling of the geophysical attribute data can be quickly completed in the modeling software, and finally the three-dimensional integrated visualization in the GIS system can be completed. Through the relative spatial position relationship between the models, the rock mass characteristics of the rock mass in front of the tunnel excavation face can be intuitively evaluated quickly, thereby providing a basis for the comprehensive analysis of various prediction methods for tunnel advanced geological prediction.

[0078] Specifically:

[0079] This embodiment specifically provides a method for three-dimensional modeling and spatial coordinate conversion of rock physical property parameters obtained by various geophysical equipment in the advanced geological prediction work during the tunnel construction stage through BIM technology. The three-dimensional geophysical model after spatial coordinate conversion is spatially integrated and visualized with the tunnel structure BIM model and the ultra-deep directional drilling three-dimensional model in the GIS system, and the spatial body with abnormal physical properties is screened, so as to predict the construction risk ahead of the tunnel excavation face.

[0080] The constructed geophysical information 3D modeling process is simple, practical and efficient. The constructed spatial conversion method can convert the railway line position (multi-line segment) in the CAD design drawing into structured data, establish the spatial conversion relationship between the engineering independent coordinate system based on mileage positioning and the geodetic WGS84 coordinate system, and finally realize the standardized 3D modeling and spatial coordinate conversion of multi-source geophysical information data (specific details are attached). Figure 1 ).

[0081] The specific steps include:

[0082] S1. Establish a relatively independent coordinate system for the tunnel (see Figure 2 shown):

[0083] Based on the starting point of tunnel excavation, a relatively independent coordinate system of the tunnel is established at the center of the tunnel entrance rail bottom with the tunnel excavation direction as the positive axis of the X axis, the horizontal right deviation direction perpendicular to the tunnel excavation direction as the positive direction of the Y axis, and the vertical direction relative to the center of the tunnel entrance rail bottom as the positive direction of the Z axis.

[0084] S2. Standard four-column multi-source geophysical 3D modeling data (X0, Y0, Z0) in the tunnel's relatively independent coordinate system:

[0085] After establishing a relatively independent coordinate system for the tunnel, the original geophysical data obtained by different geophysical methods and equipment (common ones include transient electromagnetic method, seismic wave method, ground audio frequency electromagnetic method, etc.) are collected and processed into a four-column multi-source geophysical three-dimensional modeling data standard format under the relatively independent coordinate system of the tunnel, as shown in Table 1.

[0086] Table 1

[0087]

[0088] The X0 coordinate is the tunnel mileage value. For example, the XX tunnel conducted a seismic wave geophysical exploration at the tunnel mileage DK285+900, so the starting value of the X0 coordinate is 285900; the Y0 coordinate is the distance of the obtained geophysical parameter space point offset from the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate, and the right side perpendicular to the tunnel excavation direction is a positive value; the Z0 coordinate is the altitude of the obtained geophysical parameter space point relative to the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate; the physical property parameter values ​​include shear wave wavelength, longitudinal wave wavelength, resistivity, etc.

[0089] For example, experimental research was carried out in a railway tunnel in Sichuan Province. The construction of the railway tunnel began in 2018. In order to find out the risks of rock excavation ahead, TSP (seismic wave) geophysical data collection was carried out at the PDK282+146 mileage of the tunnel guide.

[0090] An independent coordinate system for the tunnel was established, and the original data of the TSP (seismic wave) method at the PDK282+146 mileage of the tunnel guide was collected. The standard four-column multi-source geophysical 3D modeling data was obtained through the geophysical data processing software, such as Figure 3 shown.

[0091] S3. Extraction of horizontal curve features of tunnel lines from CAD drawings:

[0092] At present, railway survey and design work is basically completed using CAD software, and various disciplines conduct collaborative design on CAD software. Among them, in the CAD tunnel line design file, the tunnel line flat curve is often composed of straight line segments, transition curve segments, and arc segments. Straight line segments and transition curve segments are reflected as polylines in CAD drawings; arc segments require computer algorithms to encrypt the points on the arc before the encrypted points can be converted into polylines on the arc.

[0093] By developing a CAD plug-in, the line selected in the plan view can be quickly processed into a complete polyline (the arc segment is encrypted into a polyline, and the encryption spacing can be appropriately adjusted according to the accuracy requirements), and then the mileage before and after the break position of each broken link is interactively specified, and the program exports the horizontal curve information of the tunnel to realize the conversion from CAD line file to structured data table, and the converted tunnel line structured data information is included in the database to complete data storage. The specific extraction process is as follows: Figure 4 shown.

[0094] The process of CAD plug-in processing line position is as follows:

[0095] When the line position is a straight line segment or a transition curve segment, extract the node coordinate array;

[0096] When the line position is an arc segment, the arc coordinate array is obtained according to the accuracy through the GetPointAtDist() function based on the spacing;

[0097] Merge the node coordinate array and the arc coordinate array. The specific process is as follows: Figure 5 shown.

[0098] The standard format structure is as follows (the example includes long and short broken links) as shown in Table 2.

[0099] Table 2

[0100] Starting mileage number Starting mileage End Mileage Crown Termination Mileage geometry D1K 200 D1K 2000 …… D2K 1200 D2K 3000 …… D1K 3200 D1K 4000 …… …… …… …… …… ……

[0101] This information is in units of broken link segments. Each line records the mileage information and plane projection position information of a complete broken link segment, where geometry is the WKT of the plane projection polyline, such as "LINESTRING(475814.763358089.19,……)". The x and y in Geometry correspond to the plane coordinates in the projection coordinate system.

[0102] S4. Extraction of vertical curve records of tunnel lines from CAD drawings

[0103] In a similar way to horizontal curves, vertical curve information is extracted from the longitudinal section CAD drawing. Its format structure is consistent with that of horizontal curves, as shown in Table 3.

[0104] Table 3

[0105] Starting mileage number Starting mileage End Mileage Crown Termination Mileage geometry D1K 200 D1K 2000 …… D2K 1200 D2K 3000 …… D1K 3200 D1K 4000 …… …… …… …… …… ……

[0106] Each line records the mileage and line elevation information of a complete broken link segment. The geometry is in the WKT format of a polyline, such as "LINESTRING(355.6 562.9,……)". The x in Geometry represents the plane expansion distance from the starting mileage point in the segment, and the y represents the line elevation of the corresponding position.

[0107] S5. Space coordinate transformation relationship:

[0108] Before batch converting coordinates, the program first reads the tunnel horizontal curve section table and vertical curve section table stored in steps S3 and S4 from the database.

[0109] S5.1. Convert the coordinates of the mileage in S2 (mileage crown number, X0, Y0, Z0) into spatial coordinates (X1, Y1, Z1). Traverse the broken link table and query the broken link index i where the mileage coordinates are located. The query logic is shown in Figure 6 .

[0110] S5.2. Calculate the plane position X1, Y1

[0111] Take out the starting mileage number (sgh), starting mileage (slc), ending mileage number (egh), ending mileage (elc), and geometry from the i-th record in the flat curve broken link table.

[0112] 1) Read geometry as a polyline object;

[0113] 2) Decompose the polyline object into a single segment line array (seglines);

[0114] 3) Initialize the remaining length of the cache variable OL (overlength);

[0115] OL = X1-slc;

[0116] 4) Calculate the coordinate position X1, Y1 at the mileage position (X0) and the offset (Y0).

[0117] Traverse the single-segment line array seglines, find the corresponding interval where X1 falls into the jth interval, and take out the starting point of the jth (seglines[j]) in the single-segment line array end See the logical flow Figure 7 .

[0118] Calculate line segment vector

[0119] calculate The right normal vector of

[0120] When X se =0:X n =Y se , Y n =0;

[0121] When Y se =0:X n =0, Y n =-X se ;

[0122] When Y se >0:X n =1,

[0123] When Y se <0:X n =-1,

[0124] calculate The unit vector

[0125] calculate The unit vector

[0126] Final calculation of plane coordinates

[0127]

[0128] S5.3. Calculate the elevation Z1

[0129] Using the broken link subscript of the mileage coordinate queried by S5.1, the geometric object (geometry) with the same subscript in the vertical curve broken link table record is retrieved from the database, and the geometric object (geometry) is read as a polyline object.

[0130] The Y value of the intersection of the longitudinal extreme line geometry object (geometry) and the longitudinal extreme line on the construction mileage coordinate is the elevation.

[0131] P(X1,Y1,Y1) is the spatial position corresponding to the coordinates (mileage number, X0,Y0,Z0), thereby completing the spatial coordinate conversion of the tunnel's relatively independent coordinate system - the spatial geodetic coordinate system (WGS84).

[0132] Specifically, the horizontal curve and vertical curve features of the tunnel line are extracted using the PDK282+146 tunnel CAD spatial line file and stored in a structured manner. Figure 8It is the tunnel line horizontal curve information database; Fig. 9 This is a query diagram for the vertical curve feature database of the tunnel line.

[0133] S6. Standard modeling data file coordinate conversion

[0134] The standard four-column multi-source geophysical 3D modeling data file in the tunnel independent coordinate system obtained in step S2 is automatically converted by program using the spatial coordinate conversion relationship obtained in step S5 to obtain a standard four-column multi-source geophysical 3D modeling data file in the spatial geodetic coordinate system (WGS84).

[0135] S7. Geophysical 3D modeling and visualization

[0136] The standard four-column multi-source geophysical three-dimensional modeling data file in the spatial geodetic coordinate system (WGS84) obtained in step S6 can be used to realize three-dimensional geophysical modeling of the tunnel in a relatively independent coordinate system using a variety of three-dimensional modeling software on the market such as EVS, GOCAD, and Super Dimension Chuangxiang, to obtain a three-dimensional geophysical attribute model and realize three-dimensional visualization of the geophysical attribute model.

[0137] S8.GIS system spatial integration and visualization

[0138] The three-dimensional geophysical property model generated in step S7 is converted into a standard 3dtiles format, which can support the import and integrated visualization of common GIS systems on the market.

[0139] After establishing the spatial coordinate conversion relationship, the standard four-column multi-source geophysical 3D modeling data (X0, Y0, Z0) in the relatively independent coordinate system of the tunnel is converted into the standard four-column multi-source geophysical 3D modeling data (X1, Y1, Z1) in the spatial geodetic coordinate system. The 3D geophysical attribute modeling is performed using the common commercial 3D modeling software on the market. After conversion to the standard format 3dtiles, it is imported into the cesuim GIS platform to realize the spatial integration of the 3D geophysical attribute model and the tunnel structure BIM model (see Fig.10 ).

[0140] 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 and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A tunnel space coordinate conversion method, characterized in that: The following steps are involved: Convert the horizontal curve of the tunnel line into the horizontal curve break table data; convert the vertical curve of the tunnel line into the vertical curve break table data; wherein the horizontal curve break table data and the vertical curve break table data both include the pre-break mileage crown number and the pre-break mileage, the post-break mileage crown number and the post-break mileage of the broken chain position, and the plane coordinates of the multi-line segment; The three-dimensional coordinate data (mileage number, X0, Y0, Z0) under the relatively independent coordinates of the tunnel are respectively matched with the horizontal curve fracture table data and the vertical curve fracture table data using the mileage number; Read the plane coordinates of the multi-line segments in the corresponding horizontal curve break table data as a polyline object; disassemble the polyline object into a single-segment line array; initialize the remaining length of the cache variable; calculate the coordinate position (X1, Y1) below the offset (Y0) at the mileage position (X0); read the plane coordinates of the multi-line segments in the corresponding vertical curve break table data as a polyline object; construct the Y value of the intersection of the longitudinal extreme line geometric object and the longitudinal extreme line on the mileage coordinate as the elevation Z1; Realize the conversion of the tunnel's relatively independent coordinate system (mileage number, X0, Y0, Z0) to the spatial geodetic coordinate system (X1, Y1, Z1).

2. The tunnel space coordinate conversion method according to claim 1, characterized in that: The horizontal curve of the tunnel line is converted into the horizontal curve fracture table data, and / or the vertical curve of the tunnel line is converted into the vertical curve fracture table data, specifically using the following method: Use 2D image software to open the geological plan of the tunnel construction site; Select the railway centerline position in the 2D image software; Load the plug-in and process the selected line into a polyline; Use the plug-in to interactively read the start and end mileage and plane position of the broken chain; Export plane coordinates and start and end mileage information to complete structured data conversion.

3. The tunnel space coordinate conversion method according to claim 2, characterized in that: Load the plug-in and process the selected line into a polyline. The specific operation process is as follows: When the line position is a straight line segment or a transition curve segment, extract the node coordinate array; When the line position is an arc segment, the arc coordinate array is obtained according to the accuracy through the GetPointAtDist() function based on the spacing; Merge the node coordinate array and the arc coordinate array.

4. The tunnel space coordinate conversion method according to claim 1, characterized in that: Calculate the coordinate position (X1, Y1) under the offset (Y0) at the mileage position (X0). The specific operation process is as follows: Traverse the single segment line array, find the corresponding interval of X1, and take out the starting point of the jth (seglines[j]) in the single segment line array end Calculate line segment vector calculate The right normal vector of When X se =0:X n =Y se , Y n =0; when Y se =0:X n =0, Y n =-X se ; When Y se >0:X n =1, When Y se <0:X n =-1, calculate The unit vector calculate The unit vector Final calculation of plane coordinates 5. The tunnel space coordinate conversion method according to any one of claims 1 to 4, characterized in that: The relatively independent coordinate system of the tunnel is established based on the starting point of tunnel excavation, with the tunnel excavation direction as the positive axis of the X-axis, the horizontal right deviation direction perpendicular to the tunnel excavation direction as the positive direction of the Y-axis, and the vertical direction relative to the center of the tunnel entrance rail bottom as the positive direction of the Z-axis. X0 is the tunnel mileage value; Y0 is the distance of the obtained geophysical parameter space point offset from the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate; Z0 is the altitude of the geophysical parameter space point relative to the center of the rail bottom at the mileage corresponding to the tunnel X0 coordinate.

6. A three-dimensional modeling method for multi-source geophysical information, characterized in that: The following steps are involved: Collecting rock physical characteristic parameters of the rock mass in front of the tunnel to be excavated; the rock physical characteristic parameters include at least one of shear wave wavelength, longitudinal wave wavelength and resistivity; Processing the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel, wherein the three-dimensional modeling data file I includes (mileage number, X0, Y0, Z0) coordinate data; The three-dimensional modeling data file I is subjected to coordinate transformation using the tunnel space coordinate transformation method as described in any one of claims 1 to 5 to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system; The three-dimensional modeling data file II is three-dimensionally modeled using three-dimensional modeling software to obtain a three-dimensional geophysical property model.

7. The three-dimensional modeling method of multi-source geophysical exploration information according to claim 6, characterized in that: It also includes converting the three-dimensional geophysical property model into 3dtiles format and using the GIS system to achieve import and integrated visualization.

8. A system using the three-dimensional modeling method of multi-source geophysical exploration information as claimed in any one of claims 6 to 7, characterized in that: include: The first module collects the rock physical property parameters of the rock mass ahead of the tunnel to be excavated; The second module is used to process the rock physical characteristic parameters into a three-dimensional modeling data file I in a relatively independent coordinate system of the tunnel; The third module is used to obtain the spatial coordinate transformation relationship between the tunnel's relatively independent coordinate system and the spatial geodetic coordinate system; The fourth module is used to perform coordinate transformation on the three-dimensional modeling data file I to obtain a three-dimensional modeling data file II in a spatial geodetic coordinate system; The fifth module is used to perform three-dimensional modeling on the three-dimensional modeling data file II using three-dimensional modeling software to obtain a three-dimensional geophysical property model.

9. A computer-readable storage medium storing a computer program, characterized in that: When the program is executed by a processor, the three-dimensional modeling method of multi-source geophysical exploration information as described in any one of claims 6 to 7 is implemented.

10. A computer device comprising: A processor and a memory storing a computer program, characterized in that the processor is configured to execute the three-dimensional modeling method of multi-source geophysical exploration information as described in any one of claims 6-7 when running the computer program.

Images