Disclosure of Invention
The invention aims to solve the defects in the prior art, and provides a geological three-dimensional modeling method and system based on multi-data fusion.
In order to achieve the purpose, the invention adopts the following technical scheme that the geological three-dimensional modeling method based on multi-data fusion comprises the following steps:
S1, acquiring a node coordinate sequence according to geological profile mapping results, extracting horizontal spacing information among nodes, calculating average features by combining all node spacing to form a node spacing proportion sequence, and dividing a profile into differential density sections according to the proportion features to obtain a profile section density label;
S2, searching the elevation and horizontal coordinates of nodes at corresponding positions in the geophysical exploration section by utilizing the subdivision sections determined by the section density labels, and combining the elevation difference and the space position difference between the nodes to form a subdivision section adjustment coefficient sequence;
S3, performing directional offset on the node elevation of the high-density section according to the subdivision section adjustment coefficient sequence, determining an offset trend by combining elevation correspondence information of the nodes in geological drilling columnar data, correcting the node angle change amplitude in the low-density section, and forming a multi-layer section reconstruction primitive after integration processing;
And S4, extracting the length relation of line segments at two sides of the node and the node turning angle information aiming at the node in the multi-layer profile reconstruction primitive to form a coupling characteristic set, and comparing the coupling characteristic set with similar coupling characteristics of adjacent profile nodes to obtain multi-layer profile coupling difference characteristics.
As a further scheme of the invention, the section density label comprises a node spacing proportion sequence, a section difference density characteristic and a section density mark, the subdivision section adjustment coefficient sequence comprises an elevation difference coefficient, a space position difference coefficient and a node adjustment sequence, the multi-layer section reconstruction primitive comprises a high-density section node offset result, a low-density section angle correction result and a section reconstruction node set, and the multi-layer section coupling difference characteristic comprises a line segment length relation difference value, a turning angle difference value and an inter-section coupling characteristic difference value.
As a further scheme of the invention, the specific steps of S1 are as follows:
S101, acquiring a node coordinate sequence in geological profile mapping results, and determining the horizontal distance between nodes through a space difference relation between corresponding coordinates based on plane coordinate values of adjacent nodes to obtain a node horizontal distance sequence;
S102, comparing node spacing information in the sequence with total spacing characteristics of the sequence according to the node horizontal spacing sequence, and obtaining a node spacing proportion sequence according to a proportion relation between the node spacing and the total characteristics;
S103, calling the node spacing proportion sequence, identifying the position exceeding the proportion difference threshold according to the difference condition between adjacent proportion values, dividing the section into continuous sections with different densities, and labeling the sections with corresponding density attributes to obtain the section density label of the section.
As a further scheme of the invention, the specific steps of S2 are as follows:
s201, acquiring a subdivision section determined by the section density label, calling a corresponding position node on a geophysical exploration section based on the subdivision section, extracting an elevation value and a horizontal coordinate value recorded by the node, and generating node space distribution information;
S202, according to the node space distribution information, aiming at elevation values and horizontal coordinate values between adjacent nodes, forming associated data for describing elevation fluctuation and position transfer by comparing elevation change conditions and space position change conditions, and obtaining a characteristic value sequence of inter-node position change;
s203, calling the inter-node position change characteristic value sequence, determining a corresponding relation by applying a corresponding elevation difference and a corresponding position difference for each group of characteristic values, and collecting adjustment index sets in the subdivision section to obtain a subdivision section adjustment coefficient sequence.
As a further scheme of the invention, the specific steps of S3 are as follows:
s301, adjusting a coefficient sequence based on the subdivision section, calling the elevation value of each node in the section and combining the coefficients, sequentially deriving the offset amplitude of the node elevation value in the corresponding section, and orderly sorting the offset amplitudes according to the node positions in the same subdivision section to obtain a node elevation offset amplitude sequence;
S302, according to the node elevation offset amplitude sequence, combining elevation information of the nodes in geological drilling columnar data, sequentially comparing differences of the node offset amplitude and the columnar elevation, determining offset trends of the nodes in the sections according to positive and negative relations of the differences, and obtaining a node offset trend sequence;
S303, calling the node offset trend sequence, sequentially searching the node angle change value in the low-density section, correcting the node angle, and sequentially integrating the corrected node angle of the low-density section with the node offset amplitude and the offset trend of the high-density section to obtain the multi-layer profile reconstruction primitive.
As a further scheme of the invention, the specific steps of S4 are as follows:
S401, acquiring nodes in the multi-layer profile reconstruction primitive, determining a proportional relationship between the lengths of left side and right side line segments of the nodes based on the length information of the line segments on two sides of the nodes and the position coordinates of the nodes, extracting length relationship attributes by utilizing the combined characteristics of the nodes and the line segments on two sides, and generating a node length relationship value;
s402, according to the node length relation value, calling angle information formed between line segments on two adjacent sides of a node, determining corresponding coordination conditions between the length relation value and the angle value, forming joint description through geometric characteristics of the node, and obtaining a node coupling characteristic value;
S403, calling the coupling characteristic values of the nodes corresponding to the adjacent sections according to the node coupling characteristic values, determining the numerical difference degree according to the corresponding relation between the nodes, forming the coupling characteristic change expression of the node positions on the multi-layer sections, and obtaining the multi-layer section coupling difference value characteristics.
As a further scheme of the invention, the coupling characteristic values of the nodes corresponding to the adjacent sections refer to the coupling characteristic values extracted from the nodes on different sections in the multi-layer section data, and specifically are node coupling characteristics formed by coupling the length relation coefficient values of the line segments on two sides of each node with the included angle values;
The coupling characteristic change expression is based on the obtained coupling characteristic values of the nodes corresponding to the adjacent sections, and the obtained data expression describing the characteristic change condition of the nodes between the sections is obtained by comparing the coupling characteristic value differences among the corresponding nodes one by one.
As a further aspect of the present invention, the method further includes:
s5, forming geological three-dimensional splicing information based on the multi-layer profile coupling difference characteristic and combining the segment length relation of the nodes and the elevation offset information in the preamble step and the original coordinates of the nodes in geological profile mapping results and geophysical exploration profiles;
The geological three-dimensional splicing information comprises node segment length relations, elevation offset information and node original coordinates.
As a further scheme of the invention, the specific steps of S5 are as follows:
S501, based on the multi-layer profile coupling difference characteristic, calling original coordinates and elevation offset information of a geological profile mapping result and a node in a geophysical exploration profile, and extracting a corresponding segment length relation according to a position relation of the node in the layer profile to obtain a node segment length difference characteristic value set;
S502, calling the node segment length difference value characteristic value set, and combining elevation offset information of nodes, and aiming at the same-layer nodes, performing corresponding relation arrangement between the segment length difference value characteristic value and the elevation offset value to obtain a node offset relation sequence;
and S503, integrating the spatial correspondence of the nodes in the longitudinal direction and the transverse direction by combining the original coordinate values of the nodes based on the node offset relation sequence to generate geological three-dimensional splicing information.
A geological three-dimensional modeling system based on multi-data fusion, comprising:
The node proportion density module acquires a node coordinate sequence in geological profile mapping, calculates the horizontal distance between nodes, forms a node distance proportion sequence by utilizing the proportion relation between the node distance and the average value, and divides the profile difference density section according to the proportion sequence to obtain a profile section density label;
The section adjustment coefficient module determines the position of a subdivision section based on the section density label of the section, acquires the elevation and horizontal coordinates of nodes in the geophysical exploration section, calculates the elevation difference value and the horizontal difference value between the nodes to form a ratio sequence, and obtains the subdivision section adjustment coefficient sequence;
The node elevation offset module extracts the elevation of the high-density section node and the corresponding elevation in the geological drilling columnar data according to the subdivision section adjustment coefficient sequence, calculates a difference value to form an offset sequence, calculates the offset sequence and the subdivision section adjustment coefficient sequence, and performs ratio screening by combining the angle change amplitude of the low-density section node to obtain a multi-layer section reconstruction primitive;
The multi-layer coupling difference module extracts the length relation of line segments at two sides of a node and the turning angle of the node based on the multi-layer profile reconstruction primitive to form a coupling feature set, and performs difference calculation with the feature set of the node corresponding to the adjacent profile to obtain multi-layer profile coupling difference features;
and the three-dimensional splicing information module invokes the multi-layer section coupling difference value characteristic, and performs cumulative calculation by combining the node section length relation, the multi-layer section reconstruction primitive and the original coordinates of the nodes in the difference section to generate geological three-dimensional splicing information.
Compared with the prior art, the invention has the advantages and positive effects that:
According to the invention, the density sections are divided according to the pitch ratio of the nodes, the adjustment coefficient is formed by combining the difference of the elevation and the horizontal coordinate, the node direction deviation of the high density region and the angle correction of the low density region are accurately realized, the space form detail expression is enhanced, the coupling characteristics are formed by the lengths and the turning angles of the two sides of the nodes, the difference value is extracted in the multi-section comparison, the three-dimensional splicing is realized by combining the length of the section, the elevation deviation and the original coordinate, the geologic body structure spreading is more coherent, and the continuous processing with the density segmentation, the node fine adjustment and the coupling difference value as cores effectively releases the space association among different geological data, and the structure restoration and interpretation precision is remarkably improved.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
In the description of the present invention, it should be understood that the terms "length," "width," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like indicate orientations or positional relationships based on the orientation or positional relationships shown in the drawings, merely to facilitate describing the present invention and simplify the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and therefore should not be construed as limiting the present invention. Furthermore, in the description of the present invention, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
Referring to fig. 1, the geological three-dimensional modeling method based on multi-data fusion comprises the following steps:
S1, acquiring a node coordinate sequence according to geological profile mapping results, extracting horizontal spacing information among nodes, combining average features calculated by all node spacing to form a node spacing proportion sequence, and dividing a profile into differential density sections according to proportion features to obtain a profile section density label;
S2, searching elevation and horizontal coordinates of corresponding position nodes in the geophysical exploration section by utilizing subdivision sections determined by section density labels, and combining elevation differences and space position differences among the nodes to form subdivision section adjustment coefficient sequences;
S3, performing directional offset on the node elevation of the high-density section according to the subdivision section adjustment coefficient sequence, determining the offset trend by combining elevation corresponding information of the node in geological drilling columnar data, correcting the node angle change amplitude in the low-density section, and forming a multi-layer section reconstruction primitive after integration treatment;
S4, aiming at the nodes in the multi-layer profile reconstruction primitive, extracting the length relation of line segments at two sides of the nodes and the turning angle information of the nodes to form a coupling characteristic set, and comparing the coupling characteristic set with similar coupling characteristics of adjacent profile nodes to obtain multi-layer profile coupling difference characteristics;
and S5, forming geological three-dimensional splicing information based on the multi-layer profile coupling difference characteristic and combining the segment length relation of the nodes and the elevation offset information in the preamble step and the original coordinates of the nodes in the geological profile mapping result and the geophysical exploration profile.
The section density label comprises a node spacing proportion sequence, section difference density characteristics and section density marks, the subdivision section adjustment coefficient sequence comprises an elevation difference coefficient, a space position difference coefficient and a node adjustment sequence, the multi-layer section reconstruction primitive comprises a high-density section node offset result, a low-density section angle correction result and a section reconstruction node set, the multi-layer section coupling difference characteristics comprise a line segment length relation difference value, a turning angle difference value and an inter-section coupling characteristic difference value, and the geological three-dimensional splicing information comprises a node segment length relation, elevation offset information and node original coordinates.
Referring to fig. 2, the specific steps of s1 are:
S101, acquiring a node coordinate sequence in geological profile mapping results, and determining the horizontal distance between nodes through a space difference relation between corresponding coordinates based on plane coordinate values of adjacent nodes to obtain a node horizontal distance sequence;
When a node coordinate sequence in a geological profile mapping result is obtained, node data are firstly extracted from completed profile mapping files, the files are DWG or DXF format drawings, node coordinate lists are output through CAD tools or are directly converted into CSV files, X coordinates and Y coordinates of each node are sequentially read, the X coordinates and the Y coordinates are sequentially read from small to large according to node numbers, then horizontal distances between the current node and the next node are sequentially calculated step by step from the first node, specifically, the node 1 coordinates are 100.0 and 50.0, the node 2 coordinates are 103.0 and 54.0, the X coordinate difference is calculated first, the Y coordinate difference is 4.0, the horizontal distances can be understood as length values of 5.0 under the geometric relation of 3.0 and 4.0, the same operation is continuously carried out on all nodes, a complete node horizontal distance sequence is sequentially obtained, if a plurality of values in the sequence are remarkably larger or smaller than the current node and the next node, the absolute difference is firstly judged by a simple method, for example, the absolute difference is applied to the absolute difference is larger than the absolute difference, the absolute difference is found, the accurate distance is found in the corresponding to the node horizontal distance sequence, and the accurate distance is recorded in the final step, and the error is recorded, and the error is found in the corresponding to the node horizontal distance sequence.
S102, comparing node spacing information in a sequence with overall spacing characteristics of the sequence according to the node horizontal spacing sequence, and obtaining a node spacing proportion sequence according to a proportion relation between the node spacing and the overall characteristics;
According to the node horizontal spacing sequence, total feature statistics needs to be performed on all the pitches, including calculation of a group of average values, total distribution range and fluctuation condition of the horizontal spacing, wherein the average values are understood as that the node horizontal spacing is divided by the node spacing number after the total sum of all the node horizontal spacing is calculated, for example, if the pitches are sequentially obtained and are 5.0, 4.8, 5.3, 7.5, 4.6 and 5.1, the total average value is that the sum of the values is divided by 6 and is about 5.4, then each node spacing is compared with the average value, the ratio of each node spacing to the average value is calculated, for example, the ratio of 5.0 to 5.4 is about 0.93,7.5 to 5.4 is about 1.39, the node spacing ratio sequence is gradually obtained and is 0.93, 0.89, 0.98, 1.39, 0.86 and 0.95, and then threshold value division is performed on the ratios, for example, the ratio is considered as a normal interval between 0.8 and the ratio is more than 1.2, the ratio is judged to be significantly higher, the ratio is significantly lower, if the ratio is 1.39 falls in the significant high range, the ratio is about the average value, the ratio is 1.86 falls in the normal range, the ratio is the normal spacing is marked as the normal range, and the normal spacing is further marked by the node spacing is further marked and the normal range.
S103, calling a node spacing proportion sequence, identifying the position exceeding a proportion difference threshold according to the difference condition between adjacent proportion values, dividing the section into continuous sections with different densities, and labeling the sections with corresponding density attributes to obtain section density labels;
After the node spacing proportion sequence is called, the difference between each pair of adjacent proportions is gradually analyzed, the proportion difference value is calculated by subtracting the previous proportion from the latter proportion and taking the absolute value, for example, the difference between 0.89 and 0.93 is about 0.04,1.39 and the difference between 0.98 is about 0.41, the threshold value is continuously set to be 0.3 in the proportion difference value sequence for distinguishing normal change from obvious jump, if the difference value of 0.41 exceeds 0.3, the difference is marked as a obvious difference point, then the whole sequence is divided into a plurality of continuous sections by taking the obvious difference point as a boundary, for example, the sequence is divided into a first section from the head to a third node, only the sections with the proportions of 0.93, 0.89 and 0.98 are contained, the fourth node is singly formed into a section with the proportion of 1.39, the fifth node and the sixth node form the next section with the proportions of 0.86 and 0.95, the density attribute is judged according to the average proportion value of each section, if the average proportion of the sections is marked as medium density and is greater than 1.2 and less than 0.8 and less than 0.39 is marked as the average density, the section with the average proportion is marked as the average proportion is 1.39, the average section is marked as the average section with the average proportion and the average value is marked as the important section with the average value.
Referring to fig. 3, the specific steps of s2 are:
S201, acquiring a subdivision section determined by a section density label, calling a corresponding position node on a geophysical exploration section based on the subdivision section, extracting an elevation value and a horizontal coordinate value recorded by the node, and generating node space distribution information;
After the node space distribution information is provided, any pair of adjacent nodes in the subdivision section are processed one by one, elevation values are subtracted one by one respectively, and elevation changes among the adjacent nodes are obtained, for example, the elevation difference is-0.4 m when the node X is between 10m and 20m and the elevation difference is-0.4 m when the node H is between 121.3m and 120.9m, and the elevation drops by 0.4m from the position of 10m to the position of 20 m. The horizontal distance between the two points, Δx, here 10m, is recorded simultaneously. Continuing the processing of all nodes in the subdivided segment in pairs may form a series of elevation differences and corresponding horizontal distances, e.g., (-0.4 m,10 m), (0.3 m,10 m), (-0.2 m,10 m), etc. To reflect this change, the change in elevation per unit horizontal distance, i.e., slope, can be recalculated to a value of the difference in elevation divided by the horizontal distance, e.g., -0.4 divided by 10 to give-0.04, representing a height of Cheng Jiangdi cm per 1 m. For convenience of subsequent processing according to the gradient change degree, a reference value may be set, for example, set to ±0.02, which represents that when the absolute value of the gradient is greater than 0.02, the change is considered to be more obvious, and less than 0.02 is considered to be mild, so that the above-mentioned-0.04 is judged to be more obvious descending because of less than-0.02, and if a certain segment in the subsequent calculation is 0.03, the information and the numerical value are collected one by one to form a complete inter-node position change characteristic value list, and the information and the numerical value are sequentially arranged to be readable (from 10m to 20m: elevation difference-0.4 m, horizontal distance 10m, gradient-0.04). (from 20m to 30m: elevation difference 0.3m, horizontal distance 10m, gradient 0.03), etc., up to the end of the subdivision section.
S202, according to node space distribution information, aiming at elevation values and horizontal coordinate values between adjacent nodes, forming associated data for describing elevation fluctuation and position transfer by comparing elevation change conditions and space position change conditions, and obtaining inter-node change characteristic value sequences;
The specific calculation formula for the altitude change condition and the space position change condition is as follows:
;
calculating the composite change ratio between nodesThe method is used for obtaining a characteristic value sequence of inter-node change subsequently;
Wherein, theRepresents the firstAnd the firstDimensionless composite change ratio between adjacent nodes,Represents the firstThe elevation value (in m) of each node,Represents the firstThe elevation value (in m) of each node,Represents the firstThe horizontal coordinate values (units m) of the individual nodes,Represents the firstThe horizontal coordinate values (units m) of the individual nodes,Represents the firstAnother dimension of the individual nodes is a horizontal coordinate value (in m),Represents the firstAnother dimension of the individual nodes is a horizontal coordinate value (in m),Representing the average value (in m) of all node elevation values,Representing the average absolute deviation of the elevation values of all nodes in the node set from the average elevation value,Representing the total number of nodes,Representing the selected reference length for elevation normalization (units m);
the meaning and the acquisition mode of the parameters are as follows:
Represent the firstThe elevation values (in meters) of the individual nodes are obtained through total station GNSS elevation measurement data;
Represent the firstThe elevation values (in meters) of the individual nodes are also obtained by a high Cheng Shice value;
representing the average value of all node elevation values, and calculating the arithmetic average value of all node elevation sampling values;
Represent the firstAnd the firstThe horizontal x coordinate value (unit is meter) of each node is derived from the coordinate acquisition result of a measurement coordinate positioning instrument or a geographic information system;
Represent the firstAnd the firstThe horizontal y coordinate value (unit is meter) of each node, and the data acquisition mode is the same as that of the previous modeRepresenting the total node number for calculation, and directly obtaining the total node number through the space node layout total number;
Representing normalized reference length, selecting a principle as an average value of the horizontal distance of the maximum node in the area, and adopting a calculation mode: WhereinFor the number of all pairs of adjacent points,Is the reference point coordinates;
The actual sampling data according to a certain monitoring area is as follows:
Meter (measured by a laser scanner);
Meter (measured by a laser scanner);
In a node setThe elevation values are 113.6, 116.2, 114.9, 112.8 and 117.5 meters respectively;
Node coordinates:,,, all coordinate units are meters;
normalized reference lengthRice calculated from the average value of the maximum horizontal distance of adjacent nodes
Parameter intermediate value calculation is carried out: Rice;
;
the mean absolute deviation is: Rice;
after normalization, the method comprises the following steps:;
Substituting the molecular moiety: After normalization:;
horizontal distance section:;
Overall molecule:;
Substitution denominator part:, The product is: taking the absolute value still as 2.511;
And (3) denominator integration:;
The main formula is carried in:;
The result shows that the composite change ratio between nodes is 0.824, the larger the numerical value is, the more obvious the combined difference between elevation change and horizontal distance is, the result is used for subsequently judging the sequencing basis of the mutation degree in the inter-node change characteristic value sequence, and the elevation disturbance sensitive area can be further identified by combining a plurality of node sequences;
The method comprises the steps of summing squares of node elevation normalization difference values and horizontal Euclidean distances in a numerator part, adding root numbers, subtracting an average value of elevation deviation in a node set to form a composite change value of elevation and horizontal position coupling, dividing the normalized elevation difference by a reference length in operation to eliminate the problem of direct addition of different physical dimensions, calculating straight line distances of two points in a plane coordinate system by a structure of the squares and the root numbers, ensuring that quantized offset degrees are all non-negative representation integral deviation degrees by absolute values, multiplying the normalized elevation difference values and the horizontal coordinate difference values to obtain a combined change amplitude, adding a zero division risk and keeping a numerical smoothness, ensuring that the absolute value ensures that the direction does not influence scale judgment, and enabling local elevation fluctuation and horizontal displacement to be simultaneously included in the calculation of molecules of change factors, and finally forming a non-dimensional ratio reflecting the feature of the coupling offset between the nodes and the horizontal position by the coupling inhibition effect of the denominator;
The concrete meaning of calculating the composite change ratio between the nodes is that by fusing the high Cheng Qifu amplitude between the nodes and the horizontal position moving distance, a quantization index capable of comprehensively measuring the landform or the field change state between the two points is constructed, the ratio not only reflects the elevation difference of the adjacent nodes in the vertical direction, but also combines the space distance between the adjacent nodes on the plane position, thereby effectively describing the integral offset strength caused by the fluctuation, settlement, lifting or horizontal displacement of the terrain in the local area, and the size of the ratio can be used for identifying the potential abnormal change section in the surface form or the engineering structure and providing basis for the subsequent sorting of the difference degree between the nodes and the position with potential deformation hidden danger.
S203, calling a characteristic value sequence of inter-node position change, determining a corresponding relation by applying corresponding elevation difference and position difference for each group of characteristic values, and collecting an adjustment index set in a subdivision section to obtain a subdivision section adjustment coefficient sequence;
The inter-node change characteristic value sequence obtained through the point-by-point calculation can be used for taking out the gradient and horizontal distance values of each section group by group, and then combining with a pre-selected adjustment sensitivity parameter to form an adjustment coefficient sequence. For example, the sensitivity coefficient is 2.5, and the gradient value is multiplied by the sensitivity coefficient and then added to 1, so that the adjustment coefficient of each section is obtained. Taking a slope of-0.04 as an example, 1 plus 2.5 times-0.04, resulting in 0.9, this means that 90% adjustment is required for this segment, and if the slope is 0.03, the adjustment factor is 1.075, meaning that 107.5% adjustment is required. For batch processing, a column of gradient data can be directly filled in the electronic form, and the other column is filled in by a formula to multiply gradient values of corresponding rows by 1 plus 2.5, so that the adjustment coefficients can be automatically calculated row by row. After repeated calculation segment by segment, an adjustment coefficient sequence in the whole subdivision section is formed, for example (0.9,1.075,0.95,1.02.) is formed, and the sequence can be directly used for adjusting the earthwork volume, elevation mark or other engineering data of the corresponding section in the design drawing or calculation, so as to ensure that the subdivision section is reasonably distributed and changed according to the actually measured height Cheng Bodong.
Referring to fig. 4, the specific steps of s3 are:
s301, adjusting a coefficient sequence based on a subdivision section, calling an elevation value of each node in the section and combining the coefficients, sequentially deriving the offset amplitude of the node elevation value in the corresponding section, and orderly arranging the offset amplitudes according to the node positions in the same subdivision section to obtain a node elevation offset amplitude sequence;
Firstly, dividing the whole line range into a plurality of small sections with equal length, for example, dividing the whole line range into small sections with the length of 500m along the line, namely, dividing the small sections into 10m of each section, totaling 50 sections, then grading the small sections according to geological weathering conditions or lithology differences, for example, assigning an adjustment coefficient of 0.6 to a fully weathered section, 0.8 to a strong weathered section, 1.0 to a weak weathered section and 1.2 to an un-weathered section, forming a sequence and arranging the sequences in sequence, then calling the original elevation value of each node according to the position of the node, for example, the node 1 is 105.3m, the node 2 is 104.8m, directly multiplying the adjustment coefficient of the section of the node with the elevation of the node, for example, the node 1 is in the fully weathered section, the corresponding adjusted value is 63.2m, the corresponding value of the 2 nd node is 83.8m in the strong wind section, all nodes are continuously processed, a group of node elevation offset amplitude sequences can be obtained through sequential calculation, the node elevation and section adjustment coefficients can be input in batches to automatically generate offset amplitude values through batch operation of a spreadsheet tool or a script, the offset amplitude values are orderly arranged according to the line mileage positions of the nodes, the order arrangement from the starting point to the end point in a list is ensured, for example, the starting point 0m corresponds to 63.2m, the mileage 20m corresponds to 83.8m, the mileage 30m corresponds to 104.2m, the whole node elevation offset amplitude sequence arrangement is completed gradually to the line end point.
S302, according to the node elevation offset amplitude sequence and combining elevation information of the nodes in geological drilling columnar data, sequentially comparing differences of the node offset amplitude and the columnar elevation, determining offset trends of the nodes in the sections according to positive and negative relations of the differences, and obtaining a node offset trend sequence;
according to the sorted node elevation offset amplitude sequence, the corresponding node column elevation values in the drilling columnar geological data are further combined and compared one by one, for example, at the 1 st node position, the offset amplitude is 63.2m, the position in the drilling columnar data is 60.2m, the difference value is 3m, in this case, if the difference value is positive, the offset trend is considered to be downward, if the difference value is negative, the offset trend is considered to be upward, the offset direction can be represented by +1 downwards, -1 represents upwards, the difference values are compared point by point, for example, the 2 nd node offset amplitude is 83.8m, the columnar data is 81.0m, the difference value is 2.8m, the direction is +1, the 5 th node offset amplitude is 91.5m, the columnar data is 95.0m, the difference value is-3.5 m, the direction is-1, the judgment of the offset trend is continuously completed from the 1 st node to the last node in this way, for example, the sequence is obtained as [ +1, -1, ], the process can be represented by a column offset trend is calculated in an electronic table, the number of the offset amplitude is subtracted from the column amplitude in a simple format, the corresponding to the two directions are directly recorded in the corresponding column orientation, and the corresponding directions are recorded in the positive-phase, and the corresponding directions are clearly defined.
S303, calling a node offset trend sequence, sequentially searching node angle change values in the low-density section, correcting the node angles, and sequentially integrating the corrected node angles of the low-density section with the node offset amplitude and the offset trend of the high-density section to obtain a multi-layer profile reconstruction primitive;
The specific calculation formulas for sequentially searching the node angle change value and correcting the node angle in the low-density section are as follows:
;
calculating a node correction angle valueSequentially integrating the corrected node angles of the low-density sections with the node offset amplitude and offset trend of the high-density sections to obtain a multi-layer profile reconstruction primitive;
Wherein, theRepresents the firstThe corrected angle values for the individual low density segment nodes,Represents the firstThe individual nodeThe calculated angle change values between the associated nodes,Represents the firstThe individual nodeThe weight factor of the individual nodes is used,Represents the firstThe initial angular offset value of the individual nodes,Representing the average of all node initial angular offset values for the current low density segment,Represents the firstStandard deviation of the angle change values of the individual nodes,Representative and the firstThe total number of nodes associated with each node;
Angle correction valueIn the calculation of (a), the node number is first determinedThe node is located in the low-density section, and the total number of adjacent nodes isThe three-dimensional laser scanning system is used for sampling the structure of the section, and the following angle change values and initial angle offset values are obtained:
the angle change values between the node 4 and the five associated nodes are as follows, and are obtained by performing inverse cosine difference on the direction vector change amounts of the nodes in the continuous five-frame data frame, wherein the angle units are:
;
Weighting factorCalculating from Euclidean distance and signal to noise ratio measured value among nodes, and according to a quantization standard:
The weight is based on the signal-to-noise ratio (SNR) ratio standardized value, and the quantization mode is the SNR of the current comparison nodeCarrying out Sigmoid normalization transformation on the ratio of the signal to noise ratio average values of all the associated nodes;
The standardized weights obtained are respectively:
;
Thereby calculating a weighted sum:
,
Original angular offset value of node 4The initial angular offset values of the 5 nodes of the low density section are in order:
3.1, 2.7, 3.4, 2.8, 3.0, average:
;
Thus:
;
standard deviation ofThe calculation mode is based on eachAverage evolution of the square of the difference between the value and the mean:
Average value:;
standard deviation:;
The denominator is calculated as:;
Finally, all the components are replaced into a main formula:;
The result shows that the correction angle of the node 4 is 1.36 degrees, the original angle value has certain discreteness, the correction value can be used for reconstructing and integrating the multi-layer section node direction vectors, and the value is taken into the subsequent node position vector fitting or contour combination logic step to serve as the reference input of the final section vector field direction judgment;
The arithmetic logic of the formula is that the absolute value of the difference value of the node offset and the section average offset multiplied by the total number of the nodes is subtracted and the root number is opened to reflect the deviation degree of the consistency of the node and the whole structure and weaken the extreme deviation influence in a square root mode, the denominator part adopts the operation of dividing the standard deviation of the node angle change by the addition of the root number of the adjacent nodes, so as to increase the denominator to inhibit the abnormal correction amplitude when the fluctuation is large, ensure that the angle correction tends to be smooth under the influence of multiple adjacent nodes, ensure that the output angle correction value is non-negative through the absolute value by the whole formula, and form a set of angle correction logic framework based on the combined action of local weighting, deviation adjustment and variance smoothing;
the node correction angle value represents a numerical adjustment result of the original angle of the node by comprehensively considering the angle change relation between the node and adjacent nodes, the importance degree of each adjacent node and the offset degree of the node relative to the whole of the same section when the multi-node section structure is analyzed, and the value is directly used for replacing the original angle information of the node so as to embody the optimization direction of the node on the basis of the local structure trend and the whole offset consistency, and ensure that the subsequent data input by angles more conforming to the continuity of the real structure and the local smooth characteristic can be used as the basis when the multi-layer section profile, the node direction vector reconstruction and the node space position calculation are generated.
Referring to fig. 5, the specific steps of s4 are:
S401, acquiring nodes in a multi-layer profile reconstruction primitive, determining a proportional relationship between the lengths of left side and right side line segments of the nodes based on the length information of the line segments on two sides of the nodes and the position coordinates of the nodes, extracting length relationship attributes by utilizing the combined characteristics of the nodes and the line segments on two sides, and generating a node length relationship value;
Acquiring nodes in a multi-layer profile reconstruction primitive, firstly directly reading profile coordinate information of each layer through three-dimensional scanning or CAD profile data files, traversing each profile line point by point from beginning to end, collecting coordinates of each node in a list, for example, for a typical multi-layer profile file of a part, generally, each layer has about 150 nodes, all node data are accumulated to form a complete node point set, then sequentially selecting each node, firstly identifying the position of each node in the current profile, recording the front node on the left side of the node and the rear node on the right side, sequentially acquiring lengths of two side line segments through coordinate differences between the nodes, if the distance between two points on the left side of a certain node is 2.3mm and the distance between two points on the right side is 4.6mm, the length ratio of the left side and the right side of the node can be directly expressed by 2.3 to 4.6, the ratio value is about 0.5, then the left side and right side ratio values of all nodes are listed one by one to form a sequence, the sequence is continuously checked, the ratio value is limited between 0.4 and 2.5 in a reasonable range, if the ratio is smaller than 0.4 or larger than 2.5, the corresponding node is put into a separate sequence for marking, then the node in a normal interval is combined with the coordinates of the node and the length information on two sides to generate a length relation attribute combination, for example, the node with the ratio value of 0.5 is formed, the recordable length relation attribute is that the total length is 2.3 plus 4.6 to obtain 6.9, the difference value is 2.3 minus 4.6 to obtain-2.3, the ratio value of 0.5 is added, after all the nodes execute the process, the complete node length relation attribute sequence is formed, the ratio value is extracted separately from the complete node length relation attribute sequence, and the numerical value sequence of the node length relation is formed for subsequent calculation.
S402, according to the node length relation value, calling angle information formed between line segments on two adjacent sides of a node, determining corresponding coordination conditions between the length relation value and the angle value, and forming joint description through geometric characteristics of the node to obtain a node coupling characteristic value;
According to the numerical value of the node length relation, angle information formed by line segments on two adjacent sides of a node is further called, the position of a node is firstly taken out from a node list, then the two line segments on the left side and the right side of the node are checked by taking the node as the center, the trend of the node is directly compared in a text mode through a three-dimensional coordinate point, for example, if the left line segment extends 1 along the X direction and 2 along the Y direction, the right line segment extends 2 along the X direction and is basically unchanged, the angle is understood to be between 60 and 70 degrees through geometry, the angle value is judged by changing the coordinate position of each node one by one in the mode, then the length proportion value of the node and the corresponding angle value form a pair, for example, the length proportion value of the node is 63.5 degrees for the proportion value, the node can be recorded as (0.5,63.5 degrees), the node is continuously traversed and the corresponding sequence of the length proportion value and the angle is formed, data is provided for further analysis, and for the detection situation of a specific mechanical blade cavity, the reasonable angle range can be set to be 30 degrees to 150 degrees, for the angle value which is smaller than 30 degrees or 150 degrees, namely, the angle value is understood to be the angle value is directly calculated to be the angle value which is the angle is smaller than 30 degrees or equal to or greater than 150 degrees, the angle value is the approximate to the value, the value is approximately equal to the value and the value is approximately equal to the value of the value and is approximately and the value is the value of the value.
S403, aiming at the node coupling characteristic values, calling the coupling characteristic values of the nodes corresponding to the adjacent sections, determining the numerical difference degree through the corresponding relation among the nodes, forming the coupling characteristic variation expression of the node positions on the multi-layer section, and obtaining the multi-layer section coupling difference value characteristics;
For the coupling characteristic values of all the nodes, starting to perform comparison on the multi-layer profile, for example, sequentially finding out the nodes corresponding to each other in the upper layer and the lower layer through coordinate matching or shortest distance matching on one layer of the current profile and then taking out the coupling characteristic values of the nodes corresponding to the lower layer for comparison, calculating the difference of the coupling values, for example, the coupling value of the nodes of the current layer is 31.7, the value of the corresponding lower layer is 29.8, the difference is 1.9, thus processing all the corresponding nodes to obtain a group of coupling difference sequences, screening can be performed by setting a threshold value of 2.5 in the process of detecting workpieces, if the difference of the nodes is greater than 2.5, the nodes with larger difference are regarded as, and continuously counting the proportion of the nodes in all the nodes, for example, the statistics result shows that 90% of the nodes have the difference of 1.5 to 2.5, the rest about 10% of nodes can be concentrated in a larger difference section, the spatial position and the number of layers of the nodes need to be additionally recorded for later structural evaluation or correction, and finally the coupling characteristic of the nodes on the profile with the difference can be formed.
Referring to fig. 6, the specific steps of s5 are:
S501, based on the multi-layer profile coupling difference value characteristic, calling original coordinates and elevation offset information of a geological profile mapping result and a node in a geophysical exploration profile, and extracting a corresponding segment length relation according to the position relation of the node in the layer profile to obtain a node segment length difference value characteristic value set;
based on the demarcation point data extracted from multi-layer geological section mapping results, node numbers and corresponding coordinate information are read point by point, the node numbers and the corresponding coordinate information comprise X, Y coordinates and Z elevations, offset in Gao Chengji X, Y directions of corresponding nodes are obtained from the geophysical section results, corresponding indexes between the nodes are established, for adjacent nodes of the same layer, the spatial distance between two points is firstly calculated according to mapping original coordinates, for example, the distance between the nodes 1 and 2 is about 11.36m, the distance between the same two points is calculated according to geophysical coordinates to obtain about 12.68m, a segment length difference value 1.32m is obtained through the subtraction of the two points, the calculation method is repeatedly applied to all nodes, a group of difference value sets such as 1.32m, 0.95m and 1.10m are obtained, the group of data is subjected to partition processing according to absolute values, the threshold value is set to be 1m, all the protruding change segments are marked with the values of which are 1m, the distance is smaller than or equal to 1m, the relative stable segments are represented with 0, a group of sequences similar to the values 1, the distances between the nodes are calculated to be about 11.36m, the distances between the two points are calculated according to the geophysical coordinates, the difference values are 1.32m, the difference values are obtained, the difference values are obviously obtained in a map-like or the map-like section extension map segments are formed in a clear section, the subsequent color difference values are formed by the map section, and the difference values are obviously apparent in a map-like section extension method is formed.
S502, calling a node segment length difference value characteristic value set, and combining elevation offset information of nodes, and aiming at the same-layer nodes, performing corresponding relation arrangement between the segment length difference value characteristic value and the elevation offset value to obtain a node offset relation sequence;
After the characteristic value set of the segment length difference of the node is obtained, the data is paired with the elevation offset information of the node, for example, the segment length difference of the nodes 1, 2 and 3 is 1.32m, 0.95m and 1.10m, the elevation offset is 1m, the two are respectively combined into a node 1 (1.32,1), a node 2 (0.95,1) and a node 3 (1.10,1), then all the segment length difference values are uniformly scaled according to the maximum value, the maximum value is 1.32m as a reference, the node 2 is about 72%, the node 3 is about 83%, the elevation offset is 100% in the example, the segment length difference value accounts for 60% through the distribution weight, the elevation offset accounts for 40%, thus obtaining a node offset integrated value, the node 2 is about 83%, the node 3 is about 90%, the offset significance is divided according to the interval for the integrated value, 85% -100% is recorded as a higher section, 60% -85% is a middle section, and 60% -is a lower section, therefore, the nodes 1, 2 and 3 are all judged as a higher section, the offset relation sequence can be formed in software table, the whole trend graph can be conveniently drawn, the three-dimensional trend graph can be stored in the software, the graph can be conveniently mapped, the whole graph is ready, and the three-dimensional trend graph can be stored in the graph, and the graph can be conveniently stored.
S503, integrating the spatial correspondence of the nodes in the longitudinal direction and the transverse direction based on the node offset relation sequence and combining the original coordinate values of the nodes to generate geological three-dimensional splicing information;
According to the obtained node offset relation sequence, the space coordinates of the nodes are adjusted one by one according to the proportion, the height is Cheng Diejia in the longitudinal direction according to the offset proportion, for example, the offset of the node 1 is about 96%, the reference value for longitudinal superposition is 2m, therefore, the height of the node 1 is changed from 50m to about 52m, the node 2 is lifted from 48m to about 50m, the node 3 is lifted from 47m to about 49m, fine adjustment is carried out in the plane direction according to the same offset proportion, half of the space coordinates are respectively distributed in the X and Y directions by adopting the reference 1m, the node 1 is increased by about 0.68m in the X direction, the Y direction is also increased by about 0.68m, new coordinates are 100.68m and 200.68m, new X, Y coordinates of all the nodes are calculated by analogy, finally, the updated X, Y, Z coordinate list of all the nodes is arranged and output, and a geological three-dimensional splicing information is generated, and can be used for importing and checking the distribution state of the nodes in the three-dimensional space in the point cloud software.
Referring to fig. 2, a geological three-dimensional modeling system based on multi-data fusion includes:
The node proportion density module acquires a node coordinate sequence in geological profile mapping, calculates the horizontal distance between nodes, forms a node distance proportion sequence by utilizing the proportion relation between the node distance and the average value, and divides the profile difference density section according to the proportion sequence to obtain a profile section density label;
the section adjustment coefficient module determines the position of a subdivision section based on the density label of the section, acquires the elevation and horizontal coordinates of nodes in the geophysical exploration section, calculates the elevation difference value and the horizontal difference value between the nodes to form a ratio sequence, and obtains the subdivision section adjustment coefficient sequence;
The node elevation offset module extracts the elevation of the high-density section node and the corresponding elevation in the geological drilling columnar data according to the subdivision section adjustment coefficient sequence, calculates a difference value to form an offset sequence, calculates the offset sequence and the subdivision section adjustment coefficient sequence, and performs ratio screening by combining the angle change amplitude of the low-density section node to obtain a multi-layer section reconstruction primitive;
The multi-layer coupling difference module extracts the length relation of line segments at two sides of a node and the turning angle of the node based on the multi-layer profile reconstruction primitive to form a coupling feature set, and performs difference calculation with the feature set of the node corresponding to the adjacent profile to obtain multi-layer profile coupling difference features;
The three-dimensional splicing information module invokes the multi-layer section coupling difference characteristic, and performs cumulative calculation by combining the node section length relation, the multi-layer section reconstruction primitive and the original coordinates of the nodes in the difference section to generate geological three-dimensional splicing information.
The present invention is not limited to the above embodiments, and any equivalent embodiments which can be changed or modified by the technical disclosure described above can be applied to other fields, but any simple modification, equivalent changes and modification made to the above embodiments according to the technical matter of the present invention will still fall within the scope of the technical disclosure.