CN117805923A - Fracture plane position identification method and system based on Tanh function - Google Patents
Fracture plane position identification method and system based on Tanh function Download PDFInfo
- Publication number
- CN117805923A CN117805923A CN202410174976.4A CN202410174976A CN117805923A CN 117805923 A CN117805923 A CN 117805923A CN 202410174976 A CN202410174976 A CN 202410174976A CN 117805923 A CN117805923 A CN 117805923A
- Authority
- CN
- China
- Prior art keywords
- gravity anomaly
- total horizontal
- horizontal derivative
- gravity
- fracture plane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000005484 gravity Effects 0.000 claims abstract description 130
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 230000002159 abnormal effect Effects 0.000 claims abstract description 23
- 238000009933 burial Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 5
- 239000011707 mineral Substances 0.000 abstract description 5
- 238000011160 research Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a fracture plane position identification method and system based on a Tanh function, comprising the following steps: utilizing the Bragg gravity anomaly of the target area to obtain the total horizontal derivative gradient of the gravity anomaly; processing the gravity anomaly total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity anomaly total horizontal derivative gradient which is free of information influence outside the boundary of the geologic body; according to the gravity anomaly total horizontal derivative gradient influenced by the information except the geologic body boundary, the gravity anomaly representing the fracture plane position distribution is obtained by means of calculation of a Tanh function, and then the fracture plane position distribution characteristics of the target area are obtained. According to the invention, the abnormal information generated at the edge of the geologic body is enhanced and sharpened by using the total horizontal derivative gradient method, the influence of different burial depths on the abnormal total horizontal derivative gradient of gravity is balanced by using the Tanh function, the identification result of the fracture plane position of the target area with higher resolution can be obtained, and the invention can provide basis for basic geological research, mineral and oil-gas resource exploration and the like.
Description
Technical Field
The invention relates to the technical field of heavy and magnetic exploration, in particular to a fracture plane position identification method and system based on a Tanh function.
Background
Fracture is one of main control factors of mineral products and oil and gas reservoirs, and plays an important role in mineral product and oil and gas resource exploration. The boundary line of the geologic body such as fracture can control the two sides of the boundary line to form different construction units, the continuity of the original stratum is damaged, the density or magnetism of the two sides of the boundary line is transversely different, the abnormal step band or abnormal dislocation characteristic is shown on the heavy magnetic force abnormality, and the fracture plane position identification by utilizing the heavy magnetic field method has unique advantages.
The method for identifying the fracture plane position by using the gravity and magnetic potential fields can be divided into three categories of numerical statistics, numerical calculation and other three categories. The mathematical statistics type edge recognition method comprises two types of small subdomain filtering and normalized standard deviation; other methods mainly include a directional filtering method, an image enhancement technology, a magnetic anomaly component conversion method and a method for carrying out edge recognition by using a neural network; the numerical calculation type edge recognition method is the most studied and widely applied edge recognition method, and a plurality of numerical calculation type edge recognition methods existing at present are mainly designed based on vertical derivatives and horizontal derivatives. Evjen (1936) first uses the zero value of the vertical first derivative (VDR) for edge recognition, but this approach may cause the creation of false boundaries in complex models (Liu et al 2023); cordell (1979) proposed edge recognition using the maximum value of the Total Horizontal Derivative (THDR), rosest et al (1992) proposed a method of edge recognition using the maximum value of the resolved signal amplitude (ASA, total gradient modulus) consisting of the horizontal derivative and the vertical derivative, but the total horizontal derivative and the resolved signal amplitude edge recognition method had a problem in that the resolution was not high. The edge recognition is carried out by using the higher derivative, so that a sharper signal can be obtained at the edge of the geologic body, the resolution of the edge recognition result is improved, but the edge information of different buried depth field sources cannot be displayed at the same time.
Disclosure of Invention
In order to solve the problems, the invention aims to provide a fracture plane position identification method and system based on a Tanh function.
In order to achieve the above object, the present invention provides the following solutions:
a fracture plane position identification method based on a Tanh function comprises the following steps:
step 1: extracting the Bragg gravity anomaly of the target area according to the terrain elevation and the free space gravity anomaly data of the target area;
step 2: utilizing the Bragg gravity anomaly of the target area to obtain the total horizontal derivative gradient of the gravity anomaly;
step 3: processing the gravity anomaly total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity anomaly total horizontal derivative ladder with information influence outside the boundary of the geologic body removed;
step 4: according to the gravity anomaly total horizontal derivative gradient influenced by the information except the boundary of the geologic body, calculating by using a Tanh function to obtain gravity anomaly representing the distribution of the fracture plane position;
step 5: and carrying out fracture plane position identification on the target area by utilizing the gravity anomaly representing the fracture plane position distribution, so as to obtain the fracture plane position distribution characteristics of the target area.
Preferably, the step 2: the method for calculating the gravity anomaly total horizontal derivative gradient by utilizing the Bragg gravity anomaly of the target area comprises the following steps:
calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target area by using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents the gravity anomaly total horizontal derivative gradient at the target zone (x, y), THDR-VDR n (x, y) represents the vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents the total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location.
The invention also provides a fracture plane position identification system based on the Tanh function, which comprises the following steps:
the Bragg gravity anomaly extraction module is used for extracting Bragg gravity anomalies of the target area according to the terrain elevation and free space gravity anomaly data of the target area;
the gravity anomaly total horizontal derivative gradient solving module is used for solving the gravity anomaly total horizontal derivative gradient by utilizing the Bragg gravity anomaly of the target area;
the abnormal signal removing module is used for processing the gravity abnormal total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity abnormal total horizontal derivative ladder which is used for removing information influence outside the boundary of the geologic body;
the different burial depth anomaly information equalization module is used for calculating and obtaining gravity anomalies representing the distribution of the fracture plane positions by using a Tanh function according to the gradient of the total horizontal derivative of the gravity anomalies affected by the information except the boundary of the geologic body;
and the fracture plane position dividing module is used for identifying the fracture plane position of the target area by utilizing the gravity anomaly representing the fracture plane position distribution to obtain the fracture plane position distribution characteristic of the target area.
Preferably, the gravity anomaly total horizontal derivative gradient solving module comprises:
a derivative calculation unit for calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target region using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents the gravity anomaly total horizontal derivative gradient at the target zone (x, y), THDR-VDR n (x, y) represents the vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents the total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location.
The invention also provides a computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the steps in a fracture plane position identification method based on a Tanh function.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
compared with the prior art, the method and the system for identifying the fracture plane position based on the Tanh function enhance and sharpen abnormal information generated at the edge of the geologic body by utilizing a total horizontal derivative gradient method, and balance the influence of different burial depths on the total horizontal derivative gradient of the gravity anomaly by utilizing the Tanh function, so that a target area fracture plane position identification result with higher resolution can be obtained, and the method and the system can provide basis for basic geological research, mineral products, oil gas resource exploration and the like.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a fracture plane position identification method based on a Tanh function;
fig. 2 is a graph of recognition results of the distribution of the positions of breaking planes of islands according to the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In order to achieve the above object, the present invention provides the following solutions:
referring to fig. 1-2, a fracture plane position identification method based on a Tanh function includes:
step 1: extracting the Bragg gravity anomaly of the target area according to the terrain elevation and the free space gravity anomaly data of the target area;
step 2: utilizing the Bragg gravity anomaly of the target area to obtain the total horizontal derivative gradient of the gravity anomaly;
further, step 2 includes:
calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target area by using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents a gravity anomaly total horizontal derivative gradient at the target zone (x, y) location, THDR-VDRn (x, y) represents a vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents a total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location.
Step 3: processing the gravity anomaly total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity anomaly total horizontal derivative ladder with information influence outside the boundary of the geologic body removed;
step 4: according to the gravity anomaly total horizontal derivative gradient influenced by the information except the boundary of the geologic body, calculating by using a Tanh function to obtain gravity anomaly representing the distribution of the fracture plane position;
step 5: and carrying out fracture plane position identification on the target area by utilizing the gravity anomaly representing the fracture plane position distribution, so as to obtain the fracture plane position distribution characteristics of the target area.
Referring to fig. 1, the following description of the present invention further describes a fracture plane location recognition method based on the Tanh function according to the present invention with reference to a specific embodiment:
step 101: data is input. Mainly the data of the Bragg gravity anomaly.
In step 101, the position of the fracture plane is identified to need to input Bragg gravity anomaly data, the topographic elevation and free space gravity anomaly data are collected through a public website, the Bragg gravity anomaly is obtained by using a generalized topography correction technology (Lei-Hui ,1984; an Yulin, etc., 2010), and the data are all input in a grid file (. Grd) format.
Step 102: and (5) calculating a derivative. And obtaining the gravity anomaly total horizontal derivative gradient by using a frequency domain derivative calculation technology.
In step 102, for the Bragg gravity anomaly obtained in step 101, a frequency domain derivative calculation technology is utilized to calculate a Bragg gravity anomaly total level derivative (THDR) firstly, then calculate a vertical derivative and a total level derivative of the total level derivative to obtain a vertical derivative (THDR-VDR) of the gravity anomaly total level derivative and a total level derivative (THDR-THDR) of the gravity anomaly total level derivative respectively, and finally calculate a ratio of the vertical derivative of the gravity anomaly total level derivative to the total level derivative of the total level derivative to obtain a gravity anomaly total level derivative gradient (HG), wherein the result effectively enhances and sharpens anomaly information reflecting fracture, and improves the resolution of edge recognition results.
Step 103: and removing abnormal information outside the boundary of the geologic body. Based on the gravity anomaly total horizontal derivative gradient calculated in the step 102, the influence of anomaly information outside the boundary of the geologic body is removed by using a threshold calculation technology which is larger than zero.
In step 103, abnormal information other than the boundary of the geologic body is removed. Based on the gravity anomaly total horizontal derivative gradient calculated in the step 102, the influence of anomaly information outside the boundary of the geologic body is removed by using a threshold value calculation technology larger than zero, so that the obtained edge recognition result is simple and clear.
Step 104: different burial depths are balanced in abnormal information. The gravity anomaly total level derivative gradient which is obtained by the step 103 and is influenced by the anomaly information outside the boundary of the geologic body is removed, the total level derivative gradient information HG can clearly reflect the plane position of the fracture, but the fracture with larger burial depth or smaller scale generates smaller amplitude of the anomaly information and even can be covered by the anomaly information of other fractures; according to the invention, by utilizing the characteristic that the amplitude range of the Tanh function is 0-1 when the independent variable is greater than zero and the slope is larger when the independent variable is smaller, the abnormal information with smaller amplitude generated by inputting the fracture with larger burial depth or smaller scale can be amplified better, so that the influence of the abnormal amplitude generated by the burial depth of the edge of different geological bodies on the edge recognition result is balanced, and the abnormal information representing the fracture plane position distribution of different burial depths (scales), namely the gravity abnormal TANH-HG, is calculated.
Step 105: fracture partitioning. And (3) combining the gradient information of the gravity anomaly total horizontal derivative which is calculated in step 104 and balances the influences of different burial depths, and identifying the distribution characteristics of the fracture plane positions.
In the embodiment of the invention, according to the flow, an edge recognition method based on a Tanh function is utilized to obtain a fracture partitioning result of an island and an attached gravity anomaly TANH-HG is shown in figure 2.
From the above description, it can be seen that the following technical effects are achieved: according to the invention, the abnormal information generated at the edge of the geologic body is effectively enhanced and sharpened by using the total horizontal derivative gradient method, the influence of different burial depths on the abnormal total horizontal derivative gradient of gravity is balanced by using the Tanh function, the identification result of the fracture plane position of the target area with higher resolution can be obtained, and the basis can be provided for basic geological research, mineral products, oil gas resource exploration and the like.
The invention also provides a fracture plane position identification system based on the Tanh function, which comprises the following steps:
the Bragg gravity anomaly extraction module is used for extracting Bragg gravity anomalies of the target area according to the terrain elevation and free space gravity anomaly data of the target area;
the gravity anomaly total horizontal derivative gradient solving module is used for solving the gravity anomaly total horizontal derivative gradient by utilizing the Bragg gravity anomaly of the target area;
the abnormal signal removing module is used for processing the gravity abnormal total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity abnormal total horizontal derivative ladder which is used for removing information influence outside the boundary of the geologic body;
the different burial depth anomaly information equalization module is used for calculating and obtaining gravity anomalies representing the distribution of the fracture plane positions by using a Tanh function according to the gradient of the total horizontal derivative of the gravity anomalies affected by the information except the boundary of the geologic body;
and the fracture plane position dividing module is used for identifying the fracture plane position of the target area by utilizing the gravity anomaly representing the fracture plane position distribution to obtain the fracture plane position distribution characteristic of the target area.
Preferably, the gravity anomaly total horizontal derivative gradient solving module comprises:
a derivative calculation unit for calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target region using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents a gravity anomaly total horizontal derivative gradient at the target zone (x, y) location, THDR-VDRn (x, y) represents a vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents a total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location. Compared with the prior art, the fracture plane position identification system based on the Tanh function has the same beneficial effects as the fracture plane position identification method based on the Tanh function in the technical scheme, and is not repeated herein.
The invention also provides a computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the steps in a fracture plane position identification method based on a Tanh function. Compared with the prior art, the beneficial effects of the computer readable storage medium provided by the invention are the same as those of the fracture plane position identification method based on the Tanh function in the technical scheme, and are not repeated here.
The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In summary, the present description should not be construed as limiting the invention.
Claims (5)
1. A fracture plane position identification method based on a Tanh function is characterized by comprising the following steps:
step 1: extracting the Bragg gravity anomaly of the target area according to the terrain elevation and the free space gravity anomaly data of the target area;
step 2: utilizing the Bragg gravity anomaly of the target area to obtain the total horizontal derivative gradient of the gravity anomaly;
step 3: processing the gravity anomaly total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity anomaly total horizontal derivative ladder with information influence outside the boundary of the geologic body removed;
step 4: according to the gravity anomaly total horizontal derivative gradient influenced by the information except the boundary of the geologic body, calculating by using a Tanh function to obtain gravity anomaly representing the distribution of the fracture plane position;
step 5: and carrying out fracture plane position identification on the target area by utilizing the gravity anomaly representing the fracture plane position distribution, so as to obtain the fracture plane position distribution characteristics of the target area.
2. The method for identifying the fracture plane position based on the Tanh function according to claim 1, wherein the step 2: the method for calculating the gravity anomaly total horizontal derivative gradient by utilizing the Bragg gravity anomaly of the target area comprises the following steps:
calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target area by using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents the gravity anomaly total horizontal derivative gradient at the target zone (x, y), THDR-VDR n (x, y) represents the vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents the total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location.
3. A fracture plane location identification system based on a Tanh function, comprising:
the Bragg gravity anomaly extraction module is used for extracting Bragg gravity anomalies of the target area according to the terrain elevation and free space gravity anomaly data of the target area;
the gravity anomaly total horizontal derivative gradient solving module is used for solving the gravity anomaly total horizontal derivative gradient by utilizing the Bragg gravity anomaly of the target area;
the abnormal signal removing module is used for processing the gravity abnormal total horizontal derivative gradient by using a threshold value calculation method which is larger than zero to obtain a gravity abnormal total horizontal derivative ladder which is used for removing information influence outside the boundary of the geologic body;
the different burial depth anomaly information equalization module is used for calculating and obtaining gravity anomalies representing the distribution of the fracture plane positions by using a Tanh function according to the gradient of the total horizontal derivative of the gravity anomalies affected by the information except the boundary of the geologic body;
and the fracture plane position dividing module is used for identifying the fracture plane position of the target area by utilizing the gravity anomaly representing the fracture plane position distribution to obtain the fracture plane position distribution characteristic of the target area.
4. A fracture plane location identification system based on a Tanh function according to claim 3, wherein the gravity anomaly total horizontal derivative gradient solving module comprises:
a derivative calculation unit for calculating a gravity anomaly total horizontal derivative gradient based on the Bragg gravity anomaly of the target region using a frequency domain derivative calculation method; the calculation formula of the gravity anomaly total horizontal derivative gradient is as follows:
wherein HG (x, y) represents the gravity anomaly total horizontal derivative gradient at the target zone (x, y), THDR-VDR n (x, y) represents the vertical derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location, and THDR-THDR (x, y) represents the total horizontal derivative of the gravity anomaly total horizontal derivative at the target zone (x, y) location.
5. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of a fracture plane location identification method based on a Tanh function as claimed in any one of claims 1-2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410174976.4A CN117805923B (en) | 2024-02-07 | 2024-02-07 | Fracture plane position identification method and system based on Tanh function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410174976.4A CN117805923B (en) | 2024-02-07 | 2024-02-07 | Fracture plane position identification method and system based on Tanh function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117805923A true CN117805923A (en) | 2024-04-02 |
| CN117805923B CN117805923B (en) | 2024-07-09 |
Family
ID=90431668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410174976.4A Active CN117805923B (en) | 2024-02-07 | 2024-02-07 | Fracture plane position identification method and system based on Tanh function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117805923B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017200692A1 (en) * | 2016-05-20 | 2017-11-23 | Exxonmobil Research And Engineering Company | Shape-based geophysical parameter inversion |
| CN116861347A (en) * | 2023-05-22 | 2023-10-10 | 青岛海洋地质研究所 | Magnetic force abnormal data calculation method based on deep learning model |
-
2024
- 2024-02-07 CN CN202410174976.4A patent/CN117805923B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017200692A1 (en) * | 2016-05-20 | 2017-11-23 | Exxonmobil Research And Engineering Company | Shape-based geophysical parameter inversion |
| CN116861347A (en) * | 2023-05-22 | 2023-10-10 | 青岛海洋地质研究所 | Magnetic force abnormal data calculation method based on deep learning model |
Non-Patent Citations (3)
| Title |
|---|
| WANG WANYIN, PAN YU, QIU ZHIYUN: "A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential fi eld data", APPLIED GEOPHYSICS, vol. 6, no. 3, 22 October 2009 (2009-10-22), pages 226 - 233 * |
| 许文强;袁炳强;刘必良;姚长利;: "多种重磁位场边缘识别方法及在南黄海北部断裂构造识别中的应用研究", 物探与化探, no. 04, 6 August 2020 (2020-08-06) * |
| 魏华敬等: "位场数据处理算法的数据挖掘试验与应用", 实验室研究与探索, vol. 38, no. 01, 15 January 2019 (2019-01-15), pages 9 - 15 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117805923B (en) | 2024-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dannemiller et al. | A new method for determination of magnetization direction | |
| Carranza | Objective selection of suitable unit cell size in data-driven modeling of mineral prospectivity | |
| Gumiel et al. | Analysis of the fractal clustering of ore deposits in the Spanish Iberian Pyrite Belt | |
| CN101504659B (en) | Method for extracting geoscience spatial information based on generalized self-similarity principle | |
| Pardo‐igúzquiza et al. | Morphometric analysis of karst depressions on a Mediterranean karst massif | |
| CN107966732A (en) | The seismic properties change rate acquiring method being oriented to based on space structure | |
| CN105572736A (en) | Quantitative prediction method for slope break belt by utilizing seismic data | |
| Yang et al. | 3D gravity inversion with optimized mesh based on edge and center anomaly detection | |
| Liu et al. | Fractal and multifractal characteristics of lineaments in the Qianhe Graben and its tectonic significance using remote sensing images | |
| Zhu et al. | Normalized vertical derivatives in the edge enhancement of maximum-edge-recognition methods in potential fields | |
| Kebede et al. | Processing and interpretation of full tensor gravity anomalies of Southern Main Ethiopian Rift | |
| CN117805923B (en) | Fracture plane position identification method and system based on Tanh function | |
| CN109254320A (en) | Seismic properties optimization and sand body Overlay District prediction technique based on image processing method | |
| Florio et al. | A fractional vertical derivative technique for regional-residual separation | |
| CN116561536B (en) | Landslide hidden danger identification method, terminal equipment and medium | |
| Zhang et al. | Quantification and classification of grainflow morphology on natural dunes | |
| CN108090929B (en) | Novel method for analyzing and extracting linear abnormity of mining area | |
| Li et al. | Statistics and spectra of turbulence under different roughness in the near‐neutral atmospheric surface layer | |
| CN110989032B (en) | Gravity horizontal total gradient fracture identification method based on inclination angle | |
| Drăguţ et al. | Automated classification of topography from SRTM data using object-based image analysis | |
| CN113900141B (en) | Oil-gas distribution prediction method and device | |
| Ventura et al. | Geomorphological map of the Somma-Vesuvius volcanic complex (Italy) | |
| CN103439755A (en) | Method for identifying geologic-body boundaries through use of potential-field high-order horizontal gradient model | |
| de Souza et al. | Fourier domain vertical derivative of the nonpotential squared analytical signal of dike and step magnetic anomalies: A case of serendipity | |
| Laffan et al. | Predicting regolith properties using environmental correlation: a comparison of spatially global and spatially local approaches |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |