Ore formation simulation method and device for ore formation research area

文档序号：136335 发布日期：2021-10-22 浏览：11次中文

阅读说明：本技术 成矿研究区的成矿模拟方法及装置 (Ore formation simulation method and device for ore formation research area ) 是由王功文黄蕾蕾李瑞喜于 2021-06-07 设计创作，主要内容包括：本发明提供一种成矿研究区的成矿模拟方法及装置,所述方法包括：基于地质的成矿模式与平面、剖面地质资料建立三维地质模型；通过测试研究区的成矿期包裹体样品,获得成矿期的温度、压力等成矿临界状态；通过测试围岩与赋矿岩体的岩石力学特性,获得数值模拟的基本参数；根据成矿模式与成矿区研究,对三维地质模型赋予初始状态条件；调整参数,借助临界条件数据进行模拟,直到获得与临界条件相当的模拟模型；将临界模型进行三维后处理,获得三维成矿范围与成矿潜力区。本发明能准确成矿研究区中的矿藏分布预测。(The invention provides a method and a device for simulating mineralization of an mineralization research area, wherein the method comprises the following steps: establishing a three-dimensional geological model based on a geological mineralization mode and plane and section geological data; obtaining the mineralizing critical states of temperature, pressure and the like in the mineralizing period through testing the mineralizing period inclusion sample in the research area; obtaining basic parameters of numerical simulation by testing the rock mechanical properties of the surrounding rock and the mine-endowing rock mass; according to the mining mode and the mining area research, giving an initial state condition to the three-dimensional geological model; adjusting parameters, and simulating by means of critical condition data until a simulation model equivalent to the critical condition is obtained; and (4) carrying out three-dimensional post-processing on the critical model to obtain a three-dimensional mineralization range and a mineralization potential area. The method can accurately predict the mineral deposit distribution in the mineralization research area.)

1. A method of simulating mineralization in an mineralization study area, the method comprising:

establishing a three-dimensional geological model based on a geological mineralization mode and plane and section geological data;

obtaining at least numerical simulation parameters in the geological model by a sampling experiment testing mode, wherein the numerical simulation parameters at least comprise at least one of the following parameters: rock mechanics parameters, tectonic dynamics parameters, hydromechanics parameters, and thermodynamic parameters;

importing the three-dimensional rock mass-structure-channel geological model into a geological numerical simulation parameter assignment application, selecting and determining a numerical simulation constitutive model, and at least endowing rock mechanical parameters, structure kinetic parameters, fluid mechanical parameters and thermodynamic parameters to different medium materials in the three-dimensional rock mass-structure-channel geological model;

according to the statistics of the field geological parameters of the mineralization research area and the test inversion of the samples, obtaining the initial geological parameters in the mineralization process and the mineralization geological parameters of each rock mass medium in the final mineralization, and assigning values to the initial conditions of the mineralization research area based on the initial geological parameters; setting boundary conditions of the mineralization research area according to the mineralization geological parameters of each rock mass medium;

setting a numerical simulation step in a mode that the step length is set from small to large, carrying out dynamic mineralization simulation, generating a corresponding simulation dynamic graph, recording the current parameters when the actual mineralization characteristics are the same as certain iteration times, and readjusting the iteration times in the simulation to enable the final simulation numerical model to conform to the three-dimensional mineralization model;

after the numerical simulation is finished, deriving a three-dimensional structural stress value and three-dimensional fluid state data of the numerical simulation;

and selecting parameters meeting the mineralization requirement from the three-dimensional structural stress value and the three-dimensional fluid state data, performing Boolean operation to obtain a target area, and taking the target area as an mineralization area in an mineralization research area.

2. The method of claim 1, wherein the geological based mineralization pattern is modeled with a planar, profile geological data in three dimensions, including:

based on at least one of a geological map, a roadway engineering plan, a drilling record map, an exploration profile map and a geophysical interpretation profile map, at least carrying out vectorization, standardization and element information extraction processing, summarizing horizontal and vertical spatial information of geological information of strata, lithology and geophysical properties, carrying out data integration in three-dimensional modeling application, and obtaining the spatial position and contact relation of a geologic body;

performing near infrared spectrum analysis on rocks of a drilling sample, a tunnel sample and an earth surface sample in the geology to obtain the alteration type of the geology; when three-dimensional alteration mapping is carried out, the spatial change information of alteration lithology is obtained based on the alteration type, and the position of a main channel of hydrothermal fluid is indicated to be used as the spatial position of a hydrothermal alteration channel;

carrying out occurrence statistics on the geological fracture layer to obtain the direction and scale of the structure motion before and after the mineralization period; through geochemical test of the rocks of two fault disks, the times of the tectonic movement are obtained; obtaining the structural stress of an ore forming period by carrying out crystal form analysis on minerals in the fault, and constructing a fracture space position;

acquiring drilling core sample bit data, gallery exploration sample bit data and surface chemical exploration data of a geological ore carrying rock mass, an ore forming rock mass and surrounding rocks based on X-ray fluorescence spectrum analysis and a main trace element method, and determining the space position of the ore forming body;

and under the comprehensive geological constraints of the plate structure characteristics, ore collection area ore control factors and the ore formation mode of the ore region, the spatial position and the contact relation of the geologic body, the spatial position of the hydrothermal alteration passage, the spatial position of the fracture and the spatial position of the ore body are integrated, and a three-dimensional rock mass-structure-passage geological model is established under the constraint of three-dimensional modeling application.

3. The method of claim 1, wherein the rock-mechanics parameter comprises at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

4. The method of claim 1, wherein the formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass;

correspondingly, the obtaining of numerical simulation parameters in the geological model comprises:

obtaining the structure movement direction and scale distance before and after the mineralization period through the statistics of fracture fault occurrence of the mineralization area; through geochemical test of two disks of rocks of the fault, the times and the time of the movement of the structure of the mineralization research area are obtained; and obtaining the tectonic stress of the mineralization research area in the mineralization period by analyzing the mineral crystal form in the fault;

and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

5. The method of claim 1, further comprising:

testing the rock inclusion to obtain the fluid pressure and the fluid temperature of the fluid in the mineralization period, and carrying out a temperature penetration experiment on the surrounding rock sample to obtain the permeability and the porosity of the rock and the heat conduction data of the rock.

6. An mineralization simulation apparatus for an mineralization study area, the apparatus comprising:

the building unit is used for building a three-dimensional geological model based on geological mineralization patterns, plane and section geological data;

an obtaining unit, configured to obtain at least a numerical simulation parameter in the geological model by a sampling experiment testing manner, where the numerical simulation parameter at least includes at least one of: rock mechanics parameters, tectonic dynamics parameters, hydromechanics parameters, and thermodynamic parameters;

the evaluation unit is used for importing the three-dimensional rock mass-structure-channel geological model into a geological numerical simulation parameter evaluation application, selecting and determining a numerical simulation constitutive model, and endowing different medium materials in the three-dimensional rock mass-structure-channel geological model with at least rock mechanical parameters, structure kinetic parameters, fluid mechanical parameters and thermodynamic parameters;

the setting unit is used for obtaining initial geological parameters at the beginning in the mineralization process and mineralization geological parameters of each rock mass medium in the final mineralization process according to the statistics of the field geological parameters of the mineralization research area and the sample test inversion, and assigning values to the initial conditions of the mineralization research area based on the initial geological parameters; setting boundary conditions of the mineralization research area according to the mineralization geological parameters of each rock mass medium;

the simulation unit is used for setting a numerical simulation step in a mode that the step length is set from small to large, carrying out dynamic mineralization simulation, generating a corresponding simulation dynamic graph, recording the current parameters when the actual mineralization characteristics are the same as certain iterative times, and readjusting the iterative times in the simulation to enable the final simulation numerical model to conform to the three-dimensional mineralization model;

the derivation unit is used for deriving the three-dimensional structural stress value and the three-dimensional fluid state data of the numerical simulation after the numerical simulation is finished;

and the mineralization determining unit is used for selecting parameters meeting mineralization requirements from the three-dimensional structural stress values and the three-dimensional fluid state data, performing Boolean operation to obtain a target area, and taking the target area as an mineralization area in an mineralization research area.

7. The apparatus of claim 6, wherein the establishing unit is further configured to:

acquiring drilling core sample bit data, gallery exploration sample bit data and surface chemical exploration data of a geological ore carrying rock mass, an ore forming rock mass and surrounding rocks based on X-ray fluorescence spectrum analysis and a main trace element device so as to determine the space position of the ore forming body;

8. The apparatus of claim 6, wherein the rock mechanics parameter comprises at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

9. The apparatus of claim 6, wherein the formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass;

correspondingly, the obtaining unit is further configured to:

and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

10. The apparatus of claim 6, wherein the obtaining unit is further configured to:

Technical Field

The invention relates to a technology for determining a mineralization area, in particular to a method and a device for simulating mineralization of an mineralization research area.

Background

The mineralization and mineralization are complex dynamic processes, the basic problem of the mineral deposit formation factor, and the conclusion is the dynamic problem of the mineralization. The mining kinetics mainly researches the speed, mechanism and process of mining action, and is the core problem of mineral deposit cause, the traditional mineral deposit science mainly summarizes and summarizes the characteristics, rules and experiences of mineral deposits, explains the rules and facts, and further deduces, and the mining kinetics summarizes and explains the mining characteristics and rules in a wider area. Both are based on the fact that the most essential characteristics are that the results of ore formation are taken as the main objects of research, the main purpose is to find the formation rule of ore deposit, the research idea is to deduce the reason from the results, and the deduction is very not strict with guess and hypothetical components to a large extent. The mining kinetics theory is mainly aimed at researching the mining process, and aims to find the reason and mechanism of mineral deposit formation, and the research method mainly comprises the reason to the result.

In recent years, along with the development of computer technology, computers provide a plurality of auxiliary technical means for the mineralization dynamics, and the computer technology is used for carrying out the numerical simulation of the mineralization dynamics. The method comprehensively utilizes research results such as geology, geophysical, geochemistry, remote sensing and the like to establish a geological model which is not limited by time and space, and becomes one of important methods for modern earth science research. Since the numerical simulation technology is widely applied to the geological field in the nineties of the 20 th century, the numerical simulation of a geological system becomes a hotspot and development trend in the current geological field and is also a necessary way for the precision and quantification of geological science. The numerical simulation of the mineralization dynamics can simulate and reproduce the geological process of the mineralization and becomes one of the important methods for the modern geoscience research. It combines geological data and experimental results on the basis of detailed geological work, quantifies geological phenomena, and establishes an ore-forming geological process model through induction and deduction logic speculation methods to help people to know the geological structure evolution process. Through numerical simulation experimental study, the coupling relation among fluid, stress and temperature in the geological history period can be deeply understood; in geological evolution, the evolution process of rock ring structures, stress strain fields and temperature fields in each historical period, the transmission and conversion of energy and the accumulation process of mineral liquid of ore bodies in the forming process of ore, the forming process of ore deposits and the forming mechanism thereof and the like.

Disclosure of Invention

In view of the above, an embodiment of the present invention provides a method and an apparatus for simulating mineralization in an mineralization research area.

The invention provides a method for simulating mineralization of an mineralization research area, which comprises the following steps:

establishing a three-dimensional geological model based on a geological mineralization mode and plane and section geological data;

after the numerical simulation is finished, deriving a three-dimensional structural stress value and three-dimensional fluid state data of the numerical simulation;

Preferably, the building of the three-dimensional geological model based on the geological mineralization pattern and the plane and section geological data comprises:

Preferably, the rock-mechanical parameters include at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

Preferably, the formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass;

correspondingly, the obtaining of numerical simulation parameters in the geological model comprises:

and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

Preferably, the method further comprises:

Another aspect of the present invention provides a mineralization simulation apparatus for an mineralization study area, comprising:

the building unit is used for building a three-dimensional geological model based on geological mineralization patterns, plane and section geological data;

Preferably, the establishing unit is further configured to:

acquiring drilling core sample bit data, gallery exploration sample bit data and surface chemical exploration data of a geological ore carrying rock mass, an ore forming rock mass and surrounding rocks based on X-ray fluorescence spectrum analysis and a main trace element device so as to determine the space position of the ore forming body;

Preferably, the rock-mechanical parameters include at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

Preferably, the formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass;

correspondingly, the obtaining unit is further configured to:

and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

Preferably, the obtaining unit is further configured to:

Another aspect of the present invention provides an electronic device, including: a processor and a memory for storing processor executable instructions, wherein the processor is configured to be able to perform the steps of the mineralization simulation method of the mineralization study area upon invocation of the executable instructions in the memory.

Another aspect of the present invention provides a computer readable storage medium having computer instructions stored thereon, wherein the instructions, when executed by a processor, implement the steps of the method for simulating mineralization in a mineralization study area.

According to the method and the device for simulating the mineralization of the mineralization research area, disclosed by the embodiment of the invention, a dynamic structure fluid channel model of an mineralization key period is obtained by performing fine three-dimensional modeling on strata, rock masses, structural fractures and the like of the mineralization research area, the data states of temperature, fluid pressure, dynamic motion direction, stress magnitude and the like of the mineralization fluid are obtained by means of a traditional mineral deposit research method, the mineralization dynamics numerical simulation is performed on three aspects of the mineralization area, a deep section and a mineral deposit under the guidance of a classical mineralization theory, the spatial distribution rule and the mineral control characteristic analysis of the mineral body and various geological elements are revealed, the mineralization mechanism and the evolution process are known, and the mineralization theory is further verified and refined. The invention standardizes the complex mineralization simulation mode in a flow process, converts a plurality of complex processes into a flow modeling-simulation process, can fill in simulation adjustment parameters according to requirements after preparing corresponding data in the early stage, and can obtain a simulation result through iterative adjustment.

Drawings

Fig. 1 is a schematic flow chart of a method for simulating mineralization in an mineralization study area according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of a three-dimensional modeling process of an ore-forming related element according to an embodiment of the invention;

FIG. 3 is a schematic diagram illustrating the collection of numerical simulation parameters in a geological model according to an embodiment of the present invention;

fig. 4 is a schematic structural diagram of a composition of a mineralization simulation apparatus in an mineralization study area according to an embodiment of the present invention.

Detailed Description

The essence of the technical scheme of the invention is explained in detail in the following with the accompanying drawings.

Fig. 1 is a schematic flow chart of a method for simulating mineralization in an mineralization research area according to the present invention, and as shown in fig. 1, the method for simulating mineralization in an mineralization research area according to an embodiment of the present invention includes the following processing steps:

step 101, establishing a three-dimensional geological model based on geological mineralization patterns and plane and section geological data.

The establishment of the three-dimensional geological model in the embodiment of the invention specifically comprises the following steps:

And 102, at least obtaining numerical simulation parameters in the geological model in a sampling experiment testing mode.

Wherein the numerical simulation parameters include at least one of: rock mechanics parameters, tectonic dynamics parameters, hydromechanics parameters, and thermodynamic parameters. The rock mechanics parameters include at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

The formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass; correspondingly, the obtaining of numerical simulation parameters in the geological model comprises: obtaining the structure movement direction and scale distance before and after the mineralization period through the statistics of fracture fault occurrence of the mineralization area; through geochemical test of two disks of rocks of the fault, the times and the time of the movement of the structure of the mineralization research area are obtained; and obtaining the tectonic stress of the mineralization research area in the mineralization period by analyzing the mineral crystal form in the fault; and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

The method for acquiring the numerical simulation parameters further comprises the following steps: testing the rock inclusion to obtain the fluid pressure and the fluid temperature of the fluid in the mineralization period, and carrying out a temperature penetration experiment on the surrounding rock sample to obtain the permeability and the porosity of the rock and the heat conduction data of the rock.

103, obtaining initial geological parameters at the beginning of the mineralization process and mineralization geological parameters of each rock mass medium during final mineralization according to the statistics of the field geological parameters of the mineralization research area and sample test inversion, and assigning values to the initial conditions of the mineralization research area based on the initial geological parameters; and setting boundary conditions of the mineralization research area according to the mineralization geological parameters of each rock mass medium.

And 104, setting a numerical simulation step in a mode that the step length is set from small to large, performing dynamic mineralization simulation, generating a corresponding simulation dynamic graph, recording the current parameters when the actual mineralization characteristics are the same as certain iterative times, and readjusting the iterative times in the simulation to enable the final simulation numerical model to conform to the three-dimensional mineralization model.

And 105, after the numerical simulation is finished, deriving a three-dimensional structural stress value and three-dimensional fluid state data of the numerical simulation.

And 106, selecting parameters meeting the mineralization requirement from the three-dimensional structural stress value and the three-dimensional fluid state data, performing Boolean operation to obtain a target area, and taking the target area as an mineralization area in an mineralization research area.

The essence of the technical solution of the embodiment of the present invention is further clarified by specific examples below.

The method aims to obtain a dynamic structure fluid channel model of an ore-forming key period by performing fine three-dimensional modeling on strata, rock masses, structural fractures and the like in an ore-forming research area, obtain data states of temperature, fluid pressure, dynamic motion direction, stress magnitude and the like of an ore-forming fluid by means of traditional ore deposit research, perform ore-forming dynamic numerical simulation on three aspects of an ore-forming area, a deep section and an ore deposit under the guidance of a classical ore-forming theory, reveal the spatial distribution rule and ore control characteristic analysis of an ore body and various geological elements, understand an ore-forming mechanism and an evolution process, and further verify and refine an ore-forming theory. The post-processing is carried out on the result of the mineralization simulation, and an important reference basis can be provided for the prediction of the deep blind ore body.

The embodiment of the invention needs to determine the following preset conditions:

the ore deposit mineralization type is hydrothermal mineralization, has a definite mineralization action mode, and has definite dynamic structure movement or thermosetting coupling mineralization action formed by hot liquid fluid in the mineralization period; the stratum-rock mass space characteristics of the key period of mineralization can be established according to geological data, detailed data are provided for the space form of a fracture structure, and an accurate three-dimensional geological model is established through plane and section geological data, such as exploratory groove outcrop, geophysical data, well drilling and exploration line data, underground mining sections and the like; the main dynamic parameters of the mineralization stage such as the direction of tectonic movement, stress magnitude, the temperature, fluid pressure and flow rate of the mineralizing hydrothermal fluid and various dynamic physical properties of rocks of surrounding rocks, initial states, critical states and the like are obtained through geochemical research.

Specifically, the method for simulating the mineralization of the mineralization research area in the embodiment of the invention comprises the following processing steps:

establishing a three-dimensional geological model through an ore forming mode and plane and section geological data;

by testing an inclusion sample in an mineralization research area in an mineralization period, acquiring mineralization critical states such as temperature, pressure and the like in the mineralization period;

obtaining basic parameters of numerical simulation by testing the rock mechanical properties of the surrounding rock and the mine-endowing rock mass;

according to the mining mode and the mining area research, giving an initial state condition to the three-dimensional geological model;

adjusting parameters, and simulating by means of critical condition data until a simulation model equivalent to the critical condition is obtained;

the mineralization process can be obtained in a simulation; and (4) carrying out three-dimensional post-processing on the critical model to obtain a three-dimensional mineralization range and a mineralization potential area.

Fig. 2 is a schematic diagram of a three-dimensional modeling process of mining-related elements of the present invention, and as shown in fig. 2, the three-dimensional modeling process of the embodiment of the present invention is based on two mainstream modeling methods, namely, data-driven modeling and knowledge framework modeling, converts geological information data of a mining research area into a data model, and constructs a three-dimensional geological model under geological constraints of a deposit mining mode, and specifically includes the following steps:

first is the collection of three-dimensional geological framework information. For example, a two-dimensional geological map, a roadway engineering plan, a drilling record map, an exploration profile, a geophysical (heavy electromagnetic seismic) interpretation profile, the horizontal and vertical spatial information of geological information such as stratum, lithology and geophysical properties of an mineralization research area is summarized through vectorization, standardization and element information extraction of various software such as Autocad, Mapgis and Arcgis, and data integration is carried out in three-dimensional modeling software (such as GOCAD) to obtain the spatial position and contact relation of a geological body;

and secondly, determining the spatial positions of the ore carrying rock mass, the ore forming rock mass and the surrounding rock. The main method is to test the grade of the mineral elements by means of X-ray fluorescence spectrum analysis (XRF), main trace element analysis and other methods. The main objects of the test are the drilling core sample data, the gallery exploration sample data, the surface chemical exploration data and the like. A three-dimensional model of a known ore body can be established by spatial position information and numerical values of sampling points and by means of geostatistics and three-dimensional modeling software interpolation functions. If the mine has the three-dimensional ore body model, the mine can be used for reference;

hydrothermal channels are primarily indicative of channel location through different lithological alterations. Lithologic alterations can be traced back in underground tunnels through pulse visual interpretation, but the best effect is near infrared spectrum alteration interpretation. And (3) obtaining the alteration type by near infrared spectrum analysis of rocks of the drilling sample, the tunnel sample and the surface sample. Obtaining the spatial variation information of the altered lithology in the process of three-dimensional altered mapping, and further indicating the position of a main channel of the hydrothermal solution;

the structural movement is to obtain the structural movement direction and scale before and after the mineralization period through the statistics of fracture fault occurrence of the mineralization area; obtaining the times of tectonic movement by geochemical tests, such as years testing, of rocks of two disks of a fault; obtaining the tectonic stress of the mineral in the mineralization period by analyzing the mineral crystal form in the fault;

and finally, integrating the geological information under the comprehensive geological constraints of the plate structure characteristics, ore collection area ore control factors and ore formation modes of the ore region, and establishing a three-dimensional rock mass-structure-channel geological model under the constraint of three-dimensional modeling software (GOCAD).

Fig. 3 is a schematic diagram illustrating acquisition of numerical simulation parameters in a geological model according to an embodiment of the present invention, and as shown in fig. 3, the numerical simulation parameters according to the embodiment of the present invention mainly include rock mechanics parameters, structural dynamics parameters, fluid mechanics parameters, and thermodynamic parameters. The above parameters are tested by sampling experiment to obtain reference values.

The mechanical parameters of the rock comprise density, shear modulus, volume modulus, cohesive force, internal friction angle, expansion angle and tensile strength. The tests can respectively sample according to different lithologies in the three-dimensional model, and the tests are carried out for multiple times under a triaxial stress test instrument to obtain an average reference value.

The tectonic kinetic parameters comprise rock mass movement speed and tectonic stress. The tectonic movement of the rock mass in the mineralization period can obtain the tectonic movement direction and scale distance before and after the mineralization period through the statistics of fracture fault occurrence of the mineralization area; obtaining the times and time of the tectonic movements through geochemical tests, such as years testing, of rocks of two fault disks; obtaining the tectonic stress of the mineral in the mineralization period by analyzing the mineral crystal form in the fault; and obtaining the average motion speed of the robot by integrating the relative motion distance and the motion time.

In the embodiment of the invention, the statistical calculation of the structure movement speed-structure geological information is as follows:

the method comprises the steps of utilizing structural geological knowledge to conduct structural information survey statistics in a research area, counting structural deformation of a fractured zone in an ore deposit with zonal characteristics similar to hydrothermal alteration, and calculating a structural rose diagram after counting information of a main fractured surface and a secondary fractured surface of the fractured zone, and the reason, the (scratch) occurrence, the scale, the motion mode and the like, and carrying out inversion to obtain information of the structural motion direction, the motion speed and the like of the fracture.

The stress-quartz deformation and electron backscatter diffraction (EBSD) texture was constructed based on the foregoing information. Specifically, a quartz dynamic recrystallization sub-particle size method is adopted to calculate the corrosion change zone difference stress value of the Newcastle gold deposit, and experimental researches on metals and minerals show that the size of dynamic recrystallization particles is related to the stress when deformation reaches a steady state: the greater the stress, the finer the particles. And (3) completing statistics on the sizes of dynamic recrystallization particles of the samples of the Newcastle gold deposit in different structural deformation periods in AutoCAD by adopting a line cutting method. Using the relationship: sigma₁-σ₃＝AD^-mAnd (6) carrying out measurement and calculation. Wherein D is the dynamic recrystallization grain size (. mu.m), A is 6.1 and m is 0.68. Typically the size of the dynamically recrystallized grains is in the range of a few microns to a few tens of microns. The difference stress value change range before the new-city gold deposit is calculated by the relational expression is 85.63-106.85 MPa, and the difference stress change range during the ore forming period is 65.91-76.77 MPa.

And testing the rock inclusion to obtain the fluid pressure and the fluid temperature of the fluid in the mineralization period, and performing a temperature permeability experiment on the surrounding rock sample to obtain the permeability and the porosity of the rock and the heat conduction data of the rock.

The quartz inclusion in the ore-forming rock body is utilized, and the solid-liquid gas-phase object form change is observed through the cold and hot treatment of the thin sheet, so that the temperature during the formation is inverted.

In the embodiment of the present invention, the inclusion is set as follows: and (4) homogenizing the system. The trapped fluid is a homogeneous system when the inclusion is formed, i.e. the primary mineral is grown in a homogeneous system. And (4) closing the system. The fluid filled (retained) in the crystal defect is a main mineral closed system, and a separate closed system is formed, and no foreign substance is added and no internal substance is escaped. An isochoric system. After the inclusion is formed, the volume is basically constant, and the characteristics of an isochoric system are kept, so that various related physical and chemical phase diagrams can be utilized.

The quartz inclusion may be determined by:

1. a freezing method: refers to the process of observing the liquid-to-solid phase transition (i.e., solidification) as the inclusion cools below room temperature. The basic principle is to observe the phase change process of the inclusion by changing the temperature on a cold stage.

Obey raoult's law-for dilute concentration solutions, the freezing point depression value of the solution is independent of the type and nature of the solute, and depends only on the concentration of the solute dissolved in water (solvent); the freezing point depression temperature is the same for each solute having the same concentration.

2. The homogeneous method comprises the following steps: according to the basic assumptions and premises of the inclusion, the fluid captured by the inclusion is an originally homogeneous single phase fluid that fills the entire inclusion space. Along with the temperature reduction, the contraction coefficient of the fluid (gas or liquid) is larger than that of the solid (main mineral), the inclusion body evolves along an isochoric line until the position of a two-phase interface, if the fluid which is larger than the critical density is captured originally, a gas phase is separated, and after the gas escapes, the gas forms spherical bubbles at favorable positions due to the influence of surface tension; if a gas-rich fluid less than the critical density is originally captured, the gas condenses out of a liquid phase in the fluid, forming a two-phase inclusion with a large bubble.

If the temperature is raised in the hot and cold stages, a reversible phase change phenomenon is visible: it is first observed that the ratio of gas and liquid phases changes with increasing temperature, and when the temperature is raised to a certain temperature, phase transition occurs, i.e. two (or more) phases are changed into one phase, i.e. phase homogeneity is achieved, and the temperature is the homogeneous temperature.

In the embodiment of the invention, the principle of thermometric analysis is realized in such a way that the thermometric analysis can be widely applied to various transparent (or semitransparent) minerals with geological significance by only adding a thermometric device to an optical microscope. The method is to measure the temperature of various instantaneous phase changes by heating or freezing on the basis of observing and identifying various phases (solid phase, gas phase and liquid phase) of the fluid contained in the inclusion in detail. Compared with a linkam-TS150 type hot table applicable to a temperature measurement range of 0-1500K, the linkam THMSG-600 type cold and hot table can achieve temperature measurement analysis, and therefore the inclusion can be determined.

The following describes a method of determining the permeability.

According to Darcy' S law Q/S ═ k Δ P/η L, where Q is the flow (m)³S); s is the cross-sectional area (m) of the sample²) (ii) a L is the sample length (m); η is the fluid viscosity coefficient (Pa · s); k is the permeability (m)²) (ii) a Δ P is the pressure difference (Pa) between the upstream and downstream of the sample. Applying a stable pressure difference deltap at the upstream and downstream ends of the rock sample, and obtaining the permeability by measuring the flow Q through the sample, or measuring the pressure difference deltap at the upstream and downstream ends while maintaining a constant flow Q.

In the embodiment of the invention, a periodic oscillation method can be adopted for measuring the ultra-low permeability, and the lower limit of the measurement can reach 10^-10μ m 2. The periodic oscillation method was originally proposed by Kranz et al by reference to thermal diffusivity measurement and was applied to the measurement of permeability of low permeability rock. Specifically, the test core acts as an attenuation and damper. Firstly, a sine oscillation pressure wave with consistent amplitude and frequency is input at one end of the core, a sine wave with the same frequency and attenuated and delayed amplitude and phase is output at the other end of the core, and the amplitude attenuation and the phase delay are related to the rock permeability. The propagation process of a sinusoidal pressure wave in the core resembles one-dimensional diffusion:

in formula (1): p is the pore pressure; k is the permeability; μ is the viscosity coefficient of the fluid; beta is a_sIs the specific reservoir flow rate of the core-fluid system. The initial and boundary conditions of equation (1) are:

in formula (2): x is 0 is selected at the downstream end face; sd is the solution of the downstream vessel volume in combination with the initial and boundary conditions, formula (2) and formula (1)

As can be seen from equation (3), the periodic oscillation method can be considered as a hybrid method combining the steady state method and the pulse decay method. The response of the downstream pressure to the upstream stress consists of two parts, namely the sudden presence of the stress due to the pressure

A transient response (also decaying exponentially) that becomes sinusoidal, and a steady state response that also exists with a stable sinusoidal period as seen in the first half of equation (3), relative to the upstream pressure, the downstream response

Attenuation should occur in amplitude (attenuation factor a) and phase delay theta. The permeability of the sample can be calculated by measuring the upstream and downstream amplitude ratio alpha and the phase delay theta.

According to the findings of the Fischer study, the upstream and downstream pressure amplitude ratios α and the phase delay θ can be expressed as a function of 2 dimensionless parameters ψ and γ:

in formula (4):

in formula (5): γ, ψ are 2 dimensionless values generated in the superposition process. From the above various parameters, a permeability of

In formula (6): t is the oscillation period(s) of the pore pressure.

As can be seen from the above equation, 2 dimensionless values ψ and γ are obtained from the amplitude ratio α and the phase shift θ measured in the laboratory, and the permeability k can be obtained from the equation (6).

In the present example, the rock heat transfer derivative (thermal conductivity) was determined in the following manner.

There are many methods for measuring thermal conductivity at normal temperature and pressure, including two major methods, the steady state method and the transient state method, the heat flow meter method, the protective hot plate method, the circular tube method, the hot wire method, the flash method and the transient plane heat source method. The ideal heat conductivity measurement method should be performed in situ, but because the laboratory measurement method is simple and has high measurement precision, a large amount of rock heat conductivity measurement is performed in a laboratory. In laboratory measurements, the rock sample should be subjected to conditions of original moisture and temperature.

The steady state method refers to that the temperature distribution on a sample to be measured in the experiment is measured after being stabilized, and the starting point of the analysis is a steady-state heat conduction differential equation which can directly measure the heat conduction coefficient. The method is characterized in that the experimental formula is simple, and the heat conduction quantity (directly or indirectly) and the temperature of a plurality of points need to be measured. The disadvantages are that the test period is too long, the sample loading is complex, and the method is difficult to be applied to the experiment needing to test a large number of samples, so that the number of people using the method is small.

The steady state method corresponds to a transient method, and is a relatively common thermal conductivity measurement method. The change of the temperature of a sample along with time in the measurement process needs to be measured at a plurality of points on the sample along with time. The analysis method has the advantages that the analysis starting point is a transient heat conduction differential equation, the characteristics are that the formula is complex, but the sample assembly is simple, the test period is short, the transient method can simultaneously measure the heat conductivity and the thermal diffusion coefficient of the sample, and the requirement on the boundary condition of the sample temperature is not high.

The measurement principle is that a certain thermal interference is applied to a sample in a thermal equilibrium state, the response (temperature or heat flow change along with time) of the sample to the thermal interference is measured, and then the value of the thermophysical property parameter is determined according to a response curve.

The thermal conductivity coefficient is measured by a thermal constant analyzer, and a thermal conductivity equation is deduced and solved by assuming that a probe is positioned in an infinite sample.

When the thermal constant analyzer probe is electrically heated, the equation for the resistance rise over time is given as (1-1):

R＝R₀{1+α[ΔT_i+△T(τ)]} (1-1)

in the formula: r-sensor resistance, R₀Is the resistance of the sensor before it is heated;

α -Temperature Coefficient of Resistance (TCR);

ΔT_i-temperature differential of the nickel and the insulating thin layer;

Δ T (τ) — the sample surface temperature rise value.

The value of the temperature increase recorded by the sensor from equation (2-1) is shown in equation (1-2):

ΔT_iin effect, represents the difference between the recorded temperature rise value of the sensor and the temperature rise value of the surface of the sample.

Temperature difference Δ T_iAt a very short time Δ t_iThen becomes constant, and can be estimated as shown in (1-3) as follows:

in the formula: delta-thickness of the insulating layer, mm;

K_ithermal diffusivity, mm, of the insulating layer material²/s。

The increase in the surface temperature of the sample as a function of time is shown in (1-4):

in the formula: p is a radical of₀-sensor heating power, W;

a-radius of the probe, mm;

λ -the heat transfer coefficient of the test sample, W/(m.K);

d (τ) — a size-independent time-dependent equation, such as (1-5).

In the formula: t-the time, s, corresponding to the transient record;

Θ -characteristic time, defined as (1-6).

From the formula (1-4), it can be seen that the average temperature rise Δ T (τ) of the sample surface is linearly related to D (τ), whereThe slope of the straight line is used, and the heat conductivity coefficient of the sample can be obtained through the slope. By repeatedly testing the k value, the average temperature rise delta T (tau) and D (tau) are in strict linear correspondence, and the heat conduction coefficient can be obtained. And then the volume specific heat is obtained through the relationship among the heat conductivity coefficient, the thermal diffusion coefficient and the volume specific heat.

Before the numerical simulation is carried out on the three-dimensional model, the properties of media such as rocks are assigned, and the mechanical parameters of the rocks are obtained by taking a certain gold deposit as an example and testing according to the method shown in the following table 1:

TABLE 1

And obtaining initial geological parameters at the beginning of the mineralization process and mineralization geological parameters of each rock mass medium during final mineralization according to field statistics and sample test inversion.

The structural analysis and the microscopic structure research of the gold deposit show that before the gold deposit is formed, the northwest-southeast (NW-SE) undergoes toughness-brittleness left-line shear deformation under the extrusion condition, and the differential stress is 61.37-111.09 MPa. Finite element numerical simulation of the gold deposit in the north section of the Jiajia gold mine zone applying NW-SE is carried out, wherein the main compression stress is 100Mpa, and the NE-SW shearing force is 40 Mpa. As shown in table 2.

TABLE 2

The pressure range of the formed quartz capsule of the gold deposit is 1082-. The maximum value is taken as the simulated minimum pore pressure value of the new city, namely 337500000, 8.4e 7-3.375 e8 pa.

The test data can be used as parameters, initial conditions and boundary conditions of numerical simulation to participate in operation.

In the embodiment of the invention, a block model established by three-dimensional geological modeling software is imported into numerical simulation software, and by taking FLAC3D software as an example, an inp format block file is imported, and the code of the inp format block file is as follows:

model new

model configure thermal fluid

zone import'gjl.INP'

zone group'GJL'

zone import'll.inp'

zone group'LL'range group'GJL'not

zone import'fault.INP'

zone group'fault'range group'LL'not group'GJL'not

zone import'fault-b.INP'

zone group'fault-b'range group'fault'not group'LL'not group'GJL'not

zone import'fault-t.INP'

zone group'fault-t'range group'fault-b'not group'fault'not group'LL'not group'GJL'not

zone attach by-face。

in the embodiment of the invention, for the mineralization numerical simulation, rock is taken as a typical application object, so a typical molar coulomb model (zone cmodel assign mohr-coulomb) is selected.

In the embodiment of the invention, mechanical parameters, fluid mechanical parameters and the like are given to different medium materials in the model, and corresponding parameters are given to different media according to the obtained attribute values. Examples are as follows:

zone initialize density 2680range group'LL'

zone initialize density 2650range group'GJL'

zone initialize density 2810range group'fault-b'

zone initialize density 2810range group'fault-t'

zone initialize density 2000range group'fault'

zone prop bulk 52e9 shear 34e9 coh 5.3e7 fri 52ten 7.5e6 dil 2range group'LL'

zone prop bulk 55e9 shear 38e9 coh 5.4e7 fri 53ten 8.2e6 dil 2range group'GJL'

；zone prop bulk 52e9 shear 34e9 coh 5.3e7 fri 52ten 7.5e6 dil 2range group'fault-t'

；zone prop bulk 55e9 shear 38e9 coh 5.4e7 fri 53ten 8.2e6 dil 2range group 'fault-b'

zone prop bulk 57e9 shear 41e9 coh 4.2e7 fri 54ten 7.6e6 dil 2range group'fault-b'

zone prop bulk 57e9 shear 41e9 coh 4.2e7 fri 54ten 7.6e6 dil 2range group 'fault-t'

zone prop bulk 9e9 shear 4e9 coh 0.7e7 fri 40ten 0.9e6 dil 2range group 'fault'

zone fluid cmodel assign isotropic

zone fluid property permeability 5e-12range group'LL'

zone fluid property porosity 0.5range group'LL'

zone fluid property permeability 5e-12range group'GJL'

zone fluid property porosity 0.5range group'GJL'。

in the embodiment of the invention, initial conditions are also assigned, such as structural stress, movement speed, void pressure and fluid temperature, and corresponding objects of the model are assigned according to the obtained initial conditions, such as structural stress assignment, temperature field assignment and fluid seepage initial pressure assignment. Examples are as follows:

zone fluid cmodel assign anisotropic

zone fluid property permeability-xz 1e-6range group'fault-t'

zone fluid property porosity 0.9range group'fault-t'

zone fluid property permeability-xz 1e-2range group'fault-b'

zone fluid property porosity 0.3range group'fault-b'

zone fluid property permeability-xz 1e-4range group'fault'

zone fluid property porosity 0.9range group'fault'

zone fluid property permeability-xz 3e-10range group'LL'

zone fluid property porosity 0.1range group'LL'

zone fluid property permeability-xz 2e-15range group'GJL'

zone fluid property porosity 0.1range group'GJL'。

in the embodiment of the invention, boundary conditions are set, such as pressure-temperature, pore pressure value and the like of a mining area.

And (3) carrying out inversion according to the information such as corresponding pressure-temperature-pore pressure and the like of the rock body in the mineralization period to obtain corresponding parameters and then assigning the boundary conditions of the model. Examples are as follows:

zone fluid biot on

zone initialize fluid-density 1000

zone gridpoint initialize fluid-modulus＝2.18e9

zone gridpoint initialize fluid-tension 0

zone gridpoint initialize saturation 0

；zone face apply pore-pressure 7e8 range group'fault'position-z-1501-1499

；zone face apply pore-pressure 7e8 range group'fault-t'position-z-1501-1499

zone face apply pore-pressure 7e20 range group'fault-b'position-z-1501-1499

zone face APPLY stress-x 100000000range group'LL'position-x 663 665

zone face apply stress-x 100000000range group'GJL'position-x 5280 5282；

zone face APPLY stress-y-40000000range group'GJL'position-y 30553057；

zone face apply stress-y-40000000range group'LL'position-y-1 1；

zone face APPLY stress-y-40000000range group'fault-b'position-y 30553057；

zone face apply stress-y-40000000range group'fault-t'position-y-1 1；

zone gridpoint fix velocity range position-y-1 1

zone gridpoint fix velocity range position-y 3055 3057

zone gridpoint fix velocity range position-z-1501-1499

zone thermal cmodel isotropic

zone thermal property conductivity 1.5expansion 1e-6specific-heat 1000range group'LL'

zone thermal property conductivity 1.3expansion 2e-6specific-heat 1000range group'GJL'

zone thermal property conductivity 1expansion 3e-6specific-heat 1000range group'fault-b'

zone thermal property conductivity 1.1expansion 3e-6specific-heat 1000range group'fault-t'

zone thermal property conductivity 0.7expansion 8e-6specific-heat 1000range group'fault'；

zone gridpoint initialize temperature 20

zone face apply temperature 5000range group'fault-b'position-z-1501-1499。

after the assignment is completed, the simulation is started. Specifically, a step of numerical simulation is set, the step length is set from small to large, a dynamic mineralization simulation motion picture is realized, when the actual mineralization characteristics are the same as certain iterative times, the current parameters are recorded, and the iterative times are readjusted in the simulation, so that the final simulation numerical model accords with the three-dimensional mineralization model. Examples of the recorded results are as follows:

model fluid active on

model thermal active on

model mechanical active on；

model thermal timestep fix 6.48e3

model gravity 10

model solve ratio 1e-5

；model solve cycles 5000。

in the embodiment of the invention, after the simulation is finished, the simulated three-dimensional structural stress value and the three-dimensional fluid state data are led out for post-processing. And (3) exporting the numerical simulation result by using a list command: the Zone grid point position derives position information (XYZ coordinates) of each data, and the Zone list stress derives pressure data. Parameters such as fluid saturation, fluid pressure and the like can be derived by using list, and the derived data is in txt format. And adjusting the data into GOCAD data import, selecting parameters meeting the mineralization requirement, and performing Boolean operation to obtain a target area.

In the embodiment of the invention, the exported simulation data is imported into the three-dimensional model and compared with the three-dimensional model of the existing ore body, if the exported simulation data does not meet the requirement, the simulation iteration times are adjusted until the model meets the requirement; when the simulation model is consistent with the actual ore body model, according to the ore formation control factors, such as a structural crushing zone (structural degeneration simulation model), a fluid osmotic pressure model and the like, the simulation model is input into the three-dimensional geological model for post-processing, under the constraint of a proper threshold value, a favorable ore formation prediction area is obtained, and Boolean operation is carried out to obtain the optimal ore formation prediction area.

Fig. 4 is a schematic structural diagram of a composition of a mineralization simulator in an mineralization research area according to an embodiment of the present invention, and as shown in fig. 4, the mineralization simulator in the mineralization research area according to the embodiment of the present invention includes:

the building unit 40 is used for building a three-dimensional geological model based on geological mineralization patterns and plane and section geological data;

an obtaining unit 41, configured to obtain at least a numerical simulation parameter in the geological model by a sampling experiment testing manner, where the numerical simulation parameter at least includes at least one of: rock mechanics parameters, tectonic dynamics parameters, hydromechanics parameters, and thermodynamic parameters;

the assignment unit 42 is configured to introduce the three-dimensional rock mass-structure-channel geological model into a geological numerical simulation parameter assignment application, select and determine a numerical simulation constitutive model, and assign at least rock mechanical parameters, structure kinetic parameters, fluid mechanical parameters and thermodynamic parameters to different medium materials in the three-dimensional rock mass-structure-channel geological model;

the setting unit 43 is configured to obtain an initial geological parameter at the beginning of a mineralization process and an mineralization geological parameter of each rock medium during final mineralization according to statistics and sample test inversion of on-site geological parameters of an mineralization research area, and assign an initial condition of the mineralization research area based on the initial geological parameter; setting boundary conditions of the mineralization research area according to the mineralization geological parameters of each rock mass medium;

the simulation unit 44 is used for setting a numerical simulation step in a mode that the step length is set from small to large, performing dynamic mineralization simulation, generating a corresponding simulation dynamic graph, recording the current parameters of the simulation when the actual mineralization characteristics are the same as certain iterative times, and readjusting the iterative times in the simulation to enable the final simulation numerical model to conform to the three-dimensional mineralization model;

the deriving unit 45 is used for deriving the three-dimensional structural stress value and the three-dimensional fluid state data of the numerical simulation after the numerical simulation is finished;

and the mineralization determining unit 46 is used for selecting parameters meeting mineralization requirements from the three-dimensional structural stress values and the three-dimensional fluid state data, performing Boolean operation to obtain a target area, and taking the target area as an mineralization area in an mineralization research area.

As an implementation manner, the establishing unit 40 is further configured to:

acquiring drilling core sample bit data, gallery exploration sample bit data and surface chemical exploration data of a geological ore carrying rock mass, an ore forming rock mass and surrounding rocks based on X-ray fluorescence spectrum analysis and a main trace element device so as to determine the space position of the ore forming body;

In an embodiment of the invention, the rock mechanical parameter comprises at least one of: density, shear modulus, bulk modulus, cohesion, internal friction angle, expansion angle, and tensile strength.

The formation kinetic parameters include at least one of: the movement speed and the tectonic stress of the rock mass;

correspondingly, the obtaining unit 41 is further configured to:

and comprehensively constructing the moving direction, the scale distance and the moving time to obtain the moving average speed of the mineralization research area.

As an implementation manner, the obtaining unit 41 is further configured to:

In an exemplary embodiment, the establishing Unit 40, the obtaining Unit 41, the assigning Unit 42, the setting Unit 43, the simulating Unit 44, the deriving Unit 45, and the mineralization determining Unit 46 may be implemented by one or more Central Processing Units (CPUs), Graphics Processing Units (GPUs), Baseband Processors (BPs), Application Specific Integrated Circuits (ASICs), the DSP, the Programmable Logic Device (PLD), the Complex Programmable Logic Device (CPLD), the Field Programmable Gate Array (FPGA), the general processor, the Controller, the Microcontroller (MCU), the Microprocessor (Microprocessor), or other electronic components, and may be implemented in combination with one or more Radio Frequency (RF) antennas, so as to perform the steps of the network data collection method of the foregoing embodiments.

In the embodiment of the present disclosure, the specific manner in which each unit in the mineralization simulation apparatus of the mineralization study area shown in fig. 4 performs operations has been described in detail in the embodiment related to the method, and will not be described in detail herein.

An embodiment of the present invention further provides an electronic device, including: a processor and a memory for storing processor executable instructions, wherein the processor is configured to be capable of performing the steps of the mineralization simulation method of the mineralization study area of the aforementioned embodiments when the executable instructions in the memory are invoked.

The embodiment of the invention also describes a computer-readable storage medium, which has computer instructions stored thereon, and is characterized in that the instructions are executed by a processor to implement the steps of the mineralization simulation method for the mineralization research area of the foregoing embodiment.

In this embodiment, the at least one processor may constitute any physical device having circuitry to perform logical operations on one or more inputs. For example, at least one processor may include one or more Integrated Circuits (ICs) including an Application Specific Integrated Circuit (ASIC), a microchip, a microcontroller, a microprocessor, all or a portion of a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), a Digital Signal Processor (DSP), a Field Programmable Gate Array (FPGA), or other circuitry suitable for executing instructions or performing logical operations. The instructions executed by the at least one processor may be preloaded into a memory integrated with or embedded in the controller, for example, or may be stored in a separate memory. The memory may include Random Access Memory (RAM), Read Only Memory (ROM), hard disk, optical disk, magnetic media, flash memory, other permanent, fixed, or volatile memory, or any other mechanism capable of storing instructions. Optionally, the at least one processor may comprise more than one processor. Each processor may have a similar structure, or the processors may have different configurations that are electrically connected or disconnected from each other. For example, the processor may be a separate circuit or integrated in a single circuit. When more than one processor is used, the processors may be configured to operate independently or cooperatively. The processors may be coupled electrically, magnetically, optically, acoustically, mechanically or by other means allowing them to interact.

In the present embodiment, the non-transitory computer readable storage medium may be a ROM, a Random Access Memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, and the like.

Furthermore, the features and benefits of the present invention are described with reference to exemplary embodiments. Accordingly, the invention is expressly not limited to these exemplary embodiments illustrating some possible non-limiting combination of features which may be present alone or in other combinations of features.

Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the invention being indicated by the following claims.

24页详细技术资料下载

上一篇：一种医用注射器针头装配设备

下一篇：一种多天体飞越探测轨道设计方法

Ore formation simulation method and device for ore formation research area

相关技术

网友询问留言