41 results about "Earth crust" patented technology

The present disclosure relates to the technical field of stability control of surrounding rock of tunnelling, and provides an excavation compensation method for tunnelling in deep rock engineering, including: acquiring an engineering geological information of the tunnelling in deep rock engineering; determining an engineering hazard type based on the engineering geological information; determining an excavation compensation support strategy for the surrounding rock of the tunnelling in deep rock engineering based on the engineering hazard type; and performing a supplementary support control on the surrounding rock of the tunnelling in deep rock engineering based on the excavation compensation support strategy. Through the supplementary support strategy, the difference value between a radial stress of the surrounding rock of the tunnelling in deep rock engineering and an initial crustal stress can be reduced to within a preset approximate value range, further effectively preventing a stress concentration phenomenon occurred in a tangential stress.

The invention discloses a three-dimensional elastic block-fault motion deformation coupling model construction method and system. Comprising the following steps: acquiring geodetic survey data of earth crust movement between earthquakes and an active fault distribution map, and dividing different active blocks; the movable blocks are associated with geodetic measurement data through spatial position attributes, and motion features of different blocks are distinguished; adjacent boundaries are extracted to serve as modeling boundary faults, and discretization processing is carried out; constructing a three-dimensional elastic block-fault motion deformation coupling model including block rotation, internal strain, fault normal velocity discontinuity counteracting, inclined fault positive slippage and negative dislocation response; inverting parameters of the three-dimensional elastic block-fault motion deformation coupling model by adopting a constrained least square method; and according to the coupling coefficient, identifying the potential pregnancy and earthquake concave-convex body on the fault. According to the method, a more complete inter-earthquake three-dimensional elastic model is established, and support is provided for accurate research of active fault coupling characteristics, recognition of gestational earthquake concave-convex bodies and judgment of major earthquake dangerous areas.

The invention discloses a physical simulation method and device for seismic response certainty of an urban underground rail transit network, and the method comprises the steps: constructing a mixed seismic source model based on the active fault parameters and kinematic parameters of a target region, the seismic source model integrates a low-wavenumber deterministic slip component and a high-wavenumber random slip component on a fault fracture surface to generate a broadband seismic wave field source reflecting the physical process of a seismic source; the constructed seismic source model is used as input, a regional crustal velocity structure model is combined, propagation of seismic waves in a layered geological medium is simulated through a frequency-wavenumber domain method, and a seismic oscillation space-time field is obtained through calculation; and establishing a finite element model of the urban underground rail transit network, applying the seismic oscillation time-space field as non-consistent seismic oscillation input to the finite element model, and solving through dynamic time-history analysis to obtain the dynamic response of the whole underground rail transit network in the earthquake process. According to the invention, a complete physical chain from seismic source fracture to structural response can be reproduced.

The invention discloses a new-type barite deposit mineralization mode construction method based on a same sedimentary fault. The method belongs to the field of geological exploration and mineral resource prediction, and comprises the following steps: analyzing the earth crust cracking activity history of an earthquake-Hangjian broken basin, and recovering a basin'checkerboard 'structure prototype; by means of field geological survey and research, drilling engineering, geophysics, remote sensing geological big data and the like, north-east-east-direction and north-west-direction co-sedimentary fault networks are identified and modeled, and a deep metallogenic fluid rising and jet flow channel is determined; the metallogenic process is concluded into a hydrothermal jet flow stage and an underwater gravity flow deposition stage, deep fluid is revealed to rise, jet flow and precipitate along a fault, and an ore-rich zone of a composite seabed fan body group is formed; and a geochemical prospecting prediction model combining trace element and sulfur isotope characteristics is established. The method systematically clarifies a coupling mechanism of ore control of the same sedimentary fault, ore-forming fluid migration and gravity flow deposition, and provides a theoretical basis and a technical means for prospecting prediction of the barite deposit in the Hendan-Hangjian fault basin.

The utility model relates to the technical field of seepage heat transfer simulation, and discloses a high-efficiency seepage heat transfer simulation device in a high-temperature and high-pressure environment, which comprises a supporting mechanism for stable supporting, and a seepage sealing assembly for simulating a rock-soil layer in an earth crust is arranged at the upper end of the supporting mechanism. The upper end of the seepage sealing assembly is provided with a sealing heating monitoring assembly used for heating and monitoring fluid to be percolated, and the upper end of the sealing heating monitoring assembly is provided with a pressurizing mechanism used for pressurizing the fluid to be percolated in the sealing heating monitoring assembly. Simulation operation of seepage fluid is achieved, damage to equipment under the high-pressure condition is avoided, the driving mechanism accurately adjusts pressure through gear transmission, stable operation of a pressurization system and coordinated operation of the whole system are ensured, the precision of seepage and heat transfer simulation is improved, and the operability and safety of experiments are effectively improved.

The utility model discloses an earth crust deformation detection device convenient to install, which comprises a top frame, two sliding blocks are connected in the top frame in a sliding manner, two convex buttons are connected in the two sliding blocks in a sliding manner, and clamping holes are uniformly formed in the outer part of the top frame. The two supporting rods are separated, meanwhile, the two clamping blocks are clamped in the two bases respectively, then the inserting rods jointly penetrate through the bases and the circular holes in the clamping blocks, the clamping blocks are clamped and limited in the bases, the top frame can be erected and assembled conveniently and rapidly, meanwhile, the clamping plates are shifted left and right, and the top frame can be conveniently and rapidly erected and assembled. The two clamping plates can jointly clamp and limit the detection device body, the sliding blocks are synchronously controlled to slide and move in the top frame when the clamping plates are shifted, the sliding blocks are clamped in the clamping holes through the convex buttons in the process, the sliding blocks are limited after sliding in the top frame, the convenience of the whole device is improved, and the detection efficiency is improved. In the process, the whole device can be quickly assembled without using a wrench tool.

The invention discloses a sedimentary basin well-free area earth heat flow prediction method, device and equipment, and the method comprises the steps: determining the thickness of each stratum layer of a sedimentary basin well-free area based on seismic data, a sedimentary cover layer, a base structure and earth crust layering features; determining the rock heat generation rate of the sedimentary cap layer and the wrinkled substrate based on the rock core sample and / or natural gamma logging data of the basin mature exploration area; based on the seismic longitudinal wave velocity, determining rock heat generation rates of the lower earth crust and the middle earth crust; determining the rock heat generation rate of the upper earth crust based on the rock heat generation rate of the wrinkled substrate and the rock heat generation rate of the middle earth crust; based on the valance heat flow values of the upheaval area and the depression area of the basin mature exploration area and the structural evolution characteristics of the basin where the sedimentary basin well-free area is located, determining the valance heat flow of the sedimentary basin well-free area; according to the method, a new thought is provided for finely depicting the current ground temperature field characteristics of the sedimentary basin well drilling deficient region.

The invention discloses a metal intercalation boron-carbon superconducting material GeB2C2 under normal pressure and a preparation method of the metal intercalation boron-carbon superconducting material GeB2C2. The material has a two-dimensional Kagome lattice structure, the space group is P6 / mmm (No.191), the lattice constants are a = b = 2.732, and c = 16.784. Through calculation and prediction of a first principle, the material has a superconducting transition temperature of about 48 K under a normal pressure (0 GPa) condition, and is dynamically stable. The superconducting mechanism is derived from a strong electroacoustic coupling effect, and the coupling constant lambda is about 1.64. The material is composed of Ge, B and C elements rich in earth crust, the raw materials are easy to obtain, the cost is low, and the large-scale preparation potential is achieved. Compared with a traditional normal-pressure superconducting material, the material disclosed by the invention has a relatively high superconducting transition temperature; compared with a high-voltage superconductor, a high-voltage environment is not needed, and the preparation and application threshold is remarkably reduced. The material can be applied to the fields of high-performance superconducting filters, superconducting nanowire single-photon detectors, small superconducting magnets and the like, and has important scientific research values and good practical application prospects.

The application provides a crust deformation research and analysis system and method based on GNSS data, relates to the technical field of crust deformation monitoring and research, and comprises a hardware subsystem and a software subsystem; the hardware subsystem comprises a GNSS receiving module, a data storage module, a communication transmission module and a power supply module which cooperate with each other, and is used for realizing the collection, local storage and remote stable transmission of GNSS data; the software subsystem runs on a terminal device and comprises a data preprocessing module, a data solving module, an analysis module and a visual display module which are sequentially associated and cooperate with each other, and is used for realizing an integrated analysis process of GNSS original data preprocessing, data solving, crust deformation characteristic analysis and result visual display. The application can realize the integration of GNSS data collection, processing, analysis and display, improve the data processing efficiency, simplify the operation process, and provide convenient technical support for crust deformation research.

The invention discloses a method for multi-source data joint identification and inversion of a middle and lower earth crust high conductor, and relates to the technical field of geophysical exploration. The method for multi-source data joint identification and inversion of the middle and lower earth crust high conductor mainly comprises the following steps: carrying out potential field separation on Bouguer gravity anomaly data to obtain residual gravity anomaly, regional gravity anomaly and field source depth, and further determining earth crust residual gravity anomaly and middle and lower earth crust density disturbance characteristics; and calculating residual gravity low anomalies reflecting the middle and lower earth crust high conductors in the research area, delineating plane distribution of the middle and lower earth crust high conductors, and performing random sub-domain three-dimensional automatic inversion for the residual gravity low anomalies to obtain a three-dimensional space distribution form of the middle and lower earth crust high conductors. By implementing the multi-source data combined identification and inversion method for the mid-lower earth crust high conductor provided by the invention, the mid-lower earth crust high conductor can be rapidly and effectively identified and inverted, and the three-dimensional space distribution of the mid-lower earth crust high conductor can be obtained.

The invention relates to the technical field of geology, in particular to a geological density information stripping determination method and device, electronic equipment and a medium, and the method comprises the steps: obtaining surface observation data of a target area; the earth surface observation data are input into the stripping model, the stripping model outputs density distribution of a plurality of material layers of a target area, a lower interface of an earth surface shallow layer model of the stripping model target area starts to extend downwards, a new interface in each extension and a lower interface of a previous material layer form a new material layer, and for each material layer, the density distribution of a plurality of material layers in the stripping model target area is determined; performing density filling on the material layer according to the global crustal structure model, optimizing the density distribution of the material layer until the geoid of the material layer is matched with the benchmark of the earth surface observation data, and continuing to extend downwards until the downward extending distance reaches the target distance; and generating geological density information of the target area according to the density distribution of the plurality of material layers. Therefore, the problem that geological density information measured by related technologies is not accurate enough is solved.

The invention provides a fracture blocking rate inversion method based on GNSS earth crust deformation observation and elastic dislocation theory constraint, and the method comprises the steps: constructing a fracture three-dimensional model and an elastic dislocation constraint model, and calculating a fracture blocking rate. Through multi-period result comparison, strain field verification and earthquake activity verification, a high-risk region of pregnancy and earthquake can be identified. According to the method, the GNSS earth crust deformation observation elastic dislocation theory is adopted to constrain the negative dislocation model, and the precision and stability of fracture blocking rate inversion are remarkably improved. According to the invention, multiple periods of GNSS data are integrated to realize dynamic inversion, and the evolution trend of the locking state can be captured. According to the method, secondary blocks are finely divided, strain rate field constraints are introduced, locking distribution of a local strain enhancement area is accurately identified, the reliability of judgment of a large earthquake risk area is improved, a fracture high-locking-rate area is defined, through theoretical constraint and multi-source data fusion, the fracture locking rate inversion precision and the dynamic capture capability are improved, and the method is suitable for large-scale popularization and application. And reliable technical support is provided for earthquake risk assessment.

The application discloses a method, device and equipment for determining crust shortening amount and a storage medium, and belongs to the technical field of geophysical exploration. The method comprises the following steps: acquiring gravity data and geological data of a detection area comprising a compression structure area and a gravity background area; determining a first gravity anomaly integral amount of the detection area based on the gravity data of the detection area; determining a gravity contrast model corresponding to a specified area based on the geological data of the detection area, and further calculating a second gravity anomaly integral amount of the detection area; and calculating the crust shortening amount of the detection area based on the first gravity anomaly integral amount and the second gravity anomaly integral amount. According to the method, the gravity data is combined with the deformation of the compression structure based on the gravity anomaly of the crust compression structure, so that the crust shortening amount can be obtained simply, quickly, accurately and reliably, and reliable prior information can be provided for compression structure research and geophysical structure modeling.

This invention relates to the field of engineering earthquake and disaster prevention engineering technology, and particularly to a method and system for predicting the probability of fault deformation in engineering fault resistance. The method includes acquiring basic data of the target active fault, including the fault's geometric characteristics, kinematic characteristics, upper magnitude limit, recurrence period, elapsed time, focal depth, and deformation width; inputting the basic data as input parameters into a computational model; calculating the fault length based on the fault's geometric characteristics; and constructing a two-dimensional computational grid for along and across the fault. This invention provides a novel probabilistic prediction method for assessing the permanent displacement of the shallow crust caused by fault displacement during large earthquakes, as well as the intense deformation of the shallow crust near the fault. It can provide sub-meter level data on fault-alternating dislocation and cross-fault deformation at important engineering sites under a certain exceedance probability.

The invention relates to a method and a device for determining the mixed dyeing proportion of earth crust components in magma. The method comprises the following steps: dating a sample to be detected to obtain the crystallization age t of the magma; performing Nd and Sr isotope testing on the whole rock to obtain an epsilon Nd (t) value and an ISr value, generating a scatter diagram, and determining the epsilon Nd (t) value and the ISr value of the initial magma and the surrounding rock; performing quantitative simulation calculation on the change of the Sr-Nd isotope of the initial magma under the joint control of separation crystallization and surrounding rock mixed dyeing; and generating a scatter diagram by the epsilon Nd (t) values of the magma and ISr under different evolution processes to obtain an AFC evolution trend line of the Sr-Nd isotope of the magma, comparing the AFC evolution trend line with data points of a sample to be detected, and reading the crustal component mixed dyeing proportion according to a calculation result. According to the method, the magma Sr-Nd isotope is taken as a carrier, two magma evolution effects of mineral separation crystallization and crustal mixed dyeing are quantitatively described based on the AFC model, and quantitative estimation of the crustal mixed dyeing proportion in the AFC process is realized.

FIELD: mining industry.SUBSTANCE: invention is intended for use in selecting locations for the placement of coal waste dumps during the development of natural energy resource deposits and the disposal of accumulated solid waste from coal mining. The method for selecting locations for coal waste dumps is that in a mining region, using geodynamic zoning, the boundaries of geodynamically active blocks of the earth's crust are identified and the width of the zone of influence of the boundaries of such blocks is determined. Then the coal waste dump site is placed by placing the centre of the site equidistant from the boundaries of the zone of influence of geodynamically active blocks, which represent linear sections of tectonic destruction of the rock massif, and orienting the longitudinal axis of symmetry of the site and its long sides in the east-west direction. After which a coal waste dump is formed. The coal waste dump site is placed near the heat consumer.EFFECT: integrated use of energy resources and improved environmental condition.1 cl, 1 dwg, 1 ex

The invention discloses a non-decoupling rapid calculation method for a seismic electromagnetic field caused by a piezomagnetic effect, and belongs to the technical field of geophysics, and the method comprises the steps: S1, converting a frequency space domain piezomagnetic effect three-dimensional seismoelectric coupling equation into a frequency wavenumber domain through introducing a surface harmonic coordinate system, and carrying out the joint solving of an elastic wave kinetic equation and an electromagnetic field Maxwell equation; s2, solving the three groups of coupled equations in a frequency wavenumber domain to obtain corresponding earthquake and electromagnetic wave fields; and S3, transforming the wave field from a frequency wavenumber domain to a time-space domain through a two-time independent domain transformation method. According to the calculation method, it is proposed that the piezomagnetic effect of ferromagnetic minerals in the earth crust can cause the electromagnetic phenomenon in the earthquake fracture process, a complete three-dimensional piezomagnetic effect seismoelectric coupling equation is given, displacement, stress, an electric field and a magnetic field in a layered model can be obtained at the same time through combined solving, and finally, when conversion is carried out back to a frequency-space domain, the electromagnetic phenomenon can be generated in the earthquake fracture process. Wave number integration only needs to be carried out once.

This invention belongs to the field of seismic information processing technology and discloses a method and system for constructing a three-dimensional model of seismogenic faults based on aftershock precision location data. The method includes: acquiring seismic geological data of the study area; establishing a high-resolution three-dimensional crustal velocity model of the study area and determining precise source parameters; obtaining tectonic stress field characteristics through source mechanism solutions and GNSS observation data; and studying the seismogenic mechanism of the study area and obtaining information on potential strong earthquake hazards. This invention combines two grid intervals, 0.2°×0.2° and 0.1°×0.1°, to obtain both the velocity structure of the non-source area with a 0.2°×0.2° grid interval and the velocity structure of the source area with a 0.1°×0.1° grid interval. This is suitable for studying the velocity structure of source areas with a relatively small number of stations.

This application discloses a method for extracting the residual density field of the Earth's crust based on the separation of magnetic mineral density effects, relating to the field of geophysics. The method includes: obtaining the three-dimensional density and magnetization structure of a target region; obtaining the three-dimensional temperature structure of the target region; converting magnetization into magnetite content; and stripping away the density effect of magnetite to obtain the residual density of the Earth's crust. This invention, through quantitative estimation and separation of the magnetic mineral density effect, explicitly removes the contribution of magnetic minerals to the density field from the total density structure, avoiding interference from magnetic mineral enrichment in the interpretation of crustal density anomalies. It estimates the volume fraction of magnetite using magnetic information and combines it with the temperature-corrected density structure to establish a crustal density component decomposition method based on a volume fraction mixing model, achieving physical differentiation between magnetic and non-magnetic components in the density field. It does not rely on specific gravity or magnetic inversion algorithms and has good versatility and scalability.

This invention relates to the field of gas-sensitive materials technology, specifically to an improved location method and system for Pn-Pg cross-arrival time difference, comprising the following steps: acquiring Pg and Pn phase data of a reference earthquake and the earthquake to be measured; determining the array epicentral distance range of Pn and Pg based on crustal thickness and focal depth; calculating the theoretical arrival times of Pg and Pn at different epicentral distances using ray theory based on a velocity model; obtaining the arrival time difference between theoretical and measured data; obtaining arrival time correction parameters based on the reference earthquake; establishing an inversion objective function based on the corrected linearized inversion model; performing causal analysis on the earthquake to be inverted using the optimal solution; and performing error analysis. By using the cross-arrival time data of Pn and Pg phases to determine the earthquake focal location, not only can the calculation efficiency be greatly improved, but the reliability of the inversion results can also be improved.

By analyzing crustal deformation data, we can clearly identify precursory anomalies such as pre-slips that precede major earthquakes, and use this data to predict earthquakes and estimate the optimal epicenter location from among potential epicenters. [Solution] The crustal deformation vector calculation unit calculates the difference vector of the change in position information at each observation station at a fixed interval such as 5 minutes, and calculates this as the crustal deformation vector for each observation station. The correlation value calculation unit performs processing based on the crustal deformation vectors for each observation station output by the crustal deformation vector calculation unit. The correlation value calculation unit selects one station as the central observation station and selects several observation stations in its vicinity as surrounding observation stations. In a certain situation, the correlation value calculation unit can select any number of surrounding observation stations from among observation stations whose distance from the central observation station is within a predetermined distance.

The invention relates to the technical field of municipal engineering, in particular to a fabricated structure of a pipeline in an underground comprehensive pipe gallery in an earth crust active area, comprising a water pipe, a first flange plate and a second flange plate fixedly sleeve two ends of the water pipe respectively, bolts are arranged in a plurality of threaded ports, and a branch pipe is fixedly mounted on the surface of the water pipe; a plurality of spring rods are fixedly installed on the surface of the first flange plate in a circumferential shape, a mounting plate is fixedly installed at one ends of the spring rods, a penetrating opening is formed in the surface of the mounting plate, a hollow rod is fixedly installed in the penetrating opening, and a suction cup is fixedly installed at the open end of the hollow rod. After the ends of the multiple water pipes are spliced and assembled, the multiple suckers can be controlled to tightly abut against and attract the ends of the multiple bolts respectively, so that the effect of limiting rotation of the bolts is achieved, the situation that the bolts are affected by external vibration acting force to rotate, and consequently the splicing and assembling positions of the ends of the water pipes are loosened is avoided, and the service life of the water pipes is prolonged. And the water pipe splicing and assembling stability is improved.

The invention relates to the technical field of computers, in particular to a method, a device, equipment, a medium and a program product for determining an underwater gravity field, and the method comprises the steps: obtaining digital topographic data and seawater dynamic topographic data of a target area; generating a surface shallow three-dimensional space model of the target area according to the digital topographic data, the seawater dynamic topographic data and the reference ellipsoid continuation data, and performing density filling on the surface shallow three-dimensional space model by using a refined crustal density model; gravitational potential information generated by the target area is calculated based on the surface shallow three-dimensional space model after density filling, and an underwater gravitational field of the target area can be calculated by combining an external gravitational field model. Therefore, by comprehensively applying multi-source data and a high-precision calculation method, the calculation efficiency of the underwater gravity field and the accuracy of the result are improved.

The application discloses a kind of ultra-deep reflection imaging method and related device based on passive source seismic data, it is related to deep structure imaging technical field, the method is related to the following steps: first, for original continuous seismic waveform data, the seismic waveform data of P-wave first arrival and S-wave non-arrival are preprocessed, to obtain teleseismic event data section.Subsequently, high signal-to-noise ratio teleseismic event is identified using time-frequency-wave number domain analysis, and those data greatly affected by noise are removed, then, using multidimensional deconvolution technique, the reconstruction of reflection wave field is carried out to high-quality teleseismic data segment and deconvolution source wavelet, to obtain ultra-deep reflection wave field data.Finally, reflection imaging processing is carried out in combination with velocity analysis and stacking imaging technology, the analysis ability of deep reflection interface is optimized, and finally the deep reflection imaging profile of earth is obtained, the continuous structural interface extending from shallow crust to deep mantle can be clearly revealed, and the coupling characteristics between shallow and deep earth structure can be effectively captured.

The application provides an evaluation method for favorable exploration targets of submarine hydrogen energy in an ultra-extended tectonic environment, and belongs to the field of submarine hydrogen energy exploration and development. The application integrates an evaluation index system for inorganic serpentinization hydrogen generation mechanism in a "hydrogen source-migration-reservoir-seal-thermal structure" five-in-one manner, so as to achieve the design purpose of significantly reducing exploration cost and risk and providing an industrialization target optimization tool for submarine hydrogen energy exploration of an ultra-extended passive continental margin. Specifically, by using existing seismic data and regional geological conditions, combining with geological environment and favorable reservoir-seal combination conditions required for submarine hydrogen energy generation, and aiming at tectonic units with hydrogen generation capacity in different basins, the development range of serpentinization is determined by comprehensively analyzing tectonic environment, extension, lithospheric change, large-scale detachment fault and deep fault characteristics, identifying upper and lower crust thickness and Moho depth, and combining temperature field and seismic reflection anomaly.

The utility model relates to the field of earth crust stress strain observation, in particular to an indoor complete machine calibration device for a vertical drilling strain gauge, which comprises a mounting frame, a pressurizing assembly is arranged on one side of the mounting frame and used for changing the vertical stress of surrounding rock, and a mounting plate is fixedly connected to one side of the mounting frame and used for fixing the mounting plate. One side of the mounting plate is fixedly connected with a positioning assembly, the cylindrical surface of the positioning assembly is provided with a driving mechanism, the driving mechanism is in threaded connection with the positioning assembly, the cylindrical surface of the positioning assembly is rotatably connected with a fixing device, and the fixing device is matched with the driving mechanism. And meanwhile, indoor calibration of the whole machine is effectively performed.

The invention provides a simulation material for a lithosphere scale model and application of the simulation material. The lithosphere scale model comprises an upper crust, a lower crust, a lithosphere ground curtain and a soft flow circle ground curtain. The simulation material of the supercrust comprises the following raw materials: silica powder, quartz sand and carborundum; the simulation materials of the lower earth crust and the lithosphere mantle comprise the following raw materials: silica gel, quartz sand and carborundum; the simulation material of the soft flow ring mantle comprises any one of syrup, honey or a sodium iodide solution. The simulation material provided by the invention can quantitatively complete a lithosphere scale structure physical simulation experiment, and scans a lithosphere scale model through an industrial CT technology to obtain a structure evolution process image, so that the structure evolution characteristics are analyzed, the rheological structure and layered stretching deformation characteristics of the lithosphere can be better simulated, and the simulation material can be applied to the field of lithosphere engineering. And an effective technical means is provided for researching development and evolution of lithosphere structural deformation.

This invention provides a method for generating a construction scheme for a CM pile composite foundation, comprising: acquiring multi-source data of a target construction area, including target building data, underground rock strata data, historical weather data, crustal activity data, and man-made facility data; generating a first construction scheme based on the target building data and underground rock strata data; revising the first construction scheme based on the historical weather data, crustal activity data, and man-made facility data to obtain a second construction scheme; and iteratively optimizing the second construction scheme based on preset variables, preset constraints, and optimization objectives to obtain a target construction scheme. This method rapidly generates a second construction scheme from multi-source data and iteratively optimizes it to obtain the target construction scheme, resulting in high efficiency in scheme generation and high completeness of the target construction scheme obtained through iterative optimization.

In a logging identification method for a small fault based on stress anisotropy and crustal stress types, a magnitude and direction of a current crustal stress are determined by means of logging calculation of the crustal stress, and a transformation location of a crustal stress type is determined; a development location of the small fault is determined according to the transformation location of the crustal stress type; and an included angle between the current crustal stress and the fault is determined according to an amplitude of variation of stress anisotropy. The present disclosure provides the logging identification method for a small fault based on stress anisotropy and crustal stress types. According to the method, prediction cost is low, operability is high, and a prediction result has a reference value for guiding exploration and development of a complex oil and gas reservoir and construction of a gas storage.

The present invention addresses the problem of clearly identifying precursory anomalies such as precursory slippage (pre-slip) which are precursors of an earthquake from crustal movement data, predicting an earthquake, and inferring the most appropriate epicenter position among a plurality of epicenter candidates. [Solution] The present invention is capable of acquiring, at prescribed intervals (e.g. 5 minute intervals), crustal movement vectors relative to given reference points defined for respective monitoring stations from data obtained at reception stations by means of a GNSS satellite 102, analyzing said crustal movement vectors, and detecting or predicting an anomaly in a crustal movement vector for each value. An anomaly detection device 120 comprises a virtual epicenter position database 210, a monitoring station position data database 220, a virtual epicenter position-to-monitoring station distance calculation unit 230, a crustal movement vector calculation unit 240, a correlation value calculation unit 250, an epicenter selection unit 260, and an anomaly assessment unit 270.