[go: up one dir, main page]

CN118568446B - Comprehensive geological exploration information management system - Google Patents

Comprehensive geological exploration information management system Download PDF

Info

Publication number
CN118568446B
CN118568446B CN202411050620.6A CN202411050620A CN118568446B CN 118568446 B CN118568446 B CN 118568446B CN 202411050620 A CN202411050620 A CN 202411050620A CN 118568446 B CN118568446 B CN 118568446B
Authority
CN
China
Prior art keywords
geological
data
formula
dynamic
key
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202411050620.6A
Other languages
Chinese (zh)
Other versions
CN118568446A (en
Inventor
郭有劲
董明华
张诏飞
刘志亮
席伟
吴昱诚
王海峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway Resources Group Survey And Design Co ltd
Original Assignee
China Railway Resources Group Survey And Design Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Railway Resources Group Survey And Design Co ltd filed Critical China Railway Resources Group Survey And Design Co ltd
Priority to CN202411050620.6A priority Critical patent/CN118568446B/en
Publication of CN118568446A publication Critical patent/CN118568446A/en
Application granted granted Critical
Publication of CN118568446B publication Critical patent/CN118568446B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及地质数据管理技术领域,具体为一种综合地质勘探信息管理系统,系统包括:非线性动态数据建模模块将地层振动与压力变化的时间序列数据转换为多维相空间,通过延迟坐标映射地质内部动态结构,分析重构的多维相空间中的关键转折点,生成地质动态吸引子图。本发明中,通过非线性动态数据建模,能有效捕捉地质内部的动态变化,分析和预测地质结构的关键转折点,构建基于物理和化学特性的岩层连续性网络,使得对岩层间的连续性和潜在断层的监测更为直观和实时,进一步强化了对地质结构稳定性的评估,应用高级降维技术对地质数据进行处理,有效减少了数据处理的复杂度同时保留了关键信息,从而使得地质特征的提取更加高效和准确。

The present invention relates to the field of geological data management technology, specifically to a comprehensive geological exploration information management system, the system comprising: a nonlinear dynamic data modeling module converts the time series data of stratum vibration and pressure change into a multidimensional phase space, maps the internal dynamic structure of the geology through delayed coordinates, analyzes the key turning points in the reconstructed multidimensional phase space, and generates a geological dynamic attractor diagram. In the present invention, through nonlinear dynamic data modeling, the dynamic changes inside the geology can be effectively captured, the key turning points of the geological structure can be analyzed and predicted, and a rock layer continuity network based on physical and chemical properties can be constructed, so that the monitoring of the continuity between rock layers and potential faults is more intuitive and real-time, and the assessment of the stability of the geological structure is further strengthened. The geological data is processed using advanced dimensionality reduction technology, which effectively reduces the complexity of data processing while retaining key information, thereby making the extraction of geological features more efficient and accurate.

Description

Comprehensive geological exploration information management system
Technical Field
The invention relates to the technical field of geological data management, in particular to a comprehensive geological exploration information management system.
Background
The field of geological data management relates to the collection, storage, processing, analysis and interpretation of data related to earth structure, composition and history. Such data typically originate from geological exploration, such as seismic surveying, drilling and surface sampling. The key task in the field is to systemize and digitize a large amount of geological information to support various applications such as mineral exploration, oil and gas field development, geological disaster prevention and the like. Geological data management systems can help scientists and engineers better understand subsurface structures and material distribution through maps, three-dimensional models, and other analytical tools. In addition, the technical field also comprises visualization, sharing and updating of data, and accuracy and instantaneity of the data in scientific research and business decision are ensured.
The comprehensive geological exploration information management system refers to an integrated technical solution, and is used for managing and analyzing various data generated in a geological exploration process. The main purpose of such systems is to provide data support for geological exploration, including seismic data, formation data, sample analysis results, etc., so that geologists and exploration engineers can effectively assess the distribution and development potential of subsurface resources. The system typically includes data acquisition, storage, processing, and visualization functions that enable a user to create comprehensive geologic models and conduct complex data analysis to guide actual exploration and development activities.
The prior art often faces the problems of slow processing speed and insufficient accuracy when processing large-scale geological data. Particularly when analyzing seismic data and formation information, conventional methods may not be able to effectively identify minor variations and complex patterns in the data, thereby affecting the accuracy of geologic structure and resource assessment. In addition, the lack of effective dynamic monitoring tools results in inadequate formation continuity and identification of potential faults, which in turn affects the decision-making efficiency of mining strategies and disaster prevention. This technical limitation is particularly evident in emergency situations, such as when a geological disaster is imminent, the response time and processing power of existing systems may be insufficient to provide the necessary support, increasing the difficulty of risk management.
Disclosure of Invention
The invention aims to solve the defects in the prior art and provides a comprehensive geological exploration information management system.
In order to achieve the above purpose, the present invention adopts the following technical scheme: a comprehensive geological survey information management system comprising:
The nonlinear dynamic data modeling module converts time series data of stratum vibration and pressure change into a multidimensional phase space, maps a geological internal dynamic structure through a delay coordinate, analyzes key turning points in the reconstructed multidimensional phase space and generates a geological dynamic attractor graph;
The stratum continuity dynamic evaluation module utilizes the geological dynamic attractor graph and combines physical and chemical characteristics of rocks to construct and update nodes and edges of an undirected graph or a directed graph in real time, monitors stratum continuity and potential faults, and evaluates stability indexes of a network by calculating shortest paths among the nodes in real time and identifying strong communication components to generate a stratum continuity network graph;
The high-dimensional dimension reduction analysis module utilizes the structural features in the stratum continuity network diagram, applies local linear embedding and equidistant mapping to reduce dimension of the high-dimensional geological data, extracts key geological features and establishes a dimension reduction geological feature set;
And the geological change mode prediction module analyzes the key geological structure and the dynamic change mode based on the dimensionality reduction geological feature set, identifies the change trend of the geological structure, predicts the development trend of a geological event and generates a geological event prediction analysis result.
As a further scheme of the invention, the acquisition steps of the geology dynamic attractor graph specifically comprise:
Converting the time series data of formation vibration and pressure change into a multi-dimensional data point set, and using the formula:
Mapping the time sequence data into a high-dimensional space to generate a multi-dimensional data point set;
wherein, Representative time of dayIs a function of the data of (a),In order to provide for the time interval of time,Is the dimension number;
using the multi-dimensional data point set, adopting the formula:
Calculating Euclidean distance between data points to generate a distance matrix;
wherein, Representing data pointsAndThe distance between the two plates is set to be equal,As the number of dimensions to be used,Is a time interval;
and screening key turning points from the distance matrix, wherein the formula is adopted:
Acquiring a key turning point set;
wherein, Representing the set of data points that are selected,Representation pointsPoint of attachmentThe distance between the two electrodes is equal to the distance between the two electrodes,For each pair of pointsAndThe weight between the two is set according to the relative importance of the data points or the connection strength;
based on the set of critical turning points, the formula is used:
constructing a geological dynamic attractor graph;
wherein, A quantitative description of the geostatistical dynamic attractor graph is presented,As a function of the data points,Is a scaling parameter for adjusting the influence of the distance, making the attraction force more sensitive or more suppressed.
As a further aspect of the present invention, the step of obtaining the stratum continuity network map specifically includes:
According to the physical and chemical characteristics of the geology dynamic attractor graph and the rock, fusing the attractor graph and the characteristic data, and adopting the formula:
generating a preliminary stratum network map;
wherein, Representing a preliminary stratigraphic network map,On behalf of the node(s),The representative edge of the sheet is represented by,Representing a dynamic attractor graph of the geology,Representing the data of the physical characteristics,Representing chemical property data;
calculating the shortest path between nodes by using the preliminary stratum network diagram, and adopting the formula:
Calculating shortest paths among nodes, and generating a shortest path analysis result among the nodes;
wherein, Representative nodeAndThe shortest path between the two paths is defined,The weight of the representative edge is calculated,Representing the adjustment parameters;
and identifying strong connected components in the network by using the shortest path analysis result between the nodes, and adopting the formula:
Wherein the method comprises the steps of Generating a strong connected component analysis result as a threshold value;
wherein, Representing a strong connected component of the signal,On behalf of the node(s),AndRepresenting the nodes of the differentiation and,Representing a threshold value;
integrating the shortest path analysis result and the strong communication component analysis result among the nodes, and adopting the formula:
generating a stratum continuity network diagram;
wherein, A network map representing the continuity of the formation,Representing a strong connected component of the signal,Representing the shortest path between the nodes,Representing a network stability indicator.
As a further scheme of the invention, the step of acquiring the dimension-reducing geological feature set specifically comprises the following steps:
by utilizing the structural characteristics in the stratum continuity network diagram, a local linear embedding method is applied, and the formula is adopted:
Calculating local neighborhood embedding to generate a local embedding mapping result;
wherein, Representing the weight of the object to be weighed,Representing data points;
and combining the local embedded mapping result, using an equidistant mapping method, and adopting the following formula:
Optimizing the distance mapping of the global geological data to generate an equidistant mapping result;
wherein, Representing a core distance adjustment factor;
Integrating the local embedded mapping result and the equidistant mapping result, and using the formula:
Optimizing the dimension reduction process to generate a dimension reduction result;
wherein, Represents an adjustment coefficient;
extracting key geological features from the dimension reduction result, and using the formula:
generating a dimension-reducing geological feature set;
wherein, Representing a geologic feature.
As a further scheme of the invention, the geological event prediction analysis result obtaining step specifically comprises the following steps:
Based on the dimensionality reduction geological feature set, analyzing a key geological structure, and adopting the formula:
Calculating the significance of the key geological structure, and generating a key geological structure analysis result;
wherein, Representing the result of the analysis of the key geological structure,The correlation coefficient is represented by a correlation coefficient,The characteristic item is represented by a characteristic term,Representing the number of features;
According to the analysis result of the key geological structure, analyzing a dynamic change mode, and adopting the formula:
quantifying the geologic structure change trend to generate a dynamic change mode analysis result;
wherein, Representing the result of the analysis of the dynamic change pattern,The rate of change is indicated as being indicative of,Representing the result of the analysis of the key geological structure,Representing the number of features;
And using the dynamic change mode analysis result to adopt the formula:
Determining the most obvious change trend, and generating a geological structure change trend identification result;
wherein, Representing the result of identifying the change trend of the geological structure,Representing trend intensity;
based on the geological structure change trend identification result, the development trend of geological events is predicted, and the formula is adopted:
integrating all the variation trends to generate a geological event prediction analysis result;
wherein, Representing the result of the predictive analysis of the geological event,The prediction parameters are represented by a set of parameters,Representing the result of identifying the change trend of the geological structure,Representing the number of features.
Compared with the prior art, the invention has the advantages and positive effects that:
According to the invention, through nonlinear dynamic data modeling, dynamic changes in the geology can be effectively captured, and key turning points of the geological structure can be analyzed and predicted. And a stratum continuity network based on physical and chemical characteristics is constructed, so that the continuity between strata and the monitoring of potential faults are more visual and real-time, and the evaluation of the stability of a geological structure is further enhanced. The advanced dimension reduction technology is applied to process the geological data, so that the complexity of data processing is effectively reduced, and key information is reserved, thereby enabling the extraction of geological features to be more efficient and accurate. The processing means not only optimizes the data analysis flow, but also enhances the prediction capability of the geological event development trend, and provides more scientific data support for geological exploration and related decisions.
Drawings
FIG. 1 is a system flow diagram of the present invention;
FIG. 2 is a flowchart of the steps for obtaining a geostatistical dynamic attractor graph according to the present invention;
FIG. 3 is a flowchart of the steps for obtaining a map of formation continuity according to the present invention;
FIG. 4 is a flowchart of the steps for acquiring a reduced-dimension set of geologic features of the present invention;
FIG. 5 is a flowchart illustrating the steps for obtaining the result of the analysis of a prediction of a geological event according to the present invention.
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.
Examples
Referring to fig. 1, a comprehensive geological exploration information management system includes:
The nonlinear dynamic data modeling module converts time series data of stratum vibration and pressure change into a multidimensional phase space, maps a geological internal dynamic structure through a delay coordinate, analyzes key turning points in the reconstructed multidimensional phase space and generates a geological dynamic attractor graph;
The stratum continuity dynamic evaluation module utilizes a geological dynamic attractor graph and combines physical and chemical characteristics of rocks to construct and update nodes and edges of an undirected graph or a directed graph in real time, monitors stratum continuity and potential faults, calculates the shortest path between the nodes in real time and identifies strong communication components, evaluates stability indexes of a network, and generates a stratum continuity network graph;
The high-dimensional geological data dimension reduction analysis module utilizes structural features in the stratum continuity network diagram, applies local linear embedding and equidistant mapping to reduce dimension of the high-dimensional geological data, extracts key geological features and establishes a dimension reduction geological feature set;
The geological change mode prediction module analyzes the key geological structure and the dynamic change mode based on the dimensionality reduction geological feature set, recognizes the change trend of the geological structure, predicts the development trend of a geological event, and generates a geological event prediction analysis result.
The geological dynamic attractor graph specifically comprises key turning points, multidimensional phase space and delay coordinate mapping, the stratum continuity network graph specifically comprises nodes, edges and strata continuity of an undirected graph or a directed graph and potential faults, the dimension-reducing geological feature set specifically comprises key geological features, high-dimensional geological data and structural features, and the geological event prediction analysis result specifically comprises key geological structures, dynamic change modes and geological event development trends.
Referring to fig. 2, the steps for obtaining the geosteering subgraph specifically include:
Converting the time series data of formation vibration and pressure change into a multi-dimensional data point set, and using the formula:
Mapping the time sequence data into a high-dimensional space to generate a multi-dimensional data point set;
wherein, Representative time of dayIs a function of the data of (a),In order to provide for the time interval of time,Is the dimension number;
using a multidimensional data point set, adopting the formula:
Calculating Euclidean distance between data points to generate a distance matrix;
wherein, Representing data pointsAndThe distance between the two plates is set to be equal,As the number of dimensions to be used,Is a time interval;
and (3) selecting key turning points from the distance matrix, and adopting the formula:
Acquiring a key turning point set;
wherein, Representing the set of data points that are selected,Representation pointsPoint of attachmentThe distance between the two electrodes is equal to the distance between the two electrodes,For each pair of pointsAndThe weight between the two is set according to the relative importance of the data points or the connection strength;
based on the set of key turning points, the formula is used:
constructing a geological dynamic attractor graph;
wherein, A quantitative description of the geostatistical dynamic attractor graph is presented,As a function of the data points,Is a scaling parameter for adjusting the influence of the distance, making the attraction force more sensitive or more suppressed.
The formula:
Detailed description of parameters
: At a single momentIs assumed to be the intensity of the formation vibrations.
: Time intervals are used to define sampling intervals of data points in a time series.
: The embedding dimension, i.e. the number of delay time sequences included in the vector.
Calculation example
Assume thatIntensity value, time for seismic monitoringIs 5.0, time intervalSet to 1 hour, embed dimensionSet to 3.
The calculation flow is as follows:
In the time-course of which the first and second contact surfaces,
When it is assumed that
When it is assumed that
Thus, the first and second substrates are bonded together,
This vector reflects the slave timeSeismic intensity values for the first three consecutive time points.
The formula:
Detailed description of parameters
The starting point of the time series is differentiated.
The dimensions are embedded.
Calculation example
Continued use of dataAnd
The calculation flow is as follows:
Distance value represents a vector AndThe euclidean distance between the two time series points reflects the degree of fitting between the two time series points.
The formula:
detailed description of parameters:
: the selected set of data points.
: Data pointsAndDistance between them.
: Each pair of pointsAndThe weight between the two is set according to the relative importance or the connection strength.
Calculation example
Suppose three data points are selectedAnd assuming weights for these points1 (I.e., the relationship between all points is equally important).
From already calculatedAssume that
The calculation flow is as follows:
this value represents the set of data points Weighted distance sum between all points in the model. Selecting a set of critical turning points that minimizes the valueIn this example, because the weights are the same, the sum of all points is the same, and the key turning point set is directly obtained
The formula:
Detailed description of parameters
: Quantitative description of the geosynamic attractor graph.
: Data points.
: And the scaling parameter is used for adjusting the influence of the distance.
: Data pointsAndDistance between them.
Calculation example
Assume thatSet to 1, known asAnd assume their distance as follows:
The calculation flow is as follows:
For each point The attractor value is calculated.
Calculation ofAttractor value of (2):
Attractor value:
Calculation of Attractor value of (2):
Attractor value:
Calculation of Attractor value of (2):
Attractor value:
Calculating a geostationary subgraph :
This value isThe total attractive force of the geology dynamic attractor graph is represented and reflects the comprehensive attractive force of all key turning points.
Referring to fig. 3, the formation continuity network map is obtained by the steps of:
According to the physical and chemical characteristics of the geological dynamic attractor graph and the rock, fusing the attractor graph and the characteristic data, and adopting the formula:
generating a preliminary stratum network map;
wherein, Representing a preliminary stratigraphic network map,On behalf of the node(s),The representative edge of the sheet is represented by,Representing a dynamic attractor graph of the geology,Representing the data of the physical characteristics,Representing chemical property data;
calculating the shortest path between nodes by using the preliminary stratum network diagram, and adopting the formula:
Calculating shortest paths among nodes, and generating a shortest path analysis result among the nodes;
wherein, Representative nodeAndThe shortest path between the two paths is defined,The weight of the representative edge is calculated,Representing the adjustment parameters;
and identifying strong connected components in the network by utilizing the shortest path analysis result among the nodes, and adopting the formula:
Wherein the method comprises the steps of Generating a strong connected component analysis result as a threshold value;
wherein, Representing a strong connected component of the signal,On behalf of the node(s),AndRepresenting the nodes of the differentiation and,Representing a threshold value;
integrating the shortest path analysis result and the strong connected component analysis result among the nodes, and adopting the formula:
generating a stratum continuity network diagram;
wherein, A network map representing the continuity of the formation,Representing a strong connected component of the signal,Representing the shortest path between the nodes,Representing a network stability indicator.
The formula:
Wherein:
: representing nodes.
: Edge set, by attractor graphAnd chemical properties of rockAnd (5) deriving.
: The geosynamic attractor graph, assuming calculated from previous steps.
: A function combined with physicsAnd chemistryCharacteristics to define the weights of the edges.
The calculation process comprises the following steps:
assumed geology dynamic attractor graph Comprising 3 nodes, each representing a rock sample of a site. Connectivity between nodes is defined by physical distanceDegree of chemical fittingAnd (5) determining. The weights of the edges are calculated from the feature synthesis. For example:
physical distance between node 1 and node 2 Unit (B)Chemical similarity
Assume thatThis means that edge weights between the nearer and chemically similar nodes are larger.
Edge weights between node 1 and node 2 are calculated:
Thus, preliminary stratigraphic network map Is constructed to include nodes and edges having the weights described above.
The formula:
Wherein:
: weighting of edges.
: Adjusting coefficients of path sensitivity, assuming
The calculation process comprises the following steps:
It is assumed that in addition to the connection between node 1 and node 2 described above, there is a connection between node 2 and node 3, the weight of which 0.005. The shortest path from node 1 to node 3 needs to be captured.
Direct path:
With node 2, the weights are calculated as follows:
thus, the weight of the shortest path between node 1 and node 3 is 0.000495.
The formula:
Wherein:
: threshold value of
The calculation process comprises the following steps:
for each node, calculate the sum of the inverse of the path weights of all other nodes to this node and check if it is greater than a threshold
For node 1:
Because of Node 1 does not belong to any strongly connected component.
By similar calculations it can be determined whether other nodes belong to a strongly connected component.
The formula:
Wherein:
: network stability index, assuming that it was derived from previous analysis
The calculation process comprises the following steps:
combining information of shortest path and strong connectivity components among all nodes and network stability index And constructing a final stratum continuity network diagram. This figure can show the connection strength between nodes and the overall stability of the network.
Referring to fig. 4, the steps for acquiring the dimension-reduced geological feature set specifically include:
by utilizing the structural characteristics in the stratum continuity network diagram, a local linear embedding method is applied, and the formula is adopted:
Calculating local neighborhood embedding to generate a local embedding mapping result;
wherein, Representing the weight of the object to be weighed,Representing data points;
in combination with the local embedding mapping result, using an equidistant mapping method, the following formula is used:
Optimizing the distance mapping of the global geological data to generate an equidistant mapping result;
wherein, Representing a core distance adjustment factor;
integrating the local embedded mapping result and the equidistant mapping result, and using the formula:
Optimizing the dimension reduction process to generate a dimension reduction result;
wherein, Represents an adjustment coefficient;
extracting key geological features from the dimension reduction result, and using the formula:
generating a dimension-reducing geological feature set;
wherein, Representing a geologic feature.
The formula:
detailed description of parameters:
: weight coefficient representing data point AndIs the relative importance or the strength of the connection.
Data points in a high-dimensional space.
Computational flow and example:
assume three data points And weight coefficient. Calculation ofThe following are provided:
Substituting the value into the formula:
Results 7.6 represent the locally linear embedded total distance based on the given weights, which can be used for further data analysis and feature extraction.
The formula:
detailed description of parameters:
: a core distance adjustment factor for adjusting the point AndThe distance between them.
: And E, calculating Euclidean distance.
Computational flow and example:
using the same data points And (2) and. Calculation ofThe following are provided:
Substituting the value into the formula:
result 9 represents the global distance sum based on the given distance adjustment factor for globally optimizing the data mapping.
The formula:
Detailed description of parameters
AndIs an adjustment factor for balancing the effects of local embedding and global mapping in the dimension reduction result.
AndThe results of the above-mentioned calculated local embedding and global mapping, respectively.
Computing flow and examples
The assumption is that the results are respectivelyAndIf set upAndThen:
Substituting the value into the formula:
Results 8.44 represent a dimension-reduced dataset that combines local and global information for further analysis and feature extraction.
The formula:
Detailed description of parameters
Representing key geologic features extracted from the comprehensive reduced-dimension dataset.
Computing flow and examples
Hypothesis resultsIs used to identify key features in the dataset. The model is simplified here, assuming three key features are identified. The specific extraction method depends on the feature extraction technique employed, such as Principal Component Analysis (PCA).
It is assumed that these features represent the main direction of change or cluster center, respectively, in the dataset. For example, ifRepresenting a geological dataset after the dimension reduction process,And respectively represent key indexes such as stratum density, mineral content, geologic age and the like.
Referring to fig. 5, the steps for obtaining the result of geological event prediction analysis specifically include:
based on the dimensionality reduction geological feature set, analyzing a key geological structure, and adopting the formula:
Calculating the significance of the key geological structure, and generating a key geological structure analysis result;
wherein, Representing the result of the analysis of the key geological structure,The correlation coefficient is represented by a correlation coefficient,The characteristic item is represented by a characteristic term,Representing the number of features;
according to the analysis result of the key geological structure, the dynamic change mode is analyzed, and the formula is adopted:
quantifying the geologic structure change trend to generate a dynamic change mode analysis result;
wherein, Representing the result of the analysis of the dynamic change pattern,The rate of change is indicated as being indicative of,Representing the result of the analysis of the key geological structure,Representing the number of features;
using the dynamic change mode analysis result, adopting the formula:
Determining the most obvious change trend, and generating a geological structure change trend identification result;
wherein, Representing the result of identifying the change trend of the geological structure,Representing trend intensity;
Based on the geological structure change trend identification result, the development trend prediction of geological events is carried out, and the formula is adopted:
integrating all the variation trends to generate a geological event prediction analysis result;
wherein, Representing the result of the predictive analysis of the geological event,The prediction parameters are represented by a set of parameters,Representing the result of identifying the change trend of the geological structure,Representing the number of features.
The formula:
Parameter interpretation and derivation process:
: analysis of key geologic structures results, quantifying the significance of each structure.
: The correlation coefficients of the key geologic structures represent the importance of each feature in the geologic analysis.
: And extracting the geological features from the dimension-reducing geological feature set.
: Total number of features.
Calculating:
Assume three features Correlation coefficientCharacteristic value
The calculation process comprises the following steps:
Analysis of results:
Obtained by Representing a significance score for the key geologic structure calculated based on the weight correlation coefficients and the eigenvalues.
The formula:
Parameter interpretation and derivation process:
: and the analysis result of the dynamic change mode represents the trend quantification of the change of the geological structure.
: And (5) quantifying the change intensity of each geological structure along with time according to the change rate obtained by time sequence analysis.
: Significance scores for key geologic structures are known.
Calculating:
Using the preceding steps Value and assumed rate of change
The calculation process comprises the following steps:
Analysis of results:
representing the trend score of the key geologic structure over time.
The formula:
Parameter interpretation and derivation process:
: and the geologic structure change trend identification result shows the most obvious one of all change trends.
Representing the intensity of a single trend.
Calculating:
Assume three trend strengths
The calculation process comprises the following steps:
Analysis of results:
the most pronounced trend intensity is 5 out of the three trends.
The formula:
Parameter interpretation and derivation process:
: and (5) predicting analysis results of the geological event, and synthesizing the estimated development of all the trend predicted geological events.
: Prediction parameters, representing the prediction weights of the multiple trends.
: The most pronounced trend intensity.
Calculating:
assuming prediction parameters And before
The calculation process comprises the following steps:
Analysis of results:
total predictive score representing the geological time resulting from the integrated multi-trend.
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.

Claims (5)

1.一种综合地质勘探信息管理系统,其特征在于,所述系统包括:1. A comprehensive geological exploration information management system, characterized in that the system comprises: 非线性动态数据建模模块将地层振动与压力变化的时间序列数据转换为多维相空间,通过延迟坐标映射地质内部动态结构,分析重构的多维相空间中的关键转折点,生成地质动态吸引子图;The nonlinear dynamic data modeling module converts the time series data of formation vibration and pressure change into multi-dimensional phase space, maps the internal dynamic structure of the geology through delayed coordinates, analyzes the key turning points in the reconstructed multi-dimensional phase space, and generates a geological dynamic attractor diagram; 地层连续性动态评估模块利用所述地质动态吸引子图,结合岩石的物理和化学特性,构建并实时更新无向图或有向图的节点与边,监测岩层连续性和潜在断层,通过实时计算节点间的最短路径和识别强连通分量,评估网络的稳定性指标,生成地层连续性网络图;The dynamic assessment module of stratum continuity utilizes the geological dynamic attractor graph and combines the physical and chemical properties of rocks to construct and update the nodes and edges of undirected or directed graphs in real time, monitors stratum continuity and potential faults, and evaluates the stability index of the network by calculating the shortest path between nodes in real time and identifying strongly connected components, thereby generating a stratum continuity network graph; 高维数据降维分析模块利用所述地层连续性网络图中的结构特征,应用局部线性嵌入与等距映射对高维地质数据进行降维,提取关键地质特征,建立降维地质特征集;The high-dimensional data dimension reduction analysis module uses the structural features in the stratigraphic continuity network diagram to reduce the dimension of the high-dimensional geological data by applying local linear embedding and isometric mapping, extracts key geological features, and establishes a dimension reduction geological feature set; 地质变化模式预测模块基于所述降维地质特征集,分析关键地质结构和动态变化模式,识别地质结构变化趋势,进行地质事件的发展趋势预测,生成地质事件预测分析结果。The geological change pattern prediction module analyzes key geological structures and dynamic change patterns based on the reduced-dimensional geological feature set, identifies geological structure change trends, predicts the development trends of geological events, and generates geological event prediction analysis results. 2.根据权利要求1所述的综合地质勘探信息管理系统,其特征在于,所述地质动态吸引子图的获取步骤具体为:2. The comprehensive geological exploration information management system according to claim 1, characterized in that the step of obtaining the geological dynamic attractor graph is specifically: 将地层振动与压力变化的时间序列数据转换为多维数据点集,使用公式:Convert the time series data of formation vibration and pressure changes into a multidimensional data point set using the formula: 将时间序列数据映射为高维空间,生成多维数据点集;Map time series data into a high-dimensional space to generate a multidimensional data point set; 其中,代表时刻的数据,为时间间隔,为维度数;in, Representative moments data, is the time interval, is the number of dimensions; 利用所述多维数据点集,采用公式:Using the multidimensional data point set, the formula is adopted: 计算数据点间的欧氏距离,生成距离矩阵;Calculate the Euclidean distance between data points and generate a distance matrix; 其中,表示数据点之间的距离,为维度数,为时间间隔;in, Represents data points and The distance between is the number of dimensions, is the time interval; 从所述距离矩阵中筛选关键转折点,采用公式:The key turning points are screened from the distance matrix using the formula: 获取关键转折点集合;Get the key turning point set; 其中,表示被选中的数据点集合,表示点与点间的距离,为每对点之间的权重,根据数据点相对重要性或连接强度设定;in, represents the set of selected data points, Indicate point With point The distance between For each pair of points and The weights between them are set according to the relative importance of data points or the strength of connection; 基于所述关键转折点集合,使用公式:Based on the set of key turning points, the formula is used: 构建地质动态吸引子图;Constructing geological dynamic attractor diagrams; 其中,表示地质动态吸引子图的定量描述,为数据点,是缩放参数,用于调整距离的影响力,使吸引力更敏感或更抑制。in, Represents a quantitative description of the geological dynamic attractor diagram, is the data point, is a scaling parameter that adjusts the influence of distance, making attraction more sensitive or more suppressed. 3.根据权利要求2所述的综合地质勘探信息管理系统,其特征在于,所述地层连续性网络图的获取步骤具体为:3. The comprehensive geological exploration information management system according to claim 2, characterized in that the step of obtaining the stratigraphic continuity network diagram is specifically: 根据所述地质动态吸引子图和岩石的物理与化学特性,将吸引子图和特性数据融合,采用公式:According to the geological dynamic attractor diagram and the physical and chemical properties of rocks, the attractor diagram and characteristic data are fused and the formula is adopted: 生成初步的地层网络图;Generate preliminary stratigraphic network maps; 其中,代表初步地层网络图,代表节点,代表边,代表地质动态吸引子图,代表物理特性数据,代表化学特性数据;in, Represents a preliminary stratigraphic network diagram, Represents a node, Represents the side, Represents the geological dynamic attractor diagram, Represents physical property data, Represents chemical property data; 利用所述初步的地层网络图,计算节点间的最短路径,采用公式:Using the preliminary stratigraphic network diagram, the shortest path between nodes is calculated using the formula: 计算节点间最短路径,生成节点间最短路径分析结果;Calculate the shortest path between nodes and generate the shortest path analysis results between nodes; 其中,代表节点之间的最短路径,代表边的权重,代表调节参数;in, Representative Node and The shortest path between represents the weight of the edge, represents the adjustment parameter; 利用所述节点间最短路径分析结果,识别网络中的强连通分量,采用公式:Using the shortest path analysis results between nodes, the strongly connected components in the network are identified using the formula: 其中为阈值,生成强连通分量分析结果;in is the threshold value, generating the strongly connected component analysis result; 其中,代表强连通分量,代表节点,代表差异化节点,代表阈值;in, represents the strongly connected component, Represents a node, and represents a differentiation node, represents the threshold value; 整合所述节点间最短路径分析结果和强连通分量分析结果,采用公式:Integrate the shortest path analysis results between nodes and the strongly connected component analysis results, using the formula: 生成地层连续性网络图;Generate stratigraphic continuity network diagrams; 其中,代表地层连续性网络图,代表强连通分量,代表节点间最短路径,代表网络稳定性指标。in, Represents the stratigraphic continuity network diagram, represents the strongly connected component, represents the shortest path between nodes, Represents the network stability index. 4.根据权利要求3所述的综合地质勘探信息管理系统,其特征在于,所述降维地质特征集的获取步骤具体为:4. The comprehensive geological exploration information management system according to claim 3, characterized in that the step of obtaining the reduced-dimensional geological feature set is specifically: 利用所述地层连续性网络图中的结构特征,应用局部线性嵌入方法,采用公式:Using the structural characteristics of the stratigraphic continuity network diagram, a local linear embedding method is applied, using the formula: 计算局部邻域嵌入,生成局部嵌入映射结果;Calculate local neighborhood embedding and generate local embedding mapping results; 其中,代表权重,代表数据点;in, represents the weight, represents a data point; 结合所述局部嵌入映射结果,使用等距映射法,通过公式:Combined with the local embedding mapping results, the isometric mapping method is used, through the formula: 优化全局地质数据的距离映射,生成等距映射结果;Optimize the distance mapping of global geological data and generate isometric mapping results; 其中,代表核心距离调整因子;in, represents the core distance adjustment factor; 整合所述局部嵌入映射结果和等距映射结果,使用公式:Integrate the local embedding mapping results and the isometric mapping results, using the formula: 优化降维过程,生成降维结果;Optimize the dimensionality reduction process and generate dimensionality reduction results; 其中,代表调节系数;in, represents the adjustment coefficient; 从所述降维结果中提取关键地质特征,使用公式:From the dimensionality reduction results, key geological features are extracted using the formula: 生成降维地质特征集;Generate a reduced dimension geological feature set; 其中,代表地质特征。in, Represents geological features. 5.根据权利要求4所述的综合地质勘探信息管理系统,其特征在于,所述地质事件预测分析结果的获取步骤具体为:5. The comprehensive geological exploration information management system according to claim 4, characterized in that the step of obtaining the geological event prediction and analysis results is specifically: 基于所述降维地质特征集,分析关键地质结构,采用公式:Based on the reduced dimension geological feature set, key geological structures are analyzed using the formula: 计算关键地质结构的显著性,生成关键地质结构分析结果;Calculate the significance of key geological structures and generate key geological structure analysis results; 其中,代表关键地质结构分析结果,表示相关系数,表示特征项,代表特征数量;in, Represents key geological structure analysis results, represents the correlation coefficient, represents the feature item, represents the number of features; 根据所述关键地质结构分析结果,分析动态变化模式,采用公式:According to the analysis results of the key geological structures, the dynamic change mode is analyzed using the formula: 对地质结构变化趋势进行量化,生成动态变化模式分析结果;Quantify the trend of geological structure changes and generate dynamic change pattern analysis results; 其中,代表动态变化模式分析结果,表示变化率,代表关键地质结构分析结果,代表特征数量;in, Represents the results of the dynamic change pattern analysis, represents the rate of change, Represents key geological structure analysis results, represents the number of features; 利用所述动态变化模式分析结果,采用公式:The dynamic change pattern analysis results are used, using the formula: 确定最显著的变化趋势,生成地质结构变化趋势识别结果;Determine the most significant change trend and generate geological structure change trend identification results; 其中,代表地质结构变化趋势识别结果,表示趋势强度;in, Represents the identification result of geological structure change trend, Indicates the strength of the trend; 基于所述地质结构变化趋势识别结果,进行地质事件的发展趋势预测,采用公式:Based on the identification results of the geological structure change trend, the development trend of geological events is predicted using the formula: 整合所有变化趋势,生成地质事件预测分析结果;Integrate all change trends to generate geological event prediction and analysis results; 其中,代表地质事件预测分析结果,表示预测参数,代表地质结构变化趋势识别结果,代表特征数量。in, Represents the results of geological event prediction and analysis. represents the prediction parameter, Represents the identification result of geological structure change trend, Represents the number of features.
CN202411050620.6A 2024-08-01 2024-08-01 Comprehensive geological exploration information management system Active CN118568446B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411050620.6A CN118568446B (en) 2024-08-01 2024-08-01 Comprehensive geological exploration information management system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411050620.6A CN118568446B (en) 2024-08-01 2024-08-01 Comprehensive geological exploration information management system

Publications (2)

Publication Number Publication Date
CN118568446A CN118568446A (en) 2024-08-30
CN118568446B true CN118568446B (en) 2024-11-19

Family

ID=92467336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411050620.6A Active CN118568446B (en) 2024-08-01 2024-08-01 Comprehensive geological exploration information management system

Country Status (1)

Country Link
CN (1) CN118568446B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104699979A (en) * 2015-03-24 2015-06-10 北京工商大学 Urban lake and reservoir algal bloom chaos time sequence predication method based on complicated network
CN117724348A (en) * 2024-02-06 2024-03-19 深圳市万斯得自动化设备有限公司 Accurate pressure regulation and control system based on explosion testing machine

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050165566A1 (en) * 2002-06-03 2005-07-28 Omnigon Technologies Ltd. Method for solving waveform sequence-matching problems using multidimensional attractor tokens
CN117725537A (en) * 2024-02-01 2024-03-19 深圳凯升联合科技有限公司 Real-time metering data processing platform
CN117891857B (en) * 2024-03-13 2024-05-24 广东工业大学 Data mining method and system based on big data
CN118394873B (en) * 2024-06-28 2024-09-06 山东省国土测绘院 AI-based geographical information crowdsourcing data analysis method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104699979A (en) * 2015-03-24 2015-06-10 北京工商大学 Urban lake and reservoir algal bloom chaos time sequence predication method based on complicated network
CN117724348A (en) * 2024-02-06 2024-03-19 深圳市万斯得自动化设备有限公司 Accurate pressure regulation and control system based on explosion testing machine

Also Published As

Publication number Publication date
CN118568446A (en) 2024-08-30

Similar Documents

Publication Publication Date Title
US10921470B2 (en) Inversion for tectonic stress
CN112689778A (en) Inverse formation modeling using mixed linear and nonlinear algorithms
CA2641867A1 (en) Methods and apparatus for predicting the hydrocarbon production of a well location
CN106920176B (en) Mining area scale mineral resource estimation method and system
CA3177867A1 (en) Structured representations of subsurface features for hydrocarbon system and geological reasoning
CN116911216B (en) Reservoir oil well productivity factor assessment and prediction method
US11391856B2 (en) Stochastic dynamic time warping for automated stratigraphic correlation
Chang et al. SegLog: Geophysical logging segmentation network for lithofacies identification
US11221425B1 (en) Generating a model for seismic velocities in a subsurface region using inversion with lateral variations
CN118759600A (en) A prospecting prediction method suitable for tectonic altered rock type gold deposits
CN113052968B (en) A Knowledge Graph Construction Method for 3D Structural Geological Models
US20210404329A1 (en) Method for identifying geological and drilling patterns in a unidimensional space via periodic orthogonal functions applied to drilling parameter data
CN118568446B (en) Comprehensive geological exploration information management system
CN119443418A (en) A multi-dimensional fusion prospecting prediction method based on multi-source data
US11899150B2 (en) Velocity model for sediment-basement interface using seismic and potential fields data
Chen et al. A geospatial case‐based reasoning model for oil–gas reservoir evaluation
US12468058B2 (en) Thermal conductivity mapping from rock physics guided seismic inversion
CN117034657A (en) Geophysical data processing module based on priori information constraint
CN112213772B (en) A method, device and system for identifying underground sedimentary body morphology by combining well-seismic analysis
Sun et al. Dynamic graph convolutional networks for fluid identification of well logging data transformed through the gram angle field
Li et al. Massive Spatial Well Clustering Based on Conventional Well Log Feature Extraction for Fast Formation Heterogeneity Characterization
US20250021729A1 (en) Porosity-assisted facies modeling for hydrocarbon extraction
RU2563862C2 (en) Strain attribute in rocks
US20240210591A1 (en) Characterization of Subsurface Geological Formations
Bhowmik et al. Predicting the migration of CO2 plume using injection data and a distance-metric approach to reservoir-model selection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant