Generation method of paleogeographic coordinates of deposit evolution track based on plate motion constraint

Abstract

The application discloses a deposit evolution track paleo-geographic coordinate generation method based on plate motion constraint, relates to the technical field of geological space simulation data processing, and comprises the following steps: performing time slice matching operation on a deposit space-time attribute chain and a plate attribution judgment range, determining bound plates of each geological time slice, and serially connecting the bound plates according to the geological time sequence to obtain a deposit plate dynamic binding chain; extracting Euler rotation parameters of the corresponding bound plate from plate motion constraint data along the deposit plate dynamic binding chain, taking present deposit positioning data as a simulation starting point, obtaining initial paleo-geographic coordinates through reverse rotation simulation, and calculating adjacent coordinate rotation constraint paths in combination with the corresponding Euler rotation parameters; and correcting the initial paleo-geographic coordinates according to the adjacent coordinate rotation constraint paths and modifying constraint errors to obtain corrected paleo-geographic coordinates. The application realizes the cooperative constraint of paleo-coordinate inversion and plate rotation geometry.

Description

Technical Field

[0001] This invention relates to the field of geological spatial simulation data processing technology, and in particular to a method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints. Background Technology

[0002] With the development of deep-time digital earth, plate reconstruction, and geological spatial simulation modeling technologies, the generation of paleogeographic coordinates for ore deposit evolution trajectories has gradually evolved from relying on single-point paleocoordinate recovery and manual trajectory connection to an automated inversion method based on multi-source geological constraint data and plate movement parameters. Existing technologies typically utilize information such as current ore deposit location data, geological age data, paleomagnetic data, plate contours, and Eulerian rotation parameters to reconstruct the paleogeographic location of ore deposits at different geological periods, and analyze the spatial migration process of ore deposits caused by plate movement based on this.

[0003] However, existing technologies still have shortcomings: First, they focus on determining plate affiliation under a single geological time slice, making it difficult to jointly constrain the continuity of plate numbering, the proximity of plate outlines, and the consistency of rotation direction between adjacent geological time slices, which can easily lead to binding plate jumps in plate boundary areas. Second, paleogeographic coordinate inversion and trajectory connection are usually processed separately, making it difficult to further convert Euler rotation parameters into rotation constraint paths between adjacent coordinates, resulting in problems such as folding, error accumulation, and insufficient continuity in the generated ore deposit evolution trajectory. Summary of the Invention

[0004] In view of the aforementioned existing problems, the present invention is proposed.

[0005] Therefore, this invention provides a method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints, which solves the problems of insufficient reliability of paleogeographic coordinate inversion caused by plate instability in different geological time slices and discontinuous connection of discrete paleocoordinates.

[0006] To solve the above-mentioned technical problems, the present invention provides the following technical solution: This invention provides a method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints. The method includes: collecting ore deposit evolution constraint data of the target ore deposit and performing spatiotemporal geological normalization processing, then concatenating the data according to geological time sequence to form a ore deposit spatiotemporal attribute chain; dividing geological time slices according to the ore deposit spatiotemporal attribute chain, extracting plate tectonics constraint data and current ore deposit location data from the ore deposit evolution constraint data, then performing time slice matching and plate contour reconstruction on the plate tectonics constraint data to form the plate affiliation judgment range for each geological time slice; and performing time slice matching calculations on the ore deposit spatiotemporal attribute chain and the plate affiliation judgment range to determine... The bound plates of various geological time slices are concatenated in geological time sequence to obtain the dynamic binding chain of the ore deposit plates. Euler rotation parameters of the corresponding bound plates are extracted from the plate motion constraint data along the dynamic binding chain of the ore deposit plates. Using the current location data of the ore deposit as the simulation starting point, the initial paleogeographic coordinates are obtained through reverse rotation simulation. Then, the rotation constraint paths of adjacent coordinates are calculated in combination with the corresponding Euler rotation parameters. Based on the rotation constraint paths of adjacent coordinates, the initial paleogeographic coordinates are corrected for path continuity and constraint error to obtain the corrected paleogeographic coordinates. Finally, the paleogeographic coordinate set of the ore deposit evolution trajectory is generated in geological time sequence.

[0007] As a preferred embodiment of the method for generating paleogeographic coordinates of ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the ore deposit evolution constraint data includes current ore deposit location data, geological age constraint data, paleomagnetic constraint data, and plate tectonics constraint data.

[0008] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the specific steps for forming the spatiotemporal attribute chain of the ore deposit are as follows: Collect ore evolution constraint data of the target ore deposit, and perform spatiotemporal benchmark unification and outlier data cleaning on the ore evolution constraint data to obtain normalized ore evolution constraint data; Based on the geological time corresponding to the normalized ore deposit evolution constraint data, the normalized ore deposit evolution constraint data is grouped by time, and attribute merging, time stamp binding and node encapsulation are performed on the normalized ore deposit evolution constraint data under the same geological time to obtain spatiotemporal attribute nodes; The spatiotemporal attribute nodes are sorted according to geological time from ancient to the present, and a succession number is added to adjacent spatiotemporal attribute nodes. Then, the spatiotemporal attribute nodes are connected in series according to the succession number to form the spatiotemporal attribute chain of the ore deposit.

[0009] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the specific steps for extracting plate tectonics constraint data and current ore deposit location data from ore deposit evolution constraint data are as follows: Extract time nodes arranged in geological time sequence from the spatiotemporal attribute chain of the ore deposit, and divide the time interval between adjacent time nodes into geological time slices; Using geological time slices as search criteria, plate tectonics constraint data and current location data of ore deposits are extracted from ore deposit evolution constraint data.

[0010] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectories based on plate tectonics constraints described in this invention, the specific steps for determining the plate affiliation range of each geological time slice are as follows: Based on geological time slices, plate movement constraint data are retrieved by time to obtain the plate number, plate outline, and Euler rotation parameters corresponding to each geological time slice. The plate contours under the same plate number are reconstructed by rotating the Euler rotation parameters corresponding to each geological time piece to obtain the reconstructed plate contours. The reconstructed plate outline is closed by boundary closure, and the closed reconstructed plate outline is associated with the corresponding geological time slice and plate number to form the plate affiliation judgment range for each geological time slice.

[0011] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the specific steps for obtaining the dynamic binding chain of the ore deposit plates are as follows: By using geological time slices to analyze the spatial location attributes of the corresponding spatiotemporal attribute nodes in the spatiotemporal attribute chain of the ore deposit, the ore deposit location corresponding to each geological time slice is obtained. Then, the ore deposit location is spatially matched with the plate affiliation judgment range of the same geological time slice to determine the candidate binding plate of each geological time slice. For candidate bound plates in adjacent geological time slices, plate number comparison and plate outline proximity determination are performed to generate plate number consistency markers and outline proximity markers. The rotation direction angle is calculated based on the Euler rotation parameters corresponding to the candidate bound plates. Then, a ore deposit plate binding status map is constructed based on the plate number consistency markers, outline proximity markers and rotation direction angles. Based on the ore deposit plate binding status map, target binding paths that run through all geological time slices are selected. Then, candidate binding plates corresponding to each geological time slice in the target binding path are determined as binding plates. The binding plates are then connected in sequence according to geological time to obtain the dynamic binding chain of ore deposit plates.

[0012] As a preferred embodiment of the method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints described in this invention, the specific steps for obtaining the initial paleogeographic coordinates through reverse rotation simulation are as follows: The binding plates corresponding to each geological time piece are determined one by one along the dynamic binding chain of the ore deposit plates. The binding plates are then matched with the plate numbers corresponding to the same geological time piece in the plate motion constraint data to determine the Euler rotation parameters corresponding to the binding plates. The Euler rotation parameters corresponding to the bound plate are decomposed to obtain the Euler pole and the reverse rotation angle, and the reverse rotation matrix of the corresponding geological time slice is constructed based on the Euler pole and the reverse rotation angle. Using the current location data of the ore deposit as the starting point for simulation, the current location data of the ore deposit is converted into current unit spherical coordinates. Then, according to the geological time slices from present to past, the current unit spherical coordinates are sequentially input into the corresponding reverse rotation matrix for cumulative coordinate transformation, and the initial paleogeographic coordinates are output.

[0013] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the specific steps for calculating the adjacent coordinate rotation constraint path are as follows: Adjacent initial paleogeographic coordinates are paired in geological time sequence and converted into spherical starting and ending points. The unit rotation axis is determined based on the Euler pole, and the reverse rotation angle is inverted to generate the rotation angle. Starting from the spherical origin, the system rotates segmentally around the unit rotation axis according to the rotation angle, and uses the spherical endpoint as the termination constraint of the segmented rotation path. The segmented rotation path is then normalized and resampled to calculate the rotation constraint path between adjacent coordinates.

[0014] As a preferred embodiment of the paleogeographic coordinate generation method for ore deposit evolution trajectory based on plate tectonics constraints described in this invention, the specific steps for obtaining the corrected paleogeographic coordinates are as follows: The initial paleogeographic coordinates are converted into initial unit spherical coordinates, and the corresponding adjacent coordinate rotation constraint paths are sampled at equal intervals to obtain path sampling points. Then, the spherical angular distance between the initial unit spherical coordinates and each path sampling point is calculated, and the path sampling point with the smallest spherical angular distance is determined as the path projection point. The spherical angle between the path projection point and the initial unit spherical coordinates is measured and the tangential direction is decomposed to obtain the spherical angular distance and the path deviation direction. Then, the amplitude is limited and normalized and the amplitude is synthesized to obtain the constraint error correction amount. The initial unit spherical coordinates are corrected for path continuity based on the constraint error correction amount, and the corrected unit spherical coordinates are then converted into corrected paleogeographic coordinates.

[0015] As a preferred embodiment of the method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints described in this invention, the specific steps for generating the paleogeographic coordinate set of ore deposit evolution trajectories in geological time sequence are as follows: The corrected paleogeographic coordinates are merged with the corresponding geological time slices, bound plate numbers, adjacent coordinate rotation constraint paths, and constraint error corrections to form paleogeographic coordinate nodes. The paleogeographic coordinate nodes are sorted in geological time sequence, and adjacency markers and trajectory numbers are added to adjacent paleogeographic coordinate nodes after sorting to form a chain of mineral deposit evolution trajectory nodes. The node fields of the ore deposit evolution trajectory node chain are extracted and the coordinate records are encapsulated. Then, the time series index is arranged according to the trajectory number to obtain the paleogeographic coordinate set of the ore deposit evolution trajectory.

[0016] The beneficial effects of this invention are as follows: By constructing a ore deposit plate binding state diagram by combining the ore deposit location, plate attribution judgment range, and rotation direction angle, continuous determination of the bound plates of various geological time slices is achieved, and the attribution jump of the ore deposit near the plate boundary is avoided. This improves the stability of deep-time plate attribution and the reliability of paleogeographic coordinate inversion. By generating initial paleogeographic coordinates through an inverse rotation matrix and combining them with the rotation constraint path of adjacent coordinates for path continuity correction, the synergistic constraint of paleocoordinate inversion and plate rotation geometry is achieved, which suppresses trajectory breakpoints, reduces error accumulation, and improves the verifiability of ore deposit evolution trajectory. Attached Figure Description

[0017] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the following description of the embodiments will be briefly introduced. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 This is a flowchart of a method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints.

[0019] Figure 2 A flowchart for determining the tectonic plate attribution range for various geological timescales.

[0020] Figure 3 A flowchart for obtaining the dynamic binding chain of the mineral deposit plate.

[0021] Figure 4 A flowchart for generating a paleogeographic coordinate set of ore deposit evolution trajectories. Detailed Implementation

[0022] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0023] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0024] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in different places in this specification does not necessarily refer to the same embodiment, nor is it a single or selective embodiment that is mutually exclusive with other embodiments.

[0025] Reference Figures 1-4 This is one embodiment of the present invention, which provides a method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints, including the following steps: S1. Collect the ore deposit evolution constraint data of the target ore deposit and perform spatiotemporal geological normalization processing, and then connect them in sequence according to geological time to form the ore deposit spatiotemporal attribute chain.

[0026] Collect ore evolution constraint data of the target ore deposit, and perform spatiotemporal benchmark unification and outlier data cleaning on the ore evolution constraint data to obtain normalized ore evolution constraint data.

[0027] The specific process includes: during the target ore deposit data acquisition process, collecting current ore deposit location data through ore deposit location records, collecting geological age constraint data through geological age records, collecting paleomagnetic constraint data through paleomagnetic measurement records, and collecting plate movement constraint data through plate reconstruction records. The current ore deposit location data, geological age constraint data, paleomagnetic constraint data, and plate movement constraint data are then aggregated according to the target ore deposit to generate ore deposit evolution constraint data. The current ore deposit location data undergoes coordinate benchmark unification processing, the geological age constraint data undergoes age scale unification processing, the paleomagnetic constraint data and plate movement constraint data undergo time stamping unification processing, and the ore deposit evolution constraint data undergoes anomaly cleaning to obtain normalized ore deposit evolution constraint data.

[0028] It should be noted that the current location data of the ore deposit is used to determine the current spatial location of the target ore deposit, the geological age constraint data is used to provide an accurate age scale for time grouping, the paleomagnetic constraint data is used to verify the rationality of plate affiliation, and the plate movement constraint data is used to provide plate reconstruction parameters under different geological time periods.

[0029] Based on the geological time corresponding to the normalized ore deposit evolution constraint data, the normalized ore deposit evolution constraint data is grouped by time, and attribute merging, time stamp binding, and node encapsulation are performed on the normalized ore deposit evolution constraint data under the same geological time to obtain spatiotemporal attribute nodes.

[0030] The specific process includes: verifying the time stamp of each normalized ore deposit evolution constraint data according to the geological time corresponding to the normalized ore deposit evolution constraint data; grouping normalized ore deposit evolution constraint data with the same geological time into the same time group to ensure that the current location data, geological age constraint data, paleomagnetic constraint data, and plate movement constraint data of the ore deposit under the same geological time are consistent in time scale; organizing the normalized ore deposit evolution constraint data with the same target ore deposit orientation and the same geological time orientation within the same time group according to field correspondence; binding the normalized ore deposit evolution constraint data that has completed attribute merging with the corresponding geological time with time stamps to give the normalized ore deposit evolution constraint data that has completed attribute merging with time stamps the data with time stamps, so that the normalized ore deposit evolution constraint data that has completed attribute merging with time stamps under the same geological time into nodes according to the node format to obtain spatiotemporal attribute nodes.

[0031] The spatiotemporal attribute nodes are sorted according to geological time from ancient to the present, and a succession number is added to adjacent spatiotemporal attribute nodes. Then, the spatiotemporal attribute nodes are connected in series according to the succession number to form the spatiotemporal attribute chain of the ore deposit.

[0032] The specific process includes: verifying the time of all spatiotemporal attribute nodes according to the geological time corresponding to the spatiotemporal attribute nodes; arranging spatiotemporal attribute nodes with earlier geological time first and those with later geological time last, so that all spatiotemporal attribute nodes form an arrangement relationship from ancient to modern times; determining adjacent spatiotemporal attribute nodes in sequence according to the arrangement relationship from ancient to modern times, and writing the succession relationship between the previous spatiotemporal attribute node and the next spatiotemporal attribute node into the succession number, so that adjacent spatiotemporal attribute nodes have a front-to-back correspondence; and connecting all spatiotemporal attribute nodes in series according to the succession number, so that the previous spatiotemporal attribute node can point to the next spatiotemporal attribute node, and connecting all spatiotemporal attribute nodes continuously from ancient to modern times according to geological time, forming a spatiotemporal attribute chain of the ore deposit.

[0033] S2. Based on the spatiotemporal attribute chain of the ore deposit, geological time slices are divided, and plate movement constraint data and the current location data of the ore deposit are extracted from the ore deposit evolution constraint data. Then, time slice matching and plate outline reconstruction are performed on the plate movement constraint data to form the plate affiliation judgment range of each geological time slice.

[0034] Extract time nodes arranged in geological time sequence from the spatiotemporal attribute chain of the ore deposit, and divide the time interval between adjacent time nodes into geological time slices.

[0035] The specific process includes: using the time markers of each spatiotemporal attribute node in the spatiotemporal attribute chain of the ore deposit as the order basis, sequentially reading the geological age information corresponding to each spatiotemporal attribute node; screening and removing time data with missing or abnormal markers; and recording the remaining valid geological times sequentially according to the arrangement order of the spatiotemporal attribute nodes in the spatiotemporal attribute chain of the ore deposit to form an initial time node sequence; standardizing the initial time node sequence to convert geological times in different representations into the same time scale, and sorting the time nodes according to the order of geological time from ancient to the present; verifying the continuity and succession relationship of the start and end times of adjacent time nodes after sorting, identifying time overlap, discontinuity, and reverse order, and correcting and marking abnormal intervals to make the time nodes arranged according to geological time continuous and orderly; and finally dividing the continuous time range between adjacent time nodes into geological time slices.

[0036] Using geological time slices as search criteria, plate tectonics constraint data and current location data of ore deposits are extracted from ore deposit evolution constraint data.

[0037] The specific process includes: using the start and end times of geological time slices as the search scope, filtering plate movement constraint data whose time markers fall into the corresponding geological time slices from the ore deposit evolution constraint data, and merging the filtered plate movement constraint data according to the geological time slices; locating the current location data of the target ore deposit from the ore deposit evolution constraint data, and verifying the integrity of the longitude, latitude, and elevation fields in the current location data; establishing a correspondence between the merged plate movement constraint data and the verified current location data of the ore deposit according to the geological time slices, and obtaining the plate movement constraint data and current location data of the ore deposit corresponding to each geological time slice.

[0038] By performing time-based retrieval of plate tectonics constraint data according to geological time slices, the plate number, plate outline, and Euler rotation parameters corresponding to each geological time slice are obtained.

[0039] The specific process includes using geological time slices as time retrieval indexes, matching plate tectonics constraint data slice by slice, grouping and classifying the matched plate tectonics constraint data according to the plate number field, and reading the corresponding plate number field value as the plate number under the corresponding geological time slice; extracting plate contour data and Euler rotation parameters from the corresponding plate tectonics records based on the plate number, verifying the boundary integrity of the plate contour data, verifying the validity of the Euler rotation parameters, and then associating and organizing the plate number, plate contour, and Euler rotation parameters according to the geological time slices to obtain the plate number, plate contour, and Euler rotation parameters corresponding to each geological time slice.

[0040] The plate contours under the same plate number are reconstructed by rotating the Euler rotation parameters corresponding to each geological time piece.

[0041] The specific process includes extracting boundary coordinate points from the plate outline, converting the boundary coordinate points into unit spherical coordinates, and then performing a spherical rotation transformation on the unit spherical coordinates according to the Euler rotation parameters of the corresponding geological time slice, so that the plate outline returns to the paleogeographic location of the corresponding geological time slice; and connecting the rotated unit spherical coordinates according to the boundary order of the plate outline to obtain the reconstructed plate outline.

[0042] The reconstructed plate outline is closed by boundary closure, and the closed reconstructed plate outline is associated with the corresponding geological time slice and plate number to form the plate affiliation judgment range for each geological time slice.

[0043] The specific process includes traversing the boundary coordinate points in the reconstructed plate outline according to the order of the boundary points of the plate outline; connecting the endpoints when there is a gap between the starting boundary coordinate point and the ending boundary coordinate point; and arranging the order of the boundary coordinate points so that the reconstructed plate outline forms a closed structure; and associating the closed reconstructed plate outline with the corresponding geological time slice and plate number to form the plate affiliation judgment range for each geological time slice.

[0044] S3. Perform time-slice matching calculation on the spatiotemporal attribute chain of the ore deposit and the plate affiliation judgment range to determine the bound plates of each geological time slice and connect them in sequence according to geological time to obtain the dynamic binding chain of the ore deposit plate.

[0045] By using geological time slices to analyze the spatial location attributes of the corresponding spatiotemporal attribute nodes in the spatiotemporal attribute chain of the ore deposit, the location of the ore deposit corresponding to each geological time slice is obtained. Then, the location of the ore deposit is spatially matched with the plate affiliation judgment range of the same geological time slice to determine the candidate binding plate of each geological time slice.

[0046] The specific process includes: based on the time range of the geological time slice, sequentially locating the spatiotemporal attribute nodes corresponding to the geological time slice from the spatiotemporal attribute chain of the ore deposit, and performing spatial positioning attribute parsing on the located spatiotemporal attribute nodes to obtain the spatial location of the ore deposit in the spatiotemporal attribute nodes; unifying the coordinate format of the longitude, latitude, and elevation fields in the spatial location of the ore deposit, and verifying the completeness of the longitude, latitude, and elevation fields to ensure that the spatial location of the ore deposit is in the same space as the plate affiliation judgment range, thereby obtaining the ore deposit location corresponding to each geological time slice.

[0047] The plate affiliation judgment range under the same geological time slice is retrieved according to the geological time slice, and the location of the deposit is spatially matched with the reconstructed plate outline after closure in the plate affiliation judgment range to determine the plate affiliation judgment range corresponding to the location of the deposit. Then, the plate associated with the plate affiliation judgment range corresponding to the location of the deposit is determined as the candidate binding plate for each geological time slice.

[0048] For candidate bound plates in adjacent geological time slices, plate number comparison and plate outline proximity determination are performed to generate plate number consistency markers and outline proximity markers. The rotation direction angle is calculated based on the Euler rotation parameters corresponding to the candidate bound plates. Finally, a ore deposit plate binding status map is constructed based on the plate number consistency markers, outline proximity markers and rotation direction angles.

[0049] The specific process includes: grouping candidate bound plates from adjacent geological time slices according to geological time sequence to form adjacent candidate bound plate groups; comparing the plate numbers of all adjacent candidate bound plate groups for consistency; treating comparisons with the same plate number as consecutive plate numbers and comparisons with different plate numbers as variable plate numbers; generating a consistent plate number marker based on the consecutive and variable plate number cases; determining the proximity of reconstructed plate contours in all adjacent candidate bound plate groups; retrieving the reconstructed plate contours corresponding to candidate bound plates in adjacent geological time slices; obtaining the nearest point spherical angle distance between two reconstructed plate contours; and generating a contour proximity marker when the nearest point spherical angle distance is less than or equal to the contour proximity determination threshold.

[0050] The Euler rotation parameters corresponding to all adjacent candidate bound plate groups are decomposed and the rotation direction angle is calculated. The candidate bound plates in each geological time slice are used as state nodes. The plate number consistency mark, contour proximity mark, and rotation direction angle between adjacent geological time slices are used as connection attributes between state nodes. The state nodes and connection attributes are associated and connected according to the geological time sequence to obtain the ore deposit plate binding state map.

[0051] It should be noted that the threshold for determining the proximity of the contours is determined based on the sampling interval of the boundary coordinate points of the plate contours in the plate motion constraint data, the boundary closure error of the reconstructed plate contours, and the span of adjacent geological time slices. For example, the threshold is taken as five percent of the smaller value of the circumscribing radii of the two reconstructed plate contours.

[0052] The expression for calculating the angle between the rotation directions: ; in, Indicates the first The geological time slice and the first The angle between the rotational directions of candidate binding plates between geological timepieces. Indicating an index for a geological time slice, Indicates the first The rotation angle of each geological timepiece corresponds to the candidate bound plate. Indicates the first The rotation angle of each geological timepiece corresponds to the candidate bound plate. Indicates the first Each geological time slice corresponds to a unit rotation axis of the candidate bound plate. Indicates the first Each geological time slice corresponds to a unit rotation axis of the candidate bound plate. Indicates the first The directional symbol of the rotation angle of a geological time slice Indicates the first The directional symbol for the rotation angle of a geological time slice.

[0053] It should be noted that when When the value is greater than or equal to 0, the value is 1. When less than 0, the value is -1; when When the value is greater than or equal to 0, the value is 1. When the value is less than 0, it takes the value of -1; when the Euler rotation angle is equal to zero, the direction sign is set to 1, which is used to keep the vector norm non-zero when the unit rotation axis participates in the calculation of the rotation direction angle.

[0054] Based on the ore deposit plate binding status map, target binding paths that run through all geological time slices are selected. Then, candidate binding plates corresponding to each geological time slice in the target binding path are determined as binding plates. The binding plates are then connected in sequence according to geological time to obtain the dynamic binding chain of ore deposit plates.

[0055] The specific process includes tracing candidate binding plate connection paths that cover all geological time slices in the ore deposit plate binding status map along the geological time sequence, and determining the continuity of candidate binding plate connection paths based on the consistency of plate number markers, the proximity of contour markers, and the rotation direction angle. Candidate binding plate connection paths that do not meet the continuity determination are eliminated, and candidate binding plate connection paths with continuous plate numbers, spatially adjacent contours, and smooth rotation direction changes are retained. When there is a candidate binding plate connection path that spans all geological time slices and meets the continuity determination, the candidate binding plate connection path is determined as the target binding path. When there is no candidate binding plate connection path that spans all geological time slices and meets the continuity determination, the candidate binding plate connection paths are sorted according to the number of geological time slices covered, the number of contour proximity markers, and the size of the rotation direction angle. The candidate binding plate connection path with the most geological time slices covered, the most contour proximity markers, and the smallest rotation direction angle is selected as the second-best binding path, and the second-best binding path is determined as the target binding path.

[0056] Candidate binding plates corresponding to various geological time slices in the target binding path are identified as binding plates, and the binding plates are sequentially numbered. The binding plates are then connected in sequence according to geological time from ancient to modern. At the same time, the binding plate numbers and Euler rotation parameters corresponding to all geological time slices are extracted. The binding plate numbers and Euler rotation parameters are associated and matched with the corresponding geological time slices to obtain the dynamic binding chain of the ore deposit plates.

[0057] S4. Extract the Euler rotation parameters of the corresponding bound plates from the plate motion constraint data along the dynamic binding chain of the ore deposit plates, and use the current location data of the ore deposit as the simulation starting point. Obtain the initial paleogeographic coordinates through reverse rotation simulation, and then calculate the rotation constraint path of adjacent coordinates by combining the corresponding Euler rotation parameters.

[0058] The binding plates corresponding to each geological timepiece are determined one by one along the dynamic binding chain of the ore deposit plates. The binding plates are then matched with the corresponding plate numbers in the plate motion constraint data for the same geological timepiece to determine the Euler rotation parameters corresponding to the binding plates.

[0059] The specific process includes: extracting each geological time slice and its corresponding bound plate number according to the geological time sequence of the dynamic binding chain of the ore deposit plates; using each geological time slice and its corresponding bound plate number as a joint matching condition; filtering plate motion constraint data in the plate motion constraint data that have the same time mark as the geological time slice and the same bound plate number, and removing plate motion constraint data with inconsistent time marks and mismatched bound plate numbers; extracting the corresponding Euler rotation parameters from the filtered plate motion constraint data and performing field integrity verification; associating and binding the verified Euler rotation parameters with the geological time slice and bound plate number in the dynamic binding chain of the ore deposit plates, so that the bound plate number under each geological time slice corresponds to a unique Euler rotation parameter.

[0060] When multiple Euler rotation parameters correspond to the same geological time slice and the same bound plate number, the Euler rotation parameters are arranged according to the degree of correspondence between the time stamp and the geological time slice. Then, the completeness of the fields of Euler pole position, rotation angle, and rotation angle direction attributes is checked. Euler rotation parameters with inconsistent time stamps, missing Euler pole positions, missing rotation angles, and missing rotation angle direction attributes are removed. Among the remaining Euler rotation parameters, the Euler rotation parameter with the smallest difference between the time stamp and the center time point of the geological time slice is selected. When the difference between the time stamp and the center time point of the geological time slice is the same, the Euler rotation parameter with the highest values ​​of Euler pole position, rotation angle, and rotation angle direction attributes is selected. When the values ​​are the same, the Euler rotation parameter that ranks first in the plate motion constraint data is selected to obtain the Euler rotation parameters corresponding to the bound plates under each geological time slice.

[0061] The Euler rotation parameters corresponding to the bound plate are decomposed to obtain the Euler pole and the reverse rotation angle, and the reverse rotation matrix of the corresponding geological time slice is constructed based on the Euler pole and the reverse rotation angle. The specific process includes decomposing and validating the Euler rotation parameters corresponding to the bound plate to obtain the Euler pole and the reverse rotation angle in the Euler rotation parameters; determining the spherical rotation center based on the Euler pole, and performing matrix transformation using the reverse rotation angle as the spherical rotation angle to obtain the reverse rotation matrix of the corresponding geological time slice.

[0062] Using the current location data of the ore deposit as the starting point for simulation, the current location data of the ore deposit is converted into current unit spherical coordinates. Then, according to the geological time slices from present to past, the current unit spherical coordinates are sequentially input into the corresponding reverse rotation matrix for cumulative coordinate transformation, and the initial paleogeographic coordinates are output.

[0063] The specific process includes: using the current location data of the ore deposit as the starting point for simulation, verifying the coordinate format of the longitude and latitude in the current location data of the ore deposit, and converting the longitude and latitude into three-dimensional coordinates on a unit sphere to form the current unit sphere coordinates; according to the geological time slices from present to past, retrieving the reverse rotation matrix of each geological time slice, and inputting the current unit sphere coordinates into the reverse rotation matrix with the closest time sequence to the present for coordinate transformation, and then inputting the unit sphere coordinates after coordinate transformation into the reverse rotation matrix corresponding to the previous geological time slice, so that the unit sphere coordinates are continuously transmitted along the geological time slice sequence from present to past, and converting the unit sphere coordinates after cumulative coordinate transformation into longitude and latitude to output the initial paleogeographic coordinates.

[0064] Adjacent initial paleogeographic coordinates are paired in geological time sequence and converted into spherical starting and ending points. The unit rotation axis is determined based on the Euler pole, and the reverse rotation angle is inverted to generate the rotation angle.

[0065] The specific process includes: organizing the initial paleogeographic coordinates according to the geological time sequence, taking two initial paleogeographic coordinates corresponding to adjacent geological time slices as a set of adjacent initial paleogeographic coordinates; converting the initial paleogeographic coordinates corresponding to the earlier geological time in a set of adjacent initial paleogeographic coordinates into a spherical starting point, and converting the initial paleogeographic coordinates corresponding to the later geological time into a spherical ending point; extracting the rotation axis from the Euler poles corresponding to adjacent geological time slices, and converting the Euler poles into a unit rotation axis, so that the unit rotation axis corresponds to the spherical starting point and the spherical ending point; and inverting the reverse rotation angle to restore the reverse rotation angle to the angle direction used for the rotational connection between adjacent initial paleogeographic coordinates, thus generating the rotation angle.

[0066] Starting from the spherical origin, the system rotates segmentally around the unit rotation axis according to the rotation angle, and uses the spherical endpoint as the termination constraint of the segmented rotation path. The segmented rotation path is then normalized and resampled to calculate the rotation constraint path between adjacent coordinates.

[0067] The specific process includes: using the spherical starting point as the starting position of the segmented rotation path; determining the spherical rotation direction based on the unit rotation axis; dividing the spherical rotation trajectory corresponding to the unit rotation axis along the spherical starting point into continuous rotation segments according to the rotation angle; and determining the corresponding spherical positions along the continuous rotation segments to extend the segmented rotation path from the spherical starting point to the spherical ending point; using the spherical ending point as the termination constraint of the segmented rotation path; constraining and correcting the positional difference between the end of the segmented rotation path and the spherical ending point to ensure the continuity of the segmented rotation path between the spherical starting point and the spherical ending point; measuring the path progress of the segmented rotation path after completing the termination constraint; generating normalized segmented parameters based on the relative positions of each sampled position between the spherical starting point and the spherical ending point in the segmented rotation path; resampling the segmented rotation path according to the normalized segmented parameters; arranging the resampled path points continuously along the segmented rotation path; and calculating the adjacent coordinate rotation constraint path.

[0068] The expression for calculating the adjacent coordinate rotation constraint path: ; in, Indicates the first Adjacent initial paleogeographic coordinates in normalized piecewise parameters The corresponding rotation constraint path points, Indices representing adjacent initial paleogeographic coordinates. This represents the normalized piecewise parameters of the adjacent coordinate rotation constraint path. Indicates the first The spherical starting point in a set of adjacent initial paleogeographic coordinates. Indicates the first The unit rotation axis of the bound plate corresponds to the initial paleogeographic coordinates of adjacent groups. Indicates the first The rotation angle of the bound plate corresponding to the initial paleogeographic coordinates of the adjacent group.

[0069] It should be noted that the normalized segment parameters are obtained by measuring the path progress of the segmented rotation path after the termination constraint is completed, and by unifying each sampling position between the spherical start point and the spherical end point to the same segment scale according to the path progress. They are used to mark the relative position of each sampling position in the segmented rotation path.

[0070] S5. Based on the adjacent coordinate rotation constraint path, perform path continuity correction and constraint error correction on the initial paleogeographic coordinates to obtain the corrected paleogeographic coordinates, and then generate a set of paleogeographic coordinates of the ore deposit evolution trajectory according to the geological time sequence.

[0071] The initial paleogeographic coordinates are converted into initial unit spherical coordinates, and the corresponding adjacent coordinate rotation constraint paths are sampled at equal intervals to obtain path sampling points. Then, the spherical angular distance between the initial unit spherical coordinates and each path sampling point is calculated, and the path sampling point with the smallest spherical angular distance is determined as the path projection point.

[0072] The specific process includes: organizing the initial paleogeographic coordinates corresponding to each geological time slice into coordinate formats according to longitude and latitude; converting the formatted initial paleogeographic coordinates into initial unit spherical coordinates; then, according to the geological time slice, associating and matching the initial unit spherical coordinates with the corresponding adjacent coordinate rotation constraint paths to obtain the total spherical angle length of the adjacent coordinate rotation constraint paths from the spherical start point to the spherical end point; determining the number of sampling points based on the total spherical angle length of the path and the longitude and latitude of the initial paleogeographic coordinates; dividing the total spherical angle length of the path equally according to the number of sampling points to obtain the sampling step size; then, sampling adjacent coordinate rotation constraint paths at equal intervals according to the sampling step size to obtain path sampling points; calculating the spherical angle distance between the initial unit spherical coordinates and each path sampling point one by one; comparing and filtering according to the spherical angle distance in ascending order, and determining the path sampling point with the smallest spherical angle distance as the path projection point.

[0073] The expression for calculating the spherical angular distance between the initial unit spherical coordinates and each path sampling point is as follows: ; in, Indicates the first The initial unit spherical coordinates and the... spherical angular distance between each path sampling point Indicates the index of the initial unit spherical coordinates. Indicates the index of the path sampling point. Indicates the first An initial unit spherical coordinate system to be corrected. Indicates the first The first initial unit spherical coordinate corresponds to the adjacent coordinate rotation constraint path on the path of the first... Each path sampling point.

[0074] The spherical angle between the path projection point and the initial unit spherical coordinates is measured and the tangential direction is decomposed to obtain the spherical angular distance and the path deviation direction. Then, amplitude limiting normalization and amplitude synthesis are performed to obtain the constraint error correction amount.

[0075] The specific process includes: measuring the spherical angle between the path projection point and the initial unit spherical coordinates to obtain the spherical angular distance between the path projection point and the initial unit spherical coordinates; determining the path tangential direction based on the preceding and following path sampling points of the path projection point on the adjacent coordinate rotation constraint path, and decomposing the deviation direction of the initial unit spherical coordinates relative to the path projection point into the path tangential direction to obtain the path deviation direction; limiting and normalizing the spherical angular distance to ensure that the spherical angular distance is within the allowable range of path continuity correction; and combining the amplitude-limited and normalized spherical angular distance with the path deviation direction to obtain the constraint error correction amount.

[0076] It should be noted that the path continuity correction range is determined based on the total spherical angle length of the adjacent coordinate rotation constraint paths, the path sampling step size, and the longitude and latitude accuracy of the initial paleogeographic coordinates. The lower limit of the path continuity correction range is zero, and the upper limit of the path continuity correction range is three percent of the total spherical angle length of the adjacent coordinate rotation constraint paths. When the spherical angle distance is less than the lower limit of the path continuity correction range, the spherical angle distance is corrected to the lower limit of the path continuity correction range. When the spherical angle distance is greater than the upper limit of the path continuity correction range, the spherical angle distance is corrected to the upper limit of the path continuity correction range. When the spherical angle distance is within the path continuity correction range, the spherical angle distance is retained.

[0077] The initial unit spherical coordinates are corrected for path continuity based on the constraint error correction amount, and the corrected unit spherical coordinates are then converted into corrected paleogeographic coordinates.

[0078] The specific process includes: performing path continuity correction on the initial unit spherical coordinates based on the constraint error correction amount, so that the initial unit spherical coordinates move closer to the adjacent coordinate rotation constraint path according to the path deviation direction and the spherical angular distance after amplitude limiting and normalization; and normalizing the initial unit spherical coordinates after path continuity correction to obtain the corrected unit spherical coordinates; and converting the corrected unit spherical coordinates to longitude and latitude to obtain the corrected paleogeographic coordinates.

[0079] The corrected paleogeographic coordinates are merged with the corresponding geological time slices, bound plate numbers, adjacent coordinate rotation constraint paths, and constraint error corrections to form paleogeographic coordinate nodes.

[0080] The specific process includes: verifying the corrected paleogeographic coordinates by time stamping according to the geological time slices to ensure that the corrected paleogeographic coordinates are consistent with the corresponding geological time slices; extracting the bound plate number under the corresponding geological time slice from the dynamic binding chain of the ore deposit plate; extracting the constraint error correction amount from the path continuity correction process; and merging the corrected paleogeographic coordinates, corresponding geological time slices, bound plate numbers, adjacent coordinate rotation constraint paths, and constraint error correction amounts according to the same geological time slice to form paleogeographic coordinate nodes.

[0081] The paleogeographic coordinate nodes are sorted in geological time sequence, and adjacency markers and trajectory numbers are added to adjacent paleogeographic coordinate nodes after sorting to form a chain of mineral deposit evolution trajectory nodes.

[0082] The specific process includes: verifying and arranging the geological time of paleogeographic coordinate nodes according to the corresponding geological time slices recorded in the paleogeographic coordinate nodes, so that each paleogeographic coordinate node is arranged from ancient to modern geological time; among the sorted paleogeographic coordinate nodes, adjacent paleogeographic coordinate nodes are determined according to their arrangement positions, and adjacency markers are added to adjacent paleogeographic coordinate nodes, so that each paleogeographic coordinate node is connected to adjacent paleogeographic coordinate nodes in geological time order; then, trajectory numbers are added to the paleogeographic coordinate nodes in geological time order, so that each paleogeographic coordinate node has a corresponding arrangement position in the mineral deposit evolution trajectory; and the paleogeographic coordinate nodes are connected in series based on the adjacency markers and trajectory numbers to form a mineral deposit evolution trajectory node chain.

[0083] The node fields of the ore deposit evolution trajectory node chain are extracted and the coordinate records are encapsulated. Then, the time series index is arranged according to the trajectory number to obtain the paleogeographic coordinate set of the ore deposit evolution trajectory.

[0084] The specific process includes: extracting paleogeographic coordinate nodes item by item along the node chain of the ore deposit evolution trajectory according to the trajectory sequence number; then extracting geological time slices, corrected paleogeographic coordinates, bound plate numbers, adjacent coordinate rotation constraint paths, constraint error correction amounts, adjacency markers, and trajectory sequence numbers from the paleogeographic coordinate nodes; verifying the time stamps of the extracted geological time slices; verifying the coordinate format of the extracted corrected paleogeographic coordinates; and verifying the correspondence of the extracted bound plate numbers, adjacent coordinate rotation constraint paths, constraint error correction amounts, adjacency markers, and trajectory sequence numbers; encapsulating the verified node fields according to the coordinate record format, and using the trajectory sequence number as the basis for time series indexing, so that the coordinate records are arranged continuously in geological time order, thus obtaining the paleogeographic coordinate set of the ore deposit evolution trajectory.

[0085] In summary, this invention achieves continuous determination of the bound plates in various geological time slices and avoids attribute jumps near plate boundaries by constructing a ore deposit plate binding state map by combining the ore deposit location, plate attribution judgment range, and rotation direction angle. This improves the stability of deep-time plate attribution and the reliability of paleogeographic coordinate inversion. Furthermore, by generating initial paleogeographic coordinates through an inverse rotation matrix and performing path continuity correction by combining adjacent coordinate rotation constraint paths, this invention achieves synergistic constraints between paleocoordinate inversion and plate rotation geometry, thereby suppressing trajectory jumps, reducing error accumulation, and improving the verifiability of ore deposit evolution trajectories.

[0086] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints, characterized in that, include: Collect ore evolution constraint data of the target ore deposit and perform spatiotemporal geological normalization processing, then connect them in sequence according to geological time to form a spatiotemporal attribute chain of the ore deposit; Geological time slices are divided based on the spatiotemporal attribute chain of the ore deposit, and plate movement constraint data and current ore deposit location data are extracted from the ore deposit evolution constraint data. Then, time slice matching and plate outline reconstruction are performed on the plate movement constraint data to form the plate affiliation judgment range of each geological time slice. Time-slice matching calculations are performed on the spatiotemporal attribute chain of the ore deposit and the plate affiliation judgment range to determine the bound plates of each geological time slice and connect them in geological time sequence to obtain the dynamic binding chain of the ore deposit plate. The Euler rotation parameters of the corresponding bound plates are extracted from the plate motion constraint data along the dynamic binding chain of the ore deposit plates. The current location data of the ore deposit is used as the simulation starting point. The initial paleogeographic coordinates are obtained through reverse rotation simulation. The rotation constraint paths of adjacent coordinates are then calculated in combination with the corresponding Euler rotation parameters. The initial paleogeographic coordinates are corrected for path continuity and constraint error by rotating the adjacent coordinates to obtain the corrected paleogeographic coordinates. Then, a set of paleogeographic coordinates of the ore deposit evolution trajectory is generated in geological time sequence.

2. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The ore deposit evolution constraint data includes current ore deposit location data, geological age constraint data, paleomagnetic constraint data, and plate tectonics constraint data.

3. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps for forming the spatiotemporal attribute chain of the ore deposit are as follows: Collect ore evolution constraint data of the target ore deposit, and perform spatiotemporal benchmark unification and outlier data cleaning on the ore evolution constraint data to obtain normalized ore evolution constraint data; Based on the geological time corresponding to the normalized ore deposit evolution constraint data, the normalized ore deposit evolution constraint data is grouped by time, and attribute merging, time stamp binding and node encapsulation are performed on the normalized ore deposit evolution constraint data under the same geological time to obtain spatiotemporal attribute nodes; The spatiotemporal attribute nodes are sorted according to geological time from ancient to the present, and a succession number is added to adjacent spatiotemporal attribute nodes. Then, the spatiotemporal attribute nodes are connected in series according to the succession number to form the spatiotemporal attribute chain of the ore deposit.

4. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps for extracting plate tectonics constraint data and current ore deposit location data from ore deposit evolution constraint data are as follows: Extract time nodes arranged in geological time sequence from the spatiotemporal attribute chain of the ore deposit, and divide the time interval between adjacent time nodes into geological time slices; Using geological time slices as search criteria, plate tectonics constraint data and current location data of ore deposits are extracted from ore deposit evolution constraint data.

5. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps for determining the tectonic plate affiliation range of the various geological timepieces are as follows: Based on geological time slices, plate movement constraint data are retrieved by time to obtain the plate number, plate outline, and Euler rotation parameters corresponding to each geological time slice. The plate contours under the same plate number are reconstructed by rotating the Euler rotation parameters corresponding to each geological time piece to obtain the reconstructed plate contours. The reconstructed plate outline is closed by boundary closure, and the closed reconstructed plate outline is associated with the corresponding geological time slice and plate number to form the plate affiliation judgment range for each geological time slice.

6. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps for obtaining the dynamic binding chain of the mineral deposit segment are as follows: By using geological time slices to analyze the spatial location attributes of the corresponding spatiotemporal attribute nodes in the spatiotemporal attribute chain of the ore deposit, the ore deposit location corresponding to each geological time slice is obtained. Then, the ore deposit location is spatially matched with the plate affiliation judgment range of the same geological time slice to determine the candidate binding plate of each geological time slice. For candidate bound plates in adjacent geological time slices, plate number comparison and plate outline proximity determination are performed to generate plate number consistency markers and outline proximity markers. The rotation direction angle is calculated based on the Euler rotation parameters corresponding to the candidate bound plates. Then, a ore deposit plate binding status map is constructed based on the plate number consistency markers, outline proximity markers and rotation direction angles. Based on the ore deposit plate binding status map, target binding paths that run through all geological time slices are selected. Then, candidate binding plates corresponding to each geological time slice in the target binding path are determined as binding plates. The binding plates are then connected in sequence according to geological time to obtain the dynamic binding chain of ore deposit plates.

7. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps for obtaining the initial paleogeographic coordinates through reverse rotation simulation are as follows: The binding plates corresponding to each geological time piece are determined one by one along the dynamic binding chain of the ore deposit plates. The binding plates are then matched with the plate numbers corresponding to the same geological time piece in the plate motion constraint data to determine the Euler rotation parameters corresponding to the binding plates. The Euler rotation parameters corresponding to the bound plate are decomposed to obtain the Euler pole and the reverse rotation angle, and the reverse rotation matrix of the corresponding geological time slice is constructed based on the Euler pole and the reverse rotation angle. Using the current location data of the ore deposit as the starting point for simulation, the current location data of the ore deposit is converted into current unit spherical coordinates. Then, according to the geological time slices from present to past, the current unit spherical coordinates are sequentially input into the corresponding reverse rotation matrix for cumulative coordinate transformation, and the initial paleogeographic coordinates are output.

8. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 7, characterized in that, The specific steps for calculating the adjacent coordinate rotation constraint path are as follows: Adjacent initial paleogeographic coordinates are paired in geological time sequence and converted into spherical starting and ending points. The unit rotation axis is determined based on the Euler pole, and the reverse rotation angle is inverted to generate the rotation angle. Starting from the spherical origin, the system rotates segmentally around the unit rotation axis according to the rotation angle, and uses the spherical endpoint as the termination constraint of the segmented rotation path. The segmented rotation path is then normalized and resampled to calculate the rotation constraint path between adjacent coordinates.

9. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 1, characterized in that, The specific steps to obtain the corrected paleogeographic coordinates are as follows: The initial paleogeographic coordinates are converted into initial unit spherical coordinates, and the corresponding adjacent coordinate rotation constraint paths are sampled at equal intervals to obtain path sampling points. Then, the spherical angular distance between the initial unit spherical coordinates and each path sampling point is calculated, and the path sampling point with the smallest spherical angular distance is determined as the path projection point. The spherical angle between the path projection point and the initial unit spherical coordinates is measured and the tangential direction is decomposed to obtain the spherical angular distance and the path deviation direction. Then, the amplitude is limited and normalized and the amplitude is synthesized to obtain the constraint error correction amount. The initial unit spherical coordinates are corrected for path continuity based on the constraint error correction amount, and the corrected unit spherical coordinates are then converted into corrected paleogeographic coordinates.

10. The method for generating paleogeographic coordinates of ore deposit evolution trajectories based on plate tectonics constraints as described in claim 9, characterized in that, The specific steps for generating the paleogeographic coordinate set of ore deposit evolution trajectories according to geological time sequence are as follows: The corrected paleogeographic coordinates are merged with the corresponding geological time slices, bound plate numbers, adjacent coordinate rotation constraint paths, and constraint error corrections to form paleogeographic coordinate nodes. The paleogeographic coordinate nodes are sorted in geological time sequence, and adjacency markers and trajectory numbers are added to adjacent paleogeographic coordinate nodes after sorting to form a chain of mineral deposit evolution trajectory nodes. The node fields of the ore deposit evolution trajectory node chain are extracted and the coordinate records are encapsulated. Then, the time series index is arranged according to the trajectory number to obtain the paleogeographic coordinate set of the ore deposit evolution trajectory.

Images